Line data Source code
1 :
2 : module mo_sim_dat
3 :
4 : private
5 : public :: set_sim_dat
6 :
7 : contains
8 :
9 1536 : subroutine set_sim_dat
10 :
11 : use chem_mods, only : clscnt, cls_rxt_cnt, clsmap, permute, adv_mass, fix_mass, crb_mass
12 : use chem_mods, only : diag_map
13 : use chem_mods, only : phtcnt, rxt_tag_cnt, rxt_tag_lst, rxt_tag_map
14 : use chem_mods, only : pht_alias_lst, pht_alias_mult
15 : use chem_mods, only : extfrc_lst, inv_lst, slvd_lst
16 : use chem_mods, only : enthalpy_cnt, cph_enthalpy, cph_rid, num_rnts, rxntot
17 : use cam_abortutils,only : endrun
18 : use mo_tracname, only : solsym
19 : use chem_mods, only : frc_from_dataset
20 : use chem_mods, only : is_scalar, is_vector
21 : use shr_kind_mod, only : r8 => shr_kind_r8
22 : use cam_logfile, only : iulog
23 :
24 : implicit none
25 :
26 : !--------------------------------------------------------------
27 : ! ... local variables
28 : !--------------------------------------------------------------
29 : integer :: ios
30 :
31 1536 : is_scalar = .false.
32 1536 : is_vector = .true.
33 :
34 1536 : clscnt(:) = (/ 2, 0, 0, 317, 0 /)
35 :
36 7680 : cls_rxt_cnt(:,1) = (/ 6, 0, 0, 2 /)
37 7680 : cls_rxt_cnt(:,4) = (/ 2, 255, 647, 317 /)
38 :
39 : solsym(:319) = (/ 'ALKNIT ','ALKOOH ','AOA_NH ','APIN ','bc_a1 ', &
40 : 'bc_a4 ','BCARY ','BENZENE ','BENZOOH ','BEPOMUC ', &
41 : 'BIGALD1 ','BIGALD2 ','BIGALD3 ','BIGALD4 ','BIGALK ', &
42 : 'BIGENE ','BPIN ','BR ','BRCL ','BRO ', &
43 : 'BRONO2 ','BRY ','BZALD ','BZOOH ','C2H2 ', &
44 : 'C2H4 ','C2H5OH ','C2H5OOH ','C2H6 ','C3H6 ', &
45 : 'C3H7OOH ','C3H8 ','C6H5OOH ','CCL4 ','CF2CLBR ', &
46 : 'CF3BR ','CFC11 ','CFC113 ','CFC114 ','CFC115 ', &
47 : 'CFC12 ','CH2BR2 ','CH2O ','CH3BR ','CH3CCL3 ', &
48 : 'CH3CHO ','CH3CL ','CH3CN ','CH3COCH3 ','CH3COCHO ', &
49 : 'CH3COOH ','CH3COOOH ','CH3OH ','CH3OOH ','CH4 ', &
50 : 'CHBR3 ','CL ','CL2 ','CL2O2 ','CLO ', &
51 : 'CLONO2 ','CLY ','CO ','CO2 ','COF2 ', &
52 : 'COFCL ','CRESOL ','DHPMPAL ','DMS ','dst_a1 ', &
53 : 'dst_a2 ','dst_a3 ','E90 ','EOOH ','F ', &
54 : 'GLYALD ','GLYOXAL ','H ','H2 ','H2402 ', &
55 : 'H2O2 ','H2SO4 ','HBR ','HCFC141B ','HCFC142B ', &
56 : 'HCFC22 ','HCL ','HCN ','HCOCH2OOH ','HCOOH ', &
57 : 'HF ','HMHP ','HNO3 ','HO2NO2 ','HOBR ', &
58 : 'HOCL ','HONITR ','HPALD1 ','HPALD4 ','HPALDB1C ', &
59 : 'HPALDB4C ','HYAC ','HYDRALD ','HYPERACET ','ICHE ', &
60 : 'IEPOX ','INHEB ','INHED ','ISOP ','ISOPFDN ', &
61 : 'ISOPFDNC ','ISOPFNC ','ISOPFNP ','ISOPHFP ','ISOPN1D ', &
62 : 'ISOPN2B ','ISOPN3B ','ISOPN4D ','ISOPNBNO3 ','ISOPNOOHB ', &
63 : 'ISOPNOOHD ','ISOPOH ','ISOPOOH ','IVOC ','LIMON ', &
64 : 'MACR ','MACRN ','MACROOH ','MEK ','MEKOOH ', &
65 : 'MPAN ','MVK ','MVKN ','MVKOOH ','MYRC ', &
66 : 'N ','N2O ','N2O5 ','NC4CHO ','ncl_a1 ', &
67 : 'ncl_a2 ','ncl_a3 ','NH3 ','NH4 ','NH_5 ', &
68 : 'NH_50 ','NO ','NO2 ','NO3 ','NO3CH2CHO ', &
69 : 'NOA ','num_a1 ','num_a2 ','num_a3 ','num_a4 ', &
70 : 'num_a5 ','O ','O3 ','O3S ','OCLO ', &
71 : 'OCS ','ONITR ','PAN ','PBZNIT ','PHENO ', &
72 : 'PHENOL ','PHENOOH ','pom_a1 ','pom_a4 ','POOH ', &
73 : 'ROOH ','S ','SF6 ','SO ','SO2 ', &
74 : 'SO3 ','so4_a1 ','so4_a2 ','so4_a3 ','so4_a5 ', &
75 : 'soa1_a1 ','soa1_a2 ','soa2_a1 ','soa2_a2 ','soa3_a1 ', &
76 : 'soa3_a2 ','soa4_a1 ','soa4_a2 ','soa5_a1 ','soa5_a2 ', &
77 : 'SOAG0 ','SOAG1 ','SOAG2 ','SOAG3 ','SOAG4 ', &
78 : 'SQTN ','ST80_25 ','SVOC ','TEPOMUC ','TERP1OOH ', &
79 : 'TERP2AOOH ','TERPA ','TERPA2 ','TERPA2PAN ','TERPA3 ', &
80 : 'TERPA3PAN ','TERPACID ','TERPACID2 ','TERPACID3 ','TERPAPAN ', &
81 : 'TERPDHDP ','TERPF1 ','TERPF2 ','TERPFDN ','TERPHFN ', &
82 : 'TERPK ','TERPNPS ','TERPNPS1 ','TERPNPT ','TERPNPT1 ', &
83 : 'TERPNS ','TERPNS1 ','TERPNT ','TERPNT1 ','TERPOOH ', &
84 : 'TERPOOHL ','TOLOOH ','TOLUENE ','XYLENES ','XYLENOOH ', &
85 : 'XYLOL ','XYLOLOOH ','NHDEP ','NDEP ','ACBZO2 ', &
86 : 'ALKO2 ','APINNO3 ','APINO2 ','APINO2VBS ','BCARYNO3 ', &
87 : 'BCARYO2 ','BCARYO2VBS ','BENZO2 ','BENZO2VBS ','BPINNO3 ', &
88 : 'BPINO2 ','BPINO2VBS ','BZOO ','C2H5O2 ','C3H7O2 ', &
89 : 'C6H5O2 ','CH3CO3 ','CH3O2 ','DICARBO2 ','ENEO2 ', &
90 : 'EO ','EO2 ','HO2 ','HOCH2OO ','IEPOXOO ', &
91 : 'ISOPB1O2 ','ISOPB4O2 ','ISOPC1C ','ISOPC1T ','ISOPC4C ', &
92 : 'ISOPC4T ','ISOPED1O2 ','ISOPED4O2 ','ISOPN1DO2 ','ISOPN2BO2 ', &
93 : 'ISOPN3BO2 ','ISOPN4DO2 ','ISOPNBNO3O2 ','ISOPNO3 ','ISOPNOOHBO2 ', &
94 : 'ISOPNOOHDO2 ','ISOPO2VBS ','ISOPZD1O2 ','ISOPZD4O2 ','IVOCO2VBS ', &
95 : 'LIMONNO3 ','LIMONO2 ','LIMONO2VBS ','MACRO2 ','MALO2 ', &
96 : 'MCO3 ','MDIALO2 ','MEKO2 ','MVKO2 ','MYRCNO3 ', &
97 : 'MYRCO2 ','MYRCO2VBS ','NC4CHOO2 ','O1D ','OH ', &
98 : 'PHENO2 ','PO2 ','RO2 ','TERP1OOHO2 ','TERP2OOHO2 ', &
99 : 'TERPA1O2 ','TERPA2CO3 ','TERPA2O2 ','TERPA3CO3 ','TERPA3O2 ', &
100 : 'TERPA4O2 ','TERPACO3 ','TERPF1O2 ','TERPF2O2 ','TERPNPS1O2 ', &
101 : 'TERPNPT1O2 ','TERPNS1O2 ','TERPNT1O2 ','TOLO2 ','TOLUO2VBS ', &
102 491520 : 'XYLENO2 ','XYLEO2VBS ','XYLOLO2 ','H2O ' /)
103 :
104 : adv_mass(:319) = (/ 133.141340_r8, 104.142600_r8, 28.010400_r8, 136.228400_r8, 12.011000_r8, &
105 : 12.011000_r8, 204.342600_r8, 78.110400_r8, 160.122200_r8, 126.108600_r8, &
106 : 84.072400_r8, 98.098200_r8, 98.098200_r8, 112.124000_r8, 72.143800_r8, &
107 : 56.103200_r8, 136.228400_r8, 79.904000_r8, 115.356700_r8, 95.903400_r8, &
108 : 141.908940_r8, 99.716850_r8, 106.120800_r8, 124.135000_r8, 26.036800_r8, &
109 : 28.051600_r8, 46.065800_r8, 62.065200_r8, 30.066400_r8, 42.077400_r8, &
110 : 76.091000_r8, 44.092200_r8, 110.109200_r8, 153.821800_r8, 165.364506_r8, &
111 : 148.910210_r8, 137.367503_r8, 187.375310_r8, 170.921013_r8, 154.466716_r8, &
112 : 120.913206_r8, 173.833800_r8, 30.025200_r8, 94.937200_r8, 133.402300_r8, &
113 : 44.051000_r8, 50.485900_r8, 41.050940_r8, 58.076800_r8, 72.061400_r8, &
114 : 60.050400_r8, 76.049800_r8, 32.040000_r8, 48.039400_r8, 16.040600_r8, &
115 : 252.730400_r8, 35.452700_r8, 70.905400_r8, 102.904200_r8, 51.452100_r8, &
116 : 97.457640_r8, 100.916850_r8, 28.010400_r8, 44.009800_r8, 66.007206_r8, &
117 : 82.461503_r8, 108.135600_r8, 136.100200_r8, 62.132400_r8, 135.064039_r8, &
118 : 135.064039_r8, 135.064039_r8, 28.010400_r8, 78.064600_r8, 18.998403_r8, &
119 : 60.050400_r8, 58.035600_r8, 1.007400_r8, 2.014800_r8, 259.823613_r8, &
120 : 34.013600_r8, 98.078400_r8, 80.911400_r8, 116.948003_r8, 100.493706_r8, &
121 : 86.467906_r8, 36.460100_r8, 27.025140_r8, 76.049800_r8, 46.024600_r8, &
122 : 20.005803_r8, 64.038800_r8, 63.012340_r8, 79.011740_r8, 96.910800_r8, &
123 : 52.459500_r8, 135.114940_r8, 116.112400_r8, 116.112400_r8, 116.112400_r8, &
124 : 116.112400_r8, 74.076200_r8, 100.113000_r8, 90.075600_r8, 116.112400_r8, &
125 : 118.127200_r8, 163.125340_r8, 163.125340_r8, 68.114200_r8, 226.137680_r8, &
126 : 224.122880_r8, 195.124140_r8, 197.138940_r8, 150.126000_r8, 147.125940_r8, &
127 : 147.125940_r8, 147.125940_r8, 147.125940_r8, 147.125940_r8, 163.125340_r8, &
128 : 163.125340_r8, 102.127800_r8, 118.127200_r8, 184.350200_r8, 136.228400_r8, &
129 : 70.087800_r8, 149.099540_r8, 120.100800_r8, 72.102600_r8, 104.101400_r8, &
130 : 147.084740_r8, 70.087800_r8, 149.099540_r8, 120.100800_r8, 136.228400_r8, &
131 : 14.006740_r8, 44.012880_r8, 108.010480_r8, 145.111140_r8, 58.442468_r8, &
132 : 58.442468_r8, 58.442468_r8, 17.028940_r8, 18.036340_r8, 28.010400_r8, &
133 : 28.010400_r8, 30.006140_r8, 46.005540_r8, 62.004940_r8, 105.048540_r8, &
134 : 119.074340_r8, 1.007400_r8, 1.007400_r8, 1.007400_r8, 1.007400_r8, &
135 : 1.007400_r8, 15.999400_r8, 47.998200_r8, 47.998200_r8, 67.451500_r8, &
136 : 60.076400_r8, 133.100140_r8, 121.047940_r8, 183.117740_r8, 93.102400_r8, &
137 : 94.109800_r8, 176.121600_r8, 12.011000_r8, 12.011000_r8, 92.090400_r8, &
138 : 90.075600_r8, 32.066000_r8, 146.056419_r8, 48.065400_r8, 64.064800_r8, &
139 : 80.064200_r8, 115.107340_r8, 115.107340_r8, 115.107340_r8, 115.107340_r8, &
140 : 250.445000_r8, 250.445000_r8, 250.445000_r8, 250.445000_r8, 250.445000_r8, &
141 : 250.445000_r8, 250.445000_r8, 250.445000_r8, 250.445000_r8, 250.445000_r8, &
142 : 250.445000_r8, 250.445000_r8, 250.445000_r8, 250.445000_r8, 250.445000_r8, &
143 : 283.354340_r8, 28.010400_r8, 310.582400_r8, 140.134400_r8, 186.241400_r8, &
144 : 186.241400_r8, 168.227200_r8, 154.201400_r8, 231.198340_r8, 170.200800_r8, &
145 : 247.197740_r8, 200.226000_r8, 186.200200_r8, 202.199600_r8, 245.224140_r8, &
146 : 236.254400_r8, 168.227200_r8, 110.150400_r8, 294.251880_r8, 265.253140_r8, &
147 : 138.202000_r8, 231.239540_r8, 231.239540_r8, 231.239540_r8, 231.239540_r8, &
148 : 215.240140_r8, 215.240140_r8, 215.240140_r8, 215.240140_r8, 186.241400_r8, &
149 : 218.240200_r8, 174.148000_r8, 92.136200_r8, 106.162000_r8, 188.173800_r8, &
150 : 122.161400_r8, 204.173200_r8, 14.006740_r8, 14.006740_r8, 137.112200_r8, &
151 : 103.135200_r8, 230.232140_r8, 185.234000_r8, 185.234000_r8, 298.346340_r8, &
152 : 253.348200_r8, 253.348200_r8, 159.114800_r8, 159.114800_r8, 230.232140_r8, &
153 : 185.234000_r8, 185.234000_r8, 123.127600_r8, 61.057800_r8, 75.083600_r8, &
154 : 109.101800_r8, 75.042400_r8, 47.032000_r8, 129.089600_r8, 105.108800_r8, &
155 : 61.057800_r8, 77.057200_r8, 33.006200_r8, 63.031400_r8, 149.118600_r8, &
156 : 117.119800_r8, 117.119800_r8, 85.121000_r8, 85.121000_r8, 85.121000_r8, &
157 : 85.121000_r8, 117.119800_r8, 117.119800_r8, 196.131540_r8, 196.131540_r8, &
158 : 196.131540_r8, 196.131540_r8, 196.131540_r8, 162.117940_r8, 212.130940_r8, &
159 : 212.130940_r8, 117.119800_r8, 117.119800_r8, 117.119800_r8, 233.355800_r8, &
160 : 230.232140_r8, 185.234000_r8, 185.234000_r8, 119.093400_r8, 115.063800_r8, &
161 : 101.079200_r8, 117.078600_r8, 103.094000_r8, 119.093400_r8, 230.232140_r8, &
162 : 185.234000_r8, 185.234000_r8, 194.116740_r8, 15.999400_r8, 17.006800_r8, &
163 : 175.114200_r8, 91.083000_r8, 89.068200_r8, 235.247000_r8, 235.247000_r8, &
164 : 171.208200_r8, 185.192800_r8, 187.207600_r8, 201.192200_r8, 203.207000_r8, &
165 : 161.129600_r8, 199.218600_r8, 217.232800_r8, 159.156000_r8, 280.245140_r8, &
166 : 280.245140_r8, 264.245740_r8, 264.245740_r8, 173.140600_r8, 173.140600_r8, &
167 1536 : 187.166400_r8, 187.166400_r8, 203.165800_r8, 18.014200_r8 /)
168 :
169 : crb_mass(:319) = (/ 60.055000_r8, 60.055000_r8, 12.011000_r8, 120.110000_r8, 12.011000_r8, &
170 : 12.011000_r8, 180.165000_r8, 72.066000_r8, 72.066000_r8, 72.066000_r8, &
171 : 48.044000_r8, 60.055000_r8, 60.055000_r8, 72.066000_r8, 60.055000_r8, &
172 : 48.044000_r8, 120.110000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, &
173 : 0.000000_r8, 0.000000_r8, 84.077000_r8, 84.077000_r8, 24.022000_r8, &
174 : 24.022000_r8, 24.022000_r8, 24.022000_r8, 24.022000_r8, 36.033000_r8, &
175 : 36.033000_r8, 36.033000_r8, 72.066000_r8, 12.011000_r8, 12.011000_r8, &
176 : 12.011000_r8, 12.011000_r8, 24.022000_r8, 24.022000_r8, 24.022000_r8, &
177 : 12.011000_r8, 12.011000_r8, 12.011000_r8, 12.011000_r8, 24.022000_r8, &
178 : 24.022000_r8, 12.011000_r8, 24.022000_r8, 36.033000_r8, 36.033000_r8, &
179 : 24.022000_r8, 24.022000_r8, 12.011000_r8, 12.011000_r8, 12.011000_r8, &
180 : 12.011000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, &
181 : 0.000000_r8, 12.011000_r8, 12.011000_r8, 12.011000_r8, 12.011000_r8, &
182 : 12.011000_r8, 84.077000_r8, 48.044000_r8, 24.022000_r8, 0.000000_r8, &
183 : 0.000000_r8, 0.000000_r8, 12.011000_r8, 24.022000_r8, 0.000000_r8, &
184 : 24.022000_r8, 24.022000_r8, 0.000000_r8, 0.000000_r8, 24.022000_r8, &
185 : 0.000000_r8, 0.000000_r8, 0.000000_r8, 24.022000_r8, 24.022000_r8, &
186 : 12.011000_r8, 0.000000_r8, 12.011000_r8, 24.022000_r8, 12.011000_r8, &
187 : 0.000000_r8, 12.011000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, &
188 : 0.000000_r8, 48.044000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, &
189 : 60.055000_r8, 36.033000_r8, 60.055000_r8, 36.033000_r8, 60.055000_r8, &
190 : 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, &
191 : 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, &
192 : 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, &
193 : 60.055000_r8, 60.055000_r8, 60.055000_r8, 156.143000_r8, 120.110000_r8, &
194 : 48.044000_r8, 48.044000_r8, 48.044000_r8, 48.044000_r8, 48.044000_r8, &
195 : 48.044000_r8, 48.044000_r8, 48.044000_r8, 48.044000_r8, 120.110000_r8, &
196 : 0.000000_r8, 0.000000_r8, 0.000000_r8, 60.055000_r8, 0.000000_r8, &
197 : 0.000000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, 12.011000_r8, &
198 : 12.011000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, 24.022000_r8, &
199 : 36.033000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, &
200 : 0.000000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, &
201 : 12.011000_r8, 48.044000_r8, 24.022000_r8, 84.077000_r8, 72.066000_r8, &
202 : 72.066000_r8, 72.066000_r8, 12.011000_r8, 12.011000_r8, 36.033000_r8, &
203 : 36.033000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, &
204 : 0.000000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, 0.000000_r8, &
205 : 180.165000_r8, 180.165000_r8, 180.165000_r8, 180.165000_r8, 180.165000_r8, &
206 : 180.165000_r8, 180.165000_r8, 180.165000_r8, 180.165000_r8, 180.165000_r8, &
207 : 180.165000_r8, 180.165000_r8, 180.165000_r8, 180.165000_r8, 180.165000_r8, &
208 : 180.165000_r8, 12.011000_r8, 264.242000_r8, 84.077000_r8, 120.110000_r8, &
209 : 120.110000_r8, 120.110000_r8, 108.099000_r8, 108.099000_r8, 108.099000_r8, &
210 : 108.099000_r8, 120.110000_r8, 108.099000_r8, 108.099000_r8, 120.110000_r8, &
211 : 120.110000_r8, 120.110000_r8, 84.077000_r8, 120.110000_r8, 120.110000_r8, &
212 : 108.099000_r8, 120.110000_r8, 120.110000_r8, 120.110000_r8, 120.110000_r8, &
213 : 120.110000_r8, 120.110000_r8, 120.110000_r8, 120.110000_r8, 120.110000_r8, &
214 : 120.110000_r8, 84.077000_r8, 84.077000_r8, 96.088000_r8, 96.088000_r8, &
215 : 96.088000_r8, 96.088000_r8, 0.000000_r8, 0.000000_r8, 84.077000_r8, &
216 : 60.055000_r8, 120.110000_r8, 120.110000_r8, 120.110000_r8, 180.165000_r8, &
217 : 180.165000_r8, 180.165000_r8, 72.066000_r8, 72.066000_r8, 120.110000_r8, &
218 : 120.110000_r8, 120.110000_r8, 84.077000_r8, 24.022000_r8, 36.033000_r8, &
219 : 72.066000_r8, 24.022000_r8, 12.011000_r8, 60.055000_r8, 48.044000_r8, &
220 : 24.022000_r8, 24.022000_r8, 0.000000_r8, 12.011000_r8, 60.055000_r8, &
221 : 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, &
222 : 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, &
223 : 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, &
224 : 60.055000_r8, 60.055000_r8, 60.055000_r8, 60.055000_r8, 156.143000_r8, &
225 : 120.110000_r8, 120.110000_r8, 120.110000_r8, 48.044000_r8, 48.044000_r8, &
226 : 48.044000_r8, 48.044000_r8, 48.044000_r8, 48.044000_r8, 120.110000_r8, &
227 : 120.110000_r8, 120.110000_r8, 60.055000_r8, 0.000000_r8, 0.000000_r8, &
228 : 72.066000_r8, 36.033000_r8, 36.033000_r8, 120.110000_r8, 120.110000_r8, &
229 : 108.099000_r8, 108.099000_r8, 108.099000_r8, 108.099000_r8, 108.099000_r8, &
230 : 72.066000_r8, 120.110000_r8, 120.110000_r8, 84.077000_r8, 120.110000_r8, &
231 : 120.110000_r8, 120.110000_r8, 120.110000_r8, 84.077000_r8, 84.077000_r8, &
232 1536 : 96.088000_r8, 96.088000_r8, 96.088000_r8, 0.000000_r8 /)
233 :
234 1536 : fix_mass(: 3) = (/ 0.00000000_r8, 31.9988000_r8, 28.0134800_r8 /)
235 :
236 4608 : clsmap(: 2,1) = (/ 233, 234 /)
237 : clsmap(:317,4) = (/ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, &
238 : 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, &
239 : 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, &
240 : 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, &
241 : 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, &
242 : 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, &
243 : 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, &
244 : 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, &
245 : 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, &
246 : 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, &
247 : 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, &
248 : 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, &
249 : 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, &
250 : 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, &
251 : 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, &
252 : 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, &
253 : 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, &
254 : 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, &
255 : 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, &
256 : 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, &
257 : 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, &
258 : 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, &
259 : 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, &
260 : 231, 232, 235, 236, 237, 238, 239, 240, 241, 242, &
261 : 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, &
262 : 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, &
263 : 263, 264, 265, 266, 267, 268, 269, 270, 271, 272, &
264 : 273, 274, 275, 276, 277, 278, 279, 280, 281, 282, &
265 : 283, 284, 285, 286, 287, 288, 289, 290, 291, 292, &
266 : 293, 294, 295, 296, 297, 298, 299, 300, 301, 302, &
267 : 303, 304, 305, 306, 307, 308, 309, 310, 311, 312, &
268 488448 : 313, 314, 315, 316, 317, 318, 319 /)
269 :
270 : permute(:317,4) = (/ 168, 165, 1, 219, 2, 3, 199, 76, 124, 77, &
271 : 132, 103, 121, 110, 81, 118, 202, 276, 88, 309, &
272 : 146, 4, 89, 113, 105, 149, 99, 114, 106, 228, &
273 : 123, 63, 100, 60, 71, 72, 64, 73, 65, 74, &
274 : 66, 137, 306, 156, 67, 230, 120, 61, 280, 266, &
275 : 173, 163, 238, 129, 279, 131, 315, 75, 58, 310, &
276 : 223, 5, 271, 247, 90, 92, 84, 155, 104, 6, &
277 : 7, 8, 9, 68, 217, 254, 235, 301, 251, 62, &
278 : 237, 69, 208, 91, 93, 108, 304, 79, 182, 227, &
279 : 107, 213, 302, 128, 195, 209, 193, 153, 154, 136, &
280 : 140, 258, 263, 187, 185, 241, 138, 218, 229, 250, &
281 : 204, 260, 234, 170, 272, 159, 150, 273, 166, 212, &
282 : 245, 186, 264, 53, 201, 269, 259, 233, 145, 109, &
283 : 141, 270, 265, 239, 177, 134, 80, 102, 274, 10, &
284 : 11, 12, 59, 13, 14, 15, 303, 312, 311, 232, &
285 : 236, 16, 17, 18, 19, 20, 305, 313, 21, 111, &
286 : 119, 82, 147, 70, 135, 78, 112, 22, 23, 148, &
287 : 122, 143, 24, 246, 216, 98, 25, 26, 27, 28, &
288 : 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, &
289 : 39, 40, 41, 42, 43, 44, 45, 46, 85, 282, &
290 : 133, 284, 220, 157, 255, 158, 171, 115, 116, 125, &
291 : 117, 293, 292, 215, 192, 203, 194, 160, 207, 172, &
292 : 244, 252, 242, 253, 161, 151, 167, 83, 86, 178, &
293 : 87, 126, 144, 211, 281, 287, 47, 277, 286, 48, &
294 : 139, 49, 290, 285, 50, 127, 225, 198, 164, 300, &
295 : 308, 184, 152, 101, 179, 307, 162, 221, 267, 268, &
296 : 94, 95, 96, 97, 256, 257, 205, 210, 200, 206, &
297 : 231, 275, 222, 226, 51, 262, 261, 52, 289, 283, &
298 : 54, 243, 180, 248, 188, 169, 240, 291, 288, 55, &
299 : 249, 314, 316, 130, 189, 224, 190, 196, 294, 297, &
300 : 296, 298, 278, 295, 299, 174, 191, 214, 197, 175, &
301 488448 : 176, 181, 56, 183, 57, 142, 317 /)
302 :
303 : diag_map(:317) = (/ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, &
304 : 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, &
305 : 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, &
306 : 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, &
307 : 41, 42, 43, 44, 45, 51, 57, 63, 69, 75, &
308 : 81, 87, 89, 95, 101, 107, 113, 114, 117, 120, &
309 : 123, 126, 130, 134, 138, 142, 146, 150, 153, 156, &
310 : 159, 164, 169, 174, 179, 182, 188, 192, 197, 201, &
311 : 204, 207, 211, 218, 223, 229, 237, 242, 245, 248, &
312 : 251, 256, 261, 266, 269, 272, 275, 279, 283, 287, &
313 : 291, 295, 301, 304, 310, 316, 322, 326, 332, 337, &
314 : 342, 347, 353, 358, 363, 368, 373, 377, 385, 393, &
315 : 401, 405, 411, 417, 423, 429, 437, 443, 449, 456, &
316 : 462, 468, 471, 477, 484, 488, 495, 502, 514, 521, &
317 : 528, 536, 543, 550, 556, 561, 569, 577, 585, 593, &
318 : 601, 605, 613, 626, 639, 646, 655, 664, 673, 682, &
319 : 689, 693, 702, 711, 718, 728, 741, 752, 763, 770, &
320 : 778, 783, 791, 797, 805, 813, 827, 842, 855, 862, &
321 : 873, 883, 895, 906, 914, 919, 924, 929, 939, 951, &
322 : 960, 969, 974, 986, 995,1006,1015,1024,1041,1054, &
323 : 1071,1092,1106,1110,1121,1132,1143,1150,1159,1166, &
324 : 1182,1194,1210,1218,1228,1235,1242,1252,1272,1289, &
325 : 1298,1309,1325,1336,1350,1361,1377,1385,1410,1431, &
326 : 1440,1454,1461,1468,1476,1481,1487,1494,1498,1508, &
327 : 1524,1536,1543,1557,1572,1600,1613,1620,1640,1659, &
328 : 1673,1689,1702,1712,1718,1732,1754,1770,1778,1793, &
329 : 1817,1848,1872,1896,1916,1926,1953,1991,2016,2041, &
330 : 2057,2078,2110,2140,2180,2199,2213,2235,2259,2277, &
331 : 2296,2320,2339,2359,2380,2408,2441,2471,2505,2536, &
332 : 2570,2600,2629,2654,2677,2701,2745,2792,2840,2913, &
333 : 2931,2949,3106,3129,3161,3188,3369,3463,3488,3515, &
334 1536 : 3610,3663,3766,3808,3850,4101,4128 /)
335 :
336 : extfrc_lst(: 16) = (/ 'num_a4 ','pom_a4 ','bc_a4 ','SVOC ','SO2 ', &
337 : 'NO2 ','so4_a1 ','so4_a2 ','so4_a5 ','CO ', &
338 : 'num_a1 ','num_a2 ','num_a5 ','NO ','N ', &
339 26112 : 'OH ' /)
340 :
341 : frc_from_dataset(: 16) = (/ .true., .true., .true., .true., .true., &
342 : .true., .true., .true., .true., .true., &
343 : .true., .true., .true., .false., .false., &
344 1536 : .false. /)
345 :
346 6144 : inv_lst(: 3) = (/ 'M ', 'O2 ', 'N2 ' /)
347 :
348 : slvd_lst(: 84) = (/ 'ACBZO2 ', 'ALKO2 ', 'APINNO3 ', 'APINO2 ', 'APINO2VBS ', &
349 : 'BCARYNO3 ', 'BCARYO2 ', 'BCARYO2VBS ', 'BENZO2 ', 'BENZO2VBS ', &
350 : 'BPINNO3 ', 'BPINO2 ', 'BPINO2VBS ', 'BZOO ', 'C2H5O2 ', &
351 : 'C3H7O2 ', 'C6H5O2 ', 'CH3CO3 ', 'CH3O2 ', 'DICARBO2 ', &
352 : 'ENEO2 ', 'EO ', 'EO2 ', 'HO2 ', 'HOCH2OO ', &
353 : 'IEPOXOO ', 'ISOPB1O2 ', 'ISOPB4O2 ', 'ISOPC1C ', 'ISOPC1T ', &
354 : 'ISOPC4C ', 'ISOPC4T ', 'ISOPED1O2 ', 'ISOPED4O2 ', 'ISOPN1DO2 ', &
355 : 'ISOPN2BO2 ', 'ISOPN3BO2 ', 'ISOPN4DO2 ', 'ISOPNBNO3O2 ', 'ISOPNO3 ', &
356 : 'ISOPNOOHBO2 ', 'ISOPNOOHDO2 ', 'ISOPO2VBS ', 'ISOPZD1O2 ', 'ISOPZD4O2 ', &
357 : 'IVOCO2VBS ', 'LIMONNO3 ', 'LIMONO2 ', 'LIMONO2VBS ', 'MACRO2 ', &
358 : 'MALO2 ', 'MCO3 ', 'MDIALO2 ', 'MEKO2 ', 'MVKO2 ', &
359 : 'MYRCNO3 ', 'MYRCO2 ', 'MYRCO2VBS ', 'NC4CHOO2 ', 'O1D ', &
360 : 'OH ', 'PHENO2 ', 'PO2 ', 'RO2 ', 'TERP1OOHO2 ', &
361 : 'TERP2OOHO2 ', 'TERPA1O2 ', 'TERPA2CO3 ', 'TERPA2O2 ', 'TERPA3CO3 ', &
362 : 'TERPA3O2 ', 'TERPA4O2 ', 'TERPACO3 ', 'TERPF1O2 ', 'TERPF2O2 ', &
363 : 'TERPNPS1O2 ', 'TERPNPT1O2 ', 'TERPNS1O2 ', 'TERPNT1O2 ', 'TOLO2 ', &
364 130560 : 'TOLUO2VBS ', 'XYLENO2 ', 'XYLEO2VBS ', 'XYLOLO2 ' /)
365 :
366 1536 : if( allocated( rxt_tag_lst ) ) then
367 0 : deallocate( rxt_tag_lst )
368 : end if
369 1536 : allocate( rxt_tag_lst(rxt_tag_cnt),stat=ios )
370 1536 : if( ios /= 0 ) then
371 0 : write(iulog,*) 'set_sim_dat: failed to allocate rxt_tag_lst; error = ',ios
372 0 : call endrun
373 : end if
374 1536 : if( allocated( rxt_tag_map ) ) then
375 0 : deallocate( rxt_tag_map )
376 : end if
377 1536 : allocate( rxt_tag_map(rxt_tag_cnt),stat=ios )
378 1536 : if( ios /= 0 ) then
379 0 : write(iulog,*) 'set_sim_dat: failed to allocate rxt_tag_map; error = ',ios
380 0 : call endrun
381 : end if
382 1536 : rxt_tag_lst( 1: 200) = (/ 'jh2o_b ', 'jh2o_a ', &
383 : 'jh2o_c ', 'jh2o2 ', &
384 : 'jo2_a ', 'jo2_b ', &
385 : 'jo3_a ', 'jo3_b ', &
386 : 'jhno3 ', 'jho2no2_a ', &
387 : 'jho2no2_b ', 'jn2o ', &
388 : 'jn2o5_a ', 'jn2o5_b ', &
389 : 'jno ', 'jno2 ', &
390 : 'jno3_b ', 'jno3_a ', &
391 : 'jalknit ', 'jalkooh ', &
392 : 'jbenzooh ', 'jbepomuc ', &
393 : 'jbigald1 ', 'jbigald2 ', &
394 : 'jbigald3 ', 'jbigald4 ', &
395 : 'jbzooh ', 'jc2h5ooh ', &
396 : 'jc3h7ooh ', 'jc6h5ooh ', &
397 : 'jch2o_a ', 'jch2o_b ', &
398 : 'jch3cho ', 'jacet ', &
399 : 'jmgly ', 'jch3co3h ', &
400 : 'jch3ooh ', 'jch4_a ', &
401 : 'jch4_b ', 'jco2 ', &
402 : 'jdhpmpal ', 'jeooh ', &
403 : 'jglyald ', 'jglyoxal ', &
404 : 'jhcoch2ooh ', 'jhmhp ', &
405 : 'jhonitr ', 'jhpald1 ', &
406 : 'jhpald4 ', 'jhpaldb1c ', &
407 : 'jhpaldb4c ', 'jhyac ', &
408 : 'jhydrald_b ', 'jhydrald_a ', &
409 : 'jhyperacet_c ', 'jhyperacet_p ', &
410 : 'jinheb ', 'jinhed ', &
411 : 'jisopfdn ', 'jisopfdnc ', &
412 : 'jisopfnc ', 'jisopfnp ', &
413 : 'jisophfp ', 'jisopn1d ', &
414 : 'jisopn2b ', 'jisopn3b ', &
415 : 'jisopn4d ', 'jisopnbno3 ', &
416 : 'jisopnoohb ', 'jisopnoohd ', &
417 : 'jisopooh ', 'jmacr_b ', &
418 : 'jmacr_a ', 'jmacrn ', &
419 : 'jmacrooh ', 'jmek ', &
420 : 'jmekooh ', 'jmpan ', &
421 : 'jmvk ', 'jmvkn ', &
422 : 'jmvkooh ', 'jnc4cho ', &
423 : 'jno3ch2cho ', 'jnoa ', &
424 : 'jonitr ', 'jpan ', &
425 : 'jphenooh ', 'jpooh ', &
426 : 'jrooh ', 'jtepomuc ', &
427 : 'jterp1ooh ', 'jterp2aooh ', &
428 : 'jterpa ', 'jterpa2 ', &
429 : 'jterpa2pan ', 'jterpa3 ', &
430 : 'jterpa3pan ', 'jterpacid ', &
431 : 'jterpacid2 ', 'jterpacid3 ', &
432 : 'jterpapan ', 'jterpdhdp ', &
433 : 'jterpfdn ', 'jterphfn ', &
434 : 'jterpnps ', 'jterpnps1 ', &
435 : 'jterpnpt ', 'jterpnpt1 ', &
436 : 'jterpns ', 'jterpns1 ', &
437 : 'jterpnt ', 'jterpnt1 ', &
438 : 'jterpooh ', 'jterpoohl ', &
439 : 'jtolooh ', 'jxylenooh ', &
440 : 'jxylolooh ', 'jbrcl ', &
441 : 'jbro ', 'jbrono2_b ', &
442 : 'jbrono2_a ', 'jccl4 ', &
443 : 'jcf2clbr ', 'jcf3br ', &
444 : 'jcfcl3 ', 'jcfc113 ', &
445 : 'jcfc114 ', 'jcfc115 ', &
446 : 'jcf2cl2 ', 'jch2br2 ', &
447 : 'jch3br ', 'jch3ccl3 ', &
448 : 'jch3cl ', 'jchbr3 ', &
449 : 'jcl2 ', 'jcl2o2 ', &
450 : 'jclo ', 'jclono2_a ', &
451 : 'jclono2_b ', 'jcof2 ', &
452 : 'jcofcl ', 'jh2402 ', &
453 : 'jhbr ', 'jhcfc141b ', &
454 : 'jhcfc142b ', 'jhcfc22 ', &
455 : 'jhcl ', 'jhf ', &
456 : 'jhobr ', 'jhocl ', &
457 : 'joclo ', 'jsf6 ', &
458 : 'jh2so4 ', 'jocs ', &
459 : 'jso ', 'jso2 ', &
460 : 'jso3 ', 'jsoa1_a1 ', &
461 : 'jsoa1_a2 ', 'jsoa2_a1 ', &
462 : 'jsoa2_a2 ', 'jsoa3_a1 ', &
463 : 'jsoa3_a2 ', 'jsoa4_a1 ', &
464 : 'jsoa4_a2 ', 'jsoa5_a1 ', &
465 : 'jsoa5_a2 ', 'O1D_H2 ', &
466 : 'O1D_H2O ', 'O1D_N2 ', &
467 : 'O1D_O2ab ', 'O1D_O3 ', &
468 : 'O_O3 ', 'usr_O_O ', &
469 : 'usr_O_O2 ', 'H2_O ', &
470 : 'H2O2_O ', 'H_HO2 ', &
471 : 'H_HO2a ', 'H_HO2b ', &
472 : 'H_O2 ', 'HO2_O ', &
473 : 'HO2_O3 ', 'H_O3 ', &
474 : 'OH_H2 ', 'OH_H2O2 ', &
475 : 'OH_HO2 ', 'OH_O ', &
476 : 'OH_O3 ', 'OH_OH ', &
477 : 'OH_OH_M ', 'usr_HO2_HO2 ', &
478 : 'HO2NO2_OH ', 'N_NO ', &
479 : 'N_NO2a ', 'N_NO2b ', &
480 : 'N_NO2c ', 'N_O2 ', &
481 310272 : 'NO2_O ', 'NO2_O3 ' /)
482 1536 : rxt_tag_lst( 201: 400) = (/ 'NO2_O_M ', 'NO3_HO2 ', &
483 : 'NO3_NO ', 'NO3_O ', &
484 : 'NO3_OH ', 'N_OH ', &
485 : 'NO_HO2 ', 'NO_O3 ', &
486 : 'NO_O_M ', 'O1D_N2Oa ', &
487 : 'O1D_N2Ob ', 'tag_NO2_HO2 ', &
488 : 'tag_NO2_NO3 ', 'tag_NO2_OH ', &
489 : 'usr_HNO3_OH ', 'usr_HO2NO2_M ', &
490 : 'usr_N2O5_M ', 'CL_CH2O ', &
491 : 'CL_CH4 ', 'CL_H2 ', &
492 : 'CL_H2O2 ', 'CL_HO2a ', &
493 : 'CL_HO2b ', 'CL_O3 ', &
494 : 'CLO_CH3O2 ', 'CLO_CLOa ', &
495 : 'CLO_CLOb ', 'CLO_CLOc ', &
496 : 'CLO_HO2 ', 'CLO_NO ', &
497 : 'CLONO2_CL ', 'CLO_NO2_M ', &
498 : 'CLONO2_O ', 'CLONO2_OH ', &
499 : 'CLO_O ', 'CLO_OHa ', &
500 : 'CLO_OHb ', 'HCL_O ', &
501 : 'HCL_OH ', 'HOCL_CL ', &
502 : 'HOCL_O ', 'HOCL_OH ', &
503 : 'O1D_CCL4 ', 'O1D_CF2CLBR ', &
504 : 'O1D_CFC11 ', 'O1D_CFC113 ', &
505 : 'O1D_CFC114 ', 'O1D_CFC115 ', &
506 : 'O1D_CFC12 ', 'O1D_HCLa ', &
507 : 'O1D_HCLb ', 'tag_CLO_CLO_M ', &
508 : 'usr_CL2O2_M ', 'BR_CH2O ', &
509 : 'BR_HO2 ', 'BR_O3 ', &
510 : 'BRO_BRO ', 'BRO_CLOa ', &
511 : 'BRO_CLOb ', 'BRO_CLOc ', &
512 : 'BRO_HO2 ', 'BRO_NO ', &
513 : 'BRO_NO2_M ', 'BRONO2_O ', &
514 : 'BRO_O ', 'BRO_OH ', &
515 : 'HBR_O ', 'HBR_OH ', &
516 : 'HOBR_O ', 'O1D_CF3BR ', &
517 : 'O1D_CHBR3 ', 'O1D_H2402 ', &
518 : 'O1D_HBRa ', 'O1D_HBRb ', &
519 : 'F_CH4 ', 'F_H2 ', &
520 : 'F_H2O ', 'F_HNO3 ', &
521 : 'O1D_COF2 ', 'O1D_COFCL ', &
522 : 'CH2BR2_CL ', 'CH2BR2_OH ', &
523 : 'CH3BR_CL ', 'CH3BR_OH ', &
524 : 'CH3CCL3_OH ', 'CH3CL_CL ', &
525 : 'CH3CL_OH ', 'CHBR3_CL ', &
526 : 'CHBR3_OH ', 'HCFC141B_OH ', &
527 : 'HCFC142B_OH ', 'HCFC22_OH ', &
528 : 'O1D_CH2BR2 ', 'O1D_CH3BR ', &
529 : 'O1D_HCFC141B ', 'O1D_HCFC142B ', &
530 : 'O1D_HCFC22 ', 'CH2O_HO2 ', &
531 : 'CH2O_NO3 ', 'CH2O_O ', &
532 : 'CH2O_OH ', 'CH3O2_CH3O2a ', &
533 : 'CH3O2_CH3O2b ', 'CH3O2_HO2 ', &
534 : 'CH3O2_NO ', 'CH3OH_OH ', &
535 : 'CH3OOH_OH ', 'CH4_OH ', &
536 : 'HCN_OH ', 'HCOOH_OH ', &
537 : 'HMHP_OH ', 'HOCH2OO_HO2 ', &
538 : 'HOCH2OO_M ', 'HOCH2OO_NO ', &
539 : 'O1D_CH4a ', 'O1D_CH4b ', &
540 : 'O1D_CH4c ', 'O1D_HCN ', &
541 : 'usr_CO_OH ', 'C2H2_CL_M ', &
542 : 'C2H2_OH_M ', 'C2H4_CL_M ', &
543 : 'C2H4_O3 ', 'C2H5O2_C2H5O2 ', &
544 : 'C2H5O2_CH3O2 ', 'C2H5O2_HO2 ', &
545 : 'C2H5O2_NO ', 'C2H5OH_OH ', &
546 : 'C2H5OOH_OH ', 'C2H6_CL ', &
547 : 'C2H6_OH ', 'CH3CHO_NO3 ', &
548 : 'CH3CHO_OH ', 'CH3CN_OH ', &
549 : 'CH3CO3_CH3CO3 ', 'CH3CO3_CH3O2 ', &
550 : 'CH3CO3_HO2 ', 'CH3CO3_NO ', &
551 : 'CH3COOH_OH ', 'CH3COOOH_OH ', &
552 : 'EO2_HO2 ', 'EO2_NO ', &
553 : 'EO_M ', 'EO_O2 ', &
554 : 'GLYALD_OH ', 'GLYOXAL_OH ', &
555 : 'HCOCH2OOH_OH ', 'NO3CH2CHO_OH ', &
556 : 'PAN_OH ', 'tag_C2H4_OH ', &
557 : 'tag_CH3CO3_NO2 ', 'usr_PAN_M ', &
558 : 'C3H6_NO3 ', 'C3H6_O3 ', &
559 : 'C3H7O2_CH3O2 ', 'C3H7O2_HO2 ', &
560 : 'C3H7O2_NO ', 'C3H7OOH_OH ', &
561 : 'C3H8_OH ', 'CH3COCHO_NO3 ', &
562 : 'CH3COCHO_OH ', 'HYAC_OH ', &
563 : 'HYPERACET_OH ', 'NOA_OH ', &
564 : 'PO2_HO2 ', 'PO2_NO ', &
565 : 'POOH_OH ', 'RO2_CH3O2 ', &
566 : 'RO2_HO2 ', 'RO2_NO ', &
567 : 'ROOH_OH ', 'tag_C3H6_OH ', &
568 : 'usr_CH3COCH3_OH ', 'BIGENE_NO3 ', &
569 : 'BIGENE_OH ', 'DHPMPAL_OH ', &
570 : 'ENEO2_NO ', 'ENEO2_NOb ', &
571 : 'HONITR_OH ', 'MACRN_OH ', &
572 : 'MACRO2_CH3CO3 ', 'MACRO2_CH3O2 ', &
573 : 'MACRO2_HO2 ', 'MACRO2_isom ', &
574 : 'MACR_O3 ', 'MACR_OH ', &
575 : 'MACROOH_OH ', 'MCO3_CH3CO3 ', &
576 : 'MCO3_CH3O2 ', 'MCO3_HO2 ', &
577 : 'MCO3_MCO3 ', 'MCO3_NO ', &
578 : 'MCO3_NO3 ', 'MEKO2_HO2 ', &
579 : 'MEKO2_NO ', 'MEK_OH ', &
580 : 'MEKOOH_OH ', 'MPAN_OH_M ', &
581 310272 : 'MVKN_OH ', 'MVKO2_CH3CO3 ' /)
582 1536 : rxt_tag_lst( 401: 600) = (/ 'MVKO2_CH3O2 ', 'MVKO2_HO2 ', &
583 : 'MVK_O3 ', 'MVK_OH ', &
584 : 'MVKOOH_OH ', 'tag_MCO3_NO2 ', &
585 : 'usr_MPAN_M ', 'ALKNIT_OH ', &
586 : 'ALKO2_HO2 ', 'ALKO2_NO ', &
587 : 'ALKO2_NOb ', 'ALKOOH_OH ', &
588 : 'BIGALK_OH ', 'HPALD1_OH ', &
589 : 'HPALD4_OH ', 'HPALDB1C_OH ', &
590 : 'HPALDB4C_OH ', 'HYDRALD_OH ', &
591 : 'ICHE_OH ', 'IEPOX_OH ', &
592 : 'IEPOXOO_HO2 ', 'INHEB_OH ', &
593 : 'INHED_OH ', 'ISOPB1O2_CH3CO3 ', &
594 : 'ISOPB1O2_CH3O2 ', 'ISOPB1O2_HO2 ', &
595 : 'ISOPB1O2_I ', 'ISOPB1O2_M_C ', &
596 : 'ISOPB1O2_M_T ', 'ISOPB4O2_CH3CO3 ', &
597 : 'ISOPB4O2_CH3O2 ', 'ISOPB4O2_HO2 ', &
598 : 'ISOPB4O2_I ', 'ISOPB4O2_M_C ', &
599 : 'ISOPB4O2_M_T ', 'ISOPC1C_O2_B ', &
600 : 'ISOPC1C_O2_D ', 'ISOPC1T_O2_B ', &
601 : 'ISOPC1T_O2_D ', 'ISOPC4C_O2_B ', &
602 : 'ISOPC4C_O2_D ', 'ISOPC4T_O2_B ', &
603 : 'ISOPC4T_O2_D ', 'ISOPED1O2_CH3CO3 ', &
604 : 'ISOPED1O2_CH3O2 ', 'ISOPED1O2_HO2 ', &
605 : 'ISOPED1O2_M_C ', 'ISOPED4O2_CH3CO3 ', &
606 : 'ISOPED4O2_CH3O2 ', 'ISOPED4O2_HO2 ', &
607 : 'ISOPED4O2_M ', 'ISOPFDNC_OH ', &
608 : 'ISOPFDN_OH ', 'ISOPFNC_OH ', &
609 : 'ISOPFNP_OH ', 'ISOPHFP_OH ', &
610 : 'ISOPN1DO2_HO2 ', 'ISOPN1DO2_I ', &
611 : 'ISOPN1D_O3 ', 'ISOPN1D_OH ', &
612 : 'ISOPN2BO2_HO2 ', 'ISOPN2BO2_I ', &
613 : 'ISOPN2B_OH ', 'ISOPN3BO2_HO2 ', &
614 : 'ISOPN3BO2_I ', 'ISOPN3B_OH ', &
615 : 'ISOPN4DO2_HO2 ', 'ISOPN4DO2_I ', &
616 : 'ISOPN4D_O3 ', 'ISOPN4D_OH ', &
617 : 'ISOPNBNO3O2_HO2 ', 'ISOPNBNO3_OH ', &
618 : 'ISOP_NO3 ', 'ISOPNO3_CH3CO3 ', &
619 : 'ISOPNO3_CH3O2 ', 'ISOPNO3_HO2 ', &
620 : 'ISOPNO3_ISOPNO3 ', 'ISOPNO3_NO3 ', &
621 : 'ISOPNOOHBO2_HO2 ', 'ISOPNOOHBO2_I ', &
622 : 'ISOPNOOHB_OH ', 'ISOPNOOHDO2_HO2 ', &
623 : 'ISOPNOOHDO2_I ', 'ISOPNOOHD_O3 ', &
624 : 'ISOPNOOHD_OH ', 'ISOP_O3 ', &
625 : 'ISOP_OH ', 'ISOPOH_OH ', &
626 : 'ISOPOOH_OH_abs ', 'ISOPOOH_OH_add ', &
627 : 'ISOPZD1O2_CH3CO3 ', 'ISOPZD1O2_CH3O2 ', &
628 : 'ISOPZD1O2_HO2 ', 'ISOPZD1O2_M ', &
629 : 'ISOPZD4O2_CH3CO3 ', 'ISOPZD4O2_CH3O2 ', &
630 : 'ISOPZD4O2_HO2 ', 'ISOPZD4O2_M_C ', &
631 : 'NC4CHOO2_HO2 ', 'NC4CHOO2_isom ', &
632 : 'NC4CHO_O3 ', 'NC4CHO_OH ', &
633 : 'usr_IEPOXOO_NOa ', 'usr_IEPOXOO_NOn ', &
634 : 'usr_ISOPB1O2_NOa ', 'usr_ISOPB1O2_NOn ', &
635 : 'usr_ISOPB4O2_NOa ', 'usr_ISOPB4O2_NOn ', &
636 : 'usr_ISOPED1O2_NOa ', 'usr_ISOPED1O2_NOn ', &
637 : 'usr_ISOPED4O2_NOa ', 'usr_ISOPED4O2_NOn ', &
638 : 'usr_ISOPN1DO2_NOa ', 'usr_ISOPN1DO2_NOn ', &
639 : 'usr_ISOPN2BO2_NOa ', 'usr_ISOPN2BO2_NOn ', &
640 : 'usr_ISOPN3BO2_NOa ', 'usr_ISOPN3BO2_NOn ', &
641 : 'usr_ISOPN4DO2_NOa ', 'usr_ISOPN4DO2_NOn ', &
642 : 'usr_ISOPNBNO3O2_NOa ', 'usr_ISOPNBNO3O2_NOn ', &
643 : 'usr_ISOPNO3_NOa ', 'usr_ISOPNO3_NOn ', &
644 : 'usr_ISOPNOOHBO2_NOa ', 'usr_ISOPNOOHBO2_NOn ', &
645 : 'usr_ISOPNOOHDO2_NOa ', 'usr_ISOPNOOHDO2_NOn ', &
646 : 'usr_ISOPZD1O2 ', 'usr_ISOPZD1O2_NOa ', &
647 : 'usr_ISOPZD1O2_NOn ', 'usr_ISOPZD4O2 ', &
648 : 'usr_ISOPZD4O2_NOa ', 'usr_ISOPZD4O2_NOn ', &
649 : 'usr_MACRO2_NOa ', 'usr_MACRO2_NOn ', &
650 : 'usr_MVKO2_NOa ', 'usr_MVKO2_NOn ', &
651 : 'usr_NC4CHOO2_NOa ', 'usr_NC4CHOO2_NOn ', &
652 : 'ACBZO2_HO2 ', 'ACBZO2_NO ', &
653 : 'BENZENE_OH ', 'BENZO2_HO2 ', &
654 : 'BENZO2_NO ', 'BENZOOH_OH ', &
655 : 'BZALD_OH ', 'BZOO_HO2 ', &
656 : 'BZOOH_OH ', 'BZOO_NO ', &
657 : 'C6H5O2_HO2 ', 'C6H5O2_NO ', &
658 : 'C6H5OOH_OH ', 'CRESOL_OH ', &
659 : 'DICARBO2_HO2 ', 'DICARBO2_NO ', &
660 : 'DICARBO2_NO2 ', 'MALO2_HO2 ', &
661 : 'MALO2_NO ', 'MALO2_NO2 ', &
662 : 'MDIALO2_HO2 ', 'MDIALO2_NO ', &
663 : 'MDIALO2_NO2 ', 'PHENO2_HO2 ', &
664 : 'PHENO2_NO ', 'PHENOL_OH ', &
665 : 'PHENO_NO2 ', 'PHENO_O3 ', &
666 : 'PHENOOH_OH ', 'tag_ACBZO2_NO2 ', &
667 : 'TOLO2_HO2 ', 'TOLO2_NO ', &
668 : 'TOLOOH_OH ', 'TOLUENE_OH ', &
669 : 'usr_PBZNIT_M ', 'XYLENES_OH ', &
670 : 'XYLENO2_HO2 ', 'XYLENO2_NO ', &
671 : 'XYLENOOH_OH ', 'XYLOLO2_HO2 ', &
672 : 'XYLOLO2_NO ', 'XYLOL_OH ', &
673 : 'XYLOLOOH_OH ', 'APIN_NO3 ', &
674 : 'APINNO3_APINNO3 ', 'APINNO3_CH3CO3 ', &
675 : 'APINNO3_CH3O2 ', 'APINNO3_HO2 ', &
676 : 'APINNO3_NO ', 'APINNO3_NO3 ', &
677 : 'APINNO3_TERPA2CO3 ', 'APINNO3_TERPA3CO3 ', &
678 : 'APINNO3_TERPACO3 ', 'APINO2_CH3CO3 ', &
679 : 'APINO2_CH3O2 ', 'APINO2_HO2 ', &
680 : 'APINO2_NO ', 'APINO2_NO3 ', &
681 310272 : 'APINO2_TERPA2CO3 ', 'APINO2_TERPA3CO3 ' /)
682 1536 : rxt_tag_lst( 601: 800) = (/ 'APINO2_TERPACO3 ', 'APIN_O3 ', &
683 : 'APIN_OH ', 'BCARY_NO3 ', &
684 : 'BCARYNO3_BCARYNO3 ', 'BCARYNO3_CH3CO3 ', &
685 : 'BCARYNO3_CH3O2 ', 'BCARYNO3_HO2 ', &
686 : 'BCARYNO3_NO ', 'BCARYNO3_NO3 ', &
687 : 'BCARYNO3_TERPA2CO3 ', 'BCARYNO3_TERPA3CO3 ', &
688 : 'BCARYNO3_TERPACO3 ', 'BCARYO2_CH3CO3 ', &
689 : 'BCARYO2_CH3O2 ', 'BCARYO2_HO2 ', &
690 : 'BCARYO2_NO ', 'BCARYO2_NO3 ', &
691 : 'BCARYO2_TERPA2CO3 ', 'BCARYO2_TERPA3CO3 ', &
692 : 'BCARYO2_TERPACO3 ', 'BCARY_O3 ', &
693 : 'BCARY_OH ', 'BPIN_NO3 ', &
694 : 'BPINNO3_BPINNO3 ', 'BPINNO3_CH3CO3 ', &
695 : 'BPINNO3_CH3O2 ', 'BPINNO3_HO2 ', &
696 : 'BPINNO3_NO ', 'BPINNO3_NO3 ', &
697 : 'BPINNO3_TERPA2CO3 ', 'BPINNO3_TERPA3CO3 ', &
698 : 'BPINNO3_TERPACO3 ', 'BPINO2_CH3CO3 ', &
699 : 'BPINO2_CH3O2 ', 'BPINO2_HO2 ', &
700 : 'BPINO2_NO ', 'BPINO2_NO3 ', &
701 : 'BPINO2_TERPA2CO3 ', 'BPINO2_TERPA3CO3 ', &
702 : 'BPINO2_TERPACO3 ', 'BPIN_O3 ', &
703 : 'BPIN_OH ', 'LIMON_NO3 ', &
704 : 'LIMONNO3_CH3CO3 ', 'LIMONNO3_CH3O2 ', &
705 : 'LIMONNO3_HO2 ', 'LIMONNO3_LIMONNO3 ', &
706 : 'LIMONNO3_NO ', 'LIMONNO3_NO3 ', &
707 : 'LIMONNO3_TERPA2CO3 ', 'LIMONNO3_TERPA3CO3 ', &
708 : 'LIMONNO3_TERPACO3 ', 'LIMONO2_CH3CO3 ', &
709 : 'LIMONO2_CH3O2 ', 'LIMONO2_HO2 ', &
710 : 'LIMONO2_NO ', 'LIMONO2_NO3 ', &
711 : 'LIMONO2_TERPA2CO3 ', 'LIMONO2_TERPA3CO3 ', &
712 : 'LIMONO2_TERPACO3 ', 'LIMON_O3 ', &
713 : 'LIMON_OH ', 'MYRC_NO3 ', &
714 : 'MYRCNO3_CH3CO3 ', 'MYRCNO3_CH3O2 ', &
715 : 'MYRCNO3_HO2 ', 'MYRCNO3_MYRCNO3 ', &
716 : 'MYRCNO3_NO ', 'MYRCNO3_NO3 ', &
717 : 'MYRCNO3_TERPA2CO3 ', 'MYRCNO3_TERPA3CO3 ', &
718 : 'MYRCNO3_TERPACO3 ', 'MYRCO2_CH3CO3 ', &
719 : 'MYRCO2_CH3O2 ', 'MYRCO2_HO2 ', &
720 : 'MYRCO2_NO ', 'MYRCO2_NO3 ', &
721 : 'MYRCO2_TERPA2CO3 ', 'MYRCO2_TERPA3CO3 ', &
722 : 'MYRCO2_TERPACO3 ', 'MYRC_O3 ', &
723 : 'MYRC_OH ', 'tag_TERPA2CO3_NO2 ', &
724 : 'tag_TERPA3CO3_NO2 ', 'tag_TERPACO3_NO2 ', &
725 : 'TERP1OOHO2_HO2 ', 'TERP1OOHO2_NO ', &
726 : 'TERP1OOH_OH ', 'TERP2AOOH_OH ', &
727 : 'TERP2OOHO2_HO2 ', 'TERP2OOHO2_NO ', &
728 : 'TERPA1O2_CH3CO3 ', 'TERPA1O2_CH3O2 ', &
729 : 'TERPA1O2_HO2 ', 'TERPA1O2_NO ', &
730 : 'TERPA1O2_NO3 ', 'TERPA1O2_TERPA2CO3 ', &
731 : 'TERPA1O2_TERPA3CO3 ', 'TERPA1O2_TERPACO3 ', &
732 : 'TERPA2CO3_CH3CO3 ', 'TERPA2CO3_CH3O2 ', &
733 : 'TERPA2CO3_HO2 ', 'TERPA2CO3_NO ', &
734 : 'TERPA2CO3_NO3 ', 'TERPA2CO3_TERPA2CO3 ', &
735 : 'TERPA2CO3_TERPACO3 ', 'TERPA2_NO3 ', &
736 : 'TERPA2O2_CH3CO3 ', 'TERPA2O2_CH3O2 ', &
737 : 'TERPA2O2_HO2 ', 'TERPA2O2_NO ', &
738 : 'TERPA2O2_NO3 ', 'TERPA2O2_TERPA2CO3 ', &
739 : 'TERPA2O2_TERPA3CO3 ', 'TERPA2O2_TERPACO3 ', &
740 : 'TERPA2_OH ', 'TERPA2PAN_OH ', &
741 : 'TERPA3CO3_CH3CO3 ', 'TERPA3CO3_CH3O2 ', &
742 : 'TERPA3CO3_HO2 ', 'TERPA3CO3_NO ', &
743 : 'TERPA3CO3_NO3 ', 'TERPA3CO3_TERPA2CO3 ', &
744 : 'TERPA3CO3_TERPA3CO3 ', 'TERPA3CO3_TERPACO3 ', &
745 : 'TERPA3_NO3 ', 'TERPA3O2_CH3CO3 ', &
746 : 'TERPA3O2_CH3O2 ', 'TERPA3O2_HO2 ', &
747 : 'TERPA3O2_NO ', 'TERPA3O2_NO3 ', &
748 : 'TERPA3O2_TERPA2CO3 ', 'TERPA3O2_TERPA3CO3 ', &
749 : 'TERPA3O2_TERPACO3 ', 'TERPA3_OH ', &
750 : 'TERPA3PAN_OH ', 'TERPA4O2_CH3CO3 ', &
751 : 'TERPA4O2_CH3O2 ', 'TERPA4O2_HO2 ', &
752 : 'TERPA4O2_NO ', 'TERPA4O2_NO3 ', &
753 : 'TERPA4O2_TERPA2CO3 ', 'TERPA4O2_TERPA3CO3 ', &
754 : 'TERPA4O2_TERPACO3 ', 'TERPACID2_OH ', &
755 : 'TERPACID3_OH ', 'TERPACID_OH ', &
756 : 'TERPACO3_CH3CO3 ', 'TERPACO3_CH3O2 ', &
757 : 'TERPACO3_HO2 ', 'TERPACO3_NO ', &
758 : 'TERPACO3_NO3 ', 'TERPACO3_TERPACO3 ', &
759 : 'TERPA_NO3 ', 'TERPA_OH ', &
760 : 'TERPAPAN_OH ', 'TERPDHDP_OH ', &
761 : 'TERPF1_NO3 ', 'TERPF1O2_HO2 ', &
762 : 'TERPF1O2_NO ', 'TERPF1_O3 ', &
763 : 'TERPF1_OH ', 'TERPF2_NO3 ', &
764 : 'TERPF2O2_HO2 ', 'TERPF2O2_NO ', &
765 : 'TERPF2_O3 ', 'TERPF2_OH ', &
766 : 'TERPFDN_OH ', 'TERPHFN_OH ', &
767 : 'TERPK_OH ', 'TERPNPS1O2_HO2 ', &
768 : 'TERPNPS1O2_NO ', 'TERPNPS1_OH ', &
769 : 'TERPNPS_OH ', 'TERPNPT1O2_HO2 ', &
770 : 'TERPNPT1O2_NO ', 'TERPNPT1_OH ', &
771 : 'TERPNPT_OH ', 'TERPNS1O2_HO2 ', &
772 : 'TERPNS1O2_NO ', 'TERPNS1_OH ', &
773 : 'TERPNS_OH ', 'TERPNT1O2_HO2 ', &
774 : 'TERPNT1O2_NO ', 'TERPNT1_OH ', &
775 : 'TERPNT_OH ', 'TERPOOHL_OH ', &
776 : 'TERPOOH_OH ', 'usr_TERPA2PAN_M ', &
777 : 'usr_TERPA3PAN_M ', 'usr_TERPAPAN_M ', &
778 : 'DMS_NO3 ', 'DMS_OHa ', &
779 : 'OCS_O ', 'OCS_OH ', &
780 : 'S_O2 ', 'SO2_OH_M ', &
781 310272 : 'S_O3 ', 'SO_BRO ' /)
782 1536 : rxt_tag_lst( 801: 904) = (/ 'SO_CLO ', 'S_OH ', &
783 : 'SO_NO2 ', 'SO_O2 ', &
784 : 'SO_O3 ', 'SO_OCLO ', &
785 : 'SO_OH ', 'usr_DMS_OH ', &
786 : 'usr_SO3_H2O ', 'NH3_OH ', &
787 : 'usr_GLYOXAL_aer ', 'usr_HO2_aer ', &
788 : 'usr_HONITR_aer ', 'usr_ICHE_aer ', &
789 : 'usr_IEPOX_aer ', 'usr_INHEB_aer ', &
790 : 'usr_INHED_aer ', 'usr_INOOHD_aer ', &
791 : 'usr_ISOPFDN_aer ', 'usr_ISOPFDNC_aer ', &
792 : 'usr_ISOPFNC_aer ', 'usr_ISOPFNP_aer ', &
793 : 'usr_ISOPHFP_aer ', 'usr_ISOPN1D_aer ', &
794 : 'usr_ISOPN2B_aer ', 'usr_ISOPN4D_aer ', &
795 : 'usr_N2O5_aer ', 'usr_NC4CHO_aer ', &
796 : 'usr_NH4_strat_tau ', 'usr_NO2_aer ', &
797 : 'usr_NO3_aer ', 'usr_ONITR_aer ', &
798 : 'usr_SQTN_aer ', 'usr_TERPDHDP_aer ', &
799 : 'usr_TERPFDN_aer ', 'usr_TERPHFN_aer ', &
800 : 'usr_TERPNPT1_aer ', 'usr_TERPNPT_aer ', &
801 : 'usr_TERPNT1_aer ', 'usr_TERPNT_aer ', &
802 : 'APIN_NO3_vbs ', 'APINO2_HO2_vbs ', &
803 : 'APINO2_NO_vbs ', 'APIN_O3_vbs ', &
804 : 'APIN_OH_vbs ', 'BCARY_NO3_vbs ', &
805 : 'BCARYO2_HO2_vbs ', 'BCARYO2_NO_vbs ', &
806 : 'BCARY_O3_vbs ', 'BCARY_OH_vbs ', &
807 : 'BENZENE_OH_vbs ', 'BENZO2_HO2_vbs ', &
808 : 'BENZO2_NO_vbs ', 'BPIN_NO3_vbs ', &
809 : 'BPINO2_HO2_vbs ', 'BPINO2_NO_vbs ', &
810 : 'BPIN_O3_vbs ', 'BPIN_OH_vbs ', &
811 : 'ISOP_NO3_vbs ', 'ISOPO2_HO2_vbs ', &
812 : 'ISOPO2_NO_vbs ', 'ISOP_O3_vbs ', &
813 : 'ISOP_OH_vbs ', 'IVOCO2_HO2_vbs ', &
814 : 'IVOCO2_NO_vbs ', 'IVOC_OH_vbs ', &
815 : 'LIMON_NO3_vbs ', 'LIMONO2_HO2_vbs ', &
816 : 'LIMONO2_NO_vbs ', 'LIMON_O3_vbs ', &
817 : 'LIMON_OH_vbs ', 'MYRC_NO3_vbs ', &
818 : 'MYRCO2_HO2_vbs ', 'MYRCO2_NO_vbs ', &
819 : 'MYRC_O3_vbs ', 'MYRC_OH_vbs ', &
820 : 'SVOC_OH ', 'TOLUENE_OH_vbs ', &
821 : 'TOLUO2_HO2_vbs ', 'TOLUO2_NO_vbs ', &
822 : 'XYLENES_OH_vbs ', 'XYLEO2_HO2_vbs ', &
823 : 'XYLEO2_NO_vbs ', 'het1 ', &
824 : 'het10 ', 'het11 ', &
825 : 'het12 ', 'het13 ', &
826 : 'het14 ', 'het15 ', &
827 : 'het16 ', 'het17 ', &
828 : 'het2 ', 'het3 ', &
829 : 'het4 ', 'het5 ', &
830 : 'het6 ', 'het7 ', &
831 : 'het8 ', 'het9 ', &
832 : 'E90_tau ', 'NH_50_tau ', &
833 162816 : 'NH_5_tau ', 'ST80_25_tau ' /)
834 0 : rxt_tag_map(:rxt_tag_cnt) = (/ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, &
835 : 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, &
836 : 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, &
837 : 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, &
838 : 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, &
839 : 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, &
840 : 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, &
841 : 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, &
842 : 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, &
843 : 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, &
844 : 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, &
845 : 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, &
846 : 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, &
847 : 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, &
848 : 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, &
849 : 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, &
850 : 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, &
851 : 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, &
852 : 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, &
853 : 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, &
854 : 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, &
855 : 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, &
856 : 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, &
857 : 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, &
858 : 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, &
859 : 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, &
860 : 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, &
861 : 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, &
862 : 281, 282, 283, 284, 285, 286, 287, 288, 289, 290, &
863 : 291, 292, 293, 294, 295, 296, 297, 298, 299, 300, &
864 : 301, 302, 303, 304, 305, 306, 307, 308, 309, 310, &
865 : 311, 312, 313, 314, 315, 316, 317, 318, 319, 320, &
866 : 321, 322, 323, 324, 325, 326, 327, 328, 329, 330, &
867 : 331, 332, 333, 334, 335, 336, 337, 338, 339, 340, &
868 : 341, 342, 343, 344, 345, 346, 347, 348, 349, 350, &
869 : 351, 352, 353, 354, 355, 356, 357, 358, 359, 360, &
870 : 361, 362, 363, 364, 365, 366, 367, 368, 369, 370, &
871 : 371, 372, 373, 374, 375, 376, 377, 378, 379, 380, &
872 : 381, 382, 383, 384, 385, 386, 387, 388, 389, 390, &
873 : 391, 392, 393, 394, 395, 396, 397, 398, 399, 400, &
874 : 401, 402, 403, 404, 405, 406, 407, 408, 409, 410, &
875 : 411, 412, 413, 414, 415, 416, 417, 418, 419, 420, &
876 : 421, 422, 423, 424, 425, 426, 427, 428, 429, 430, &
877 : 431, 432, 433, 434, 435, 436, 437, 438, 439, 440, &
878 : 441, 442, 443, 444, 445, 446, 447, 448, 449, 450, &
879 : 451, 452, 453, 454, 455, 456, 457, 458, 459, 460, &
880 : 461, 462, 463, 464, 465, 466, 467, 468, 469, 470, &
881 : 471, 472, 473, 474, 475, 476, 477, 478, 479, 480, &
882 : 481, 482, 483, 484, 485, 486, 487, 488, 489, 490, &
883 : 491, 492, 493, 494, 495, 496, 497, 498, 499, 500, &
884 : 501, 502, 503, 504, 505, 506, 507, 508, 509, 510, &
885 : 511, 512, 513, 514, 515, 516, 517, 518, 519, 520, &
886 : 521, 522, 523, 524, 525, 526, 527, 528, 529, 530, &
887 : 531, 532, 533, 534, 535, 536, 537, 538, 539, 540, &
888 : 541, 542, 543, 544, 545, 546, 547, 548, 549, 550, &
889 : 551, 552, 553, 554, 555, 556, 557, 558, 559, 560, &
890 : 561, 562, 563, 564, 565, 566, 567, 568, 569, 570, &
891 : 571, 572, 573, 574, 575, 576, 577, 578, 579, 580, &
892 : 581, 582, 583, 584, 585, 586, 587, 588, 589, 590, &
893 : 591, 592, 593, 594, 595, 596, 597, 598, 599, 600, &
894 : 601, 602, 603, 604, 605, 606, 607, 608, 609, 610, &
895 : 611, 612, 613, 614, 615, 616, 617, 618, 619, 620, &
896 : 621, 622, 623, 624, 625, 626, 627, 628, 629, 630, &
897 : 631, 632, 633, 634, 635, 636, 637, 638, 639, 640, &
898 : 641, 642, 643, 644, 645, 646, 647, 648, 649, 650, &
899 : 651, 652, 653, 654, 655, 656, 657, 658, 659, 660, &
900 : 661, 662, 663, 664, 665, 666, 667, 668, 669, 670, &
901 : 671, 672, 673, 674, 675, 676, 677, 678, 679, 680, &
902 : 681, 682, 683, 684, 685, 686, 687, 688, 689, 690, &
903 : 691, 692, 693, 694, 695, 696, 697, 698, 699, 700, &
904 : 701, 702, 703, 704, 705, 706, 707, 708, 709, 710, &
905 : 711, 712, 713, 714, 715, 716, 717, 718, 719, 720, &
906 : 721, 722, 723, 724, 725, 726, 727, 728, 729, 730, &
907 : 731, 732, 733, 734, 735, 736, 737, 738, 739, 740, &
908 : 741, 742, 743, 744, 745, 746, 747, 748, 749, 750, &
909 : 751, 752, 753, 754, 755, 756, 757, 758, 759, 760, &
910 : 761, 762, 763, 764, 765, 766, 767, 768, 769, 770, &
911 : 771, 772, 773, 774, 775, 776, 777, 778, 779, 780, &
912 : 781, 782, 783, 784, 785, 786, 787, 788, 789, 790, &
913 : 791, 792, 793, 794, 795, 796, 797, 798, 799, 800, &
914 : 801, 802, 803, 804, 805, 806, 807, 808, 809, 810, &
915 : 811, 812, 813, 814, 815, 816, 817, 818, 819, 820, &
916 : 821, 822, 823, 824, 825, 826, 827, 828, 829, 830, &
917 : 831, 832, 833, 834, 835, 836, 837, 838, 839, 840, &
918 : 841, 842, 843, 844, 845, 846, 847, 848, 849, 850, &
919 : 851, 852, 853, 854, 855, 856, 857, 858, 859, 860, &
920 : 861, 862, 863, 864, 865, 866, 867, 868, 869, 870, &
921 : 871, 872, 873, 874, 875, 876, 877, 878, 879, 880, &
922 : 881, 882, 883, 884, 885, 886, 887, 888, 889, 890, &
923 : 891, 892, 893, 894, 895, 896, 897, 898, 899, 900, &
924 1391616 : 901, 902, 903, 904 /)
925 1536 : if( allocated( pht_alias_lst ) ) then
926 0 : deallocate( pht_alias_lst )
927 : end if
928 1536 : allocate( pht_alias_lst(phtcnt,2),stat=ios )
929 1536 : if( ios /= 0 ) then
930 0 : write(iulog,*) 'set_sim_dat: failed to allocate pht_alias_lst; error = ',ios
931 0 : call endrun
932 : end if
933 1536 : if( allocated( pht_alias_mult ) ) then
934 0 : deallocate( pht_alias_mult )
935 : end if
936 1536 : allocate( pht_alias_mult(phtcnt,2),stat=ios )
937 1536 : if( ios /= 0 ) then
938 0 : write(iulog,*) 'set_sim_dat: failed to allocate pht_alias_mult; error = ',ios
939 0 : call endrun
940 : end if
941 0 : pht_alias_lst(:,1) = (/ ' ', ' ', ' ', ' ', &
942 : 'userdefined ', 'userdefined ', ' ', ' ', &
943 : ' ', ' ', ' ', ' ', &
944 : ' ', ' ', 'userdefined ', ' ', &
945 : ' ', ' ', ' ', ' ', &
946 : ' ', ' ', ' ', ' ', &
947 : ' ', ' ', ' ', ' ', &
948 : ' ', ' ', ' ', ' ', &
949 : ' ', ' ', ' ', ' ', &
950 : ' ', ' ', ' ', ' ', &
951 : ' ', ' ', ' ', ' ', &
952 : ' ', ' ', ' ', ' ', &
953 : ' ', ' ', ' ', ' ', &
954 : ' ', ' ', ' ', ' ', &
955 : ' ', ' ', ' ', ' ', &
956 : ' ', ' ', ' ', ' ', &
957 : ' ', ' ', ' ', ' ', &
958 : ' ', ' ', ' ', ' ', &
959 : ' ', ' ', ' ', ' ', &
960 : ' ', ' ', ' ', ' ', &
961 : ' ', ' ', ' ', ' ', &
962 : ' ', ' ', ' ', ' ', &
963 : ' ', ' ', ' ', ' ', &
964 : ' ', ' ', ' ', ' ', &
965 : ' ', ' ', ' ', ' ', &
966 : ' ', ' ', ' ', ' ', &
967 : ' ', ' ', ' ', ' ', &
968 : ' ', ' ', ' ', ' ', &
969 : ' ', ' ', ' ', ' ', &
970 : ' ', ' ', ' ', ' ', &
971 : ' ', ' ', ' ', ' ', &
972 : ' ', ' ', ' ', ' ', &
973 : ' ', ' ', ' ', ' ', &
974 : ' ', ' ', ' ', ' ', &
975 : ' ', ' ', ' ', ' ', &
976 : ' ', ' ', ' ', ' ', &
977 : ' ', ' ', ' ', ' ', &
978 : ' ', ' ', ' ', ' ', &
979 : ' ', ' ', ' ', ' ', &
980 : ' ', ' ', ' ', ' ', &
981 : ' ', ' ', ' ', ' ', &
982 258048 : ' ', ' ', ' ' /)
983 0 : pht_alias_lst(:,2) = (/ ' ', ' ', ' ', ' ', &
984 : ' ', ' ', ' ', ' ', &
985 : ' ', ' ', ' ', ' ', &
986 : ' ', ' ', ' ', ' ', &
987 : ' ', ' ', 'jch3ooh ', 'jch3ooh ', &
988 : 'jch3ooh ', 'jno2 ', 'jno2 ', 'jno2 ', &
989 : 'jno2 ', 'jno2 ', 'jch3ooh ', 'jch3ooh ', &
990 : 'jch3ooh ', 'jch3ooh ', ' ', ' ', &
991 : ' ', ' ', ' ', 'jh2o2 ', &
992 : ' ', ' ', ' ', ' ', &
993 : 'jch3ooh ', 'jch3ooh ', ' ', 'jmgly ', &
994 : 'jch3ooh ', 'jch3ooh ', 'jch2o_a ', 'jmacr_a ', &
995 : 'jmacr_a ', 'jch3ooh ', 'jch3ooh ', ' ', &
996 : 'jmacr_b ', 'jmacr_a ', 'jacet ', 'jch3ooh ', &
997 : 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', 'jch2o_a ', &
998 : 'jch2o_a ', 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', &
999 : 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', &
1000 : 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', ' ', &
1001 : ' ', 'jch2o_a ', 'jch3ooh ', 'jacet ', &
1002 : 'jch3ooh ', 'jpan ', ' ', 'jch2o_a ', &
1003 : 'jch3ooh ', 'jch2o_a ', 'jch2o_a ', 'jch2o_a ', &
1004 : 'jch3cho ', ' ', 'jch3ooh ', 'jch3ooh ', &
1005 : 'jch3ooh ', 'jno2 ', 'jch3ooh ', 'jch3ooh ', &
1006 : 'jch3cho ', 'jch3cho ', 'jpan ', 'jch3cho ', &
1007 : 'jpan ', 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', &
1008 : 'jpan ', 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', &
1009 : 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', &
1010 : 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', &
1011 : 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', 'jch3ooh ', &
1012 : 'jch3ooh ', ' ', ' ', ' ', &
1013 : ' ', ' ', ' ', ' ', &
1014 : ' ', ' ', ' ', ' ', &
1015 : ' ', ' ', ' ', ' ', &
1016 : ' ', ' ', ' ', ' ', &
1017 : ' ', ' ', ' ', ' ', &
1018 : ' ', ' ', ' ', ' ', &
1019 : ' ', ' ', ' ', ' ', &
1020 : ' ', ' ', ' ', ' ', &
1021 : ' ', ' ', ' ', ' ', &
1022 : ' ', 'jno2 ', 'jno2 ', 'jno2 ', &
1023 : 'jno2 ', 'jno2 ', 'jno2 ', 'jno2 ', &
1024 258048 : 'jno2 ', 'jno2 ', 'jno2 ' /)
1025 0 : pht_alias_mult(:,1) = (/ 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1026 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1027 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1028 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1029 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1030 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1031 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1032 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1033 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1034 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1035 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1036 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1037 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1038 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1039 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1040 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1041 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1042 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1043 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1044 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1045 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1046 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1047 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1048 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1049 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1050 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1051 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1052 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1053 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1054 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1055 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1056 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1057 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1058 258048 : 1._r8, 1._r8 /)
1059 0 : pht_alias_mult(:,2) = (/ 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1060 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1061 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1062 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1063 : 1._r8, .10_r8, .14_r8, .20_r8, .20_r8, &
1064 : .006_r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1065 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1066 : 0.28_r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1067 : 4.62_r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1068 : 0.75_r8, 1._r8, 110.0_r8, 110.0_r8, 4.62_r8, &
1069 : 4.62_r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1070 : 1._r8, 1._r8, 1._r8, 1._r8, 10.0_r8, &
1071 : 10.0_r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1072 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1073 : 1._r8, 1._r8, 1._r8, 5.8_r8, 1._r8, &
1074 : 1._r8, 1._r8, 1._r8, 1._r8, 1.26_r8, &
1075 : 1._r8, 9.2_r8, 4.3_r8, 1._r8, 1._r8, &
1076 : 1._r8, 1._r8, 1._r8, 1._r8, .10_r8, &
1077 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1078 : 1._r8, 1._r8, 0.71_r8, 0.71_r8, 0.71_r8, &
1079 : 1._r8, 2.0_r8, 1._r8, 1._r8, 1._r8, &
1080 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1081 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1082 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1083 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1084 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1085 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1086 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1087 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1088 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1089 : 1._r8, 1._r8, 1._r8, 1._r8, 1._r8, &
1090 : 1._r8, 1._r8, .0004_r8, .0004_r8, .0004_r8, &
1091 : .0004_r8, .0004_r8, .0004_r8, .0004_r8, .0004_r8, &
1092 258048 : .0004_r8, .0004_r8 /)
1093 1536 : allocate( cph_enthalpy(enthalpy_cnt),stat=ios )
1094 1536 : if( ios /= 0 ) then
1095 0 : write(iulog,*) 'set_sim_dat: failed to allocate cph_enthalpy; error = ',ios
1096 0 : call endrun
1097 : end if
1098 1536 : allocate( cph_rid(enthalpy_cnt),stat=ios )
1099 1536 : if( ios /= 0 ) then
1100 0 : write(iulog,*) 'set_sim_dat: failed to allocate cph_rid; error = ',ios
1101 0 : call endrun
1102 : end if
1103 : cph_rid(:) = (/ 170, 173, 174, 175, 178, &
1104 : 181, 182, 183, 184, 187, &
1105 : 188, 189, 192, 194, 198, &
1106 29184 : 199, 207, 208 /)
1107 : cph_enthalpy(:) = (/ 189.810000_r8, 392.190000_r8, 493.580000_r8, 101.390000_r8, 232.590000_r8, &
1108 : 203.400000_r8, 226.580000_r8, 120.100000_r8, 194.710000_r8, 293.620000_r8, &
1109 : 67.670000_r8, 165.300000_r8, 165.510000_r8, 313.750000_r8, 133.750000_r8, &
1110 29184 : 193.020000_r8, 34.470000_r8, 199.170000_r8 /)
1111 1536 : allocate( num_rnts(rxntot-phtcnt),stat=ios )
1112 1536 : if( ios /= 0 ) then
1113 0 : write(iulog,*) 'set_sim_dat: failed to allocate num_rnts; error = ',ios
1114 0 : call endrun
1115 : end if
1116 : num_rnts(:) = (/ 2, 2, 2, 2, 2, 2, 3, 3, 2, 2, &
1117 : 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, &
1118 : 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, &
1119 : 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, &
1120 : 2, 3, 2, 2, 3, 3, 3, 2, 2, 2, &
1121 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1122 : 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, &
1123 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1124 : 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, &
1125 : 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, &
1126 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1127 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1128 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1129 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1130 : 2, 3, 2, 2, 2, 1, 2, 2, 2, 2, &
1131 : 2, 2, 3, 3, 3, 2, 2, 2, 2, 2, &
1132 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1133 : 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, &
1134 : 2, 2, 3, 3, 2, 2, 2, 2, 2, 2, &
1135 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1136 : 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, &
1137 : 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, &
1138 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1139 : 3, 2, 2, 2, 2, 2, 2, 2, 3, 2, &
1140 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1141 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, &
1142 : 1, 1, 2, 2, 2, 1, 1, 1, 2, 2, &
1143 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, &
1144 : 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, &
1145 : 1, 2, 2, 2, 1, 2, 2, 1, 2, 2, &
1146 : 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1147 : 2, 2, 1, 2, 2, 1, 2, 2, 2, 2, &
1148 : 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, &
1149 : 1, 2, 1, 2, 2, 2, 2, 2, 2, 2, &
1150 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1151 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1152 : 2, 1, 2, 2, 1, 2, 2, 2, 2, 2, &
1153 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1154 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, &
1155 : 2, 2, 3, 2, 2, 3, 2, 2, 2, 2, &
1156 : 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, &
1157 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1158 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1159 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1160 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1161 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1162 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1163 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1164 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1165 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1166 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1167 : 2, 2, 2, 2, 2, 2, 3, 3, 3, 2, &
1168 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1169 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1170 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1171 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1172 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1173 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1174 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1175 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1176 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1177 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1178 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1179 : 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1180 : 2, 2, 2, 1, 1, 1, 1, 1, 1, 1, &
1181 : 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, &
1182 : 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, &
1183 : 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, &
1184 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1185 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1186 : 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, &
1187 : 2, 2, 2, 2, 2, 2, 1, 2, 1, 1, &
1188 : 1, 1, 2, 2, 2, 1, 1, 2, 2, 2, &
1189 1133568 : 1, 1, 2, 1, 1, 1, 1 /)
1190 :
1191 1536 : end subroutine set_sim_dat
1192 :
1193 : end module mo_sim_dat
|