Time integrated flux product yield

The photodissociation rate coefficients are included as source and sink terms in a system of time-dependent continuity equations for the atmosphere. Modem values for vertical (eddy) diffusion and solar photon flux are utilized. The system of 2nd-order ordinary differential equations is solved by integration, and yields chemical species abundances as a function of time and altitude. The isotope atmospheric chemistry includes only SO2 isotopologue photodissociation reactions and production of SO isotopologues. Additional isotopic reactions such as SO2 oxidation by OH, SO photolysis, SO disproportionation during self-reaction, and SO dimmer formation, have been neglected. My objective here is to focus only on SO2 photolysis as a S-MIF mechanism. 

The method of zonation was applied to the energy and material conservation equations. Based on centered finite difTerence approximations, this method can transform three partial differential equations in radial distance and time to ordinary differential equations in time only. Following this, the ordinary differential equations were solved by using Crank-Nicholson algorithm. On the basis of this, the volumetric fluxes of those tar-phase and total volatile phase components were integrated with time by using in roved Euler method to evaluate overall pyrolysis product yields, and afterwards the gas yield can be deduced. 

The product yield in the desired channel is related to the time-integrated flux... 

---