Hello,

The switch ICNTRL(15) allows us to call Update_RCONST() several times
during the KPP sub-timesteps. Some of the rate constants depend on other
concentrations. In GEOS-Chem for example, k(IONO+BrSALA) depends on
C(ind_IONO).

Why is Update_RCONST() using concentrations from the C array, which
still has the old values before the integration started? Don't we have
to use the intermediate values from the YNEW arrray inside
Update_RCONST()?