| Line | |
|---|
| 1 | #include saprcnov.spc |
|---|
| 2 | #include saprcnov.eqn |
|---|
| 3 | |
|---|
| 4 | {#LANGUAGE Fortran77 |
|---|
| 5 | #INTEGRATOR rodas3 |
|---|
| 6 | #DRIVER general} |
|---|
| 7 | |
|---|
| 8 | #LOOKATALL |
|---|
| 9 | #MONITOR O3; NO; NO2; SO2; |
|---|
| 10 | |
|---|
| 11 | #INITVALUES |
|---|
| 12 | |
|---|
| 13 | CFACTOR = 2.4476E+13; |
|---|
| 14 | |
|---|
| 15 | ALl_SPEC = 0.0E0; |
|---|
| 16 | NO = 1.0E-1; |
|---|
| 17 | NO2 = 5.0E-2; |
|---|
| 18 | HONO = 1.E-3; |
|---|
| 19 | SO2 = 5.E-2; |
|---|
| 20 | ISOPRENE = 4.33E-4; |
|---|
| 21 | TERP = 8.20E-4; |
|---|
| 22 | ETHENE = 1.89E-2; |
|---|
| 23 | MEOH = 5.89E-3; |
|---|
| 24 | HCOOH = 6.77E-4; |
|---|
| 25 | CCO_OH = 1.16E-3; |
|---|
| 26 | RCO_OH = 3.92E-4; |
|---|
| 27 | HCHO = 1.12E-2; |
|---|
| 28 | CCHO = 2.32E-3; |
|---|
| 29 | RCHO = 1.72E-3; |
|---|
| 30 | GLY = 1.21E-4; |
|---|
| 31 | MGLY = 8.37E-5; |
|---|
| 32 | METHACRO = 1.30E-3; |
|---|
| 33 | ISOPROD = 8.93E-5; |
|---|
| 34 | BALD = 7.51E-5; |
|---|
| 35 | ACET = 5.07E-3; |
|---|
| 36 | MEK = 3.26E-3; |
|---|
| 37 | PROD2 = 1.93E-3; |
|---|
| 38 | PHEN = 6.06E-4; |
|---|
| 39 | CRES = 5.60E-4; |
|---|
| 40 | ALK3 = 4.69E-2; |
|---|
| 41 | ALK4 = 4.17E-2; |
|---|
| 42 | ALK5 = 3.06E-2; |
|---|
| 43 | ARO1 = 1.18E-2; |
|---|
| 44 | ARO2 = 8.74E-3; |
|---|
| 45 | OLE1 = 1.04E-2; |
|---|
| 46 | OLE2 = 7.97E-3; |
|---|
| 47 | CH4 = 0.07E0; |
|---|
| 48 | H2O = 2.0E+04; |
|---|
| 49 | O2 = 2.09E+5; |
|---|
| 50 | AIR = 1.0E6; |
|---|
| 51 | |
|---|
| 52 | #INLINE F77_INIT |
|---|
| 53 | TSTART = 0.0D0 |
|---|
| 54 | TEND = TSTART + 2*24*3600.0D0 ! 2160.0D0*60.0D0 |
|---|
| 55 | DT = 3600.D0 |
|---|
| 56 | TEMP = 300.0D0 |
|---|
| 57 | #ENDINLINE |
|---|
| 58 | |
|---|
| 59 | #INLINE F90_INIT |
|---|
| 60 | TSTART = 0.0D0 |
|---|
| 61 | TEND = TSTART + 2*24*3600.0D0 ! 2160.0D0*60.0D0 |
|---|
| 62 | DT = 3600.D0 |
|---|
| 63 | TEMP = 300.0D0 |
|---|
| 64 | #ENDINLINE |
|---|
| 65 | |
|---|
| 66 | #INLINE MATLAB_INIT |
|---|
| 67 | global TSTART TEND DT TEMP |
|---|
| 68 | TSTART = 0.0; |
|---|
| 69 | TEND = TSTART + 2*24*3600.0D0; |
|---|
| 70 | DT = 3600.0; |
|---|
| 71 | TEMP = 300.00; |
|---|
| 72 | #ENDINLINE |
|---|
| 73 | |
|---|
| 74 | #INLINE C_INIT |
|---|
| 75 | TSTART = 0.0; |
|---|
| 76 | TEND = TSTART + 2*24*3600.0D0; /*+ 2160.0*60.0;*/ |
|---|
| 77 | DT = 3600.0; |
|---|
| 78 | TEMP = 300.00; |
|---|
| 79 | #ENDINLINE |
|---|
Note: See
TracBrowser
for help on using the repository browser.
