I changed each diagnostic to correspond to a single term in the parametrization, rather than one-and-a-half and then the sum of all four.
That is, the first diagnostic of the old version was the first term of the parametrization plus the part of the second term that does not uses variables at the future time in the implicit update step (the $\left.\frac{d(P-L)}{dr}\right|_0 (r - r_0)$ term is implicit in r; the first diagnostic included the $r_0$ part and excluded the $r(t+dt)$ part), and the second looked like the total photochemical tendency from the scheme, which the driver could accumulate on its own. In the new version, each diagnostic in run_o3prog_2015 corresponds to a term of the equation in the documentation.
I changed each diagnostic to correspond to a single term in the parametrization, rather than one-and-a-half and then the sum of all four.
That is, the first diagnostic of the old version was the first term of the parametrization plus the part of the second term that does not uses variables at the future time in the implicit update step (the $\left.\frac{d(P-L)}{dr}\right|_0 (r - r_0)$ term is implicit in
r
; the first diagnostic included the $r_0$ part and excluded the $r(t+dt)$ part), and the second looked like the total photochemical tendency from the scheme, which the driver could accumulate on its own. In the new version, each diagnostic inrun_o3prog_2015
corresponds to a term of the equation in the documentation.