|
Last change
on this file since 4115 was
2696,
checked in by kanani, 7 years ago
|
|
Merge of branch palm4u into trunk
|
|
File size:
1.8 KB
|
| Line | |
|---|
| 1 | #include strato.spc |
|---|
| 2 | #include strato.eqn |
|---|
| 3 | |
|---|
| 4 | {#LANGUAGE Fortran77 |
|---|
| 5 | #INTEGRATOR rodas3 |
|---|
| 6 | #DRIVER general} |
|---|
| 7 | |
|---|
| 8 | #LOOKATALL |
|---|
| 9 | |
|---|
| 10 | #MONITOR O; |
|---|
| 11 | |
|---|
| 12 | #INITVALUES |
|---|
| 13 | |
|---|
| 14 | CFACTOR = 8.12E+16 ; |
|---|
| 15 | |
|---|
| 16 | O = 8.158574E-09 ; |
|---|
| 17 | O1D = 2.773328E-15 ; |
|---|
| 18 | O3 = 6.560000E-06 ; |
|---|
| 19 | NO = 1.074516E-08 ; |
|---|
| 20 | NO2 = 2.758840E-09 ; |
|---|
| 21 | NO3 = 5.667175E-14 ; |
|---|
| 22 | N2O5 = 6.905712E-12 ; |
|---|
| 23 | HNO3 = 3.486344E-10 ; |
|---|
| 24 | HNO4 = 3.827485E-12 ; |
|---|
| 25 | H2O = 6.100000E-06 ; |
|---|
| 26 | H = 4.190815E-14 ; |
|---|
| 27 | OH = 2.041165E-10 ; |
|---|
| 28 | HO2 = 1.439703E-10 ; |
|---|
| 29 | H2O2 = 7.227409E-11 ; |
|---|
| 30 | H2 = 3.700000E-07 ; |
|---|
| 31 | CH4 = 4.900000E-07 ; |
|---|
| 32 | CH3 = 3.984427E-17 ; |
|---|
| 33 | CH3O = 2.496504E-15 ; |
|---|
| 34 | CH3O2 = 5.130799E-12 ; |
|---|
| 35 | CH3OOH = 4.587317E-12 ; |
|---|
| 36 | HCO = 3.010774E-19 ; |
|---|
| 37 | CH2O = 9.913988E-11 ; |
|---|
| 38 | CO = 2.000000E-08 ; |
|---|
| 39 | Cl = 8.756027E-12 ; |
|---|
| 40 | Cl2 = 8.665679E-14 ; |
|---|
| 41 | ClO = 1.011874E-09 ; |
|---|
| 42 | ClOO = 2.659306E-15 ; |
|---|
| 43 | OClO = 7.193558E-14 ; |
|---|
| 44 | Cl2O2 = 7.103997E-14 ; |
|---|
| 45 | HCl = 2.157657E-09 ; |
|---|
| 46 | HOCl = 2.196369E-10 ; |
|---|
| 47 | ClONO2 = 2.166934E-11 ; |
|---|
| 48 | Br = 4.025093E-12 ; |
|---|
| 49 | BrCl = 5.843046E-14 ; |
|---|
| 50 | BrO = 1.132088E-11 ; |
|---|
| 51 | HBr = 4.349534E-13 ; |
|---|
| 52 | HOBr = 5.066596E-12 ; |
|---|
| 53 | BrONO2 = 9.405408E-14 ; |
|---|
| 54 | O2 = 0.209; |
|---|
| 55 | N2 = 0.781; |
|---|
| 56 | |
|---|
| 57 | #INLINE F77_INIT |
|---|
| 58 | TSTART = 12.D0*3600.D0 |
|---|
| 59 | TEND = TSTART + 24.D0*3600.D0 * 3 |
|---|
| 60 | DT = 3600.D0 |
|---|
| 61 | TEMP = 241.42 |
|---|
| 62 | #ENDINLINE |
|---|
| 63 | |
|---|
| 64 | #INLINE F90_INIT |
|---|
| 65 | TSTART = 12.D0*3600.D0 |
|---|
| 66 | TEND = TSTART + 24.D0*3600.D0 * 3 |
|---|
| 67 | DT = 3600.D0 |
|---|
| 68 | TEMP = 241.42 |
|---|
| 69 | #ENDINLINE |
|---|
| 70 | |
|---|
| 71 | #INLINE MATLAB_INIT |
|---|
| 72 | global TSTART TEND DT TEMP |
|---|
| 73 | TSTART = 12.0*3600.0; |
|---|
| 74 | TEND = TSTART + 24.0*3600.0 * 3; |
|---|
| 75 | DT = 3600.0; |
|---|
| 76 | TEMP = 241.42; |
|---|
| 77 | #ENDINLINE |
|---|
| 78 | |
|---|
| 79 | #INLINE C_INIT |
|---|
| 80 | TSTART = 12.0*3600.0; |
|---|
| 81 | TEND = TSTART + 24.0*3600.0 * 3; |
|---|
| 82 | DT = 3600.0; |
|---|
| 83 | TEMP = 241.42; |
|---|
| 84 | #ENDINLINE |
|---|
Note: See
TracBrowser
for help on using the repository browser.
