|
Last change
on this file since 3477 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.
