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