Last change
on this file since 3775 was
2696,
checked in by kanani, 7 years ago
|
Merge of branch palm4u into trunk
|
File size:
1.6 KB
|
Rev | Line | |
---|
[2696] | 1 | #include saprc99.spc |
---|
| 2 | #include saprc99.eqn |
---|
| 3 | |
---|
| 4 | {#LANGUAGE Fortran77 |
---|
| 5 | #INTEGRATOR rodas3 |
---|
| 6 | #DRIVER matlab |
---|
| 7 | #HESSIAN ON} |
---|
| 8 | |
---|
| 9 | #LOOKATALL |
---|
| 10 | |
---|
| 11 | #MONITOR O3; NO; NO2; ETHENE; |
---|
| 12 | |
---|
| 13 | #INITVALUES |
---|
| 14 | |
---|
| 15 | CFACTOR = 2.4476e+13; |
---|
| 16 | |
---|
| 17 | ALL_SPEC = 0.0e0; |
---|
| 18 | NO = 1.0e-1; |
---|
| 19 | NO2 = 5.0e-2; |
---|
| 20 | HONO = 1.e-3; |
---|
| 21 | SO2 = 5.e-2; |
---|
| 22 | HCHO = 1.121e-2; |
---|
| 23 | CCHO = 2.316e-3; |
---|
| 24 | RCHO = 1.72e-3; |
---|
| 25 | ACET = 5.07e-3; |
---|
| 26 | MEK = 3.26e-3; |
---|
| 27 | MEOH = 5.89e-3; |
---|
| 28 | GLY = 1.21e-4; |
---|
| 29 | MGLY = 8.37e-5; |
---|
| 30 | PHEN = 6.06e-4; |
---|
| 31 | CRES = 5.60e-4; |
---|
| 32 | BALD = 7.51e-5; |
---|
| 33 | METHACRO = 1.30e-3; |
---|
| 34 | ISOPROD = 8.93e-5; |
---|
| 35 | PROD2 = 1.93e-3; |
---|
| 36 | ETHENE = 1.89e-2; |
---|
| 37 | ISOPRENE = 4.33e-4; |
---|
| 38 | ALK1 = 1.167e-2; |
---|
| 39 | ALK2 = 1.88e-2; |
---|
| 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 | TERP = 8.20e-4; |
---|
| 48 | XC = 0.2E0; |
---|
| 49 | CCO_OH = 1.16e-3; |
---|
| 50 | RCO_OH = 3.92e-4; |
---|
| 51 | HCOOH = 6.77e-4; |
---|
| 52 | O3P = 7.843e-9; |
---|
| 53 | H2O = 2.0e+04; |
---|
| 54 | O2 = 2.09e+5; |
---|
| 55 | AIR = 1.0e+6; |
---|
| 56 | CH4 = 1.0e0; |
---|
| 57 | |
---|
| 58 | #INLINE F77_INIT |
---|
| 59 | TSTART = 12.0D0*3600.0D0 |
---|
| 60 | TEND = TSTART + 120.0D0*3600.0D0 |
---|
| 61 | DT = 3600.D0 |
---|
| 62 | TEMP = 300.0D0 |
---|
| 63 | #ENDINLINE |
---|
| 64 | |
---|
| 65 | #INLINE F90_INIT |
---|
| 66 | TSTART = 12.0d0*3600.0d0 |
---|
| 67 | TEND = TSTART + 120.0d0*3600.0d0 |
---|
| 68 | DT = 3600.d0 |
---|
| 69 | TEMP = 300.0d0 |
---|
| 70 | #ENDINLINE |
---|
| 71 | |
---|
| 72 | #INLINE MATLAB_INIT |
---|
| 73 | global TSTART TEND DT TEMP |
---|
| 74 | TSTART = 12*3600; |
---|
| 75 | TEND = TSTART + 120*3600; |
---|
| 76 | DT = 3600; |
---|
| 77 | TEMP = 300; |
---|
| 78 | #ENDINLINE |
---|
| 79 | |
---|
| 80 | #INLINE C_INIT |
---|
| 81 | TSTART = 12.0*3600.0; |
---|
| 82 | TEND = TSTART + 120.0*3600.0; |
---|
| 83 | DT = 3600.0; |
---|
| 84 | TEMP = 300.0; |
---|
| 85 | #ENDINLINE |
---|
| 86 | |
---|
| 87 | |
---|
Note: See
TracBrowser
for help on using the repository browser.