1 |
jscott |
1.1 |
|
2 |
|
|
! ============================================================ |
3 |
|
|
! |
4 |
|
|
! CHEM_META: An include file consists of COMMONs for the |
5 |
|
|
! META sub model. |
6 |
|
|
! |
7 |
|
|
! ------------------------------------------------------------ |
8 |
|
|
! Author: Chien Wang |
9 |
|
|
! MIT Joint Program on Science and Policy |
10 |
|
|
! of Global Change |
11 |
|
|
! Created: January 29, 1999 |
12 |
|
|
! |
13 |
|
|
! Revision: |
14 |
|
|
! 062599 Chien Wang using timec data |
15 |
|
|
! 112800 Chien Wang use meta_nvar |
16 |
|
|
! |
17 |
|
|
! ============================================================ |
18 |
|
|
|
19 |
|
|
c---------------------------------------------------------- |
20 |
|
|
c Brief descriptions of variables: |
21 |
|
|
! -------------------------------------------- |
22 |
|
|
! |
23 |
|
|
! --- functions of latitudes: |
24 |
|
|
! rlati: |lat| in degree (0:65) |
25 |
|
|
! rtemp: deviation from normal temperature (-10:10) |
26 |
|
|
! rclou: tenth of mean cloud coverage (0:10) |
27 |
|
|
! rmixi: maximum mixing layer top in meter (200:2000) |
28 |
|
|
! rtime: residence time 200km/U in second (21600:259200) |
29 |
|
|
! |
30 |
|
|
! --- functions of latitudes and pollution types: |
31 |
|
|
! rso2: emission of SO2 in kg/km^2/day |
32 |
|
|
! rco: emission of SO2 in kg/km^2/day |
33 |
|
|
! rnox: emission of SO2 in kg/km^2/day |
34 |
|
|
! rvoc: deviation from total daily VOC emission (-0.7:0.7) |
35 |
|
|
! raqino |
36 |
|
|
! raqinoo |
37 |
|
|
! raqivoc |
38 |
|
|
! raqiozo |
39 |
|
|
! raqisox |
40 |
|
|
! results(58) |
41 |
|
|
! |
42 |
|
|
! --- Order of results and others: |
43 |
|
|
! --- |
44 |
|
|
! --- 1: NO mass |
45 |
|
|
! --- 2: NO flux |
46 |
|
|
! --- 3: NO2 mass |
47 |
|
|
! --- 4: NO2 flux |
48 |
|
|
! --- 5: O3 mass |
49 |
|
|
! --- 6: O3 flux |
50 |
|
|
! --- 7: HONO flux |
51 |
|
|
! --- 8: HNO3 flux |
52 |
|
|
! --- 9: N2O5 flux |
53 |
|
|
! --- 10: CO mass |
54 |
|
|
! --- 11: CO flux |
55 |
|
|
! --- 12: HCHO mass |
56 |
|
|
! --- 13: HCHO flux |
57 |
|
|
! --- 14: ALD2 mass |
58 |
|
|
! --- 15: MEK mass |
59 |
|
|
! --- 16: PAN flux |
60 |
|
|
! --- 17: ALKA mass |
61 |
|
|
! --- 18: ETHE mass |
62 |
|
|
! --- 19: ALKE mass |
63 |
|
|
! --- 20: TOLU mass |
64 |
|
|
! --- 21: AROM mass |
65 |
|
|
! --- 22: SO2 mass |
66 |
|
|
! --- 23: SO2 flux |
67 |
|
|
! --- 24: SO3 mass |
68 |
|
|
! --- 25: SO3 flux |
69 |
|
|
! --- 26: O3 max |
70 |
|
|
! --- 27: O3 peak hour |
71 |
|
|
! --- Daily averaged mass |
72 |
|
|
! --- 28: NO |
73 |
|
|
! --- 29: NO2 |
74 |
|
|
! --- 30: O3 |
75 |
|
|
! --- 31: HONO |
76 |
|
|
! --- 32: HNO3 |
77 |
|
|
! --- 33: HNO4 |
78 |
|
|
! --- 34: N2O5 |
79 |
|
|
! --- 35: NO3 |
80 |
|
|
! --- 36: HO2 |
81 |
|
|
! --- 37: CO |
82 |
|
|
! --- 38: HCHO |
83 |
|
|
! --- 39: ALD2 |
84 |
|
|
! --- 40: MEK |
85 |
|
|
! --- 41: MGLY |
86 |
|
|
! --- 42: PAN |
87 |
|
|
! --- 43: RO2 |
88 |
|
|
! --- 44: MCO3 |
89 |
|
|
! --- 45: ALKN |
90 |
|
|
! --- 46: ALKA |
91 |
|
|
! --- 47: ETHE |
92 |
|
|
! --- 48: ALKE |
93 |
|
|
! --- 49: TOLU |
94 |
|
|
! --- 50: AROM |
95 |
|
|
! --- 51: DIAL |
96 |
|
|
! --- 52: CRES |
97 |
|
|
! --- 53: NPHE |
98 |
|
|
! --- 54: H2O2 |
99 |
|
|
! --- 55: MEOH |
100 |
|
|
! --- 56: ETOH |
101 |
|
|
! --- 57: SO2 |
102 |
|
|
! --- 58: SO3 |
103 |
|
|
! ---------------------------------------------------- |
104 |
|
|
|
105 |
|
|
real rlati (nlat), rtemp (nlat), rclou (nlat), |
106 |
|
|
& rmixi (nlat), rtime (nlat) |
107 |
|
|
|
108 |
|
|
real raqino (3,nlat), raqinoo(3,nlat), |
109 |
|
|
& raqivoc(3,nlat), raqiozo(3,nlat), raqisox(3,nlat) |
110 |
|
|
|
111 |
|
|
real results_meta(meta_nvar,3,nlat), |
112 |
|
|
& results_mon(meta_nvar,3,nlat) |
113 |
|
|
|
114 |
|
|
real xc_meta(meta_nvar),w_meta(meta_nvar), |
115 |
|
|
& a_meta(meta_nvar),c_meta(meta_nvar) |
116 |
|
|
|
117 |
|
|
real alpha_nox(nlat,nchemyr),alpha_so2 (nlat), |
118 |
|
|
& alpha_co (nlat,nchemyr),alpha_vocs(nlat,nchemyr) |
119 |
|
|
|
120 |
|
|
real urban_beta_1(3), |
121 |
|
|
& urban_beta_2(3), |
122 |
|
|
& urban_beta_3(3) |
123 |
|
|
|
124 |
|
|
integer n_total_urban (nlat,nchemyr), |
125 |
|
|
& n_urban (3,nlat,nchemyr) |
126 |
|
|
|
127 |
|
|
integer nstep_meta |
128 |
|
|
|
129 |
|
|
common /chem_meta_comm/rlati, |
130 |
|
|
& rtemp, |
131 |
|
|
& rclou, |
132 |
|
|
& rmixi, |
133 |
|
|
& rtime, |
134 |
|
|
& raqino, |
135 |
|
|
& raqinoo, |
136 |
|
|
& raqivoc, |
137 |
|
|
& raqiozo, |
138 |
|
|
& raqisox, |
139 |
|
|
& results_meta, |
140 |
|
|
& results_mon |
141 |
|
|
|
142 |
|
|
! ---------------------------------------------- |
143 |
|
|
! Misc coefficients: |
144 |
|
|
! |
145 |
|
|
common /chem_meta_para0/nstep_meta |
146 |
|
|
|
147 |
|
|
! ---------------------------------------------------- |
148 |
|
|
! Monika's sin fitting coefficients: |
149 |
|
|
! |
150 |
|
|
common /chem_meta_para1/xc_meta, |
151 |
|
|
& w_meta, |
152 |
|
|
& a_meta, |
153 |
|
|
& c_meta |
154 |
|
|
|
155 |
|
|
! ---------------------------------------------------- |
156 |
|
|
! Urban area numbers of low, midium, and high pollution |
157 |
|
|
! |
158 |
|
|
common /chem_meta_para2/n_total_urban, |
159 |
|
|
& n_urban |
160 |
|
|
|
161 |
|
|
! ---------------------------------------------------- |
162 |
|
|
! Fraction of urban shares in total emissions (0-1) |
163 |
|
|
! |
164 |
|
|
common /chem_meta_para3/alpha_nox,alpha_so2, |
165 |
|
|
& alpha_co,alpha_vocs, |
166 |
|
|
& urban_beta_1, |
167 |
|
|
& urban_beta_2, |
168 |
|
|
& urban_beta_3 |
169 |
|
|
|