Last change
on this file since 3914 was
2696,
checked in by kanani, 7 years ago
|
Merge of branch palm4u into trunk
|
File size:
3.6 KB
|
Line | |
---|
1 | #include texas.spc |
---|
2 | #include texas.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 = 1. ; {Conversion Factor} |
---|
16 | NO2 = 46.0; |
---|
17 | NO = 30.0; |
---|
18 | O3P = 16.0; |
---|
19 | O3 = 48.0; |
---|
20 | NO3 = 62.0; |
---|
21 | N2O5 = 108.0; |
---|
22 | HNO3 = 63.0; |
---|
23 | O1D2 = 16.0; |
---|
24 | HO = 17.0; |
---|
25 | HONO = 47.0; |
---|
26 | HO2 = 33.0; |
---|
27 | CO = 28.0; |
---|
28 | HNO4 = 79.0; |
---|
29 | HO2H = 34.0; |
---|
30 | SO2 = 64.0; |
---|
31 | SULF = 98.0; |
---|
32 | C_O2 = 47.0; |
---|
33 | HCHO = 30.0; |
---|
34 | COOH = 48.0; |
---|
35 | MEOH = 32.0; |
---|
36 | RO2_R = 1.0; |
---|
37 | ROOH = 62.0; |
---|
38 | R2O2 = 1.0; |
---|
39 | RO2_N = 1.0; |
---|
40 | RNO3 = 133.0; |
---|
41 | MEK = 72.0; |
---|
42 | PROD2 = 100.0; |
---|
43 | CCO_O2 = 75.0; |
---|
44 | PAN = 121.0; |
---|
45 | CCO_OOH = 76.0; |
---|
46 | CCO_OH = 60.0; |
---|
47 | RCO_O2 = 89.0; |
---|
48 | PAN2 = 135.0; |
---|
49 | CCHO = 44.0; |
---|
50 | RCO_OOH = 90.0; |
---|
51 | RCO_OH = 74.0; |
---|
52 | BZCO_O2 = 137.0; |
---|
53 | PBZN = 183.0; |
---|
54 | BZ_O = 93.0; |
---|
55 | MA_RCO3 = 101.0; |
---|
56 | MA_PAN = 147.0; |
---|
57 | TBU_O = 73.0; |
---|
58 | ACET = 58.0; |
---|
59 | NPHE = 139.0; |
---|
60 | PHEN = 94.0; |
---|
61 | BZNO2_O = 139.0; |
---|
62 | HOCOO = 63.0; |
---|
63 | HCOOH = 46.0; |
---|
64 | RCHO = 58.0; |
---|
65 | GLY = 58.0; |
---|
66 | MGLY = 72.0; |
---|
67 | BACL = 86.0; |
---|
68 | CRES = 108.0; |
---|
69 | BALD = 106.0; |
---|
70 | METHACRO = 70.0; |
---|
71 | MVK = 70.0; |
---|
72 | ISOPROD = 70.0; |
---|
73 | DCB1 = 58.0; |
---|
74 | DCB2 = 72.0; |
---|
75 | DCB3 = 72.0; |
---|
76 | ETHENE = 28.0; |
---|
77 | ISOPRENE = 68.0; |
---|
78 | TRP1 = 136.0; |
---|
79 | ALK1 = 30.1; |
---|
80 | ALK2 = 36.7; |
---|
81 | ALK3 = 58.6; |
---|
82 | ALK4 = 77.6; |
---|
83 | ALK5 = 118.9; |
---|
84 | ARO1 = 98.6; |
---|
85 | ARO2 = 118.7; |
---|
86 | OLE1 = 72.3; |
---|
87 | OLE2 = 75.8; |
---|
88 | |
---|
89 | #INLINE F77_INIT |
---|
90 | TSTART = 12.0D0*3600.0D0 |
---|
91 | TEND = TSTART + 120.0D0*3600.0D0 |
---|
92 | DT = 3600.D0 |
---|
93 | TEMP = 300.0D0 |
---|
94 | #ENDINLINE |
---|
95 | |
---|
96 | #INLINE F90_INIT |
---|
97 | TSTART = 12.0d0*3600.0d0 |
---|
98 | TEND = TSTART + 120.0d0*3600.0d0 |
---|
99 | DT = 3600.d0 |
---|
100 | TEMP = 300.0d0 |
---|
101 | #ENDINLINE |
---|
102 | |
---|
103 | #INLINE MATLAB_INIT |
---|
104 | global TSTART TEND DT TEMP |
---|
105 | TSTART = 12*3600; |
---|
106 | TEND = TSTART + 120*3600; |
---|
107 | DT = 3600; |
---|
108 | TEMP = 300; |
---|
109 | #ENDINLINE |
---|
110 | |
---|
111 | #INLINE C_INIT |
---|
112 | TSTART = 12.0*3600.0; |
---|
113 | TEND = TSTART + 120.0*3600.0; |
---|
114 | DT = 3600.0; |
---|
115 | TEMP = 300.0; |
---|
116 | #ENDINLINE |
---|
Note: See
TracBrowser
for help on using the repository browser.