Appears in CMR's model where the steady state of the auxiliary variables
is not correct if set_auxiliary_variables routine is only called once.

Added a debug mode (false by default) to check that the mapping defined
by set_auxiliary_variables routine is invariant.