Last change
on this file since 4270 was
2696,
checked in by kanani, 7 years ago
|
Merge of branch palm4u into trunk
|
File size:
1.4 KB
|
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.