source: palm/trunk/UTIL/chemistry/gasphase_preproc/kpp/models/wet.def @ 3359

Last change on this file since 3359 was 2696, checked in by kanani, 7 years ago

Merge of branch palm4u into trunk

File size: 2.4 KB
Line 
1#include wet.spc
2#include wet.eqn
3
4{#LANGUAGE Fortran77
5#INTEGRATOR rodas3
6#DRIVER general}
7
8#LOOKATALL
9
10#MONITOR O3;
11
12#INITVALUES
13  CFACTOR = 2.55E+10; {ppb-to-mcm}
14  ALL_SPEC= 1.0E-10;
15{Variable species}
16  NO2        = 0.5 ;
17  NO         = 0.2 ;
18  O3         = 60.0 ;
19  HNO3       = 1.6594E-06 ;
20  HNO2       = 8.4898E-03 ;
21  H2O2       = 1.5 ;
22  CO         = 200 ;
23  HCHO       = 1.0 ;
24  ALD2       = 0.1 {3.2067E-02} ;
25  PAN        = 0.1 {1.5022E-02} ;
26  MGLY       = 7.3500E-03 ;
27  PAR        = 1.2 ;
28  OLE        = 6.6095E-03 ;
29  ETH        = 1.0 {2.4286E-02} ;
30  TOL        = 3.0184E-02 ;
31  CRES       = 6.4171E-03 ;
32  XYL        = 2.0E-02 ;
33  ROOH       = 7.0E-02 ;
34  SO2        = 3.3371E-02 ;
35  PNA        = 1.6736E-02 ;
36  OPEN       = 1.0E-03 ;
37  ISOP       = 1.0 ;
38  HONO       = 1.0E-3 ;
39  NO3        = 5.0E-05 ;
40  OH         = 6.0E-04 ;
41  HO2        = 2.8900E-02 ;
42  N2O5       = 4.6513E-05 ;
43  C2O3       = 6.1457E-05 ;
44  XO2        = 6.2020E-03 ;
45  XO2N       = 4.0E-3 ;
46  ROR        = 4.2461E-09 ;
47  TO2        = 8.6120E-07 ;
48  NH3        = 1.3637E-02 ;
49  CRO        = 3.1320E-06 ;
50  H2O        = 1.35e+7 ;
51  O2         = 2.1E+8 ;
52  CH4        = 1.7E+3 ;
53  CO2aq      = 3.0E-1 ;
54  HO2aq      = 4.2400E-04 ;
55  H2O2aq     = 9.6683E-01 ;
56  O3aq       = 4.8688E-06 ;
57  HCHOaq     = 2.0245E-02 ;
58  HCOOH      = 0.2 ;
59  HCOOHaq    = 2.0180E-03 ;
60  CH3O2Haq   = 1.1101E-04 ;
61  CH3O2aq    = 9.0994E-05 ;
62  HNO3aq     = 2.5563E-06 ;
63  NO3aq      = 5.2533E-09 ;
64  NH3aq      = 5.8021E-06 ;
65  SO2aq      = 3.0071E-07 ;
66  OHaq       = 4.0141E-05 ;
67  H2Oaq      = 4.3638E+05 ;
68  HPLUSaq    = 1.9654E-01 ;
69  OHaq       = 3.1392E-06 ;
70  HCOOMINaq  = 3.6325E-03 ;
71  O2MINaq    = 6.7840E-05 ;
72  HSO3MINaq  = 5.1188E-05 ;
73  SO32MINaq  = 3.4183E-08 ;
74  NH4aq      = 9.8636E-01 ;
75  NO3aq      = 3.9367E-01 ;
76  FE3PLUSaq  = 7.8548E-03 ;
77  FE2PLUSaq  = 3.1419E-02 ;
78  SO42MINaq  = 2.3564E+00 ;
79  CU2PLUSaq  = 2.1601E-04 ;
80
81#INLINE F77_INIT
82        TSTART = 12.D0*3600.D0
83        TEND = TSTART + 24.D0*3600.D0 * 5
84        DT = 3600.D0
85        TEMP = 288.15
86#ENDINLINE
87
88#INLINE F90_INIT
89        TSTART = 12.D0*3600.D0
90        TEND = TSTART + 24.D0*3600.D0 * 5
91        DT = 3600.D0
92        TEMP = 288.15
93#ENDINLINE
94
95#INLINE MATLAB_INIT
96  global TSTART TEND DT TEMP
97  TSTART = 12.0*3600.0;
98  TEND = TSTART + 24.0*3600.0*5;
99  DT = 3600.0;
100  TEMP = 288.15;
101#ENDINLINE
102
103#INLINE C_INIT
104  TSTART = 12.0*3600.0;
105  TEND = TSTART + 24.0*3600.0*5;
106  DT = 3600.0;
107  TEMP = 288.15;
108#ENDINLINE
109
Note: See TracBrowser for help on using the repository browser.