| 1 |
|
| 2 |
#include "ctrparam.h" |
| 3 |
|
| 4 |
! ============================================================ |
| 5 |
! |
| 6 |
! CHEMINIT.F: Subroutine for initialization of tracers |
| 7 |
! in MIT Global Chemistry Model |
| 8 |
! |
| 9 |
! ------------------------------------------------------------ |
| 10 |
! |
| 11 |
! Author: Chien Wang |
| 12 |
! MIT Joint Program on Science and Policy |
| 13 |
! of Global Change |
| 14 |
! |
| 15 |
! ---------------------------------------------------------- |
| 16 |
! |
| 17 |
! Revision History: |
| 18 |
! |
| 19 |
! When Who What |
| 20 |
! ---- ---------- ------- |
| 21 |
! 052200 Chien Wang rev. |
| 22 |
! 080200 Chien Wang repack based on CliChem3 & add cpp |
| 23 |
! 101800 Chien Wang replaced if_3gases with cpp |
| 24 |
! 112800 Chien Wang add extra edaily for meta particles |
| 25 |
! 113000 Chien Wang add read n_total_urban & n_urban |
| 26 |
! 121800 Chien Wang move initialization of xnoxltnt here |
| 27 |
! 092001 Chien Wang add bc and oc |
| 28 |
! 051604 Chien Wang add back four missing reading lines |
| 29 |
! 062304 Chien Wang add drydep 46 data & delete block data |
| 30 |
! |
| 31 |
! ========================================================== |
| 32 |
|
| 33 |
! |
| 34 |
subroutine cheminit(ISTRT1,pt,qv) |
| 35 |
! ================================= |
| 36 |
|
| 37 |
#include "chem_para" |
| 38 |
#include "chem_com" |
| 39 |
#include "chem_meta" |
| 40 |
|
| 41 |
dimension pt(nlon,nlat,nlev) |
| 42 |
dimension qv(nlon,nlat,nlev) |
| 43 |
|
| 44 |
#if ( defined CPL_CHEM ) |
| 45 |
#if ( N_LAT == 24 ) |
| 46 |
ddepno (1,1:nlat) = (/ |
| 47 |
& -0.2300E-08, -0.2300E-08, -0.2558E-08, -0.3597E-08, |
| 48 |
& -0.3272E-08, -0.3715E-08, -0.3993E-08, -0.5103E-08, |
| 49 |
& -0.6727E-08, -0.7102E-08, -0.6618E-08, -0.7047E-08, |
| 50 |
& -0.6805E-08, -0.7037E-08, -0.8250E-08, -0.9527E-08, |
| 51 |
& -0.9802E-08, -0.1062E-07, -0.9616E-08, -0.7989E-08, |
| 52 |
& -0.5587E-08, -0.3989E-08, -0.2755E-08, -0.2300E-08 |
| 53 |
& /) |
| 54 |
ddepno2 (1,1:nlat) = (/ |
| 55 |
& -0.1150E-07, -0.1150E-07, -0.1408E-07, -0.2137E-07, |
| 56 |
& -0.2119E-07, -0.2463E-07, -0.2634E-07, -0.3318E-07, |
| 57 |
& -0.4317E-07, -0.4548E-07, -0.4249E-07, -0.4514E-07, |
| 58 |
& -0.4364E-07, -0.4507E-07, -0.5254E-07, -0.6040E-07, |
| 59 |
& -0.6209E-07, -0.6715E-07, -0.6025E-07, -0.4973E-07, |
| 60 |
& -0.3368E-07, -0.2381E-07, -0.1502E-07, -0.1150E-07 |
| 61 |
& /) |
| 62 |
ddephno3 (1,1:nlat) = (/ |
| 63 |
& -0.5750E-06, -0.5750E-06, -0.7042E-06, -0.1006E-05, |
| 64 |
& -0.1059E-05, -0.1211E-05, -0.1275E-05, -0.1532E-05, |
| 65 |
& -0.1906E-05, -0.1993E-05, -0.1881E-05, -0.1980E-05, |
| 66 |
& -0.1924E-05, -0.1978E-05, -0.2258E-05, -0.2552E-05, |
| 67 |
& -0.2616E-05, -0.2806E-05, -0.2524E-05, -0.2113E-05, |
| 68 |
& -0.1471E-05, -0.1099E-05, -0.7304E-06, -0.5750E-06 |
| 69 |
& /) |
| 70 |
ddepn2o5 (1,1:nlat) = (/ |
| 71 |
& -0.5750E-06, -0.5750E-06, -0.7042E-06, -0.1006E-05, |
| 72 |
& -0.1059E-05, -0.1211E-05, -0.1275E-05, -0.1532E-05, |
| 73 |
& -0.1906E-05, -0.1993E-05, -0.1881E-05, -0.1980E-05, |
| 74 |
& -0.1924E-05, -0.1978E-05, -0.2258E-05, -0.2552E-05, |
| 75 |
& -0.2616E-05, -0.2806E-05, -0.2524E-05, -0.2113E-05, |
| 76 |
& -0.1471E-05, -0.1099E-05, -0.7304E-06, -0.5750E-06 |
| 77 |
& /) |
| 78 |
ddepo3 (1,1:nlat) = (/ |
| 79 |
& -0.8050E-07, -0.8050E-07, -0.8050E-07, -0.1010E-06, |
| 80 |
& -0.8064E-07, -0.8722E-07, -0.9428E-07, -0.1225E-06, |
| 81 |
& -0.1637E-06, -0.1732E-06, -0.1609E-06, -0.1718E-06, |
| 82 |
& -0.1657E-06, -0.1715E-06, -0.2024E-06, -0.2348E-06, |
| 83 |
& -0.2417E-06, -0.2626E-06, -0.2416E-06, -0.2037E-06, |
| 84 |
& -0.1510E-06, -0.1107E-06, -0.8731E-07, -0.8050E-07 |
| 85 |
& /) |
| 86 |
ddeph2o2 (1,1:nlat) = (/ |
| 87 |
& -0.3680E-06, -0.3680E-06, -0.5437E-06, -0.6702E-06, |
| 88 |
& -0.1025E-05, -0.1140E-05, -0.1129E-05, -0.1086E-05, |
| 89 |
& -0.1024E-05, -0.1010E-05, -0.1028E-05, -0.1012E-05, |
| 90 |
& -0.1021E-05, -0.1012E-05, -0.9654E-06, -0.9163E-06, |
| 91 |
& -0.9057E-06, -0.8741E-06, -0.7831E-06, -0.7509E-06, |
| 92 |
& -0.6073E-06, -0.6612E-06, -0.4848E-06, -0.3680E-06 |
| 93 |
& /) |
| 94 |
ddepch3o2h (1,1:nlat) = (/ |
| 95 |
& -0.1150E-07, -0.1150E-07, -0.7351E-07, -0.1198E-06, |
| 96 |
& -0.2433E-06, -0.2844E-06, -0.2812E-06, -0.2684E-06, |
| 97 |
& -0.2497E-06, -0.2454E-06, -0.2509E-06, -0.2460E-06, |
| 98 |
& -0.2488E-06, -0.2461E-06, -0.2321E-06, -0.2174E-06, |
| 99 |
& -0.2142E-06, -0.2047E-06, -0.1709E-06, -0.1565E-06, |
| 100 |
& -0.1016E-06, -0.1174E-06, -0.5326E-07, -0.1150E-07 |
| 101 |
& /) |
| 102 |
ddepbc (1,1:nlat) = (/ |
| 103 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 104 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 105 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 106 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 107 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 108 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06 |
| 109 |
& /) |
| 110 |
ddepoc (1,1:nlat) = (/ |
| 111 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 112 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 113 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 114 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 115 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 116 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06 |
| 117 |
& /) |
| 118 |
ddepref (1,1:nlat) = (/ |
| 119 |
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, |
| 120 |
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, |
| 121 |
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 |
| 122 |
& /) |
| 123 |
#endif |
| 124 |
|
| 125 |
#if ( N_LAT == 46 ) |
| 126 |
ddepno (1,1:nlat) = (/ |
| 127 |
& -0.2300E-08, -0.2300E-08, -0.2306E-08, -0.2438E-08, |
| 128 |
& -0.2604E-08, -0.3135E-08, -0.3575E-08, -0.3409E-08, |
| 129 |
& -0.3311E-08, -0.3538E-08, -0.3746E-08, -0.3888E-08, |
| 130 |
& -0.4141E-08, -0.4708E-08, -0.5356E-08, -0.6186E-08, |
| 131 |
& -0.6794E-08, -0.6985E-08, -0.7005E-08, -0.6758E-08, |
| 132 |
& -0.6713E-08, -0.6933E-08, -0.6988E-08, -0.6864E-08, |
| 133 |
& -0.6867E-08, -0.6985E-08, -0.7387E-08, -0.8007E-08, |
| 134 |
& -0.8647E-08, -0.9300E-08, -0.9619E-08, -0.9759E-08, |
| 135 |
& -0.1009E-07, -0.1051E-07, -0.1024E-07, -0.9728E-08, |
| 136 |
& -0.8965E-08, -0.8134E-08, -0.6975E-08, -0.5747E-08, |
| 137 |
& -0.4877E-08, -0.4060E-08, -0.3413E-08, -0.2782E-08, |
| 138 |
& -0.2533E-08, -0.2300E-08 |
| 139 |
& /) |
| 140 |
|
| 141 |
ddepno2 (1,1:nlat) = (/ |
| 142 |
& -0.1150E-07, -0.1150E-07, -0.1156E-07, -0.1288E-07, |
| 143 |
& -0.1440E-07, -0.1813E-07, -0.2136E-07, -0.2127E-07, |
| 144 |
& -0.2150E-07, -0.2325E-07, -0.2482E-07, -0.2569E-07, |
| 145 |
& -0.2725E-07, -0.3075E-07, -0.3473E-07, -0.3984E-07, |
| 146 |
& -0.4358E-07, -0.4476E-07, -0.4488E-07, -0.4335E-07, |
| 147 |
& -0.4308E-07, -0.4443E-07, -0.4477E-07, -0.4401E-07, |
| 148 |
& -0.4402E-07, -0.4475E-07, -0.4723E-07, -0.5105E-07, |
| 149 |
& -0.5499E-07, -0.5900E-07, -0.6096E-07, -0.6183E-07, |
| 150 |
& -0.6389E-07, -0.6648E-07, -0.6454E-07, -0.6102E-07, |
| 151 |
& -0.5604E-07, -0.5067E-07, -0.4295E-07, -0.3475E-07, |
| 152 |
& -0.2929E-07, -0.2425E-07, -0.1971E-07, -0.1522E-07, |
| 153 |
& -0.1330E-07, -0.1150E-07 |
| 154 |
& /) |
| 155 |
|
| 156 |
ddephno3 (1,1:nlat) = (/ |
| 157 |
& -0.5750E-06, -0.5750E-06, -0.5779E-06, -0.6439E-06, |
| 158 |
& -0.7176E-06, -0.8719E-06, -0.1010E-05, -0.1037E-05, |
| 159 |
& -0.1073E-05, -0.1150E-05, -0.1218E-05, -0.1251E-05, |
| 160 |
& -0.1309E-05, -0.1441E-05, -0.1590E-05, -0.1781E-05, |
| 161 |
& -0.1921E-05, -0.1966E-05, -0.1971E-05, -0.1913E-05, |
| 162 |
& -0.1903E-05, -0.1954E-05, -0.1966E-05, -0.1938E-05, |
| 163 |
& -0.1938E-05, -0.1966E-05, -0.2059E-05, -0.2202E-05, |
| 164 |
& -0.2349E-05, -0.2500E-05, -0.2573E-05, -0.2606E-05, |
| 165 |
& -0.2684E-05, -0.2781E-05, -0.2699E-05, -0.2555E-05, |
| 166 |
& -0.2360E-05, -0.2150E-05, -0.1842E-05, -0.1514E-05, |
| 167 |
& -0.1306E-05, -0.1116E-05, -0.9270E-06, -0.7386E-06, |
| 168 |
& -0.6544E-06, -0.5750E-06 |
| 169 |
& /) |
| 170 |
|
| 171 |
ddepn2o5 (1,1:nlat) = (/ |
| 172 |
& -0.5750E-06, -0.5750E-06, -0.5779E-06, -0.6439E-06, |
| 173 |
& -0.7176E-06, -0.8719E-06, -0.1010E-05, -0.1037E-05, |
| 174 |
& -0.1073E-05, -0.1150E-05, -0.1218E-05, -0.1251E-05, |
| 175 |
& -0.1309E-05, -0.1441E-05, -0.1590E-05, -0.1781E-05, |
| 176 |
& -0.1921E-05, -0.1966E-05, -0.1971E-05, -0.1913E-05, |
| 177 |
& -0.1903E-05, -0.1954E-05, -0.1966E-05, -0.1938E-05, |
| 178 |
& -0.1938E-05, -0.1966E-05, -0.2059E-05, -0.2202E-05, |
| 179 |
& -0.2349E-05, -0.2500E-05, -0.2573E-05, -0.2606E-05, |
| 180 |
& -0.2684E-05, -0.2781E-05, -0.2699E-05, -0.2555E-05, |
| 181 |
& -0.2360E-05, -0.2150E-05, -0.1842E-05, -0.1514E-05, |
| 182 |
& -0.1306E-05, -0.1116E-05, -0.9270E-06, -0.7386E-06, |
| 183 |
& -0.6544E-06, -0.5750E-06 |
| 184 |
& /) |
| 185 |
|
| 186 |
ddepo3 (1,1:nlat) = (/ |
| 187 |
& -0.8050E-07, -0.8050E-07, -0.8050E-07, -0.8050E-07, |
| 188 |
& -0.8141E-07, -0.9189E-07, -0.9964E-07, -0.8924E-07, |
| 189 |
& -0.8122E-07, -0.8459E-07, -0.8800E-07, -0.9161E-07, |
| 190 |
& -0.9804E-07, -0.1125E-06, -0.1289E-06, -0.1500E-06, |
| 191 |
& -0.1654E-06, -0.1702E-06, -0.1707E-06, -0.1645E-06, |
| 192 |
& -0.1633E-06, -0.1689E-06, -0.1703E-06, -0.1672E-06, |
| 193 |
& -0.1672E-06, -0.1702E-06, -0.1804E-06, -0.1962E-06, |
| 194 |
& -0.2125E-06, -0.2290E-06, -0.2371E-06, -0.2406E-06, |
| 195 |
& -0.2491E-06, -0.2598E-06, -0.2547E-06, -0.2439E-06, |
| 196 |
& -0.2264E-06, -0.2071E-06, -0.1814E-06, -0.1545E-06, |
| 197 |
& -0.1331E-06, -0.1125E-06, -0.9978E-07, -0.8783E-07, |
| 198 |
& -0.8398E-07, -0.8050E-07 |
| 199 |
& /) |
| 200 |
|
| 201 |
ddeph2o2 (1,1:nlat) = (/ |
| 202 |
& -0.3680E-06, -0.3680E-06, -0.3719E-06, -0.4617E-06, |
| 203 |
& -0.5493E-06, -0.6140E-06, -0.6939E-06, -0.8752E-06, |
| 204 |
& -0.1035E-05, -0.1094E-05, -0.1139E-05, -0.1133E-05, |
| 205 |
& -0.1123E-05, -0.1101E-05, -0.1076E-05, -0.1045E-05, |
| 206 |
& -0.1022E-05, -0.1014E-05, -0.1014E-05, -0.1023E-05, |
| 207 |
& -0.1024E-05, -0.1016E-05, -0.1014E-05, -0.1019E-05, |
| 208 |
& -0.1019E-05, -0.1014E-05, -0.9985E-06, -0.9747E-06, |
| 209 |
& -0.9501E-06, -0.9250E-06, -0.9128E-06, -0.9073E-06, |
| 210 |
& -0.8945E-06, -0.8783E-06, -0.8397E-06, -0.7932E-06, |
| 211 |
& -0.7702E-06, -0.7538E-06, -0.6903E-06, -0.6169E-06, |
| 212 |
& -0.6313E-06, -0.6588E-06, -0.5789E-06, -0.4887E-06, |
| 213 |
& -0.4277E-06, -0.3680E-06 |
| 214 |
& /) |
| 215 |
|
| 216 |
ddepch3o2h (1,1:nlat) = (/ |
| 217 |
& -0.1150E-07, -0.1150E-07, -0.1288E-07, -0.4457E-07, |
| 218 |
& -0.7557E-07, -0.9923E-07, -0.1280E-06, -0.1912E-06, |
| 219 |
& -0.2470E-06, -0.2680E-06, -0.2840E-06, -0.2824E-06, |
| 220 |
& -0.2795E-06, -0.2730E-06, -0.2655E-06, -0.2559E-06, |
| 221 |
& -0.2489E-06, -0.2467E-06, -0.2465E-06, -0.2493E-06, |
| 222 |
& -0.2498E-06, -0.2473E-06, -0.2467E-06, -0.2481E-06, |
| 223 |
& -0.2481E-06, -0.2467E-06, -0.2421E-06, -0.2349E-06, |
| 224 |
& -0.2275E-06, -0.2200E-06, -0.2163E-06, -0.2147E-06, |
| 225 |
& -0.2108E-06, -0.2060E-06, -0.1919E-06, -0.1747E-06, |
| 226 |
& -0.1651E-06, -0.1578E-06, -0.1333E-06, -0.1053E-06, |
| 227 |
& -0.1086E-06, -0.1167E-06, -0.8747E-07, -0.5469E-07, |
| 228 |
& -0.3284E-07, -0.1150E-07 |
| 229 |
& /) |
| 230 |
|
| 231 |
ddepbc (1,1:nlat) = (/ |
| 232 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 233 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 234 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 235 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 236 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 237 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 238 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 239 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 240 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 241 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 242 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 243 |
& -0.1150E-06, -0.1150E-06 |
| 244 |
& /) |
| 245 |
ddepoc (1,1:nlat) = (/ |
| 246 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 247 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 248 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 249 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 250 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 251 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 252 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 253 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 254 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 255 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 256 |
& -0.1150E-06, -0.1150E-06, -0.1150E-06, -0.1150E-06, |
| 257 |
& -0.1150E-06, -0.1150E-06 |
| 258 |
& /) |
| 259 |
ddepref (1,1:nlat) = (/ |
| 260 |
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, |
| 261 |
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, |
| 262 |
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, |
| 263 |
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, |
| 264 |
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, |
| 265 |
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 |
| 266 |
& /) |
| 267 |
#endif |
| 268 |
#endif |
| 269 |
|
| 270 |
namelist /chem_name_list/bio_uptake, |
| 271 |
& svi_intensity, |
| 272 |
& nhr_for_chem |
| 273 |
|
| 274 |
! ----------------------------------------------------- |
| 275 |
|
| 276 |
#if ( defined CPL_CHEM ) |
| 277 |
|
| 278 |
c-------------------------------------------------- |
| 279 |
c Read control parameters from namelist: |
| 280 |
c |
| 281 |
c open(120,file='INCLUDE/chem.nml',status='old') |
| 282 |
C open statment moved to assign.inc |
| 283 |
read(120,nml=chem_name_list) |
| 284 |
close(120) |
| 285 |
|
| 286 |
c-------------------------------------------------- |
| 287 |
c Pressure at each vertical layer in hPa for |
| 288 |
c photochemistry: |
| 289 |
|
| 290 |
#if ( N_LEV == 9 ) |
| 291 |
airpress( 1) = 959.0 |
| 292 |
airpress( 2) = 894.0 |
| 293 |
airpress( 3) = 786.0 |
| 294 |
airpress( 4) = 634.0 |
| 295 |
airpress( 5) = 468.0 |
| 296 |
airpress( 6) = 321.0 |
| 297 |
airpress( 7) = 201.0 |
| 298 |
airpress( 8) = 103.0 |
| 299 |
airpress( 9) = 027.0 |
| 300 |
#endif |
| 301 |
|
| 302 |
#if ( N_LEV == 11 ) |
| 303 |
airpress( 1) = 959.0 |
| 304 |
airpress( 2) = 894.0 |
| 305 |
airpress( 3) = 786.0 |
| 306 |
airpress( 4) = 634.0 |
| 307 |
airpress( 5) = 468.0 |
| 308 |
airpress( 6) = 321.0 |
| 309 |
airpress( 7) = 201.0 |
| 310 |
airpress( 8) = 122.0 |
| 311 |
airpress( 9) = 78.0 |
| 312 |
airpress(10) = 42.0 |
| 313 |
airpress(11) = 17.0 |
| 314 |
#endif |
| 315 |
|
| 316 |
c-------------------------------------------------- |
| 317 |
c Read solar radiative flux spectrum data: |
| 318 |
c |
| 319 |
c read(163,101)specdata |
| 320 |
c101 format(6e11.3) |
| 321 |
|
| 322 |
c-------------------------------------------------- |
| 323 |
print *,'cheminit emissions for ',nchemyr,'years' |
| 324 |
c Read integrated daily emission data: |
| 325 |
c |
| 326 |
! read(164)edailyf11, |
| 327 |
! & edailyf12, |
| 328 |
! & edailyn2o, |
| 329 |
! & edailyco, |
| 330 |
! & edailynox, |
| 331 |
! & edailych4, |
| 332 |
! & edailyso2, |
| 333 |
! & edailyco2, |
| 334 |
! & edailyhfc134a, |
| 335 |
! & edailypfc, |
| 336 |
! & edailysf6, |
| 337 |
! & edailyuco, |
| 338 |
! & edailyunmv, |
| 339 |
! & edailyunox, |
| 340 |
! & edailyusox, |
| 341 |
! & edailybc, |
| 342 |
! & edailynh3, |
| 343 |
! & edailyoc, |
| 344 |
! & edailyubc, |
| 345 |
! & edailyunh3, |
| 346 |
! & edailyuoc, |
| 347 |
! & n_total_urban, |
| 348 |
! & n_urban |
| 349 |
|
| 350 |
! Changed by Chien on October 27 of 2008 |
| 351 |
! bc and oc emissions are scaled based on Bond |
| 352 |
! edailybc=edailybc*8.6/13.3 |
| 353 |
! edailyoc=edailyoc*33.4/54. |
| 354 |
|
| 355 |
c-------------------------------------------------- |
| 356 |
c Read rk table: |
| 357 |
c |
| 358 |
read(165)rktable1 |
| 359 |
read(165)rktable2 |
| 360 |
|
| 361 |
! |
| 362 |
! --- 121800 |
| 363 |
! --- moved from chememission.F |
| 364 |
! |
| 365 |
! Lightning produced NO emission data, 062795: |
| 366 |
! based on Kumar et al., 1995: |
| 367 |
! |
| 368 |
xnoxltnt = 1.072e+4 |
| 369 |
|
| 370 |
c 092295 |
| 371 |
c ===== For start a new run: |
| 372 |
c |
| 373 |
if(ISTRT1.eq.0)then |
| 374 |
c-------------------------------------------------- |
| 375 |
c Reading data in unit of ppb(m) = 1.e-9 kg/kg |
| 376 |
c from a given data set to initialize |
| 377 |
c simulation: |
| 378 |
|
| 379 |
read(163)cfc11 |
| 380 |
read(163)cfc12 |
| 381 |
read(163)xn2o |
| 382 |
read(163)o3 |
| 383 |
read(163)co |
| 384 |
read(163)zco2 |
| 385 |
read(163)ch4 |
| 386 |
read(163)bcarbon |
| 387 |
read(163)ocarbon |
| 388 |
|
| 389 |
c |
| 390 |
c 022796 |
| 391 |
c |
| 392 |
read(160)ho |
| 393 |
read(160)ho2 |
| 394 |
read(160)h2o2 |
| 395 |
read(160)xno |
| 396 |
read(160)xno2 |
| 397 |
read(160)ch2o |
| 398 |
read(160)so2 |
| 399 |
|
| 400 |
c |
| 401 |
c 051598 |
| 402 |
c |
| 403 |
read(162)o3top |
| 404 |
close(162) |
| 405 |
|
| 406 |
do 2 i=1,n3d |
| 407 |
cfc110(i,1,1) = cfc11(i,1,1) |
| 408 |
2 continue |
| 409 |
|
| 410 |
do 3 j=1,nlat |
| 411 |
do 3 i=1,nlon |
| 412 |
|
| 413 |
c 061695 add o3 top array: |
| 414 |
c o3top (i,j) = o3(i,j,nlev) |
| 415 |
|
| 416 |
do 4 k=1,nlev |
| 417 |
atomo (i,j,k) = 0.0 |
| 418 |
o1d (i,j,k) = 0.0 |
| 419 |
atomh (i,j,k) = 0.0 |
| 420 |
xno3 (i,j,k) = 0.0 |
| 421 |
xn2o5 (i,j,k) = 0.0 |
| 422 |
hno3 (i,j,k) = 0.0 |
| 423 |
ch3 (i,j,k) = 0.0 |
| 424 |
cho (i,j,k) = 0.0 |
| 425 |
ch3o (i,j,k) = 0.0 |
| 426 |
ch3o2 (i,j,k) = 0.0 |
| 427 |
ch3o2h(i,j,k) = 0.0 |
| 428 |
hoso2 (i,j,k) = 0.0 |
| 429 |
so3 (i,j,k) = 0.0 |
| 430 |
h2so4 (i,j,k) = 0.0 |
| 431 |
|
| 432 |
hfc134a (i,j,k) = 0.0 |
| 433 |
sf6 (i,j,k) = 0.0 |
| 434 |
|
| 435 |
! === if hfc, pfc, and sf6 are included: |
| 436 |
#if ( defined INC_3GASES ) |
| 437 |
pfc(i,j,k) = cfc11(i,j,k)/127.0/137.3675 |
| 438 |
& *55.0*88.0 !use f11 to derive CF4 |
| 439 |
!initial value (55pptv at 1977) |
| 440 |
|
| 441 |
!pfc(i,j,k) = xn2o(i,j,k)/279.0/44.0 * 40.0 * 88.0 / 1e3 |
| 442 |
! use n2o to derive CF4 initial value: 40 pptv at 1892 |
| 443 |
|
| 444 |
#else |
| 445 |
pfc(i,j,k) = 0.0 |
| 446 |
#endif |
| 447 |
|
| 448 |
4 continue |
| 449 |
3 continue |
| 450 |
|
| 451 |
#if ( defined INC_3GASES ) |
| 452 |
! 10/28/2008 |
| 453 |
! reading hfc, pfc, and sf6 from initfile for Dec 1990 |
| 454 |
print *, 'reading hfc, pfc, and sf6 from initfile ' |
| 455 |
read(161)hfc134a |
| 456 |
read(161)pfc |
| 457 |
read(161)sf6 |
| 458 |
close(161) |
| 459 |
#endif |
| 460 |
|
| 461 |
do i=1,n3d |
| 462 |
sviod(i,1,1) = 0.0 |
| 463 |
enddo |
| 464 |
|
| 465 |
c---------------------------------------------- |
| 466 |
c Calculate initial chemical field: |
| 467 |
c |
| 468 |
c first give a average flux to drive the photochemistry: |
| 469 |
c |
| 470 |
do i=1,n3d |
| 471 |
solarflux(i,1,1) = 600.0 |
| 472 |
enddo |
| 473 |
|
| 474 |
c call ss and solver to derive the initial values: |
| 475 |
|
| 476 |
c call chemtrop0(1, pT,qv,3600.0,1) |
| 477 |
|
| 478 |
do i=1,n3d |
| 479 |
solarflux(i,1,1) = 0.0 |
| 480 |
enddo |
| 481 |
|
| 482 |
endif |
| 483 |
|
| 484 |
c 092295 |
| 485 |
c ===== For renew a run: |
| 486 |
c |
| 487 |
if(ISTRT1.eq.1)then |
| 488 |
|
| 489 |
c open(179,file='DUMP/renewchem.dat', |
| 490 |
c & form='unformatted', |
| 491 |
c & status='unknown') |
| 492 |
C open statment moved to assign.inc |
| 493 |
C file number changed to 279 |
| 494 |
|
| 495 |
read(279)iyearchem,mymonth,airmass, |
| 496 |
& cfc11 ,cfc110, |
| 497 |
& cfc11m , |
| 498 |
& cfc11sd, |
| 499 |
& cfc12 ,cfc12m, |
| 500 |
& cfc12sd, |
| 501 |
& xn2o ,xn2om , |
| 502 |
& xn2osd , |
| 503 |
& hfc134a,hfc134am, |
| 504 |
& pfc ,pfcm , |
| 505 |
& sf6 ,sf6m , |
| 506 |
& bcarbon,bcm , |
| 507 |
& ocarbon,ocm , |
| 508 |
& atomo , |
| 509 |
& o1d , |
| 510 |
& o3 ,o3m , |
| 511 |
& co ,com , |
| 512 |
& zco2 ,zco2m, |
| 513 |
& atomh , |
| 514 |
& ho , |
| 515 |
& ho2 ,hoxm , |
| 516 |
& h2o2 , |
| 517 |
& xno , |
| 518 |
& xno2 ,xnoxm, |
| 519 |
& xno3 , |
| 520 |
& xn2o5 ,xnoym, |
| 521 |
& hno3 , |
| 522 |
& ch4 ,ch4m , |
| 523 |
& ch3 , |
| 524 |
& cho , |
| 525 |
& ch2o , |
| 526 |
& ch3o , |
| 527 |
& ch3o2 , |
| 528 |
& ch3o2h , |
| 529 |
& so2 ,so2m , |
| 530 |
& hoso2 , |
| 531 |
& so3 , |
| 532 |
& h2so4 ,h2so4m, |
| 533 |
& sviod ,sviodm |
| 534 |
|
| 535 |
close(279) |
| 536 |
|
| 537 |
c mymonth = mymonth + 1 |
| 538 |
c if(mymonth.gt.12)then |
| 539 |
c myyear = myyear + 1 |
| 540 |
c mymonth= 1 |
| 541 |
c endif |
| 542 |
|
| 543 |
endif |
| 544 |
|
| 545 |
c---------------------------------------------- |
| 546 |
c Set monthstep = 0 and start monthly average |
| 547 |
c calculation: |
| 548 |
c |
| 549 |
monthstep = 0 |
| 550 |
do i=1,n3d |
| 551 |
cfc11m (i,1,1) = 0.0 |
| 552 |
cfc12m (i,1,1) = 0.0 |
| 553 |
xn2om (i,1,1) = 0.0 |
| 554 |
o3m (i,1,1) = 0.0 |
| 555 |
com (i,1,1) = 0.0 |
| 556 |
zco2m (i,1,1) = 0.0 |
| 557 |
hoxm (i,1,1) = 0.0 |
| 558 |
xnoxm (i,1,1) = 0.0 |
| 559 |
xnoym (i,1,1) = 0.0 |
| 560 |
ch4m (i,1,1) = 0.0 |
| 561 |
so2m (i,1,1) = 0.0 |
| 562 |
h2so4m (i,1,1) = 0.0 |
| 563 |
sviodm (i,1,1) = 0.0 |
| 564 |
|
| 565 |
hfc134am(i,1,1) = 0.0 |
| 566 |
pfcm (i,1,1) = 0.0 |
| 567 |
sf6m (i,1,1) = 0.0 |
| 568 |
|
| 569 |
bcm (i,1,1) = 0.0 |
| 570 |
ocm (i,1,1) = 0.0 |
| 571 |
bcodm (i,1,1) = 0.0 |
| 572 |
ocodm (i,1,1) = 0.0 |
| 573 |
enddo |
| 574 |
|
| 575 |
c------------------------------------------ |
| 576 |
c Write initial field into files: |
| 577 |
c |
| 578 |
! write(6,199)cfc11 |
| 579 |
! write(6,199)cfc12 |
| 580 |
! write(6,199)xn2o |
| 581 |
! write(6,199)o3 |
| 582 |
! write(6,199)co |
| 583 |
! write(6,199)zco2 |
| 584 |
! write(6,199)ch4 |
| 585 |
199 format(6e11.3) |
| 586 |
|
| 587 |
! write(166)cfc11 |
| 588 |
! write(166)cfc12 |
| 589 |
! write(166)xn2o |
| 590 |
|
| 591 |
! write(171)o3 |
| 592 |
! write(171)o1d |
| 593 |
! write(171)ho |
| 594 |
! write(171)ho2 |
| 595 |
! write(171)h2o2 |
| 596 |
|
| 597 |
! write(172)xno |
| 598 |
! write(172)xno2 |
| 599 |
! write(172)xno3 |
| 600 |
! write(172)xn2o5 |
| 601 |
! write(172)hno3 |
| 602 |
|
| 603 |
! write(173)co |
| 604 |
! write(173)zco2 |
| 605 |
! write(173)ch4 |
| 606 |
! write(173)ch2o |
| 607 |
|
| 608 |
! write(174)so2 |
| 609 |
! write(174)hoso2 |
| 610 |
! write(174)so3 |
| 611 |
! write(174)h2so4 |
| 612 |
|
| 613 |
#endif |
| 614 |
|
| 615 |
return |
| 616 |
end |
| 617 |
|
Parent Directory
| 
Revision Log
| 
Revision Graph