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