FREE ENRG HSEN ASEN / in the event of a Newton failure, take 100 timesteps of 1.E-6 TIME 100 1.00E-6 / begin on a uniform mesh of 6 points NPTS 16 / definition of the computational interval XEND 10.0 /XEND 0.5 /XCEN 0.35 /WMIX 10.0 / pressure and inlet mass flow rate PRES 1.0 (atmospheres) FLRT 0.250 (g/cm**2-sec) / adaptive mesh criteria /GRAD 0.99 /CURV 0.99 /GRAD 0.2 /CURV 0.5 GRAD 0.1 CURV 0.1 / unreacted mole fractions TFIX 400. MOLE REAC N2 0.79 REAC O2 0.21 REAC H2 0.42 REAC H2O 0.00 / estimated products PROD H2O 0.42 PROD N2 0.79 PROD H2 0.00 PROD O2 0.00 / estimated intermediate mole fractions INTM H2O2 .00001 INTM HO2 .001 INTM H .0001 INTM O .0001 / tolerances for the Newton iteration ATOL 1.E-10 RTOL 1.E-4 / tolerances for the time step Newton iteration ATIM 1.E-5 RTIM 1.E-5 / print control PRNT 1 / given temperature profile TEMP 0. 300.0 TEMP .1250 484.5 TEMP .250 583.7 TEMP .375 1672.2 TEMP .5 1753.5 TEMP .75 1901.4 TEMP 1.0 2027. TEMP 1.25 2120. TEMP 1.5 2184. TEMP 2.0 2260. TEMP 3.0 2348. TEMP 6.0 2475. TEMP 10.0 2524. MULT TDIF RSTR END