BUGSrad
 All Classes Files Functions Variables
gases_ckd.f90
Go to the documentation of this file.
1 
2 
3 ! CVS: $Id: gases_ckd.F90,v 1.1 2003/11/07 22:18:38 norm Exp $
4 ! CVS: $Name: $
5 
6 module gases_ckd
7 !Correlated-k coefficients and routines for calculating gaseous absorption per Fu, 1991
8 !The hk* variables are the weights for the k-distribution intervals for each band
9 !The variable fk1o3 is the ozone absorption coefficient in (atm-cm)^-1
10 !The various c2h2o, c3h2o, c10ch4, c10n2o etc. variables are coefficients
11 !required to calculate the absorption coefficient in (atm-cm)^-1.
12 
13  use kinds, only: int_kind, dbl_kind
14  implicit none
15  private
16 
17  public :: &
18  gases & ! Subroutine for calculation gas transmission
19  ,pscale
20 
21 !This include file contains all the correlated-k coefficient
22 !values and associated parameters, (including some public ones)
23 !CVS: $Id: gases_ckd_data.h,v 1.4 2006/11/16 18:43:30 norm Exp $
24 !CVS: $Name: $
25 
26  public ::&
27  mb & ! Total number of spectral intervals
28  ,mbs & ! Number of shortwave spectral intervals
29  ,mbir & ! Number of longwave spectral intervals
30  ,kg & ! Number of k samples per spectral interval
31  ,stanps & ! Reference pressures for shortwave correlated-k tabulation
32  ,stanpir ! Reference pressures for longwave correlated-k tabulation
33 
34  ! Number of bands
35  integer(kind=int_kind), parameter :: &
36  MB = 18 & ! Total number of spectral bands
37  ,MBS = 6 & ! Number of shortwave spectral bands
38  ,MBIR = 12 ! Number of longwave spectral bands
39 
40  ! Number of k samples per band
41  integer(kind=int_kind), parameter, dimension(MB) :: &
42  KG = (/ 10, 8, 12, 7, 12, 5, & ! Number of k samples per shortwave band
43  2, 3, 4, 4, 3, 5, & ! Number of k samples per longwave band
44  2, 10, 12, 7, 7, 8 /)
45 
46  ! Pressure and temperature resolution for k tabulations
47  integer(kind=int_kind), parameter :: &
48  NUMPS = 11 & ! Number of pressures in the shortwave k tabulations
49  ,NUMTS = 3 & ! Number of temperatures in the shortwave k tabulations
50  ,NUMPIR = 19 & ! Number of pressures in the longwave k tabulations
51  ,NUMTIR = 3 & ! Number of temperatures in the longwave k tabulations
52  ,T_1 = 1 & !
53  ,T_2 = 2 & !
54  ,T_3 = 3
55 
56  ! Reference pressures in hPa for correlated-k tabulations
57  real(kind=dbl_kind), parameter, dimension(NUMPS) :: &
58  STANPS = (/ 10.0, 15.8, 25.1, 39.8, 63.1, 100.0, &
59  158.0, 251.0, 398.0, 631.0, 1000.0 /)
60 
61  real(kind=dbl_kind), parameter, dimension(NUMPIR) :: &
62  STANPIR = (/ 0.251, 0.398, 0.631, 1.000, 1.58, 2.51, &
63  3.98, 6.31, 10.0, 15.8, 25.1, 39.8, 63.1, &
64  100.0, 158.0, 251.0, 398.0, 631.0, 1000.0 /)
65 
66 
67  integer (kind=int_kind):: &
68  i,j,k
69 
70  real (kind=dbl_kind), dimension(KG(1)) :: hk1 ! 50000 - 14500 cm^-1
71  real (kind=dbl_kind), dimension(KG(2)) :: hk2 ! 14500 - 7700 cm^-1
72  real (kind=dbl_kind), dimension(KG(3)) :: hk3 ! 7700 - 5250 cm^-1
73  real (kind=dbl_kind), dimension(KG(4)) :: hk4 ! 5250 - 4000 cm^-1
74  real (kind=dbl_kind), dimension(KG(5)) :: hk5 ! 4000 - 2850 cm^-1
75  real (kind=dbl_kind), dimension(KG(6)) :: hk6 ! 2850 - 2500 cm^-1
76  real (kind=dbl_kind), dimension(KG(7)) :: hk7 ! 2200 - 1900 cm^-1
77  real (kind=dbl_kind), dimension(KG(8)) :: hk8 ! 1900 - 1700 cm^-1
78  real (kind=dbl_kind), dimension(KG(9)) :: hk9 ! 1700 - 1400 cm^-1
79  real (kind=dbl_kind), dimension(KG(10)) :: hk10 ! 1400 - 1250 cm^-1
80  real (kind=dbl_kind), dimension(KG(11)) :: hk11 ! 1250 - 1100 cm^-1
81  real (kind=dbl_kind), dimension(KG(12)) :: hk12 ! 1100 - 980 cm^-1
82  real (kind=dbl_kind), dimension(KG(13)) :: hk13 ! 980 - 800 cm^-1
83  real (kind=dbl_kind), dimension(KG(14)) :: hk14 ! 800 - 670 cm^-1
84  real (kind=dbl_kind), dimension(KG(15)) :: hk15 ! 670 - 540 cm^-1
85  real (kind=dbl_kind), dimension(KG(16)) :: hk16 ! 540 - 400 cm^-1
86  real (kind=dbl_kind), dimension(KG(17)) :: hk17 ! 400 - 280 cm^-1
87  real (kind=dbl_kind), dimension(KG(18)) :: hk18 ! 280 - 0 cm^-1
88 
89 
90  real (kind=dbl_kind), dimension(KG(1)) :: fk1o3
91  real (kind=dbl_kind), dimension(NUMTS, NUMPS, KG(2)) :: c2h2o
92  real (kind=dbl_kind), dimension(NUMTS, NUMPS, KG(3)) :: c3h2o
93  real (kind=dbl_kind), dimension(NUMTS, NUMPS, KG(4)) :: c4h2o
94  real (kind=dbl_kind), dimension(NUMTS, NUMPS, KG(5)) :: c5h2o
95  real (kind=dbl_kind), dimension(NUMTS, NUMPS, KG(6)) :: c6h2o
96 
97  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(7)) :: c7h2o
98  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(8)) :: c8h2o
99  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(9)) :: c9h2o
100  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(10)) :: c10h2o
101  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR) :: c10ch4
102  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR) :: c10n2o
103  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(11)) :: c11h2o
104  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR) :: c11ch4
105  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR) :: c11n2o
106  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(12)) :: c12o3
107  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR) :: c12h2o
108  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(13)) :: c13h2o
109  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(14)) :: c14hca
110  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(14)) :: c14hcb
111  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(15)) :: c15hca
112  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(15)) :: c15hcb
113  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(16)) :: c16h2o
114  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(17)) :: c17h2o
115  real (kind=dbl_kind), dimension(NUMTIR, NUMPIR, KG(18)) :: c18h2o
116 
117  data hk1 /.24,.16,.24,.28,.03,.016,.01,.008,.008,.008 /
118  data hk2 /.71,.11,.06,.06,.04,.016,.0034,.0006/
119  data hk3 /.34,.11,.1 ,.09,.12 ,.1 , .06,.04,.026,.01,.0035,.0005/
120  data hk4 /.52,.21,.11,.1,.04,.015,.005/
121  data hk5 /.13,.14,.13 ,.16 ,.18 ,.14 , .07,.02,.016,.008,.004,.002 /
122  data hk6 /.3,.2,.2,.2,.1/
123  data hk7 /0.7,0.3/
124  data hk8 /0.2,0.7,0.1 /
125  data hk9 /0.22,0.51,0.22,0.05/
126  data hk10 /0.28,0.42,0.25,0.05 /
127  data hk11 /0.80,0.15,0.05 /
128  data hk12 /0.45,0.30,0.2,0.04,0.01/
129  data hk13 / 0.95, 0.05 /
130  data hk14 / .3,.3,.2,.12,.06,.012,.004,.0025,.0011,.0004 /
131  data hk15 /.24,.36,.18,.1,.05,.02,.016,.012,.01,.006,.0039,.0021/
132  data hk16 /.12,.24,.24,.20,.12,.06,.02/
133  data hk17 /.12,.26,.22,.20,.10,.085,.015/
134  data hk18 /.10,.1,.2,.25,.2,.1,.03,.02/
135 
136  data fk1o3 /.2204e-08,.1207e-01,.4537e-01,.1032e+00,.1740e+00, &
137  .1210e+01,.7367e+01,.2050e+02,.8100e+02,.2410e+03/
138 
139 
140  data((c2h2o(t_1,j,i),i=1,kg(2)),j=1,numps) / &
141  -0.1735e+02,-0.1407e+02,-0.1268e+02,-0.1131e+02,-0.9261e+01,-0.6666e+01,-0.3937e+01,-0.5448e+00,-0.1690e+02,-0.1365e+02, &
142  -0.1232e+02,-0.1101e+02,-0.9058e+01,-0.6574e+01,-0.3914e+01,-0.5529e+00,-0.1643e+02,-0.1323e+02,-0.1195e+02,-0.1068e+02, &
143  -0.8840e+01,-0.6475e+01,-0.3889e+01,-0.6143e+00,-0.1598e+02,-0.1282e+02,-0.1157e+02,-0.1035e+02,-0.8598e+01,-0.6339e+01, &
144  -0.3848e+01,-0.6636e+00,-0.1551e+02,-0.1241e+02,-0.1119e+02,-0.1001e+02,-0.8342e+01,-0.6178e+01,-0.3788e+01,-0.8181e+00, &
145  -0.1506e+02,-0.1201e+02,-0.1082e+02,-0.9692e+01,-0.8073e+01,-0.6017e+01,-0.3703e+01,-0.9003e+00,-0.1446e+02,-0.1154e+02, &
146  -0.1042e+02,-0.9332e+01,-0.7810e+01,-0.5846e+01,-0.3576e+01,-0.1083e+01,-0.1394e+02,-0.1112e+02,-0.1005e+02,-0.8992e+01, &
147  -0.7548e+01,-0.5674e+01,-0.3477e+01,-0.1266e+01,-0.1351e+02,-0.1076e+02,-0.9722e+01,-0.8702e+01,-0.7334e+01,-0.5531e+01, &
148  -0.3401e+01,-0.1524e+01,-0.1311e+02,-0.1044e+02,-0.9422e+01,-0.8423e+01,-0.7117e+01,-0.5383e+01,-0.3410e+01,-0.1785e+01, &
149  -0.1274e+02,-0.1015e+02,-0.9162e+01,-0.8190e+01,-0.6949e+01,-0.5236e+01,-0.3477e+01,-0.2082e+01 /
150 
151  data((c2h2o(t_2,j,i),i=1,kg(2)),j=1,numps) / &
152  0.2407e-02, 0.2847e-02, 0.3768e-02, 0.4626e-02, 0.5631e-02, 0.4542e-02, 0.3475e-02,-0.3085e-02, 0.2428e-02, 0.2805e-02, &
153  0.3412e-02, 0.3893e-02, 0.4773e-02, 0.3998e-02, 0.2742e-02,-0.2556e-02, 0.2428e-02, 0.2721e-02, 0.3077e-02, 0.3161e-02, &
154  0.4019e-02, 0.3224e-02, 0.2512e-02,-0.1884e-02, 0.2449e-02, 0.2617e-02, 0.2763e-02, 0.2658e-02, 0.3286e-02, 0.2617e-02, &
155  0.1989e-02,-0.1740e-02, 0.2512e-02, 0.2470e-02, 0.2470e-02, 0.2282e-02, 0.2512e-02, 0.1926e-02, 0.1465e-02,-0.2612e-02, &
156  0.2554e-02, 0.2303e-02, 0.2303e-02, 0.1842e-02, 0.2030e-02, 0.1340e-02, 0.1068e-02,-0.1413e-02, 0.2449e-02, 0.2198e-02, &
157  0.2030e-02, 0.1465e-02, 0.1528e-02, 0.9838e-03, 0.1005e-02,-0.1099e-02, 0.2868e-02, 0.2198e-02, 0.1968e-02, 0.1382e-02, &
158  0.1172e-02, 0.5652e-03, 0.6070e-03,-0.1662e-02, 0.3077e-02, 0.2219e-02, 0.1800e-02, 0.1277e-02, 0.1005e-02, 0.3349e-03, &
159  0.2512e-03,-0.1195e-02, 0.3182e-02, 0.2219e-02, 0.1758e-02, 0.1172e-02, 0.7326e-03, 0.4815e-03, 0.6280e-04,-0.1880e-02, &
160  0.3265e-02, 0.2114e-02, 0.1696e-02, 0.1298e-02, 0.4187e-03, 0.4187e-03,-0.3768e-03,-0.1467e-02 /
161 
162  data((c2h2o(t_3,j,i),i=1,kg(2)),j=1,numps) / &
163  -0.1180e-04,-0.1294e-04,-0.1142e-04,-0.7232e-05,-0.8754e-05,-0.1484e-04,-0.8373e-05, 0.1028e-04,-0.1218e-04,-0.1142e-04, &
164  -0.9515e-05,-0.1522e-05,-0.9134e-05,-0.1484e-04,-0.3425e-05, 0.1142e-06,-0.1294e-04,-0.9895e-05,-0.7231e-05,-0.4187e-05, &
165  -0.7612e-05,-0.3806e-05, 0.1522e-05,-0.3882e-05,-0.1256e-04,-0.8754e-05,-0.7612e-05,-0.6470e-05,-0.4948e-05,-0.3425e-05, &
166  0.4948e-05,-0.1054e-04,-0.1370e-04,-0.6089e-05,-0.8373e-05,-0.5709e-05,-0.3045e-05,-0.3806e-05, 0.5328e-05, 0.8678e-05, &
167  -0.1370e-04,-0.6851e-05,-0.8373e-05,-0.1522e-05,-0.3425e-05, 0.0000e+00, 0.1256e-04,-0.1572e-04,-0.1484e-04,-0.7231e-05, &
168  -0.7992e-05,-0.4567e-05,-0.2664e-05,-0.3807e-06,-0.1522e-05, 0.2169e-05,-0.1713e-04,-0.9515e-05,-0.6089e-05,-0.6851e-05, &
169  -0.3045e-05,-0.1142e-05, 0.1903e-05, 0.9363e-05,-0.1560e-04,-0.9134e-05,-0.5328e-05,-0.4948e-05, 0.0000e+00, 0.7611e-06, &
170  -0.6851e-05, 0.1252e-04,-0.1522e-04,-0.8373e-05,-0.6089e-05,-0.6089e-05,-0.3805e-06,-0.1142e-05,-0.3807e-06, 0.2512e-05, &
171  -0.1599e-04,-0.7231e-05,-0.5709e-05,-0.4567e-05, 0.1522e-05,-0.2284e-05,-0.3941e-10, 0.5290e-05 /
172 
173 
174  data((c3h2o(t_1,j,i),i=1,kg(3)),j=1,numps) / &
175  -0.1900e+02,-0.1515e+02,-0.1344e+02,-0.1224e+02,-0.1081e+02,-0.9337e+01,-0.7965e+01,-0.6585e+01,-0.4578e+01,-0.2247e+01, &
176  0.1747e+00, 0.3083e+01,-0.1854e+02,-0.1471e+02,-0.1300e+02,-0.1181e+02,-0.1039e+02,-0.8927e+01,-0.7576e+01,-0.6238e+01, &
177  -0.4317e+01,-0.2119e+01, 0.1888e+00, 0.3033e+01,-0.1808e+02,-0.1426e+02,-0.1257e+02,-0.1137e+02,-0.9966e+01,-0.8513e+01, &
178  -0.7177e+01,-0.5885e+01,-0.4053e+01,-0.1977e+01, 0.2245e+00, 0.3005e+01,-0.1763e+02,-0.1381e+02,-0.1213e+02,-0.1094e+02, &
179  -0.9542e+01,-0.8094e+01,-0.6779e+01,-0.5524e+01,-0.3788e+01,-0.1796e+01, 0.2961e+00, 0.2828e+01,-0.1716e+02,-0.1337e+02, &
180  -0.1170e+02,-0.1051e+02,-0.9116e+01,-0.7677e+01,-0.6381e+01,-0.5153e+01,-0.3493e+01,-0.1607e+01, 0.3850e+00, 0.2660e+01, &
181  -0.1670e+02,-0.1295e+02,-0.1127e+02,-0.1008e+02,-0.8690e+01,-0.7265e+01,-0.5991e+01,-0.4799e+01,-0.3212e+01,-0.1438e+01, &
182  0.4582e+00, 0.2588e+01,-0.1596e+02,-0.1231e+02,-0.1067e+02,-0.9501e+01,-0.8151e+01,-0.6793e+01,-0.5588e+01,-0.4458e+01, &
183  -0.2940e+01,-0.1257e+01, 0.4888e+00, 0.2260e+01,-0.1530e+02,-0.1184e+02,-0.1017e+02,-0.8992e+01,-0.7661e+01,-0.6369e+01, &
184  -0.5213e+01,-0.4145e+01,-0.2701e+01,-0.1108e+01, 0.4239e+00, 0.1974e+01,-0.1481e+02,-0.1144e+02,-0.9756e+01,-0.8573e+01, &
185  -0.7255e+01,-0.5994e+01,-0.4868e+01,-0.3829e+01,-0.2485e+01,-0.9738e+00, 0.3343e+00, 0.1667e+01,-0.1439e+02,-0.1108e+02, &
186  -0.9360e+01,-0.8183e+01,-0.6885e+01,-0.5646e+01,-0.4559e+01,-0.3555e+01,-0.2314e+01,-0.8904e+00, 0.2169e+00, 0.1289e+01, &
187  -0.1402e+02,-0.1073e+02,-0.8987e+01,-0.7817e+01,-0.6551e+01,-0.5335e+01,-0.4278e+01,-0.3316e+01,-0.2147e+01,-0.8695e+00, &
188  0.1587e-01, 0.8658e+00 /
189 
190  data((c3h2o(t_2,j,i),i=1,kg(3)),j=1,numps) / &
191  0.1132e-01, 0.8855e-02, 0.6698e-02, 0.5296e-02, 0.4396e-02, 0.3370e-02, 0.3245e-02, 0.4145e-02, 0.4731e-02, 0.4756e-02, &
192  0.3116e-02,-0.2763e-02, 0.1135e-01, 0.8917e-02, 0.6657e-02, 0.5170e-02, 0.4207e-02, 0.3056e-02, 0.2868e-02, 0.3433e-02, &
193  0.3726e-02, 0.4109e-02, 0.2836e-02,-0.3119e-02, 0.1135e-01, 0.8980e-02, 0.6615e-02, 0.5045e-02, 0.4061e-02, 0.2847e-02, &
194  0.2491e-02, 0.2847e-02, 0.2910e-02, 0.2671e-02, 0.2396e-02,-0.3245e-02, 0.1135e-01, 0.9043e-02, 0.6594e-02, 0.4940e-02, &
195  0.3914e-02, 0.2638e-02, 0.2156e-02, 0.2261e-02, 0.2051e-02, 0.1978e-02, 0.1566e-02,-0.3203e-02, 0.1139e-01, 0.9085e-02, &
196  0.6531e-02, 0.4835e-02, 0.3768e-02, 0.2428e-02, 0.1842e-02, 0.1612e-02, 0.1591e-02, 0.1279e-02, 0.7201e-03,-0.2763e-02, &
197  0.1143e-01, 0.9085e-02, 0.6447e-02, 0.4752e-02, 0.3684e-02, 0.2261e-02, 0.1570e-02, 0.1235e-02, 0.1151e-02, 0.7243e-03, &
198  0.6489e-04,-0.2240e-02, 0.1135e-01, 0.9001e-02, 0.5694e-02, 0.4438e-02, 0.3412e-02, 0.1968e-02, 0.1235e-02, 0.9420e-03, &
199  0.8792e-03, 0.5045e-03,-0.1821e-03,-0.1936e-02, 0.1174e-01, 0.9273e-02, 0.5882e-02, 0.4689e-02, 0.3454e-02, 0.1947e-02, &
200  0.1151e-02, 0.6070e-03, 0.6698e-03, 0.9420e-04,-0.6740e-03,-0.2707e-02, 0.1218e-01, 0.9336e-02, 0.6050e-02, 0.4731e-02, &
201  0.3475e-02, 0.1863e-02, 0.1151e-02, 0.4605e-03, 0.3768e-03,-0.1214e-03,-0.4396e-03,-0.1903e-02, 0.1235e-01, 0.9294e-02, &
202  0.6029e-02, 0.4584e-02, 0.3370e-02, 0.1800e-02, 0.1068e-02, 0.2303e-03, 0.1675e-03,-0.4501e-03,-0.7571e-03,-0.1149e-02, &
203  0.1233e-01, 0.9315e-02, 0.6029e-02, 0.4438e-02, 0.3203e-02, 0.1842e-02, 0.9629e-03, 0.0000e+00,-0.2198e-03,-0.5338e-03, &
204  -0.9721e-03,-0.7661e-03 /
205 
206  data((c3h2o(t_3,j,i),i=1,kg(3)),j=1,numps) / &
207  -0.3692e-04,-0.3844e-04,-0.2588e-04,-0.1180e-04,-0.1066e-04,-0.3426e-05,-0.2664e-05, 0.7611e-06, 0.6089e-05,-0.4568e-06, &
208  -0.2077e-04,-0.1142e-04,-0.3730e-04,-0.3806e-04,-0.2360e-04,-0.1256e-04,-0.1180e-04,-0.4567e-05,-0.3425e-05,-0.2284e-05, &
209  -0.1522e-05,-0.4225e-05,-0.9940e-05,-0.4187e-05,-0.3501e-04,-0.3844e-04,-0.2131e-04,-0.1256e-04,-0.9896e-05,-0.3806e-05, &
210  -0.4186e-05, 0.7612e-06,-0.1903e-05, 0.4110e-05, 0.1789e-05,-0.2169e-04,-0.3425e-04,-0.3882e-04,-0.1941e-04,-0.1294e-04, &
211  -0.9515e-05,-0.4567e-05,-0.4186e-05, 0.1522e-05,-0.4187e-10, 0.4605e-05,-0.2588e-05, 0.6470e-05,-0.3501e-04,-0.3730e-04, &
212  -0.1751e-04,-0.1332e-04,-0.1066e-04,-0.3806e-05,-0.4567e-05,-0.1142e-05,-0.3045e-05, 0.1104e-05,-0.1058e-04, 0.2816e-04, &
213  -0.3578e-04,-0.3501e-04,-0.1751e-04,-0.1332e-04,-0.1218e-04,-0.3806e-05,-0.3425e-05,-0.3806e-06,-0.4187e-05,-0.6090e-06, &
214  -0.6965e-05,-0.3463e-04,-0.3578e-04,-0.3349e-04,-0.1675e-04,-0.9895e-05,-0.9515e-05,-0.6090e-05,-0.6470e-05,-0.3807e-06, &
215  -0.5328e-05,-0.4186e-06,-0.3996e-05, 0.2074e-04,-0.3540e-04,-0.3083e-04,-0.1789e-04,-0.9896e-05,-0.1104e-04,-0.6470e-05, &
216  -0.5709e-05, 0.3425e-05,-0.4567e-05, 0.3463e-05, 0.5633e-05,-0.3159e-05,-0.3730e-04,-0.2740e-04,-0.1484e-04,-0.1066e-04, &
217  -0.1142e-04,-0.6470e-05,-0.6470e-05, 0.1522e-05,-0.1522e-05,-0.3045e-05, 0.3197e-05,-0.1039e-04,-0.3425e-04,-0.2284e-04, &
218  -0.1370e-04,-0.1028e-04,-0.1104e-04,-0.8373e-05,-0.4948e-05, 0.1903e-05,-0.7612e-06,-0.1104e-05, 0.2455e-05,-0.3805e-07, &
219  -0.3235e-04,-0.2093e-04,-0.1294e-04,-0.1142e-04,-0.1180e-04,-0.6851e-05,-0.3045e-05,-0.7611e-06, 0.1256e-05,-0.7231e-06, &
220  0.9924e-05, 0.3578e-05 /
221 
222 
223  data((c4h2o(t_1,j,i),i=1,kg(4)),j=1,numps) / &
224  -0.1722e+02,-0.1402e+02,-0.1202e+02,-0.1001e+02,-0.7702e+01,-0.5273e+01,-0.6530e+00,-0.1677e+02,-0.1359e+02,-0.1164e+02, &
225  -0.9662e+01,-0.7419e+01,-0.5001e+01,-0.6040e+00,-0.1630e+02,-0.1316e+02,-0.1125e+02,-0.9303e+01,-0.7092e+01,-0.4750e+01, &
226  -0.5715e+00,-0.1584e+02,-0.1274e+02,-0.1086e+02,-0.8939e+01,-0.6751e+01,-0.4458e+01,-0.4928e+00,-0.1538e+02,-0.1232e+02, &
227  -0.1048e+02,-0.8579e+01,-0.6399e+01,-0.4191e+01,-0.4683e+00,-0.1493e+02,-0.1192e+02,-0.1011e+02,-0.8241e+01,-0.6065e+01, &
228  -0.3910e+01,-0.4310e+00,-0.1440e+02,-0.1145e+02,-0.9643e+01,-0.7873e+01,-0.5710e+01,-0.3668e+01,-0.3304e+00,-0.1391e+02, &
229  -0.1104e+02,-0.9238e+01,-0.7479e+01,-0.5367e+01,-0.3387e+01,-0.3604e+00,-0.1348e+02,-0.1069e+02,-0.8918e+01,-0.7122e+01, &
230  -0.5086e+01,-0.3152e+01,-0.3030e+00,-0.1310e+02,-0.1037e+02,-0.8626e+01,-0.6790e+01,-0.4815e+01,-0.2945e+01,-0.4789e+00, &
231  -0.1275e+02,-0.1011e+02,-0.8347e+01,-0.6484e+01,-0.4584e+01,-0.2788e+01,-0.5807e+00 /
232 
233  data((c4h2o(t_2,j,i),i=1,kg(4)),j=1,numps) / &
234  0.7934e-02, 0.9231e-02, 0.1005e-01, 0.9043e-02, 0.8164e-02, 0.8980e-02, 0.6403e-02, 0.7954e-02, 0.9169e-02, 0.9797e-02, &
235  0.8687e-02, 0.7724e-02, 0.7954e-02, 0.6652e-02, 0.7954e-02, 0.9043e-02, 0.9608e-02, 0.8499e-02, 0.7347e-02, 0.7473e-02, &
236  0.6382e-02, 0.7996e-02, 0.8980e-02, 0.9378e-02, 0.8289e-02, 0.7264e-02, 0.6594e-02, 0.6674e-02, 0.8059e-02, 0.8938e-02, &
237  0.9294e-02, 0.8227e-02, 0.7201e-02, 0.6678e-02, 0.7032e-02, 0.8122e-02, 0.8896e-02, 0.9189e-02, 0.8038e-02, 0.7033e-02, &
238  0.5987e-02, 0.5475e-02, 0.8268e-02, 0.9064e-02, 0.8792e-02, 0.7975e-02, 0.6573e-02, 0.5087e-02, 0.4657e-02, 0.8541e-02, &
239  0.8980e-02, 0.9085e-02, 0.7996e-02, 0.6133e-02, 0.4501e-02, 0.3860e-02, 0.8813e-02, 0.9043e-02, 0.9294e-02, 0.8122e-02, &
240  0.5861e-02, 0.4354e-02, 0.3964e-02, 0.8875e-02, 0.8834e-02, 0.9797e-02, 0.8164e-02, 0.5463e-02, 0.4417e-02, 0.3270e-02, &
241  0.8938e-02, 0.8771e-02, 0.1005e-01, 0.8247e-02, 0.5589e-02, 0.4835e-02, 0.3033e-02 /
242 
243  data((c4h2o(t_3,j,i),i=1,kg(4)),j=1,numps) / &
244  -0.1484e-04,-0.2169e-04,-0.2436e-04,-0.2588e-04,-0.1142e-04,-0.1142e-05,-0.1519e-04,-0.1522e-04,-0.2055e-04,-0.2131e-04, &
245  -0.2398e-04,-0.4948e-05,-0.1675e-04,-0.3593e-04,-0.1522e-04,-0.2055e-04,-0.1865e-04,-0.2207e-04,-0.4948e-05,-0.1180e-04, &
246  -0.1237e-04,-0.1598e-04,-0.2017e-04,-0.1903e-04,-0.2284e-04,-0.1028e-04,-0.1865e-04,-0.2381e-04,-0.1713e-04,-0.2017e-04, &
247  -0.1827e-04,-0.2169e-04,-0.1218e-04,-0.9515e-05,-0.2415e-04,-0.1827e-04,-0.2093e-04,-0.1637e-04,-0.1827e-04,-0.9134e-05, &
248  -0.8373e-05,-0.1243e-04,-0.1560e-04,-0.1865e-04,-0.1599e-04,-0.1256e-04,-0.1066e-04,-0.1142e-05,-0.2181e-04,-0.1675e-04, &
249  -0.1560e-04,-0.1522e-04,-0.1675e-04,-0.1865e-04,-0.1865e-04,-0.9522e-05,-0.1332e-04,-0.1370e-04,-0.1446e-04,-0.2055e-04, &
250  -0.1142e-04,-0.2512e-04,-0.3343e-04,-0.1294e-04,-0.1294e-04,-0.1751e-04,-0.2512e-04,-0.1560e-04,-0.2854e-04,-0.7003e-05, &
251  -0.8753e-05,-0.1028e-04,-0.1751e-04,-0.2512e-04,-0.1713e-04,-0.1713e-04,-0.1245e-04 /
252 
253 
254  data((c5h2o(t_1,j,i),i=1,kg(5)),j=1,numps) / &
255  -0.1499e+02,-0.1267e+02,-0.1118e+02,-0.9696e+01,-0.7992e+01,-0.6323e+01,-0.4414e+01,-0.2961e+01,-0.1715e+01,-0.1406e+00, &
256  0.1612e+01, 0.3689e+01,-0.1454e+02,-0.1223e+02,-0.1075e+02,-0.9277e+01,-0.7576e+01,-0.5915e+01,-0.4043e+01,-0.2630e+01, &
257  -0.1449e+01, 0.2314e-01, 0.1708e+01, 0.3744e+01,-0.1408e+02,-0.1178e+02,-0.1031e+02,-0.8851e+01,-0.7154e+01,-0.5503e+01, &
258  -0.3666e+01,-0.2288e+01,-0.1141e+01, 0.2772e+00, 0.1819e+01, 0.3788e+01,-0.1363e+02,-0.1134e+02,-0.9876e+01,-0.8423e+01, &
259  -0.6733e+01,-0.5091e+01,-0.3286e+01,-0.1938e+01,-0.8649e+00, 0.5349e+00, 0.1969e+01, 0.3795e+01,-0.1318e+02,-0.1091e+02, &
260  -0.9452e+01,-0.8004e+01,-0.6309e+01,-0.4677e+01,-0.2904e+01,-0.1595e+01,-0.5641e+00, 0.7592e+00, 0.2109e+01, 0.3783e+01, &
261  -0.1275e+02,-0.1048e+02,-0.9028e+01,-0.7585e+01,-0.5892e+01,-0.4267e+01,-0.2524e+01,-0.1274e+01,-0.2782e+00, 0.9376e+00, &
262  0.2257e+01, 0.3714e+01,-0.1180e+02,-0.9887e+01,-0.8492e+01,-0.7014e+01,-0.5390e+01,-0.3834e+01,-0.2156e+01,-0.9775e+00, &
263  -0.3129e-01, 0.1151e+01, 0.2330e+01, 0.3592e+01,-0.1114e+02,-0.9367e+01,-0.8002e+01,-0.6514e+01,-0.4928e+01,-0.3435e+01, &
264  -0.1835e+01,-0.7064e+00, 0.2153e+00, 0.1309e+01, 0.2422e+01, 0.3488e+01,-0.1074e+02,-0.8941e+01,-0.7582e+01,-0.6116e+01, &
265  -0.4536e+01,-0.3072e+01,-0.1521e+01,-0.4651e+00, 0.4053e+00, 0.1465e+01, 0.2374e+01, 0.3260e+01,-0.1041e+02,-0.8545e+01, &
266  -0.7180e+01,-0.5745e+01,-0.4177e+01,-0.2735e+01,-0.1245e+01,-0.2356e+00, 0.5786e+00, 0.1516e+01, 0.2263e+01, 0.3074e+01, &
267  -0.1008e+02,-0.8149e+01,-0.6804e+01,-0.5409e+01,-0.3855e+01,-0.2427e+01,-0.9857e+00,-0.4939e-01, 0.7060e+00, 0.1483e+01, &
268  0.2159e+01, 0.2745e+01 /
269 
270  data((c5h2o(t_2,j,i),i=1,kg(5)),j=1,numps) / &
271  0.9985e-02, 0.8373e-02, 0.7431e-02, 0.6866e-02, 0.4584e-02, 0.2952e-02, 0.3098e-02, 0.3768e-02, 0.4013e-02, 0.3960e-02, &
272  0.3228e-02, 0.3203e-02, 0.1007e-01, 0.8436e-02, 0.7368e-02, 0.6657e-02, 0.4375e-02, 0.2617e-02, 0.2742e-02, 0.3286e-02, &
273  0.3192e-02, 0.2992e-02, 0.2612e-02, 0.1968e-02, 0.1019e-01, 0.8457e-02, 0.7264e-02, 0.6426e-02, 0.4187e-02, 0.2365e-02, &
274  0.2324e-02, 0.2614e-02, 0.2736e-02, 0.2068e-02, 0.2085e-02, 0.1005e-02, 0.1028e-01, 0.8478e-02, 0.7138e-02, 0.6259e-02, &
275  0.3998e-02, 0.2156e-02, 0.1926e-02, 0.1953e-02, 0.2250e-02, 0.1844e-02, 0.1869e-02,-0.6489e-03, 0.1030e-01, 0.8478e-02, &
276  0.7033e-02, 0.6112e-02, 0.3852e-02, 0.1989e-02, 0.1716e-02, 0.1763e-02, 0.1432e-02, 0.1193e-02, 0.1306e-02,-0.5861e-03, &
277  0.1042e-01, 0.8499e-02, 0.6887e-02, 0.5987e-02, 0.3768e-02, 0.1800e-02, 0.1549e-02, 0.1712e-02, 0.1287e-02, 0.7389e-03, &
278  0.7222e-03,-0.1130e-02, 0.8227e-02, 0.7201e-02, 0.6866e-02, 0.5903e-02, 0.3412e-02, 0.1591e-02, 0.1402e-02, 0.1346e-02, &
279  0.1041e-02, 0.8185e-03, 0.3349e-03,-0.4815e-03, 0.8268e-02, 0.6992e-02, 0.7159e-02, 0.6384e-02, 0.3286e-02, 0.1591e-02, &
280  0.1271e-02, 0.1202e-02, 0.9187e-03, 0.6531e-03,-0.4187e-03,-0.7954e-03, 0.8478e-02, 0.7159e-02, 0.7117e-02, 0.6447e-02, &
281  0.3349e-02, 0.1528e-02, 0.9964e-03, 0.9210e-03, 0.6112e-03, 0.6259e-03,-0.3768e-03,-0.1298e-02, 0.8520e-02, 0.7075e-02, &
282  0.7096e-02, 0.6405e-02, 0.3245e-02, 0.1528e-02, 0.1011e-02, 0.7877e-03, 0.7536e-03, 0.9001e-04,-0.6719e-03,-0.1026e-02, &
283  0.8561e-02, 0.6950e-02, 0.7033e-02, 0.6280e-02, 0.2993e-02, 0.1528e-02, 0.6698e-03, 0.5847e-03, 0.2847e-03,-0.6280e-04, &
284  -0.9420e-03,-0.1444e-02 /
285 
286  data((c5h2o(t_3,j,i),i=1,kg(5)),j=1,numps) / &
287  -0.1408e-04,-0.2664e-04,-0.1180e-04,-0.1903e-04,-0.9515e-05, 0.3806e-06,-0.6851e-05,-0.3806e-05,-0.4834e-05,-0.3239e-05, &
288  -0.2284e-05,-0.1028e-04,-0.1484e-04,-0.2550e-04,-0.1142e-04,-0.1827e-04,-0.9515e-05, 0.3805e-06,-0.4948e-05, 0.3806e-06, &
289  -0.2664e-06, 0.1058e-04,-0.1012e-04,-0.1142e-04,-0.1560e-04,-0.2512e-04,-0.1256e-04,-0.1865e-04,-0.9134e-05, 0.1142e-05, &
290  -0.3425e-05, 0.2474e-05,-0.9781e-05,-0.1519e-05,-0.7916e-05,-0.1294e-04,-0.1560e-04,-0.2474e-04,-0.1180e-04,-0.2017e-04, &
291  -0.7992e-05, 0.3805e-06,-0.2283e-05,-0.4453e-05,-0.1180e-05,-0.5138e-05,-0.4453e-05,-0.3425e-05,-0.1522e-04,-0.2550e-04, &
292  -0.9896e-05,-0.1903e-04,-0.9134e-05,-0.1142e-05,-0.7611e-06,-0.5252e-05,-0.4567e-06,-0.4643e-05,-0.4567e-06,-0.4567e-05, &
293  -0.1294e-04,-0.2512e-04,-0.1028e-04,-0.2055e-04,-0.9896e-05,-0.4567e-05,-0.2284e-05,-0.5100e-05,-0.4339e-06,-0.9515e-06, &
294  -0.1252e-04,-0.7612e-06,-0.2246e-04,-0.1370e-04,-0.1066e-04,-0.1598e-04,-0.8754e-05,-0.5328e-05,-0.6622e-05,-0.5138e-05, &
295  -0.8754e-07,-0.9515e-06, 0.6090e-05, 0.4187e-05,-0.3463e-04,-0.1599e-04,-0.1218e-04,-0.2093e-04,-0.9515e-05,-0.4567e-05, &
296  -0.1104e-05,-0.1903e-05,-0.1488e-05,-0.3730e-05,-0.4567e-05, 0.3045e-05,-0.3463e-04,-0.1675e-04,-0.1294e-04,-0.1979e-04, &
297  -0.1066e-04,-0.4187e-05,-0.4034e-05,-0.2893e-05,-0.2588e-05,-0.9401e-05, 0.2284e-05, 0.3045e-05,-0.2778e-04,-0.1522e-04, &
298  -0.1560e-04,-0.1751e-04,-0.1256e-04,-0.5709e-05,-0.2474e-05,-0.2577e-05,-0.2284e-05,-0.4187e-06, 0.7650e-05,-0.3425e-05, &
299  -0.3083e-04,-0.1827e-04,-0.1370e-04,-0.1751e-04,-0.1104e-04,-0.9515e-05,-0.6318e-05,-0.4358e-05,-0.7613e-07, 0.4643e-05, &
300  0.4415e-05, 0.1028e-04 /
301 
302 
303  data((c6h2o(t_1,j,i),i=1,kg(6)),j=1,numps) / &
304  -0.1905e+02,-0.1602e+02,-0.1472e+02,-0.1307e+02,-0.1024e+02,-0.1823e+02,-0.1555e+02,-0.1427e+02,-0.1266e+02,-0.9938e+01, &
305  -0.1749e+02,-0.1508e+02,-0.1381e+02,-0.1225e+02,-0.9641e+01,-0.1684e+02,-0.1462e+02,-0.1337e+02,-0.1185e+02,-0.9367e+01, &
306  -0.1630e+02,-0.1417e+02,-0.1294e+02,-0.1145e+02,-0.9123e+01,-0.1578e+02,-0.1373e+02,-0.1251e+02,-0.1108e+02,-0.8881e+01, &
307  -0.1517e+02,-0.1327e+02,-0.1209e+02,-0.1072e+02,-0.8653e+01,-0.1463e+02,-0.1284e+02,-0.1169e+02,-0.1040e+02,-0.8453e+01, &
308  -0.1421e+02,-0.1244e+02,-0.1133e+02,-0.1014e+02,-0.8312e+01,-0.1382e+02,-0.1207e+02,-0.1100e+02,-0.9887e+01,-0.8220e+01, &
309  -0.1348e+02,-0.1173e+02,-0.1071e+02,-0.9685e+01,-0.8220e+01 /
310 
311  data((c6h2o(t_2,j,i),i=1,kg(6)),j=1,numps) / &
312  0.1024e-01, 0.1842e-02, 0.6908e-03, 0.1737e-02, 0.3517e-02, 0.8394e-02, 0.2072e-02, 0.8164e-03, 0.1716e-02, 0.2805e-02, &
313  0.8143e-02, 0.2240e-02, 0.9001e-03, 0.1570e-02, 0.1800e-02, 0.8227e-02, 0.2386e-02, 0.9420e-03, 0.1486e-02, 0.1068e-02, &
314  0.8373e-02, 0.2533e-02, 0.9210e-03, 0.1319e-02, 0.9420e-03, 0.8394e-02, 0.2700e-02, 0.9629e-03, 0.1026e-02, 0.5233e-03, &
315  0.8917e-02, 0.2575e-02, 0.8792e-03, 0.7536e-03, 0.4187e-03, 0.9378e-02, 0.2617e-02, 0.7955e-03, 0.6070e-03, 0.4815e-03, &
316  0.9797e-02, 0.2638e-02, 0.6908e-03, 0.5233e-03, 0.6280e-03, 0.1009e-01, 0.2638e-02, 0.4815e-03, 0.2931e-03, 0.4815e-03, &
317  0.1036e-01, 0.2428e-02, 0.3140e-03, 0.3977e-03, 0.2093e-03 /
318 
319  data((c6h2o(t_3,j,i),i=1,kg(6)),j=1,numps) / &
320  -0.5366e-04,-0.1522e-04,-0.5709e-05,-0.2664e-05, 0.3806e-05,-0.4301e-04,-0.1484e-04,-0.4948e-05,-0.7610e-06, 0.7610e-06, &
321  -0.3920e-04,-0.1484e-04,-0.4948e-05, 0.3804e-06,-0.3806e-05,-0.3920e-04,-0.1522e-04,-0.4948e-05, 0.3425e-05, 0.1903e-05, &
322  -0.3806e-04,-0.1484e-04,-0.3045e-05, 0.2664e-05, 0.7993e-05,-0.4148e-04,-0.1408e-04,-0.3806e-05, 0.4187e-05, 0.7993e-05, &
323  -0.5481e-04,-0.1180e-04,-0.3045e-05, 0.3045e-05, 0.2284e-05,-0.5709e-04,-0.1104e-04,-0.2283e-05,-0.2664e-05,-0.1142e-05, &
324  -0.6090e-04,-0.1218e-04,-0.2664e-05, 0.3804e-06, 0.3045e-05,-0.6698e-04,-0.1218e-04,-0.2664e-05, 0.1523e-05,-0.1142e-05, &
325  -0.6508e-04,-0.1218e-04,-0.3425e-05, 0.1903e-05, 0.7612e-06 /
326 
327 
328  data((c7h2o(t_1,j,i),i=1,kg(7)),j=1,numpir) / &
329  -0.2008e+02,-0.1467e+02,-0.2004e+02,-0.1426e+02,-0.2001e+02,-0.1386e+02,-0.1998e+02,-0.1345e+02,-0.1995e+02,-0.1304e+02, &
330  -0.1992e+02,-0.1263e+02,-0.1989e+02,-0.1223e+02,-0.1986e+02,-0.1183e+02,-0.1984e+02,-0.1143e+02,-0.1758e+02,-0.1038e+02, &
331  -0.1602e+02,-0.9480e+01,-0.1469e+02,-0.8752e+01,-0.1349e+02,-0.8218e+01,-0.1255e+02,-0.7677e+01,-0.1174e+02,-0.7184e+01, &
332  -0.1110e+02,-0.6735e+01,-0.1056e+02,-0.6332e+01,-0.1019e+02,-0.5975e+01,-0.9874e+01,-0.5644e+01 /
333 
334  data((c7h2o(t_2,j,i),i=1,kg(7)),j=1,numpir) / &
335  0.2533e-02, 0.2269e-01, 0.2575e-02, 0.2263e-01, 0.2554e-02, 0.2267e-01, 0.2491e-02, 0.2250e-01, 0.2449e-02, 0.2244e-01, &
336  0.2344e-02, 0.2234e-01, 0.2219e-02, 0.2208e-01, 0.5694e-02, 0.2190e-01, 0.9650e-02, 0.2162e-01, 0.3286e-01, 0.1848e-01, &
337  0.2987e-01, 0.1578e-01, 0.2527e-01, 0.1465e-01, 0.2175e-01, 0.1386e-01, 0.2056e-01, 0.1235e-01, 0.1963e-01, 0.1116e-01, &
338  0.1926e-01, 0.1040e-01, 0.2014e-01, 0.1040e-01, 0.2024e-01, 0.1042e-01, 0.1972e-01, 0.1080e-01 /
339 
340  data((c7h2o(t_3,j,i),i=1,kg(7)),j=1,numpir) / &
341  -0.8754e-05,-0.6698e-04,-0.1104e-04,-0.6432e-04,-0.1142e-04,-0.6051e-04,-0.1180e-04,-0.6128e-04,-0.1180e-04,-0.6242e-04, &
342  -0.1218e-04,-0.6280e-04,-0.1218e-04,-0.6204e-04, 0.5328e-04,-0.5709e-04, 0.1275e-03,-0.5214e-04,-0.1370e-03,-0.4148e-04, &
343  -0.1100e-03,-0.3045e-04,-0.9248e-04,-0.3197e-04,-0.7346e-04,-0.2436e-04,-0.5100e-04,-0.2131e-04,-0.5861e-04,-0.2550e-04, &
344  -0.5328e-04,-0.3311e-04,-0.6090e-04,-0.4225e-04,-0.5443e-04,-0.4415e-04,-0.4034e-04,-0.4339e-04 /
345 
346 
347  data((c8h2o(t_1,j,i),i=1,kg(8)),j=1,numpir) / &
348  -0.2283e+02,-0.1639e+02,-0.6155e+01,-0.2237e+02,-0.1595e+02,-0.5775e+01,-0.2191e+02,-0.1551e+02,-0.5381e+01,-0.2145e+02, &
349  -0.1507e+02,-0.5004e+01,-0.2099e+02,-0.1463e+02,-0.4617e+01,-0.2053e+02,-0.1419e+02,-0.4218e+01,-0.2025e+02,-0.1375e+02, &
350  -0.3806e+01,-0.2021e+02,-0.1330e+02,-0.3403e+01,-0.2018e+02,-0.1287e+02,-0.2993e+01,-0.1998e+02,-0.1091e+02,-0.2586e+01, &
351  -0.1744e+02,-0.9171e+01,-0.2162e+01,-0.1490e+02,-0.7642e+01,-0.1763e+01,-0.1303e+02,-0.6526e+01,-0.1373e+01,-0.1113e+02, &
352  -0.5846e+01,-0.9699e+00,-0.9814e+01,-0.5280e+01,-0.5955e+00,-0.8582e+01,-0.4787e+01,-0.2510e+00,-0.8020e+01,-0.4350e+01, &
353  0.2770e-01,-0.7571e+01,-0.3942e+01, 0.2406e+00,-0.7140e+01,-0.3537e+01, 0.3567e+00 /
354 
355  data((c8h2o(t_2,j,i),i=1,kg(8)),j=1,numpir) / &
356  0.3722e-01, 0.1505e-01, 0.6615e-02, 0.3722e-01, 0.1518e-01, 0.5840e-02, 0.3720e-01, 0.1526e-01, 0.5170e-02, 0.3399e-01, &
357  0.1530e-01, 0.4773e-02, 0.3012e-01, 0.1551e-01, 0.4333e-02, 0.2625e-01, 0.1553e-01, 0.3956e-02, 0.2240e-01, 0.1562e-01, &
358  0.3454e-02, 0.1846e-01, 0.1574e-01, 0.3161e-02, 0.1446e-01, 0.1572e-01, 0.3098e-02, 0.5924e-02, 0.8875e-02, 0.2658e-02, &
359  0.2204e-01, 0.7096e-02, 0.2504e-02, 0.1591e-01, 0.5233e-02, 0.2292e-02, 0.8855e-02, 0.4249e-02, 0.2190e-02, 0.5422e-02, &
360  0.3496e-02, 0.2041e-02, 0.4919e-02, 0.3621e-02, 0.2200e-02, 0.6657e-02, 0.3663e-02, 0.2248e-02, 0.8645e-02, 0.3852e-02, &
361  0.2118e-02, 0.8771e-02, 0.3873e-02, 0.2176e-02, 0.9043e-02, 0.3747e-02, 0.2079e-02 /
362 
363  data((c8h2o(t_3,j,i),i=1,kg(8)),j=1,numpir) / &
364  -0.1568e-03,-0.4681e-04, 0.4567e-05,-0.1568e-03,-0.4605e-04,-0.3425e-05,-0.1572e-03,-0.4605e-04,-0.1104e-04,-0.2154e-03, &
365  -0.4453e-04,-0.6851e-05,-0.2843e-03,-0.4225e-04,-0.7231e-05,-0.3562e-03,-0.4110e-04,-0.7231e-05,-0.3692e-03,-0.4110e-04, &
366  -0.1028e-04,-0.3007e-03,-0.4263e-04,-0.6470e-05,-0.2325e-03,-0.3996e-04,-0.8373e-05,-0.5290e-04,-0.7612e-05,-0.4948e-05, &
367  -0.7422e-04,-0.1256e-04,-0.8449e-05,-0.3501e-04,-0.1446e-04,-0.4834e-05, 0.4529e-04,-0.2246e-04,-0.2893e-05, 0.6470e-05, &
368  -0.1789e-04,-0.7498e-05,-0.4948e-05,-0.1713e-04,-0.8183e-05,-0.5481e-04,-0.1713e-04,-0.1447e-04,-0.4986e-04,-0.1903e-04, &
369  -0.1353e-04,-0.5138e-04,-0.1484e-04,-0.1147e-04,-0.5328e-04,-0.1560e-04,-0.6588e-05 /
370 
371 
372  data((c9h2o(t_1,j,i),i=1,kg(9)),j=1,numpir) / &
373  -0.2066e+02,-0.1464e+02,-0.8301e+01,-0.3548e+01,-0.2025e+02,-0.1419e+02,-0.7905e+01,-0.3260e+01,-0.2019e+02,-0.1374e+02, &
374  -0.7495e+01,-0.2927e+01,-0.2013e+02,-0.1329e+02,-0.7078e+01,-0.2584e+01,-0.2007e+02,-0.1284e+02,-0.6675e+01,-0.2247e+01, &
375  -0.2001e+02,-0.1239e+02,-0.6268e+01,-0.1890e+01,-0.1996e+02,-0.1194e+02,-0.5853e+01,-0.1530e+01,-0.1991e+02,-0.1150e+02, &
376  -0.5441e+01,-0.1133e+01,-0.1987e+02,-0.1105e+02,-0.5022e+01,-0.7447e+00,-0.1575e+02,-0.9657e+01,-0.4191e+01,-0.3728e+00, &
377  -0.1329e+02,-0.8133e+01,-0.3638e+01, 0.1616e-01,-0.1181e+02,-0.6675e+01,-0.3178e+01, 0.4083e+00,-0.1036e+02,-0.5655e+01, &
378  -0.2731e+01, 0.7953e+00,-0.8628e+01,-0.4990e+01,-0.2303e+01, 0.1153e+01,-0.7223e+01,-0.4453e+01,-0.1877e+01, 0.1454e+01, &
379  -0.6567e+01,-0.3974e+01,-0.1461e+01, 0.1663e+01,-0.6077e+01,-0.3551e+01,-0.1071e+01, 0.1800e+01,-0.5651e+01,-0.3136e+01, &
380  -0.7005e+00, 0.1809e+01,-0.5241e+01,-0.2726e+01,-0.3859e+00, 0.1781e+01 /
381 
382  data((c9h2o(t_2,j,i),i=1,kg(9)),j=1,numpir) / &
383  0.1315e-01, 0.4542e-02, 0.3496e-02, 0.4877e-02, 0.9650e-02, 0.4542e-02, 0.3098e-02, 0.3956e-02, 0.6154e-02, 0.4626e-02, &
384  0.2763e-02, 0.3077e-02, 0.2658e-02, 0.4626e-02, 0.2512e-02, 0.2261e-02, 0.2658e-02, 0.4689e-02, 0.2219e-02, 0.1405e-02, &
385  0.2700e-02, 0.4752e-02, 0.1926e-02, 0.7473e-03, 0.2658e-02, 0.4773e-02, 0.1737e-02, 0.5066e-03, 0.4668e-02, 0.4815e-02, &
386  0.1507e-02, 0.1842e-03, 0.8541e-02, 0.4794e-02, 0.1382e-02,-0.2156e-03, 0.1022e-01, 0.2198e-02, 0.3977e-03,-0.2910e-03, &
387  0.5484e-02, 0.6698e-03, 0.0000e+00,-0.2339e-03, 0.3349e-02, 0.1068e-02,-0.2512e-03,-0.4228e-03, 0.1884e-02, 0.2093e-03, &
388  -0.3977e-03,-0.6405e-03,-0.8373e-04,-0.5233e-03,-0.4124e-03,-0.5945e-03, 0.7536e-03,-0.6698e-03,-0.4919e-03,-0.4794e-03, &
389  0.3600e-02,-0.4605e-03,-0.4375e-03,-0.3517e-03, 0.3873e-02,-0.5861e-03,-0.3203e-03,-0.4689e-03, 0.3935e-02,-0.7326e-03, &
390  -0.2072e-03,-0.4228e-03, 0.4124e-02,-0.8582e-03,-0.4187e-04,-0.5945e-03 /
391 
392  data((c9h2o(t_3,j,i),i=1,kg(9)),j=1,numpir) / &
393  -0.8525e-04, 0.1865e-04,-0.1142e-05, 0.2664e-05,-0.1313e-03, 0.1865e-04, 0.0000e+00, 0.1256e-04,-0.6470e-04, 0.1865e-04, &
394  -0.3045e-05, 0.8754e-05, 0.3805e-06, 0.1789e-04,-0.6851e-05, 0.5328e-05, 0.1142e-05, 0.1827e-04,-0.6090e-05, 0.4148e-05, &
395  0.1142e-05, 0.1865e-04,-0.3806e-05,-0.3768e-05,-0.1903e-05, 0.1751e-04,-0.4948e-05, 0.3121e-05, 0.3159e-04, 0.1979e-04, &
396  -0.3045e-05,-0.9896e-06, 0.1005e-03, 0.1789e-04,-0.6089e-05,-0.1865e-05,-0.2207e-04, 0.1941e-04, 0.1903e-05, 0.2322e-05, &
397  -0.1675e-04, 0.6090e-05,-0.7611e-06, 0.4397e-05, 0.3425e-04, 0.3806e-06, 0.1522e-05, 0.3806e-05, 0.4796e-04, 0.1522e-05, &
398  -0.3806e-06, 0.3654e-05,-0.6851e-05, 0.2664e-05,-0.3920e-05,-0.6850e-06,-0.1370e-04, 0.5328e-05,-0.6584e-05,-0.8716e-05, &
399  -0.8374e-10, 0.1522e-05,-0.6356e-05, 0.1294e-05,-0.9515e-05, 0.7612e-06,-0.3235e-05,-0.1066e-05,-0.7612e-05, 0.1142e-05, &
400  -0.4529e-05, 0.3730e-05,-0.2664e-05,-0.3806e-06,-0.3501e-05,-0.5328e-06 /
401 
402 
403  data((c10h2o(t_1,j,i),i=1,kg(10)),j=1,numpir) / &
404  -0.2023e+02,-0.1641e+02,-0.1171e+02,-0.6090e+01,-0.2016e+02,-0.1595e+02,-0.1133e+02,-0.5867e+01,-0.2011e+02,-0.1550e+02, &
405  -0.1095e+02,-0.5660e+01,-0.2005e+02,-0.1504e+02,-0.1055e+02,-0.5407e+01,-0.2001e+02,-0.1459e+02,-0.1015e+02,-0.5137e+01, &
406  -0.1997e+02,-0.1413e+02,-0.9749e+01,-0.4852e+01,-0.1993e+02,-0.1367e+02,-0.9337e+01,-0.4534e+01,-0.1990e+02,-0.1321e+02, &
407  -0.8920e+01,-0.4211e+01,-0.1987e+02,-0.1276e+02,-0.8506e+01,-0.3889e+01,-0.1645e+02,-0.1179e+02,-0.7711e+01,-0.3613e+01, &
408  -0.1442e+02,-0.1081e+02,-0.6942e+01,-0.3316e+01,-0.1308e+02,-0.9950e+01,-0.6344e+01,-0.2950e+01,-0.1212e+02,-0.9217e+01, &
409  -0.5904e+01,-0.2577e+01,-0.1131e+02,-0.8559e+01,-0.5519e+01,-0.2256e+01,-0.1064e+02,-0.7962e+01,-0.5183e+01,-0.1929e+01, &
410  -0.1013e+02,-0.7447e+01,-0.4833e+01,-0.1643e+01,-0.9712e+01,-0.7071e+01,-0.4485e+01,-0.1410e+01,-0.9305e+01,-0.6760e+01, &
411  -0.4145e+01,-0.1249e+01,-0.8966e+01,-0.6477e+01,-0.3820e+01,-0.1114e+01 /
412 
413  data((c10h2o(t_2,j,i),i=1,kg(10)),j=1,numpir) / &
414  0.7913e-02, 0.8206e-02, 0.1509e-01, 0.1869e-01, 0.4228e-02, 0.8247e-02, 0.1467e-01, 0.1783e-01, 0.2010e-02, 0.8227e-02, &
415  0.1442e-01, 0.1687e-01, 0.1947e-02, 0.8289e-02, 0.1394e-01, 0.1568e-01, 0.1863e-02, 0.8289e-02, 0.1346e-01, 0.1484e-01, &
416  0.1842e-02, 0.8415e-02, 0.1310e-01, 0.1400e-01, 0.1800e-02, 0.8457e-02, 0.1275e-01, 0.1377e-01, 0.1696e-02, 0.8478e-02, &
417  0.1220e-01, 0.1321e-01, 0.1842e-02, 0.8478e-02, 0.1189e-01, 0.1250e-01, 0.1409e-01, 0.8624e-02, 0.1254e-01, 0.1214e-01, &
418  0.9043e-02, 0.1045e-01, 0.1225e-01, 0.1260e-01, 0.8561e-02, 0.1202e-01, 0.1181e-01, 0.1296e-01, 0.1114e-01, 0.1235e-01, &
419  0.1191e-01, 0.1330e-01, 0.1199e-01, 0.1271e-01, 0.1195e-01, 0.1371e-01, 0.1415e-01, 0.1315e-01, 0.1218e-01, 0.1361e-01, &
420  0.1478e-01, 0.1338e-01, 0.1296e-01, 0.1306e-01, 0.1518e-01, 0.1375e-01, 0.1365e-01, 0.1334e-01, 0.1530e-01, 0.1411e-01, &
421  0.1392e-01, 0.1327e-01, 0.1547e-01, 0.1507e-01, 0.1390e-01, 0.1264e-01 /
422 
423  data((c10h2o(t_3,j,i),i=1,kg(10)),j=1,numpir) / &
424  -0.1089e-03,-0.2740e-04,-0.2017e-04,-0.5519e-04,-0.4491e-04,-0.2740e-04,-0.1408e-04,-0.5937e-04,-0.6090e-05,-0.2702e-04, &
425  -0.6470e-05,-0.4719e-04,-0.7232e-05,-0.2740e-04,-0.6089e-05,-0.4910e-04,-0.7231e-05,-0.2969e-04,-0.4186e-05,-0.5366e-04, &
426  -0.6090e-05,-0.3045e-04,-0.2284e-05,-0.4986e-04,-0.4568e-05,-0.3121e-04,-0.4948e-05,-0.5100e-04,-0.3426e-05,-0.3007e-04, &
427  -0.7993e-05,-0.4910e-04, 0.1522e-05,-0.2931e-04,-0.9896e-05,-0.5366e-04,-0.5823e-04,-0.1599e-04,-0.1713e-04,-0.4110e-04, &
428  -0.3121e-04,-0.1713e-04,-0.3159e-04,-0.3578e-04,-0.3996e-04,-0.1598e-04,-0.3958e-04,-0.4605e-04,-0.3349e-04,-0.1751e-04, &
429  -0.3844e-04,-0.5576e-04,-0.2626e-04,-0.2474e-04,-0.3920e-04,-0.4464e-04,-0.1979e-04,-0.3045e-04,-0.3958e-04,-0.5336e-04, &
430  -0.2893e-04,-0.3616e-04,-0.3996e-04,-0.4754e-04,-0.2398e-04,-0.3083e-04,-0.4415e-04,-0.5119e-04,-0.2702e-04,-0.2664e-04, &
431  -0.4605e-04,-0.4038e-04,-0.2398e-04,-0.2360e-04,-0.4948e-04,-0.5149e-04 /
432 
433 
434  data((c10ch4(k,j),j=1,numpir),k=1,numtir) / &
435  -0.8909e+01,-0.8464e+01,-0.8018e+01,-0.7573e+01,-0.7133e+01,-0.6687e+01,-0.6240e+01,-0.5803e+01,-0.5377e+01,-0.4534e+01, &
436  -0.3983e+01,-0.3502e+01,-0.3062e+01,-0.2648e+01,-0.2265e+01,-0.1896e+01,-0.1568e+01,-0.1234e+01,-0.9298e+00, 0.9629e-03, &
437  0.9838e-03, 0.1088e-02, 0.1172e-02, 0.1256e-02, 0.1402e-02, 0.1528e-02, 0.1633e-02, 0.1716e-02, 0.4815e-03,-0.3977e-03, &
438  -0.5652e-03,-0.5024e-03,-0.4605e-03,-0.4563e-03,-0.4438e-03,-0.4521e-03,-0.4312e-03,-0.3789e-03,-0.1294e-04,-0.1408e-04, &
439  -0.1522e-04,-0.1675e-04,-0.1751e-04,-0.1941e-04,-0.2246e-04,-0.2207e-04,-0.1827e-04,-0.1256e-04,-0.9515e-05,-0.6470e-05, &
440  -0.3045e-05,-0.3806e-05,-0.2055e-05,-0.3730e-05,-0.7612e-06,-0.3806e-05, 0.1256e-05 /
441 
442 
443  data((c10n2o(k,j),j=1,numpir),k=1,numtir) / &
444  -0.7863e+01,-0.7412e+01,-0.6963e+01,-0.6514e+01,-0.6065e+01,-0.5611e+01,-0.5167e+01,-0.4720e+01,-0.4283e+01,-0.3454e+01, &
445  -0.2858e+01,-0.2404e+01,-0.1922e+01,-0.1491e+01,-0.1097e+01,-0.7177e+00,-0.3548e+00, 0.1218e-01, 0.3088e+00, 0.4459e-02, &
446  0.4542e-02, 0.4668e-02, 0.4752e-02, 0.4815e-02, 0.4919e-02, 0.5087e-02, 0.5254e-02, 0.5296e-02, 0.2324e-02, 0.2093e-02, &
447  0.2294e-02, 0.2125e-02, 0.2058e-02, 0.1920e-02, 0.1786e-02, 0.1689e-02, 0.1788e-02, 0.2144e-02,-0.7231e-05,-0.7231e-05, &
448  -0.7231e-05,-0.6470e-05,-0.6851e-05,-0.7231e-05,-0.5709e-05,-0.6470e-05,-0.4186e-05, 0.8754e-05,-0.7612e-05,-0.9134e-06, &
449  -0.8640e-05,-0.8487e-05,-0.8259e-05,-0.9553e-05,-0.8107e-05,-0.1654e-04,-0.1858e-04 /
450 
451 
452  data((c11h2o(t_1,j,i),i=1,kg(11)),j=1,numpir) / &
453  -0.2005e+02,-0.1548e+02,-0.1021e+02,-0.2001e+02,-0.1504e+02,-0.1001e+02,-0.1997e+02,-0.1459e+02,-0.9814e+01,-0.1993e+02, &
454  -0.1416e+02,-0.9595e+01,-0.1989e+02,-0.1373e+02,-0.9349e+01,-0.1985e+02,-0.1328e+02,-0.9072e+01,-0.1982e+02,-0.1286e+02, &
455  -0.8833e+01,-0.1957e+02,-0.1243e+02,-0.8566e+01,-0.1911e+02,-0.1200e+02,-0.8276e+01,-0.1743e+02,-0.1134e+02,-0.7958e+01, &
456  -0.1625e+02,-0.1078e+02,-0.7629e+01,-0.1524e+02,-0.1036e+02,-0.7334e+01,-0.1429e+02,-0.9970e+01,-0.7051e+01,-0.1348e+02, &
457  -0.9620e+01,-0.6749e+01,-0.1282e+02,-0.9270e+01,-0.6505e+01,-0.1229e+02,-0.8932e+01,-0.6277e+01,-0.1186e+02,-0.8628e+01, &
458  -0.6120e+01,-0.1148e+02,-0.8345e+01,-0.6049e+01,-0.1112e+02,-0.8066e+01,-0.5906e+01 /
459 
460  data((c11h2o(t_2,j,i),i=1,kg(11)),j=1,numpir) / &
461  0.1842e-02, 0.2131e-01, 0.3033e-01, 0.1905e-02, 0.2137e-01, 0.2841e-01, 0.1926e-02, 0.2135e-01, 0.2696e-01, 0.1926e-02, &
462  0.2133e-01, 0.2514e-01, 0.1884e-02, 0.2154e-01, 0.2401e-01, 0.5589e-02, 0.2156e-01, 0.2321e-01, 0.9483e-02, 0.2156e-01, &
463  0.2210e-01, 0.1333e-01, 0.2150e-01, 0.2133e-01, 0.1725e-01, 0.2154e-01, 0.2074e-01, 0.2254e-01, 0.1999e-01, 0.2005e-01, &
464  0.2118e-01, 0.1926e-01, 0.1978e-01, 0.1936e-01, 0.1920e-01, 0.1963e-01, 0.1905e-01, 0.1911e-01, 0.1934e-01, 0.1909e-01, &
465  0.1903e-01, 0.1920e-01, 0.1922e-01, 0.1901e-01, 0.1899e-01, 0.1934e-01, 0.1930e-01, 0.1974e-01, 0.1966e-01, 0.1909e-01, &
466  0.2014e-01, 0.1976e-01, 0.1905e-01, 0.1984e-01, 0.1963e-01, 0.1940e-01, 0.1897e-01 /
467 
468  data((c11h2o(t_3,j,i),i=1,kg(11)),j=1,numpir) / &
469  -0.1522e-05,-0.6013e-04,-0.5062e-04,-0.2665e-05,-0.6204e-04,-0.5519e-04,-0.3806e-05,-0.6394e-04,-0.5633e-04,-0.4567e-05, &
470  -0.6280e-04,-0.5214e-04,-0.6090e-05,-0.6128e-04,-0.5290e-04, 0.6051e-04,-0.6242e-04,-0.5823e-04, 0.1313e-03,-0.6013e-04, &
471  -0.5176e-04, 0.1336e-03,-0.5747e-04,-0.4072e-04, 0.6318e-04,-0.5671e-04,-0.3996e-04,-0.5595e-04,-0.3996e-04,-0.4263e-04, &
472  -0.3958e-04,-0.4719e-04,-0.4453e-04,-0.3387e-04,-0.5138e-04,-0.5100e-04,-0.5252e-04,-0.4986e-04,-0.4491e-04,-0.5100e-04, &
473  -0.4453e-04,-0.4529e-04,-0.5176e-04,-0.4795e-04,-0.4453e-04,-0.5557e-04,-0.5176e-04,-0.5062e-04,-0.5747e-04,-0.4795e-04, &
474  -0.5633e-04,-0.5709e-04,-0.4643e-04,-0.3806e-04,-0.5481e-04,-0.5671e-04,-0.4948e-04 /
475 
476 
477  data((c11ch4(k,j),j=1,numpir),k=1,numtir) / &
478  -0.1207e+02,-0.1162e+02,-0.1116e+02,-0.1070e+02,-0.1024e+02,-0.9777e+01,-0.9319e+01,-0.8858e+01,-0.8398e+01,-0.7384e+01, &
479  -0.6643e+01,-0.6081e+01,-0.5602e+01,-0.5188e+01,-0.4822e+01,-0.4479e+01,-0.4184e+01,-0.3884e+01,-0.3627e+01, 0.1036e-01, &
480  0.1036e-01, 0.1040e-01, 0.1040e-01, 0.1045e-01, 0.1047e-01, 0.1049e-01, 0.1055e-01, 0.1059e-01, 0.1059e-01, 0.1026e-01, &
481  0.1011e-01, 0.1024e-01, 0.1049e-01, 0.1072e-01, 0.1089e-01, 0.1109e-01, 0.1153e-01, 0.1191e-01,-0.4910e-04,-0.4834e-04, &
482  -0.4910e-04,-0.4910e-04,-0.4910e-04,-0.4872e-04,-0.4834e-04,-0.4948e-04,-0.5100e-04,-0.5633e-04,-0.6166e-04,-0.5595e-04, &
483  -0.5366e-04,-0.5366e-04,-0.5328e-04,-0.5328e-04,-0.4948e-04,-0.5519e-04,-0.5595e-04 /
484 
485 
486  data((c11n2o(k,j),j=1,numpir),k=1,numtir) / &
487  -0.9461e+01,-0.9003e+01,-0.8543e+01,-0.8084e+01,-0.7629e+01,-0.7166e+01,-0.6707e+01,-0.6249e+01,-0.5793e+01,-0.5312e+01, &
488  -0.4847e+01,-0.4393e+01,-0.3974e+01,-0.3587e+01,-0.3231e+01,-0.2885e+01,-0.2602e+01,-0.2358e+01,-0.2108e+01, 0.4710e-02, &
489  0.4752e-02, 0.4773e-02, 0.4773e-02, 0.4815e-02, 0.4877e-02, 0.4898e-02, 0.4982e-02, 0.5066e-02, 0.5296e-02, 0.5149e-02, &
490  0.5129e-02, 0.5024e-02, 0.4752e-02, 0.4501e-02, 0.4270e-02, 0.4019e-02, 0.3646e-02, 0.2759e-02,-0.1484e-04,-0.1408e-04, &
491  -0.1446e-04,-0.1446e-04,-0.1522e-04,-0.1560e-04,-0.1522e-04,-0.1522e-04,-0.1598e-04,-0.1484e-04,-0.9895e-05,-0.1028e-04, &
492  -0.7612e-05,-0.1903e-05, 0.1903e-05, 0.0000e+00, 0.2283e-05, 0.6166e-05,-0.2740e-05 /
493 
494 
495  data((c12h2o(k,j),j=1,numpir),k=1,numtir) / &
496  -0.1984e+02,-0.1983e+02,-0.1982e+02,-0.1981e+02,-0.1963e+02,-0.1917e+02,-0.1871e+02,-0.1825e+02,-0.1779e+02,-0.1639e+02, &
497  -0.1545e+02,-0.1484e+02,-0.1433e+02,-0.1387e+02,-0.1345e+02,-0.1305e+02,-0.1268e+02,-0.1231e+02,-0.1196e+02, 0.6071e-03, &
498  0.2072e-02, 0.6196e-02, 0.1030e-01, 0.1436e-01, 0.1846e-01, 0.2259e-01, 0.2667e-01, 0.2993e-01, 0.2878e-01, 0.2803e-01, &
499  0.2851e-01, 0.2864e-01, 0.2874e-01, 0.2862e-01, 0.2859e-01, 0.2853e-01, 0.2868e-01, 0.2887e-01,-0.3808e-06, 0.2474e-04, &
500  0.9895e-04, 0.1728e-03, 0.1911e-03, 0.1165e-03, 0.4225e-04,-0.3121e-04,-0.8982e-04,-0.9553e-04,-0.9705e-04,-0.9591e-04, &
501  -0.9287e-04,-0.9172e-04,-0.9096e-04,-0.9134e-04,-0.9248e-04,-0.1050e-03,-0.1031e-03 /
502 
503 
504  data((c12o3(t_1,j,i),i=1,kg(12)),j=1,numpir) / &
505  -0.6590e+01,-0.3912e+01,-0.8513e+00, 0.2731e+01, 0.5515e+01,-0.6157e+01,-0.3583e+01,-0.7292e+00, 0.2740e+01, 0.5508e+01, &
506  -0.5731e+01,-0.3242e+01,-0.5800e+00, 0.2782e+01, 0.5485e+01,-0.5301e+01,-0.2901e+01,-0.4131e+00, 0.2805e+01, 0.5455e+01, &
507  -0.4879e+01,-0.2551e+01,-0.2288e+00, 0.2878e+01, 0.5416e+01,-0.4449e+01,-0.2201e+01,-0.2228e-01, 0.3000e+01, 0.5374e+01, &
508  -0.4018e+01,-0.1843e+01, 0.2055e+00, 0.3143e+01, 0.5342e+01,-0.3615e+01,-0.1502e+01, 0.4561e+00, 0.3288e+01, 0.5204e+01, &
509  -0.3228e+01,-0.1172e+01, 0.7099e+00, 0.3396e+01, 0.5077e+01,-0.2828e+01,-0.8499e+00, 0.9664e+00, 0.3463e+01, 0.4893e+01, &
510  -0.2480e+01,-0.5393e+00, 0.1229e+01, 0.3493e+01, 0.4656e+01,-0.2181e+01,-0.2653e+00, 0.1504e+01, 0.3456e+01, 0.4398e+01, &
511  -0.1950e+01,-0.1469e-01, 0.1735e+01, 0.3387e+01, 0.4115e+01,-0.1788e+01, 0.2517e+00, 0.1919e+01, 0.3251e+01, 0.3832e+01, &
512  -0.1677e+01, 0.5027e+00, 0.2032e+01, 0.3088e+01, 0.3581e+01,-0.1637e+01, 0.7373e+00, 0.2100e+01, 0.2910e+01, 0.3364e+01, &
513  -0.1650e+01, 0.9383e+00, 0.2123e+01, 0.2793e+01, 0.3150e+01,-0.1658e+01, 0.1091e+01, 0.2112e+01, 0.2683e+01, 0.3021e+01, &
514  -0.1654e+01, 0.1163e+01, 0.2099e+01, 0.2602e+01, 0.2871e+01 /
515 
516  data((c12o3(t_2,j,i),i=1,kg(12)),j=1,numpir) / &
517  0.9498e-02, 0.8894e-02, 0.1161e-01, 0.8828e-02,-0.1669e-02, 0.9613e-02, 0.8347e-02, 0.1053e-01, 0.8462e-02,-0.1612e-02, &
518  0.9700e-02, 0.7829e-02, 0.9101e-02, 0.7915e-02,-0.1439e-02, 0.9815e-02, 0.7167e-02, 0.7981e-02, 0.7282e-02,-0.1094e-02, &
519  0.9671e-02, 0.6764e-02, 0.6930e-02, 0.5613e-02,-0.8347e-03, 0.9613e-02, 0.6312e-02, 0.6225e-02, 0.4145e-02,-0.1295e-02, &
520  0.9728e-02, 0.6099e-02, 0.5293e-02, 0.2965e-02,-0.1756e-02, 0.9844e-02, 0.5915e-02, 0.4496e-02, 0.1871e-02,-0.2044e-02, &
521  0.9930e-02, 0.5817e-02, 0.3509e-02, 0.1324e-02,-0.2044e-02, 0.9988e-02, 0.5535e-02, 0.2711e-02, 0.6620e-03,-0.1813e-02, &
522  0.1034e-01, 0.5247e-02, 0.1926e-02,-0.2303e-03,-0.1842e-02, 0.1058e-01, 0.4795e-02, 0.1197e-02,-0.9498e-03,-0.2216e-02, &
523  0.1084e-01, 0.4414e-02, 0.6188e-03,-0.1123e-02,-0.2303e-02, 0.1079e-01, 0.3926e-02, 0.1756e-03,-0.1497e-02,-0.2274e-02, &
524  0.1039e-01, 0.3425e-02,-0.1900e-03,-0.1353e-02,-0.2389e-02, 0.9815e-02, 0.2769e-02,-0.6620e-03,-0.1756e-02,-0.1785e-02, &
525  0.9818e-02, 0.2444e-02,-0.1016e-02,-0.1410e-02,-0.1698e-02, 0.1074e-01, 0.3218e-02,-0.1235e-02,-0.1900e-02,-0.2533e-02, &
526  0.1145e-01, 0.3684e-02,-0.1364e-02,-0.1353e-02,-0.1957e-02 /
527 
528  data((c12o3(t_3,j,i),i=1,kg(12)),j=1,numpir) / &
529  -0.4030e-04,-0.2375e-04,-0.3814e-05,-0.4943e-04,-0.3166e-04,-0.3742e-04,-0.1871e-04,-0.1137e-04,-0.4317e-04,-0.2878e-04, &
530  -0.3526e-04,-0.2015e-04,-0.1295e-04,-0.4821e-04,-0.2303e-04,-0.3382e-04,-0.2087e-04,-0.1519e-04,-0.2231e-04,-0.1871e-04, &
531  -0.3454e-04,-0.2087e-04,-0.8109e-05,-0.6476e-05,-0.1511e-04,-0.3454e-04,-0.1820e-04,-0.1269e-05,-0.1439e-04,-0.5037e-05, &
532  -0.4173e-04,-0.2598e-04, 0.6645e-05,-0.1943e-04,-0.2087e-04,-0.3454e-04,-0.2267e-04, 0.2159e-05,-0.2231e-04,-0.2159e-05, &
533  -0.2950e-04,-0.2080e-04, 0.2159e-06,-0.4317e-05, 0.1799e-04,-0.3670e-04,-0.1590e-04,-0.4461e-05,-0.9354e-05,-0.3598e-05, &
534  -0.3216e-04,-0.1475e-04,-0.2231e-05,-0.1295e-04,-0.2878e-05,-0.3576e-04,-0.7347e-05,-0.1022e-04,-0.2159e-05,-0.7915e-05, &
535  -0.3015e-04,-0.5230e-05,-0.5109e-05,-0.6476e-05,-0.7196e-05,-0.2331e-04,-0.1079e-04,-0.4102e-05, 0.1439e-05,-0.1223e-04, &
536  -0.2216e-04,-0.1094e-04,-0.5325e-05,-0.7196e-06,-0.1655e-04,-0.1036e-04,-0.7627e-05,-0.2878e-05, 0.5037e-05,-0.1295e-04, &
537  0.1029e-04,-0.1346e-04,-0.4821e-05,-0.7915e-05, 0.7915e-05, 0.2835e-04,-0.2893e-04,-0.1367e-05,-0.7196e-05,-0.1871e-04, &
538  0.3965e-04,-0.3310e-04,-0.3310e-05,-0.7195e-06, 0.2303e-04 /
539 
540 
541  data((c13h2o(t_1,j,i),i=1,kg(13)),j=1,numpir) / &
542  -0.1992e+02,-0.1446e+02,-0.1992e+02,-0.1405e+02,-0.1991e+02,-0.1363e+02,-0.1990e+02,-0.1322e+02,-0.1989e+02,-0.1282e+02, &
543  -0.1989e+02,-0.1242e+02,-0.1988e+02,-0.1201e+02,-0.1987e+02,-0.1159e+02,-0.1986e+02,-0.1119e+02,-0.1982e+02,-0.1079e+02, &
544  -0.1817e+02,-0.1039e+02,-0.1659e+02,-0.1000e+02,-0.1537e+02,-0.9623e+01,-0.1460e+02,-0.9266e+01,-0.1406e+02,-0.8959e+01, &
545  -0.1354e+02,-0.8676e+01,-0.1309e+02,-0.8411e+01,-0.1267e+02,-0.8232e+01,-0.1229e+02,-0.8094e+01 /
546 
547  data((c13h2o(t_2,j,i),i=1,kg(13)),j=1,numpir) / &
548  0.5024e-03, 0.3199e-01, 0.5652e-03, 0.3199e-01, 0.6071e-03, 0.3211e-01, 0.6489e-03, 0.3199e-01, 0.6699e-03, 0.3178e-01, &
549  0.6908e-03, 0.3157e-01, 0.6908e-03, 0.3109e-01, 0.6698e-03, 0.3075e-01, 0.6698e-03, 0.3054e-01, 0.1474e-01, 0.3000e-01, &
550  0.3085e-01, 0.2960e-01, 0.3659e-01, 0.2935e-01, 0.3016e-01, 0.2920e-01, 0.2834e-01, 0.2895e-01, 0.2780e-01, 0.2870e-01, &
551  0.2753e-01, 0.2843e-01, 0.2755e-01, 0.2820e-01, 0.2765e-01, 0.2732e-01, 0.2769e-01, 0.2705e-01 /
552 
553  data((c13h2o(t_3,j,i),i=1,kg(13)),j=1,numpir) / &
554  0.6299e-09,-0.7993e-04,-0.3802e-06,-0.7992e-04,-0.3802e-06,-0.8525e-04,-0.3808e-06,-0.8449e-04,-0.7610e-06,-0.7764e-04, &
555  -0.1142e-05,-0.7231e-04,-0.1142e-05,-0.7345e-04,-0.2284e-05,-0.8259e-04,-0.2284e-05,-0.8031e-04, 0.2436e-03,-0.7878e-04, &
556  0.7612e-05,-0.8525e-04,-0.1248e-03,-0.9439e-04,-0.9477e-04,-0.9172e-04,-0.8982e-04,-0.8640e-04,-0.7916e-04,-0.6813e-04, &
557  -0.7574e-04,-0.6090e-04,-0.7612e-04,-0.7117e-04,-0.7498e-04,-0.7041e-04,-0.7269e-04,-0.7992e-04 /
558 
559 
560  data((c14hca(t_1,j,i),i=1,kg(14)),j=1,numpir) / &
561  -0.1847e+02,-0.1399e+02,-0.1106e+02,-0.8539e+01,-0.5852e+01,-0.3295e+01,-0.1208e+01,-0.6272e-01, 0.2055e+01, 0.6071e+01, &
562  -0.1801e+02,-0.1357e+02,-0.1067e+02,-0.8171e+01,-0.5562e+01,-0.3071e+01,-0.1073e+01, 0.1033e+00, 0.2055e+01, 0.6071e+01, &
563  -0.1755e+02,-0.1314e+02,-0.1027e+02,-0.7798e+01,-0.5224e+01,-0.2823e+01,-0.9280e+00, 0.2723e+00, 0.2165e+01, 0.5969e+01, &
564  -0.1709e+02,-0.1272e+02,-0.9868e+01,-0.7404e+01,-0.4880e+01,-0.2569e+01,-0.6908e+00, 0.4453e+00, 0.2241e+01, 0.5969e+01, &
565  -0.1663e+02,-0.1230e+02,-0.9467e+01,-0.7013e+01,-0.4535e+01,-0.2297e+01,-0.4408e+00, 0.6353e+00, 0.2359e+01, 0.5969e+01, &
566  -0.1617e+02,-0.1188e+02,-0.9050e+01,-0.6619e+01,-0.4160e+01,-0.1967e+01,-0.1687e+00, 0.8213e+00, 0.2421e+01, 0.5969e+01, &
567  -0.1571e+02,-0.1147e+02,-0.8629e+01,-0.6230e+01,-0.3771e+01,-0.1648e+01, 0.1573e+00, 0.1019e+01, 0.2511e+01, 0.5884e+01, &
568  -0.1525e+02,-0.1106e+02,-0.8215e+01,-0.5841e+01,-0.3393e+01,-0.1331e+01, 0.4013e+00, 0.1198e+01, 0.2654e+01, 0.5794e+01, &
569  -0.1480e+02,-0.1066e+02,-0.7800e+01,-0.5454e+01,-0.3032e+01,-0.9870e+00, 0.6323e+00, 0.1373e+01, 0.2905e+01, 0.5647e+01, &
570  -0.1402e+02,-0.9693e+01,-0.7206e+01,-0.4846e+01,-0.2656e+01,-0.6540e+00, 0.8323e+00, 0.1530e+01, 0.3211e+01, 0.5355e+01, &
571  -0.1343e+02,-0.9060e+01,-0.6596e+01,-0.4399e+01,-0.2294e+01,-0.3519e+00, 0.9823e+00, 0.1673e+01, 0.3420e+01, 0.5083e+01, &
572  -0.1279e+02,-0.8611e+01,-0.5785e+01,-0.4010e+01,-0.1936e+01,-0.1177e+00, 0.1134e+01, 0.1974e+01, 0.3591e+01, 0.4770e+01, &
573  -0.1230e+02,-0.8174e+01,-0.5298e+01,-0.3611e+01,-0.1607e+01, 0.3636e-01, 0.1433e+01, 0.2260e+01, 0.3539e+01, 0.4439e+01, &
574  -0.1192e+02,-0.7763e+01,-0.4946e+01,-0.3228e+01,-0.1321e+01, 0.1991e+00, 0.1720e+01, 0.2420e+01, 0.3383e+01, 0.4041e+01, &
575  -0.1154e+02,-0.7377e+01,-0.4576e+01,-0.2851e+01,-0.1093e+01, 0.4430e+00, 0.1896e+01, 0.2462e+01, 0.3122e+01, 0.3620e+01, &
576  -0.1118e+02,-0.7003e+01,-0.4210e+01,-0.2524e+01,-0.8973e+00, 0.7490e+00, 0.1966e+01, 0.2363e+01, 0.2818e+01, 0.3182e+01, &
577  -0.1080e+02,-0.6677e+01,-0.3872e+01,-0.2264e+01,-0.6846e+00, 0.9392e+00, 0.1867e+01, 0.2138e+01, 0.2505e+01, 0.2738e+01, &
578  -0.1031e+02,-0.6353e+01,-0.3596e+01,-0.1938e+01,-0.4537e+00, 0.1015e+01, 0.1659e+01, 0.1830e+01, 0.2142e+01, 0.2287e+01, &
579  -0.9695e+01,-0.5977e+01,-0.3427e+01,-0.1596e+01,-0.1979e+00, 0.9458e+00, 0.1363e+01, 0.1545e+01, 0.1743e+01, 0.1832e+01 /
580 
581  data((c14hca(t_2,j,i),i=1,kg(14)),j=1,numpir) / &
582  0.3628e-01, 0.2728e-01, 0.2213e-01, 0.1656e-01, 0.1507e-01, 0.1564e-01, 0.1623e-01, 0.1419e-01, 0.1455e-01, 0.1089e-02, &
583  0.3632e-01, 0.2740e-01, 0.2164e-01, 0.1606e-01, 0.1369e-01, 0.1418e-01, 0.1444e-01, 0.1275e-01, 0.1331e-01, 0.9210e-03, &
584  0.3636e-01, 0.2746e-01, 0.2114e-01, 0.1557e-01, 0.1239e-01, 0.1285e-01, 0.1237e-01, 0.1141e-01, 0.1141e-01, 0.9210e-03, &
585  0.3640e-01, 0.2748e-01, 0.2064e-01, 0.1516e-01, 0.1141e-01, 0.1125e-01, 0.1092e-01, 0.1026e-01, 0.1011e-01,-0.5652e-03, &
586  0.3646e-01, 0.2746e-01, 0.2024e-01, 0.1478e-01, 0.1036e-01, 0.9688e-02, 0.9610e-02, 0.9305e-02, 0.9399e-02,-0.6489e-03, &
587  0.3651e-01, 0.2734e-01, 0.1984e-01, 0.1438e-01, 0.9436e-02, 0.8486e-02, 0.8214e-02, 0.8995e-02, 0.7892e-02,-0.8582e-03, &
588  0.3655e-01, 0.2723e-01, 0.1951e-01, 0.1402e-01, 0.8716e-02, 0.7433e-02, 0.7169e-02, 0.8072e-02, 0.5443e-02,-0.1172e-02, &
589  0.3659e-01, 0.2709e-01, 0.1911e-01, 0.1379e-01, 0.8107e-02, 0.6818e-02, 0.6818e-02, 0.7033e-02, 0.3056e-02,-0.1047e-02, &
590  0.3670e-01, 0.2698e-01, 0.1890e-01, 0.1363e-01, 0.7502e-02, 0.6371e-02, 0.6558e-02, 0.6489e-02,-0.5652e-03,-0.1340e-02, &
591  0.3592e-01, 0.2238e-01, 0.1804e-01, 0.1007e-01, 0.6730e-02, 0.5512e-02, 0.6194e-02, 0.4375e-02,-0.1109e-02,-0.3559e-03, &
592  0.3609e-01, 0.2242e-01, 0.1526e-01, 0.8582e-02, 0.6284e-02, 0.5809e-02, 0.4501e-02, 0.9420e-03,-0.9001e-03,-0.1005e-02, &
593  0.3703e-01, 0.2196e-01, 0.1281e-01, 0.7860e-02, 0.5861e-02, 0.5842e-02, 0.1800e-02,-0.1591e-02,-0.1235e-02,-0.9420e-03, &
594  0.3728e-01, 0.2114e-01, 0.1347e-01, 0.6678e-02, 0.5449e-02, 0.4837e-02,-0.1084e-02,-0.1361e-02,-0.6699e-03,-0.1256e-03, &
595  0.3683e-01, 0.2061e-01, 0.1350e-01, 0.6133e-02, 0.5449e-02, 0.2111e-02,-0.1386e-02,-0.1235e-02,-0.5652e-03,-0.8373e-04, &
596  0.3656e-01, 0.1988e-01, 0.1348e-01, 0.5441e-02, 0.5149e-02,-0.8813e-03,-0.1116e-02,-0.8373e-03,-0.3140e-03,-0.6280e-04, &
597  0.3669e-01, 0.1934e-01, 0.1363e-01, 0.5035e-02, 0.3585e-02,-0.1250e-02,-0.9357e-03,-0.8227e-03,-0.3140e-03,-0.4187e-04, &
598  0.3618e-01, 0.1856e-01, 0.1390e-01, 0.3836e-02, 0.1470e-02,-0.1096e-02,-0.8080e-03,-0.4480e-03,-0.2093e-03,-0.2093e-04, &
599  0.3416e-01, 0.1741e-01, 0.1431e-01, 0.1951e-02,-0.2923e-04,-0.9422e-03,-0.4576e-03,-0.2395e-03,-0.1565e-03,-0.2799e-04, &
600  0.3219e-01, 0.1674e-01, 0.1516e-01, 0.6652e-03,-0.5051e-03,-0.7052e-03,-0.2002e-03,-0.2135e-03,-0.7633e-04,-0.7300e-04 /
601 
602  data((c14hca(t_3,j,i),i=1,kg(14)),j=1,numpir) / &
603  -0.1290e-03,-0.9934e-04,-0.5595e-04,-0.3996e-04, 0.1294e-04,-0.9134e-05, 0.1294e-05,-0.3121e-05,-0.4757e-04,-0.1979e-04, &
604  -0.1305e-03,-0.9629e-04,-0.5481e-04,-0.4301e-04, 0.1827e-04,-0.9363e-05, 0.1777e-04,-0.2185e-04,-0.1903e-04,-0.1675e-04, &
605  -0.1313e-03,-0.9439e-04,-0.5404e-04,-0.4263e-04, 0.9134e-05,-0.1020e-04, 0.3524e-04,-0.2599e-04,-0.2093e-04, 0.1675e-04, &
606  -0.1313e-03,-0.9172e-04,-0.5252e-04,-0.4567e-04, 0.4186e-05,-0.3920e-05, 0.2552e-04,-0.2059e-04,-0.2246e-04,-0.1028e-04, &
607  -0.1324e-03,-0.9210e-04,-0.5138e-04,-0.4491e-04, 0.6470e-05,-0.2131e-05, 0.1496e-04,-0.1572e-04,-0.3311e-04,-0.8754e-05, &
608  -0.1324e-03,-0.9058e-04,-0.5328e-04,-0.4225e-04, 0.1827e-05,-0.8411e-06, 0.4719e-05,-0.6813e-05,-0.2474e-04,-0.1256e-04, &
609  -0.1340e-03,-0.8868e-04,-0.5633e-04,-0.4187e-04,-0.4415e-05, 0.6055e-05,-0.1648e-04,-0.1507e-04, 0.1979e-04,-0.2131e-04, &
610  -0.1340e-03,-0.8373e-04,-0.5899e-04,-0.3920e-04,-0.4072e-05, 0.1491e-04,-0.9781e-05,-0.5328e-05, 0.3578e-04,-0.1979e-04, &
611  -0.1321e-03,-0.7954e-04,-0.5899e-04,-0.4072e-04, 0.1066e-05, 0.5728e-05,-0.5138e-05,-0.8373e-05, 0.2626e-04,-0.2436e-04, &
612  -0.1363e-03,-0.6432e-04,-0.5176e-04,-0.3083e-04, 0.2169e-05,-0.8944e-05, 0.3159e-05, 0.6470e-05,-0.4187e-05, 0.4948e-05, &
613  -0.1302e-03,-0.7802e-04,-0.3311e-04,-0.1903e-04, 0.5328e-05,-0.1884e-04, 0.1408e-04, 0.3311e-04, 0.1142e-05,-0.7613e-06, &
614  -0.1473e-03,-0.6737e-04,-0.7536e-04,-0.1085e-04,-0.1903e-05,-0.1458e-04, 0.4034e-04,-0.3941e-10,-0.7992e-05, 0.2664e-05, &
615  -0.1361e-03,-0.5709e-04,-0.8550e-04,-0.5709e-05,-0.8640e-05, 0.6523e-05, 0.1903e-05,-0.8221e-05,-0.3045e-05,-0.9134e-05, &
616  -0.1329e-03,-0.5529e-04,-0.7107e-04, 0.2664e-05,-0.9020e-05, 0.3320e-04,-0.2131e-05,-0.4187e-05,-0.7231e-05,-0.3806e-05, &
617  -0.1278e-03,-0.5247e-04,-0.6465e-04, 0.3806e-05,-0.6091e-05, 0.1245e-04,-0.3844e-05,-0.6090e-05,-0.8754e-05,-0.2664e-05, &
618  -0.1321e-03,-0.5632e-04,-0.5897e-04, 0.1012e-04, 0.1168e-04,-0.4196e-06,-0.8411e-05,-0.8868e-05,-0.1484e-04,-0.1522e-05, &
619  -0.1252e-03,-0.4907e-04,-0.5932e-04, 0.3245e-04, 0.1996e-04,-0.3325e-05,-0.5785e-05,-0.6394e-05,-0.6851e-05,-0.1142e-05, &
620  -0.1093e-03,-0.4731e-04,-0.6761e-04, 0.1808e-04, 0.1754e-04,-0.5079e-05,-0.5809e-05,-0.5649e-05,-0.3988e-05,-0.5849e-06, &
621  -0.1151e-03,-0.4965e-04,-0.7163e-04, 0.7839e-05, 0.5505e-05,-0.6084e-05,-0.3344e-05,-0.3894e-05,-0.1391e-05,-0.1327e-05 /
622 
623 
624  data((c14hcb(t_1,j,i),i=1,kg(14)),j=1,numpir) / &
625  -0.9398e+01,-0.5678e+01,-0.3606e+01,-0.2192e+01, 0.2104e+01, 0.3044e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
626  -0.9094e+01,-0.5422e+01,-0.3448e+01,-0.1650e+01, 0.2046e+01, 0.2749e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
627  -0.8760e+01,-0.5270e+01,-0.3329e+01,-0.1147e+01, 0.2112e+01, 0.2709e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
628  -0.8537e+01,-0.5152e+01,-0.3129e+01,-0.9544e+00, 0.2254e+01, 0.2771e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
629  -0.8176e+01,-0.4936e+01,-0.2680e+01,-0.9259e+00, 0.2247e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
630  -0.7836e+01,-0.4676e+01,-0.2378e+01,-0.3550e+00, 0.1396e+01, 0.1976e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
631  -0.7419e+01,-0.4122e+01,-0.2407e+01,-0.1204e-01, 0.1744e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
632  -0.7124e+01,-0.3727e+01,-0.2160e+01, 0.6158e+00, 0.1953e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
633  -0.6823e+01,-0.3324e+01,-0.1748e+01,-0.9806e-01, 0.2319e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
634  -0.5957e+01,-0.3017e+01,-0.1647e+01, 0.1398e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
635  -0.5115e+01,-0.2290e+01,-0.5273e+00, 0.5662e+00, 0.1459e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
636  -0.4162e+01,-0.1453e+01, 0.1116e+00,-0.4587e+02, 0.9569e+00,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
637  -0.3611e+01,-0.9744e+00,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
638  -0.3075e+01,-0.4176e+00,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
639  -0.3469e+01,-0.9395e+00, 0.5092e+00, 0.6200e+00,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
640  -0.3808e+01,-0.1505e+01, 0.3901e+00, 0.6264e+00,-0.1155e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
641  -0.4058e+01,-0.1818e+01, 0.2693e+00, 0.7087e+00, 0.3820e+00,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
642  -0.4262e+01,-0.2097e+01,-0.5711e-01, 0.5681e+00, 0.1310e+01, 0.7371e+00,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
643  -0.3997e+01,-0.1784e+01, 0.4388e-01, 0.5167e+00, 0.6930e+00,-0.6906e+00,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02 /
644 
645  data((c14hcb(t_2,j,i),i=1,kg(14)),j=1,numpir) / &
646  0.2944e-01, 0.2723e-01, 0.1854e-01, 0.2023e-01, 0.2254e-01, 0.3059e-02, 0.4788e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
647  0.3080e-01, 0.2549e-01, 0.1547e-01, 0.2225e-01, 0.2107e-01, 0.3059e-02, 0.4737e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
648  0.3269e-01, 0.2656e-01, 0.2125e-01, 0.2179e-01, 0.2162e-01, 0.4589e+00, 0.4643e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
649  0.3322e-01, 0.2476e-01, 0.2075e-01, 0.2139e-01, 0.1907e-01, 0.4501e+00, 0.4441e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
650  0.3387e-01, 0.2182e-01, 0.2665e-01, 0.1841e-01, 0.2506e-01, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
651  0.3532e-01, 0.2091e-01, 0.1995e-01, 0.2067e-01, 0.1949e-01, 0.4491e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
652  0.3468e-01, 0.2075e-01, 0.2587e-01, 0.1401e-01, 0.8646e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
653  0.3666e-01, 0.2430e-01, 0.1919e-01, 0.2007e-01, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
654  0.3613e-01, 0.2147e-01, 0.1892e-01, 0.1361e-01, 0.3059e-02, 0.4506e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
655  0.3129e-01, 0.1954e-01, 0.2442e-01, 0.1011e-01, 0.4420e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
656  0.3177e-01, 0.2101e-01, 0.1526e-01, 0.4376e+00, 0.4379e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
657  0.2887e-01, 0.2044e-01, 0.1285e-01, 0.3059e-02,-0.4862e-03, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
658  0.2759e-01, 0.2114e-01, 0.4303e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
659  0.2880e-01, 0.1690e-01,-0.4187e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
660  0.2852e-01, 0.2255e-01, 0.2184e-01, 0.4334e+00, 0.4217e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
661  0.2840e-01, 0.2136e-01, 0.1644e-01, 0.2812e-01, 0.4358e+00, 0.4288e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
662  0.2809e-01, 0.2173e-01, 0.1708e-01, 0.3346e-01, 0.4225e-01, 0.4419e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
663  0.2702e-01, 0.2260e-01, 0.1607e-01, 0.2720e-01, 0.3982e-01, 0.4452e+00, 0.4365e+00, 0.4345e+00, 0.4432e+00, 0.4623e+00, &
664  0.2684e-01, 0.2328e-01, 0.2099e-01, 0.3040e-01, 0.3867e-01, 0.4389e+00, 0.3132e-01, 0.3158e-01, 0.4083e-01, 0.4580e+00 /
665 
666  data((c14hcb(t_3,j,i),i=1,kg(14)),j=1,numpir) / &
667  -0.1581e-03,-0.9707e-04,-0.1250e-03, 0.2580e-03, 0.7378e-04,-0.1617e-01, 0.8646e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
668  -0.1319e-03,-0.9528e-04,-0.1710e-03, 0.7118e-04, 0.2076e-04,-0.1608e-01, 0.8552e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
669  -0.1721e-03,-0.4680e-04,-0.5522e-04,-0.6242e-04, 0.4517e-04,-0.7777e-02, 0.8382e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
670  -0.1482e-03,-0.4208e-04,-0.5216e-04,-0.6514e-04,-0.8378e-04,-0.7956e-02, 0.8013e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
671  -0.1501e-03,-0.4002e-04,-0.1664e-03, 0.2272e-04,-0.1888e-03,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
672  -0.1201e-03,-0.4709e-04,-0.5371e-04,-0.1574e-03, 0.1854e-03,-0.7712e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
673  -0.1333e-03,-0.1062e-03, 0.5785e-04,-0.4150e-04,-0.5717e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
674  -0.1212e-03,-0.8524e-04,-0.5895e-04,-0.2884e-03,-0.1581e-01,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
675  -0.8148e-04,-0.9361e-04,-0.2873e-03, 0.1883e-03,-0.1594e-01, 0.8133e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
676  -0.1221e-03,-0.1430e-04, 0.6335e-04,-0.2581e-03, 0.7977e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
677  -0.9257e-04,-0.5008e-04, 0.6389e-04,-0.7455e-02,-0.7745e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
678  -0.1186e-03,-0.9037e-04,-0.7461e-04,-0.4656e-05, 0.1168e-03,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
679  -0.8513e-04,-0.5708e-04, 0.7763e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
680  -0.1124e-03,-0.1228e-03, 0.7663e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
681  -0.1015e-03,-0.8369e-04,-0.2167e-03,-0.7548e-02, 0.7608e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
682  -0.1049e-03,-0.6414e-04,-0.1384e-03,-0.1644e-03,-0.6919e-02, 0.7736e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
683  -0.1008e-03,-0.7047e-04,-0.1276e-03,-0.2445e-03,-0.1860e-03, 0.7975e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
684  -0.9629e-04,-0.1007e-03,-0.1127e-03,-0.1527e-03,-0.3238e-03,-0.7373e-02, 0.7877e-02, 0.7840e-02, 0.7997e-02, 0.8345e-02, &
685  -0.8800e-04,-0.1072e-03,-0.1046e-03,-0.1777e-03,-0.2146e-03,-0.7016e-02, 0.1516e-01, 0.1532e-01, 0.1509e-01, 0.8268e-02 /
686 
687 
688  data((c15hca(t_1,j,i),i=1,kg(15)),j=1,numpir) / &
689  -0.1921e+02,-0.1363e+02,-0.1080e+02,-0.8392e+01,-0.6776e+01,-0.5696e+01,-0.4572e+01,-0.3752e+01,-0.2382e+01,-0.1110e+01, &
690  0.6803e+00, 0.3259e+01,-0.1875e+02,-0.1321e+02,-0.1040e+02,-0.8026e+01,-0.6449e+01,-0.5401e+01,-0.4316e+01,-0.3498e+01, &
691  -0.2141e+01,-0.9439e+00, 0.8103e+00, 0.3314e+01,-0.1829e+02,-0.1278e+02,-0.1000e+02,-0.7646e+01,-0.6089e+01,-0.5085e+01, &
692  -0.4047e+01,-0.3217e+01,-0.1872e+01,-0.7106e+00, 0.9573e+00, 0.3390e+01,-0.1783e+02,-0.1236e+02,-0.9596e+01,-0.7264e+01, &
693  -0.5735e+01,-0.4740e+01,-0.3743e+01,-0.2882e+01,-0.1587e+01,-0.4714e+00, 0.1120e+01, 0.3425e+01,-0.1737e+02,-0.1195e+02, &
694  -0.9193e+01,-0.6877e+01,-0.5371e+01,-0.4404e+01,-0.3405e+01,-0.2574e+01,-0.1298e+01,-0.1747e+00, 0.1327e+01, 0.3547e+01, &
695  -0.1691e+02,-0.1153e+02,-0.8776e+01,-0.6490e+01,-0.4993e+01,-0.4049e+01,-0.3039e+01,-0.2256e+01,-0.1012e+01, 0.1103e+00, &
696  0.1530e+01, 0.3651e+01,-0.1644e+02,-0.1112e+02,-0.8360e+01,-0.6105e+01,-0.4623e+01,-0.3688e+01,-0.2694e+01,-0.1915e+01, &
697  -0.6855e+00, 0.3993e+00, 0.1714e+01, 0.3950e+01,-0.1598e+02,-0.1073e+02,-0.7943e+01,-0.5723e+01,-0.4236e+01,-0.3314e+01, &
698  -0.2338e+01,-0.1596e+01,-0.3583e+00, 0.6963e+00, 0.1868e+01, 0.4127e+01,-0.1553e+02,-0.1034e+02,-0.7542e+01,-0.5357e+01, &
699  -0.3856e+01,-0.2942e+01,-0.1986e+01,-0.1299e+01,-0.5472e-01, 0.9443e+00, 0.2149e+01, 0.4261e+01,-0.1485e+02,-0.9661e+01, &
700  -0.7008e+01,-0.4830e+01,-0.3458e+01,-0.2566e+01,-0.1658e+01,-0.9639e+00, 0.2083e+00, 0.1182e+01, 0.2458e+01, 0.4452e+01, &
701  -0.1427e+02,-0.9166e+01,-0.6373e+01,-0.4404e+01,-0.3073e+01,-0.2209e+01,-0.1349e+01,-0.6648e+00, 0.4023e+00, 0.1452e+01, &
702  0.2739e+01, 0.4466e+01,-0.1380e+02,-0.8726e+01,-0.5772e+01,-0.3982e+01,-0.2732e+01,-0.1874e+01,-0.1052e+01,-0.4403e+00, &
703  0.5763e+00, 0.1792e+01, 0.2999e+01, 0.4335e+01,-0.1305e+02,-0.8270e+01,-0.5304e+01,-0.3586e+01,-0.2392e+01,-0.1568e+01, &
704  -0.8299e+00,-0.2650e+00, 0.8584e+00, 0.2062e+01, 0.3141e+01, 0.4168e+01,-0.1269e+02,-0.7900e+01,-0.4956e+01,-0.3205e+01, &
705  -0.2065e+01,-0.1332e+01,-0.6415e+00,-0.7921e-01, 0.1170e+01, 0.2269e+01, 0.3198e+01, 0.4066e+01,-0.1227e+02,-0.7536e+01, &
706  -0.4576e+01,-0.2859e+01,-0.1815e+01,-0.1139e+01,-0.4520e+00, 0.2272e+00, 0.1371e+01, 0.2351e+01, 0.3150e+01, 0.3935e+01, &
707  -0.1186e+02,-0.7159e+01,-0.4223e+01,-0.2538e+01,-0.1619e+01,-0.9324e+00,-0.1566e+00, 0.5151e+00, 0.1520e+01, 0.2339e+01, &
708  0.3132e+01, 0.3880e+01,-0.1120e+02,-0.6777e+01,-0.3919e+01,-0.2330e+01,-0.1387e+01,-0.6737e+00, 0.1108e+00, 0.6991e+00, &
709  0.1531e+01, 0.2163e+01, 0.3150e+01, 0.3767e+01,-0.9973e+01,-0.6279e+01,-0.3638e+01,-0.2048e+01,-0.1098e+01,-0.4407e+00, &
710  0.3043e+00, 0.7797e+00, 0.1424e+01, 0.2002e+01, 0.3122e+01, 0.3611e+01,-0.8483e+01,-0.5607e+01,-0.3357e+01,-0.1744e+01, &
711  -0.8884e+00,-0.2264e+00, 0.3800e+00, 0.7504e+00, 0.1245e+01, 0.2032e+01, 0.3097e+01, 0.3546e+01 /
712 
713  data((c15hca(t_2,j,i),i=1,kg(15)),j=1,numpir) / &
714  0.3762e-01, 0.2372e-01, 0.1643e-01, 0.1208e-01, 0.1170e-01, 0.1164e-01, 0.1214e-01, 0.1161e-01, 0.1028e-01, 0.9185e-02, &
715  0.7712e-02, 0.1001e-01, 0.3762e-01, 0.2382e-01, 0.1593e-01, 0.1145e-01, 0.1059e-01, 0.1049e-01, 0.1080e-01, 0.1057e-01, &
716  0.8894e-02, 0.7807e-02, 0.7132e-02, 0.1032e-01, 0.3764e-01, 0.2386e-01, 0.1555e-01, 0.1080e-01, 0.9692e-02, 0.9231e-02, &
717  0.9585e-02, 0.9644e-02, 0.7711e-02, 0.6443e-02, 0.6223e-02, 0.9922e-02, 0.3764e-01, 0.2395e-01, 0.1516e-01, 0.1028e-01, &
718  0.8917e-02, 0.8415e-02, 0.8457e-02, 0.8777e-02, 0.6436e-02, 0.5428e-02, 0.5499e-02, 0.8017e-02, 0.3768e-01, 0.2399e-01, &
719  0.1482e-01, 0.9692e-02, 0.8247e-02, 0.7640e-02, 0.7582e-02, 0.7783e-02, 0.5432e-02, 0.4482e-02, 0.4919e-02, 0.5903e-02, &
720  0.3770e-01, 0.2401e-01, 0.1449e-01, 0.9252e-02, 0.7620e-02, 0.6678e-02, 0.6845e-02, 0.6925e-02, 0.4939e-02, 0.3471e-02, &
721  0.4124e-02, 0.3873e-02, 0.3776e-01, 0.2395e-01, 0.1419e-01, 0.8959e-02, 0.7096e-02, 0.6184e-02, 0.6110e-02, 0.6075e-02, &
722  0.4419e-02, 0.2891e-02, 0.3056e-02, 0.1214e-02, 0.3780e-01, 0.2391e-01, 0.1392e-01, 0.8687e-02, 0.6573e-02, 0.5733e-02, &
723  0.5359e-02, 0.5009e-02, 0.4034e-02, 0.2755e-02, 0.1968e-02,-0.4187e-04, 0.3791e-01, 0.2382e-01, 0.1373e-01, 0.8561e-02, &
724  0.6060e-02, 0.5120e-02, 0.4618e-02, 0.4713e-02, 0.3965e-02, 0.2481e-02, 0.8164e-03,-0.1088e-02, 0.3843e-01, 0.2148e-01, &
725  0.1302e-01, 0.6384e-02, 0.5256e-02, 0.4260e-02, 0.4077e-02, 0.4181e-02, 0.4132e-02, 0.2135e-02,-0.2931e-03,-0.1151e-02, &
726  0.3896e-01, 0.2081e-01, 0.1097e-01, 0.5568e-02, 0.4475e-02, 0.3795e-02, 0.3828e-02, 0.3996e-02, 0.3766e-02, 0.1193e-02, &
727  -0.1089e-02,-0.9420e-03, 0.3973e-01, 0.2024e-01, 0.9943e-02, 0.4815e-02, 0.3820e-02, 0.3663e-02, 0.3568e-02, 0.3881e-02, &
728  0.2859e-02, 0.6698e-03,-0.1549e-02,-0.6280e-03, 0.3635e-01, 0.1963e-01, 0.1061e-01, 0.3812e-02, 0.3509e-02, 0.3429e-02, &
729  0.3693e-02, 0.3316e-02, 0.1120e-02, 0.6552e-03,-0.1193e-02,-0.1109e-02, 0.3631e-01, 0.1893e-01, 0.1056e-01, 0.3172e-02, &
730  0.3378e-02, 0.3164e-02, 0.2751e-02, 0.1722e-02, 0.1112e-02, 0.4354e-03,-0.7327e-03,-0.1319e-02, 0.3500e-01, 0.1828e-01, &
731  0.1050e-01, 0.2831e-02, 0.2784e-02, 0.2564e-02, 0.1469e-02, 0.7739e-03, 0.1209e-02, 0.7913e-03,-0.2512e-03,-0.1758e-02, &
732  0.3352e-01, 0.1763e-01, 0.1045e-01, 0.2401e-02, 0.1928e-02, 0.1340e-02, 0.3753e-03, 0.5794e-03, 0.9060e-03, 0.1042e-02, &
733  0.1465e-03,-0.2533e-02, 0.2880e-01, 0.1729e-01, 0.1077e-01, 0.1347e-02, 0.1194e-02,-0.1191e-03, 0.2828e-03, 0.6606e-03, &
734  0.9743e-03, 0.1002e-02, 0.0000e+00,-0.3140e-02, 0.2040e-01, 0.1585e-01, 0.1165e-01, 0.3871e-05, 0.1509e-04,-0.1046e-02, &
735  0.2444e-03, 0.4359e-03, 0.1041e-02, 0.2429e-02,-0.1721e-03,-0.2786e-02, 0.1737e-01, 0.1560e-01, 0.1240e-01,-0.2139e-03, &
736  -0.1025e-02,-0.1248e-02,-0.6934e-04, 0.1649e-03, 0.4062e-03, 0.1554e-02,-0.4179e-03,-0.7795e-03 /
737 
738  data((c15hca(t_3,j,i),i=1,kg(15)),j=1,numpir) / &
739  -0.1488e-03,-0.9248e-04,-0.2322e-04,-0.4187e-05, 0.1104e-04, 0.9895e-05,-0.2283e-05, 0.2512e-05,-0.9058e-05, 0.8449e-05, &
740  0.8297e-05,-0.3882e-04,-0.1488e-03,-0.9058e-04,-0.2398e-04,-0.5709e-05, 0.1218e-04, 0.1180e-04, 0.1522e-05, 0.6927e-05, &
741  -0.1161e-04, 0.1714e-04,-0.4948e-06,-0.3540e-04,-0.1500e-03,-0.8830e-04,-0.2474e-04,-0.8373e-05, 0.6470e-05, 0.7992e-05, &
742  0.9096e-05, 0.6737e-05,-0.1485e-04, 0.1873e-04,-0.4948e-06,-0.4491e-04,-0.1500e-03,-0.8601e-04,-0.2664e-04,-0.1028e-04, &
743  0.6851e-05, 0.6851e-05, 0.1294e-04,-0.2550e-05,-0.1520e-04, 0.2310e-04, 0.4948e-06,-0.2017e-04,-0.1507e-03,-0.8373e-04, &
744  -0.2664e-04,-0.1256e-04, 0.4567e-05, 0.1028e-04, 0.9210e-05,-0.2131e-05,-0.6995e-05, 0.7498e-05,-0.1104e-04,-0.2284e-05, &
745  -0.1519e-03,-0.8183e-04,-0.2816e-04,-0.1142e-04, 0.7611e-06, 0.7231e-05, 0.1751e-05,-0.7612e-06, 0.8312e-05, 0.2436e-05, &
746  -0.7231e-05, 0.2398e-04,-0.1530e-03,-0.7992e-04,-0.2893e-04,-0.9896e-05, 0.3806e-06, 0.8906e-05, 0.3159e-05,-0.5328e-05, &
747  0.3692e-05,-0.2093e-05,-0.6851e-05,-0.3045e-05,-0.1538e-03,-0.7536e-04,-0.3007e-04,-0.8754e-05,-0.3045e-05, 0.5138e-05, &
748  0.9134e-06,-0.1979e-06, 0.1560e-05,-0.1507e-04, 0.2284e-04, 0.9895e-05,-0.1541e-03,-0.7688e-04,-0.2969e-04,-0.5709e-05, &
749  -0.3996e-05, 0.1142e-05,-0.8373e-06, 0.1235e-04,-0.7079e-05,-0.6737e-05, 0.1028e-04, 0.3578e-04,-0.1560e-03,-0.6851e-04, &
750  -0.1903e-04,-0.4187e-05,-0.4605e-05,-0.1142e-06, 0.3878e-05, 0.3597e-05,-0.9591e-05, 0.5328e-05, 0.7612e-05,-0.4948e-05, &
751  -0.1587e-03,-0.6546e-04,-0.2740e-04,-0.7612e-06,-0.3578e-05, 0.1713e-05, 0.6064e-05,-0.9781e-05, 0.1408e-05, 0.5709e-05, &
752  0.8373e-05,-0.1256e-04,-0.1484e-03,-0.5823e-04,-0.4301e-04,-0.1522e-05, 0.7498e-05,-0.5328e-06,-0.7855e-05,-0.1599e-05, &
753  0.1964e-04,-0.2284e-05, 0.7882e-10, 0.5328e-05,-0.1238e-03,-0.5700e-04,-0.5266e-04, 0.3286e-05, 0.4910e-05,-0.8602e-05, &
754  0.6090e-06, 0.8454e-05, 0.1256e-05,-0.4072e-05,-0.1903e-05, 0.6470e-05,-0.1155e-03,-0.5231e-04,-0.4396e-04, 0.3626e-05, &
755  -0.7051e-05,-0.1743e-05, 0.9667e-05, 0.2064e-04,-0.2778e-05,-0.6546e-05,-0.4948e-05, 0.1903e-05,-0.1024e-03,-0.5129e-04, &
756  -0.4506e-04, 0.7943e-06, 0.3074e-06, 0.3243e-05, 0.2754e-04,-0.1479e-05, 0.1661e-05,-0.2969e-05,-0.1066e-04, 0.7612e-06, &
757  -0.8473e-04,-0.5418e-04,-0.4674e-04,-0.3418e-05, 0.9460e-05, 0.1151e-04, 0.5714e-05,-0.1069e-04,-0.2022e-05,-0.9061e-05, &
758  -0.1104e-04,-0.3083e-04,-0.4283e-04,-0.5037e-04,-0.4476e-04, 0.1951e-04, 0.8922e-05, 0.1296e-04,-0.4053e-05,-0.4355e-05, &
759  -0.2355e-05,-0.5004e-05,-0.1218e-04,-0.1522e-04, 0.6411e-05,-0.5937e-04,-0.5331e-04, 0.1934e-04, 0.5284e-05, 0.1129e-04, &
760  -0.2166e-05,-0.1484e-06,-0.5407e-05,-0.1364e-04,-0.3115e-05, 0.3004e-04,-0.5074e-04,-0.6256e-04,-0.5097e-04, 0.2218e-04, &
761  0.1228e-04,-0.1160e-05,-0.1105e-05, 0.1618e-06,-0.6089e-05,-0.4216e-06,-0.5314e-05, 0.7903e-05 /
762 
763 
764  data((c15hcb(t_1,j,i),i=1,kg(15)),j=1,numpir) / &
765  -0.9593e+01,-0.4078e+01,-0.2812e+01,-0.6506e+00,-0.4123e+00, 0.2055e+01, 0.4097e+01, 0.4671e+01, 0.4639e+01,-0.4587e+02, &
766  -0.4587e+02,-0.4587e+02,-0.9276e+01,-0.3757e+01,-0.2467e+01,-0.5784e+00, 0.8833e-01, 0.2232e+01, 0.3826e+01, 0.4723e+01, &
767  0.4942e+01, 0.5135e+01,-0.4587e+02,-0.4587e+02,-0.8968e+01,-0.3508e+01,-0.2116e+01,-0.1363e+00, 0.1662e+00, 0.2424e+01, &
768  0.4220e+01, 0.4513e+01, 0.1375e+01, 0.4601e+01,-0.4587e+02,-0.4587e+02,-0.8662e+01,-0.3164e+01,-0.1722e+01, 0.5178e-01, &
769  0.7288e+00, 0.2411e+01, 0.3805e+01, 0.4766e+01, 0.4342e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.8292e+01,-0.2799e+01, &
770  -0.1359e+01, 0.3271e+00, 0.1650e+01, 0.2395e+01, 0.4192e+01, 0.4758e+01, 0.2470e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
771  -0.7812e+01,-0.2404e+01,-0.1085e+01, 0.7167e+00, 0.2202e+01, 0.2922e+01, 0.4322e+01, 0.4591e+01, 0.4186e+01,-0.4587e+02, &
772  -0.4587e+02,-0.4587e+02,-0.7441e+01,-0.2066e+01,-0.7142e+00, 0.1057e+01, 0.2524e+01, 0.2946e+01, 0.4220e+01, 0.3607e+01, &
773  -0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.7191e+01,-0.1745e+01,-0.3487e+00, 0.1453e+01, 0.2739e+01, 0.3660e+01, &
774  0.4114e+01, 0.3245e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.6895e+01,-0.1326e+01,-0.3500e+00, 0.1647e+01, &
775  0.2899e+01, 0.4023e+01, 0.3361e+01, 0.3360e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.5876e+01,-0.9573e+00, &
776  0.2014e+00, 0.2130e+01, 0.3493e+01, 0.4088e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
777  -0.4429e+01,-0.3417e+00, 0.1204e+01, 0.2780e+01, 0.3843e+01, 0.3099e+01,-0.4587e+02, 0.3605e+01,-0.4587e+02,-0.4587e+02, &
778  -0.4587e+02,-0.4587e+02,-0.3122e+01, 0.2697e+00, 0.1866e+01, 0.3526e+01, 0.3569e+01, 0.1025e+01,-0.4587e+02,-0.4587e+02, &
779  -0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.2284e+01, 0.8186e+00, 0.2754e+01, 0.3206e+01, 0.3704e+01,-0.4587e+02, &
780  -0.4587e+02, 0.4625e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.1711e+01, 0.1220e+01, 0.3248e+01,-0.4587e+02, &
781  0.2565e+01, 0.3297e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.1758e+01, 0.7970e+00, &
782  0.2758e+01, 0.2926e+01, 0.2613e+01, 0.1974e+01,-0.4587e+02, 0.2310e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.4587e+02, &
783  -0.1737e+01, 0.3499e+00, 0.2246e+01, 0.2673e+01, 0.3308e+01, 0.3463e+01, 0.3103e+01, 0.2611e+01, 0.2178e+01,-0.4587e+02, &
784  -0.4587e+02,-0.4587e+02,-0.1559e+01, 0.2215e+00, 0.1875e+01, 0.2500e+01, 0.3346e+01, 0.3585e+01, 0.3946e+01, 0.3533e+01, &
785  0.3205e+01,-0.4587e+02,-0.4587e+02,-0.4587e+02,-0.1601e+01, 0.5060e-01, 0.1275e+01, 0.2176e+01, 0.3081e+01, 0.3649e+01, &
786  0.3940e+01, 0.4106e+01, 0.4112e+01, 0.4349e+01, 0.2292e+01,-0.4587e+02,-0.1222e+01, 0.3199e+00, 0.1642e+01, 0.2380e+01, &
787  0.3254e+01, 0.3534e+01, 0.3687e+01, 0.3717e+01, 0.3402e+01, 0.3868e+01,-0.4587e+02,-0.4587e+02 /
788 
789  data((c15hcb(t_2,j,i),i=1,kg(15)),j=1,numpir) / &
790  0.2967e-01, 0.1697e-01, 0.1795e-01, 0.1387e-01, 0.2032e-01, 0.1187e-01, 0.2560e-01, 0.1044e-01,-0.4560e+00, 0.3059e-02, &
791  0.3059e-02, 0.3059e-02, 0.2998e-01, 0.1586e-01, 0.1786e-01, 0.1521e-01, 0.1710e-01, 0.1061e-01, 0.2030e-01, 0.1158e-01, &
792  0.4452e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.2993e-01, 0.1551e-01, 0.1481e-01, 0.9846e-02, 0.2443e-01, 0.1150e-01, &
793  0.1865e-01, 0.1376e-01, 0.4617e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3035e-01, 0.1417e-01, 0.1438e-01, 0.1511e-01, &
794  0.1901e-01, 0.8582e-02, 0.1746e-01, 0.1450e-01, 0.4523e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.2970e-01, 0.1347e-01, &
795  0.1322e-01, 0.1252e-01, 0.1665e-01, 0.1037e-01, 0.1320e-01, 0.1199e-01, 0.4436e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
796  0.2949e-01, 0.1291e-01, 0.1671e-01, 0.1111e-01, 0.1400e-01, 0.1318e-01, 0.1060e-01, 0.1046e-01, 0.3059e-02, 0.3059e-02, &
797  0.3059e-02, 0.3059e-02, 0.3004e-01, 0.1300e-01, 0.1413e-01, 0.9085e-02, 0.9764e-02, 0.2260e-01, 0.9778e-02, 0.4671e+00, &
798  0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3086e-01, 0.1436e-01, 0.1205e-01, 0.1081e-01, 0.4681e-02, 0.1479e-01, &
799  0.1888e-01, 0.3494e-01, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3094e-01, 0.1500e-01, 0.1457e-01, 0.1060e-01, &
800  0.8319e-02, 0.8983e-02, 0.3791e-01, 0.2232e-01, 0.4631e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3158e-01, 0.1585e-01, &
801  0.1292e-01, 0.6531e-02, 0.1383e-01, 0.4605e+00, 0.4662e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
802  0.3182e-01, 0.1586e-01, 0.8724e-02, 0.5798e-02, 0.2454e-01, 0.4607e+00, 0.4560e+00, 0.4511e+00, 0.3059e-02, 0.3059e-02, &
803  0.3059e-02, 0.3059e-02, 0.2369e-01, 0.1606e-01, 0.5477e-02, 0.1228e-01, 0.4579e+00, 0.4561e+00, 0.4497e+00, 0.3059e-02, &
804  0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.2190e-01, 0.1779e-01, 0.6267e-02, 0.4535e+00, 0.4533e+00, 0.3059e-02, &
805  0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.2100e-01, 0.1653e-01, 0.7449e-02, 0.4543e+00, &
806  0.4472e+00, 0.4439e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.3059e-02, 0.1864e-01, 0.1771e-01, &
807  0.7040e-02, 0.2877e-01, 0.3381e-01, 0.2691e-01, 0.4466e+00, 0.3059e-02, 0.4613e+00, 0.3059e-02, 0.3059e-02, 0.3059e-02, &
808  0.1637e-01, 0.1641e-01, 0.8424e-02, 0.1318e-01, 0.2060e-01, 0.3426e-01, 0.4122e-01, 0.4621e+00, 0.4555e+00, 0.4525e+00, &
809  0.3059e-02, 0.3059e-02, 0.1607e-01, 0.1452e-01, 0.8013e-02, 0.1213e-01, 0.1482e-01, 0.2125e-01, 0.3379e-01, 0.3562e-01, &
810  0.4619e+00, 0.4569e+00, 0.3059e-02, 0.3059e-02, 0.1698e-01, 0.1538e-01, 0.6616e-02, 0.1147e-01, 0.1217e-01, 0.1696e-01, &
811  0.1871e-01, 0.2273e-01, 0.4513e-01, 0.4702e+00, 0.4617e+00, 0.4553e+00, 0.1700e-01, 0.1547e-01, 0.6456e-02, 0.1324e-01, &
812  0.1502e-01, 0.2095e-01, 0.2547e-01, 0.2823e-01, 0.4107e-01, 0.4676e+00, 0.4583e+00, 0.4498e+00 /
813 
814  data((c15hcb(t_3,j,i),i=1,kg(15)),j=1,numpir) / &
815  -0.6747e-05,-0.2483e-04, 0.6575e-04, 0.1026e-03, 0.3888e-03,-0.8519e-04,-0.1629e-03,-0.1808e-04,-0.8355e-02,-0.4656e-05, &
816  -0.4656e-05,-0.4656e-05,-0.2270e-04,-0.3427e-04, 0.5118e-04, 0.1218e-03, 0.1245e-03,-0.1245e-03, 0.3841e-05,-0.4151e-04, &
817  -0.8763e-02,-0.1687e-01,-0.4656e-05,-0.4656e-05,-0.4557e-04,-0.3023e-04, 0.2286e-04, 0.5656e-04, 0.4113e-04,-0.1407e-03, &
818  -0.1301e-03, 0.8503e-04,-0.7284e-02,-0.1669e-01,-0.4656e-05,-0.4656e-05,-0.5325e-04,-0.5309e-04,-0.1246e-04, 0.2244e-04, &
819  0.5136e-04,-0.1272e-03, 0.4217e-04,-0.1749e-04,-0.8435e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.6857e-04,-0.7217e-04, &
820  0.1740e-05, 0.3653e-04,-0.1490e-03,-0.4090e-04,-0.2376e-04, 0.2047e-04,-0.7974e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
821  -0.1232e-03,-0.9826e-04,-0.2849e-04, 0.1703e-04,-0.1895e-03,-0.3363e-03, 0.7102e-04,-0.1838e-05,-0.1655e-01,-0.4656e-05, &
822  -0.4656e-05,-0.4656e-05,-0.9896e-04,-0.5127e-04,-0.2704e-04,-0.1218e-04,-0.1207e-03,-0.5883e-04, 0.6893e-04,-0.7924e-02, &
823  -0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.7837e-04,-0.4980e-04, 0.6902e-05,-0.1072e-03,-0.4051e-04,-0.1991e-05, &
824  -0.1173e-03,-0.5195e-04,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.8136e-04,-0.8102e-04, 0.1254e-03,-0.4658e-04, &
825  0.3173e-04,-0.4461e-05,-0.1558e-03,-0.2036e-03, 0.8360e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.2232e-04,-0.6411e-04, &
826  0.9486e-04,-0.2322e-03,-0.8282e-04,-0.8202e-02, 0.8416e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
827  -0.1398e-03,-0.7165e-04,-0.4258e-04,-0.3970e-04,-0.2839e-03,-0.7873e-02, 0.8231e-02,-0.8213e-02,-0.4656e-05,-0.4656e-05, &
828  -0.4656e-05,-0.4656e-05,-0.6754e-04,-0.7469e-04,-0.6898e-04,-0.1702e-03,-0.8079e-02,-0.7270e-02, 0.8116e-02,-0.4656e-05, &
829  -0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.2396e-04,-0.2361e-04,-0.8664e-04,-0.8038e-02,-0.8207e-02,-0.4656e-05, &
830  -0.4656e-05,-0.1670e-01,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.5479e-04,-0.7593e-04,-0.1005e-03, 0.8199e-02, &
831  -0.7942e-02,-0.8244e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.4656e-05,-0.3806e-04,-0.5825e-04, &
832  -0.1003e-03,-0.2925e-03,-0.1506e-03, 0.3148e-04, 0.8060e-02,-0.1593e-01, 0.8327e-02,-0.4656e-05,-0.4656e-05,-0.4656e-05, &
833  -0.4706e-04,-0.3630e-04,-0.7811e-04,-0.6881e-04,-0.1822e-03,-0.3091e-03,-0.3033e-03,-0.7684e-02,-0.7663e-02, 0.8167e-02, &
834  -0.4656e-05,-0.4656e-05,-0.7669e-04,-0.4610e-04,-0.8063e-04,-0.7250e-04,-0.1094e-03,-0.1241e-03,-0.2944e-03,-0.1736e-03, &
835  -0.7886e-02, 0.8248e-02,-0.4656e-05,-0.4656e-05,-0.7138e-04,-0.4545e-04,-0.3653e-04,-0.6075e-04,-0.4528e-04,-0.1077e-03, &
836  -0.1119e-03,-0.1657e-03,-0.4695e-03,-0.8112e-02,-0.7587e-02, 0.8217e-02,-0.6812e-04,-0.4558e-04,-0.6739e-04,-0.8861e-04, &
837  -0.9386e-04,-0.1334e-03,-0.2007e-03,-0.2179e-03,-0.1650e-03,-0.8001e-02, 0.8273e-02, 0.8118e-02 /
838 
839 
840  data((c16h2o(t_1,j,i),i=1,kg(16)),j=1,numpir) / &
841  -0.2344e+02,-0.2016e+02,-0.1986e+02,-0.1655e+02,-0.1243e+02,-0.8437e+01,-0.4858e+01,-0.2298e+02,-0.2014e+02,-0.1984e+02, &
842  -0.1609e+02,-0.1198e+02,-0.8020e+01,-0.4548e+01,-0.2252e+02,-0.2012e+02,-0.1981e+02,-0.1564e+02,-0.1153e+02,-0.7596e+01, &
843  -0.4239e+01,-0.2206e+02,-0.2009e+02,-0.1957e+02,-0.1517e+02,-0.1111e+02,-0.7161e+01,-0.3871e+01,-0.2160e+02,-0.2007e+02, &
844  -0.1911e+02,-0.1472e+02,-0.1065e+02,-0.6721e+01,-0.3479e+01,-0.2113e+02,-0.2005e+02,-0.1865e+02,-0.1426e+02,-0.1021e+02, &
845  -0.6302e+01,-0.3081e+01,-0.2067e+02,-0.2003e+02,-0.1819e+02,-0.1379e+02,-0.9765e+01,-0.5883e+01,-0.2678e+01,-0.2026e+02, &
846  -0.2001e+02,-0.1773e+02,-0.1333e+02,-0.9332e+01,-0.5443e+01,-0.2253e+01,-0.2024e+02,-0.1999e+02,-0.1727e+02,-0.1288e+02, &
847  -0.8897e+01,-0.5029e+01,-0.1858e+01,-0.2026e+02,-0.1959e+02,-0.1481e+02,-0.1147e+02,-0.7477e+01,-0.4555e+01,-0.1464e+01, &
848  -0.2022e+02,-0.1632e+02,-0.1305e+02,-0.9885e+01,-0.6689e+01,-0.4108e+01,-0.1068e+01,-0.1936e+02,-0.1438e+02,-0.1163e+02, &
849  -0.8499e+01,-0.6146e+01,-0.3673e+01,-0.6816e+00,-0.1675e+02,-0.1281e+02,-0.1020e+02,-0.7716e+01,-0.5678e+01,-0.3256e+01, &
850  -0.3125e+00,-0.1510e+02,-0.1124e+02,-0.8821e+01,-0.7140e+01,-0.5243e+01,-0.2851e+01,-0.2560e-01,-0.1334e+02,-0.9708e+01, &
851  -0.8061e+01,-0.6611e+01,-0.4842e+01,-0.2459e+01, 0.1711e+00,-0.1155e+02,-0.8798e+01,-0.7440e+01,-0.6123e+01,-0.4439e+01, &
852  -0.2089e+01, 0.2480e+00,-0.1020e+02,-0.8154e+01,-0.6945e+01,-0.5681e+01,-0.4055e+01,-0.1737e+01, 0.2390e+00,-0.9464e+01, &
853  -0.7677e+01,-0.6512e+01,-0.5284e+01,-0.3707e+01,-0.1453e+01, 0.2015e+00,-0.9033e+01,-0.7246e+01,-0.6093e+01,-0.4882e+01, &
854  -0.3346e+01,-0.1264e+01, 0.1033e+00 /
855 
856  data((c16h2o(t_2,j,i),i=1,kg(16)),j=1,numpir) / &
857  0.4658e-01, 0.5840e-02, 0.4626e-02, 0.2688e-01, 0.2395e-01, 0.1804e-01, 0.2074e-01, 0.4660e-01, 0.1884e-02, 0.8561e-02, &
858  0.2690e-01, 0.2403e-01, 0.1788e-01, 0.1934e-01, 0.4660e-01, 0.1800e-02, 0.1252e-01, 0.2694e-01, 0.2393e-01, 0.1786e-01, &
859  0.1825e-01, 0.4660e-01, 0.1779e-02, 0.1649e-01, 0.2696e-01, 0.2397e-01, 0.1779e-01, 0.1765e-01, 0.4348e-01, 0.1758e-02, &
860  0.2043e-01, 0.2696e-01, 0.2393e-01, 0.1748e-01, 0.1675e-01, 0.3944e-01, 0.1737e-02, 0.2445e-01, 0.2698e-01, 0.2384e-01, &
861  0.1752e-01, 0.1549e-01, 0.3538e-01, 0.1654e-02, 0.2847e-01, 0.2702e-01, 0.2384e-01, 0.1714e-01, 0.1565e-01, 0.3127e-01, &
862  0.1570e-02, 0.3245e-01, 0.2705e-01, 0.2374e-01, 0.1712e-01, 0.1514e-01, 0.2715e-01, 0.1444e-02, 0.3540e-01, 0.2711e-01, &
863  0.2363e-01, 0.1702e-01, 0.1446e-01, 0.2960e-01, 0.1760e-01, 0.2977e-01, 0.2397e-01, 0.2087e-01, 0.1618e-01, 0.1445e-01, &
864  0.2466e-01, 0.3039e-01, 0.2428e-01, 0.2217e-01, 0.1821e-01, 0.1593e-01, 0.1463e-01, 0.2640e-01, 0.2545e-01, 0.2231e-01, &
865  0.2060e-01, 0.1773e-01, 0.1555e-01, 0.1473e-01, 0.3456e-01, 0.2135e-01, 0.2030e-01, 0.1844e-01, 0.1740e-01, 0.1559e-01, &
866  0.1428e-01, 0.3203e-01, 0.2047e-01, 0.1809e-01, 0.1760e-01, 0.1725e-01, 0.1545e-01, 0.1541e-01, 0.2137e-01, 0.1857e-01, &
867  0.1616e-01, 0.1698e-01, 0.1700e-01, 0.1537e-01, 0.1636e-01, 0.1338e-01, 0.1518e-01, 0.1580e-01, 0.1658e-01, 0.1710e-01, &
868  0.1518e-01, 0.1513e-01, 0.1570e-01, 0.1614e-01, 0.1603e-01, 0.1673e-01, 0.1706e-01, 0.1497e-01, 0.1439e-01, 0.1987e-01, &
869  0.1731e-01, 0.1601e-01, 0.1675e-01, 0.1681e-01, 0.1535e-01, 0.1425e-01, 0.2018e-01, 0.1723e-01, 0.1597e-01, 0.1691e-01, &
870  0.1666e-01, 0.1509e-01, 0.1446e-01 /
871 
872  data((c16h2o(t_3,j,i),i=1,kg(16)),j=1,numpir) / &
873  -0.2873e-03,-0.8031e-04, 0.4225e-04,-0.9287e-04,-0.6013e-04,-0.4339e-04,-0.2474e-04,-0.2862e-03,-0.8372e-05, 0.1146e-03, &
874  -0.9248e-04,-0.6166e-04,-0.3882e-04,-0.1827e-04,-0.2870e-03,-0.6851e-05, 0.1865e-03,-0.9172e-04,-0.6128e-04,-0.3616e-04, &
875  -0.7612e-05,-0.2877e-03,-0.7231e-05, 0.1880e-03,-0.9287e-04,-0.5671e-04,-0.4110e-04,-0.1104e-04,-0.3429e-03,-0.7612e-05, &
876  0.1149e-03,-0.9287e-04,-0.6356e-04,-0.4529e-04,-0.2436e-04,-0.4187e-03,-0.7992e-05, 0.4339e-04,-0.9325e-04,-0.6280e-04, &
877  -0.4225e-04,-0.3197e-04,-0.4925e-03,-0.8754e-05,-0.2740e-04,-0.9477e-04,-0.6432e-04,-0.3768e-04,-0.3361e-04,-0.5511e-03, &
878  -0.8753e-05,-0.9972e-04,-0.9515e-04,-0.6394e-04,-0.3806e-04,-0.3787e-04,-0.4792e-03,-0.1028e-04,-0.1534e-03,-0.9477e-04, &
879  -0.6356e-04,-0.3616e-04,-0.2923e-04,-0.5070e-03, 0.1922e-03,-0.1028e-03,-0.5823e-04,-0.7954e-04,-0.2550e-04,-0.3893e-04, &
880  -0.3776e-03,-0.1043e-03,-0.7993e-04,-0.7422e-04,-0.4948e-04,-0.3007e-04,-0.3863e-04, 0.8335e-04,-0.5709e-04,-0.6090e-04, &
881  -0.7840e-04,-0.3692e-04,-0.3007e-04,-0.4251e-04,-0.6204e-04,-0.4872e-04,-0.3806e-04,-0.4681e-04,-0.3463e-04,-0.3007e-04, &
882  -0.4312e-04,-0.1142e-04,-0.5176e-04,-0.5024e-04,-0.3007e-04,-0.3730e-04,-0.3037e-04,-0.3888e-04, 0.2550e-04,-0.6508e-04, &
883  -0.2512e-04,-0.3083e-04,-0.3197e-04,-0.3041e-04,-0.3750e-04, 0.1484e-04,-0.1941e-04,-0.2626e-04,-0.3349e-04,-0.3463e-04, &
884  -0.2896e-04,-0.1716e-04,-0.7231e-04,-0.3920e-04,-0.2893e-04,-0.3540e-04,-0.3311e-04,-0.3734e-04,-0.2550e-05,-0.7650e-04, &
885  -0.3159e-04,-0.2778e-04,-0.3121e-04,-0.2169e-04,-0.4365e-04,-0.1546e-04,-0.7916e-04,-0.2931e-04,-0.2854e-04,-0.3654e-04, &
886  -0.1979e-04,-0.4811e-04,-0.1435e-04 /
887 
888 
889  data((c17h2o(t_1,j,i),i=1,kg(17)),j=1,numpir) / &
890  -0.2255e+02,-0.2000e+02,-0.1703e+02,-0.1282e+02,-0.9215e+01,-0.5938e+01,-0.2009e+01,-0.2209e+02,-0.1997e+02,-0.1657e+02, &
891  -0.1236e+02,-0.8764e+01,-0.5499e+01,-0.1582e+01,-0.2163e+02,-0.1993e+02,-0.1611e+02,-0.1191e+02,-0.8324e+01,-0.5061e+01, &
892  -0.1170e+01,-0.2117e+02,-0.1990e+02,-0.1565e+02,-0.1146e+02,-0.7889e+01,-0.4631e+01,-0.7737e+00,-0.2071e+02,-0.1987e+02, &
893  -0.1519e+02,-0.1100e+02,-0.7440e+01,-0.4179e+01,-0.3719e+00,-0.2026e+02,-0.1985e+02,-0.1473e+02,-0.1054e+02,-0.6995e+01, &
894  -0.3721e+01, 0.0000e+00,-0.2024e+02,-0.1982e+02,-0.1426e+02,-0.1009e+02,-0.6549e+01,-0.3284e+01, 0.4053e+00,-0.2022e+02, &
895  -0.1980e+02,-0.1381e+02,-0.9639e+01,-0.6097e+01,-0.2821e+01, 0.8375e+00,-0.2021e+02,-0.1933e+02,-0.1335e+02,-0.9187e+01, &
896  -0.5653e+01,-0.2379e+01, 0.1272e+01,-0.2010e+02,-0.1503e+02,-0.1125e+02,-0.7665e+01,-0.4492e+01,-0.1893e+01, 0.1642e+01, &
897  -0.1747e+02,-0.1278e+02,-0.9547e+01,-0.6120e+01,-0.3756e+01,-0.1443e+01, 0.1995e+01,-0.1529e+02,-0.1095e+02,-0.8107e+01, &
898  -0.5036e+01,-0.3182e+01,-0.1032e+01, 0.2429e+01,-0.1370e+02,-0.9303e+01,-0.6691e+01,-0.4357e+01,-0.2683e+01,-0.6173e+00, &
899  0.2805e+01,-0.1150e+02,-0.7859e+01,-0.5618e+01,-0.3843e+01,-0.2234e+01,-0.2171e+00, 0.2973e+01,-0.9590e+01,-0.6537e+01, &
900  -0.4886e+01,-0.3355e+01,-0.1805e+01, 0.1615e+00, 0.3157e+01,-0.7530e+01,-0.5699e+01,-0.4306e+01,-0.2892e+01,-0.1388e+01, &
901  0.5448e+00, 0.3155e+01,-0.6758e+01,-0.5112e+01,-0.3809e+01,-0.2464e+01,-0.9947e+00, 0.8713e+00, 0.3203e+01,-0.6245e+01, &
902  -0.4610e+01,-0.3376e+01,-0.2058e+01,-0.6166e+00, 0.1073e+01, 0.3109e+01,-0.5777e+01,-0.4175e+01,-0.2963e+01,-0.1671e+01, &
903  -0.2556e+00, 0.1241e+01, 0.3014e+01 /
904 
905  data((c17h2o(t_2,j,i),i=1,kg(17)),j=1,numpir) / &
906  0.4264e-01, 0.1968e-02, 0.1863e-01, 0.1436e-01, 0.1101e-01, 0.1055e-01, 0.1281e-01, 0.4264e-01, 0.1989e-02, 0.1861e-01, &
907  0.1438e-01, 0.1095e-01, 0.1030e-01, 0.1211e-01, 0.3996e-01, 0.1968e-02, 0.1861e-01, 0.1434e-01, 0.1103e-01, 0.1019e-01, &
908  0.1160e-01, 0.3600e-01, 0.1947e-02, 0.1861e-01, 0.1442e-01, 0.1086e-01, 0.1003e-01, 0.1157e-01, 0.3203e-01, 0.5756e-02, &
909  0.1861e-01, 0.1444e-01, 0.1080e-01, 0.9922e-02, 0.1151e-01, 0.2801e-01, 0.9713e-02, 0.1859e-01, 0.1446e-01, 0.1070e-01, &
910  0.9880e-02, 0.1066e-01, 0.2393e-01, 0.1369e-01, 0.1859e-01, 0.1451e-01, 0.1057e-01, 0.9880e-02, 0.1072e-01, 0.1987e-01, &
911  0.1767e-01, 0.1863e-01, 0.1451e-01, 0.1040e-01, 0.9880e-02, 0.1057e-01, 0.1572e-01, 0.2169e-01, 0.1863e-01, 0.1442e-01, &
912  0.1022e-01, 0.9742e-02, 0.1036e-01, 0.3391e-02, 0.1884e-01, 0.1566e-01, 0.1105e-01, 0.1011e-01, 0.1001e-01, 0.1017e-01, &
913  0.1982e-01, 0.1444e-01, 0.1189e-01, 0.1030e-01, 0.9859e-02, 0.9861e-02, 0.1038e-01, 0.1748e-01, 0.1321e-01, 0.9922e-02, &
914  0.1068e-01, 0.1013e-01, 0.9937e-02, 0.9958e-02, 0.1346e-01, 0.9943e-02, 0.9566e-02, 0.1097e-01, 0.9815e-02, 0.9964e-02, &
915  0.1059e-01, 0.9817e-02, 0.7159e-02, 0.8687e-02, 0.1114e-01, 0.1007e-01, 0.1014e-01, 0.1058e-01, 0.3370e-02, 0.7264e-02, &
916  0.9378e-02, 0.1112e-01, 0.9767e-02, 0.1016e-01, 0.1101e-01, 0.2993e-02, 0.8017e-02, 0.9566e-02, 0.1116e-01, 0.9738e-02, &
917  0.1025e-01, 0.1086e-01, 0.8331e-02, 0.8771e-02, 0.1001e-01, 0.1117e-01, 0.9847e-02, 0.1076e-01, 0.1084e-01, 0.7850e-02, &
918  0.9378e-02, 0.1001e-01, 0.1105e-01, 0.9964e-02, 0.1113e-01, 0.1168e-01, 0.8038e-02, 0.9336e-02, 0.9817e-02, 0.1096e-01, &
919  0.1024e-01, 0.1175e-01, 0.1107e-01 /
920 
921  data((c17h2o(t_3,j,i),i=1,kg(17)),j=1,numpir) / &
922  -0.2188e-03,-0.2283e-05,-0.8069e-04,-0.4415e-04,-0.2284e-04,-0.4491e-04,-0.4518e-04,-0.2196e-03,-0.2665e-05,-0.8107e-04, &
923  -0.4301e-04,-0.2398e-04,-0.4795e-04,-0.4693e-04,-0.2683e-03,-0.3045e-05,-0.8107e-04,-0.4301e-04,-0.2246e-04,-0.4757e-04, &
924  -0.4152e-04,-0.3403e-03,-0.4187e-05,-0.8031e-04,-0.3996e-04,-0.1865e-04,-0.4301e-04,-0.4350e-04,-0.4118e-03, 0.6584e-04, &
925  -0.8107e-04,-0.4034e-04,-0.1903e-04,-0.4643e-04,-0.4834e-04,-0.4803e-03, 0.1378e-03,-0.8069e-04,-0.4072e-04,-0.1713e-04, &
926  -0.5176e-04,-0.3460e-04,-0.4099e-03, 0.2101e-03,-0.8069e-04,-0.3920e-04,-0.1713e-04,-0.5024e-04,-0.3524e-04,-0.3391e-03, &
927  0.2809e-03,-0.7992e-04,-0.3616e-04,-0.2017e-04,-0.5633e-04,-0.4886e-04,-0.2668e-03, 0.2078e-03,-0.8069e-04,-0.3768e-04, &
928  -0.2131e-04,-0.5580e-04,-0.5454e-04,-0.2207e-04,-0.8601e-04,-0.4643e-04,-0.2436e-04,-0.4148e-04,-0.5458e-04,-0.4579e-04, &
929  -0.5138e-04,-0.2893e-04,-0.3273e-04,-0.3882e-04,-0.3920e-04,-0.5035e-04,-0.3170e-04,-0.2169e-04,-0.3007e-04,-0.2740e-04, &
930  -0.5328e-04,-0.4491e-04,-0.4403e-04,-0.6383e-04, 0.4834e-04,-0.2702e-04,-0.4453e-04,-0.4339e-04,-0.4457e-04,-0.4551e-04, &
931  -0.8133e-04, 0.3768e-04,-0.7611e-06,-0.2626e-04,-0.4643e-04,-0.4305e-04,-0.4840e-04,-0.5149e-04, 0.7193e-04,-0.2169e-04, &
932  -0.4491e-04,-0.3996e-04,-0.4483e-04,-0.4487e-04,-0.6698e-04,-0.4834e-04,-0.3463e-04,-0.4986e-04,-0.4377e-04,-0.4514e-04, &
933  -0.5377e-04,-0.2626e-04,-0.4187e-04,-0.3692e-04,-0.5100e-04,-0.4651e-04,-0.4392e-04,-0.5386e-04,-0.4643e-04,-0.4301e-04, &
934  -0.3578e-04,-0.5176e-04,-0.4594e-04,-0.4551e-04,-0.3920e-04,-0.3425e-04,-0.4491e-04,-0.3654e-04,-0.5138e-04,-0.4377e-04, &
935  -0.5614e-04,-0.5758e-04,-0.3600e-04 /
936 
937 
938  data((c18h2o(t_1,j,i),i=1,kg(18)),j=1,numpir) / &
939  -0.2121e+02,-0.2002e+02,-0.1676e+02,-0.1274e+02,-0.8780e+01,-0.5167e+01,-0.2692e+01,-0.6275e+00,-0.2075e+02,-0.1996e+02, &
940  -0.1630e+02,-0.1228e+02,-0.8324e+01,-0.4718e+01,-0.2260e+01,-0.2303e+00,-0.2029e+02,-0.1990e+02,-0.1584e+02,-0.1182e+02, &
941  -0.7868e+01,-0.4269e+01,-0.1806e+01, 0.1645e+00,-0.2022e+02,-0.1985e+02,-0.1538e+02,-0.1136e+02,-0.7417e+01,-0.3820e+01, &
942  -0.1373e+01, 0.5657e+00,-0.2018e+02,-0.1981e+02,-0.1492e+02,-0.1090e+02,-0.6965e+01,-0.3369e+01,-0.9319e+00, 0.9577e+00, &
943  -0.2013e+02,-0.1937e+02,-0.1446e+02,-0.1044e+02,-0.6512e+01,-0.2917e+01,-0.4928e+00, 0.1376e+01,-0.2009e+02,-0.1891e+02, &
944  -0.1400e+02,-0.9984e+01,-0.6063e+01,-0.2466e+01,-0.6887e-01, 0.1768e+01,-0.2006e+02,-0.1845e+02,-0.1354e+02,-0.9530e+01, &
945  -0.5618e+01,-0.2024e+01, 0.3615e+00, 0.2196e+01,-0.2003e+02,-0.1800e+02,-0.1308e+02,-0.9075e+01,-0.5174e+01,-0.1593e+01, &
946  0.7820e+00, 0.2600e+01,-0.1827e+02,-0.1464e+02,-0.1097e+02,-0.7525e+01,-0.3733e+01,-0.1077e+01, 0.1204e+01, 0.3014e+01, &
947  -0.1525e+02,-0.1210e+02,-0.9275e+01,-0.5876e+01,-0.2768e+01,-0.6286e+00, 0.1622e+01, 0.3394e+01,-0.1298e+02,-0.1060e+02, &
948  -0.7764e+01,-0.4462e+01,-0.2154e+01,-0.2001e+00, 0.2034e+01, 0.3756e+01,-0.1157e+02,-0.8941e+01,-0.5984e+01,-0.3509e+01, &
949  -0.1651e+01, 0.2279e+00, 0.2422e+01, 0.4066e+01,-0.9986e+01,-0.7062e+01,-0.4794e+01,-0.2818e+01,-0.1196e+01, 0.6394e+00, &
950  0.2791e+01, 0.4283e+01,-0.8064e+01,-0.5512e+01,-0.3933e+01,-0.2274e+01,-0.7559e+00, 0.1036e+01, 0.3085e+01, 0.4444e+01, &
951  -0.6440e+01,-0.4863e+01,-0.3219e+01,-0.1791e+01,-0.3279e+00, 0.1427e+01, 0.3304e+01, 0.4527e+01,-0.5902e+01,-0.4207e+01, &
952  -0.2756e+01,-0.1350e+01, 0.7686e-01, 0.1776e+01, 0.3475e+01, 0.4550e+01,-0.5439e+01,-0.3739e+01,-0.2330e+01,-0.9233e+00, &
953  0.4612e+00, 0.2066e+01, 0.3564e+01, 0.4502e+01,-0.5006e+01,-0.3316e+01,-0.1906e+01,-0.5066e+00, 0.8352e+00, 0.2272e+01, &
954  0.3587e+01, 0.4419e+01 /
955 
956  data((c18h2o(t_2,j,i),i=1,kg(18)),j=1,numpir) / &
957  0.2338e-01, 0.1968e-02, 0.9503e-02, 0.3412e-02, 0.6280e-03,-0.1109e-02,-0.1089e-02,-0.1026e-02, 0.1972e-01, 0.2093e-02, &
958  0.9503e-02, 0.3391e-02, 0.6489e-03,-0.1172e-02,-0.1164e-02,-0.1158e-02, 0.1603e-01, 0.3328e-02, 0.9524e-02, 0.3391e-02, &
959  0.6489e-03,-0.1277e-02,-0.1229e-02,-0.1296e-02, 0.1229e-01, 0.7138e-02, 0.9524e-02, 0.3370e-02, 0.6070e-03,-0.1319e-02, &
960  -0.1264e-02,-0.1610e-02, 0.8478e-02, 0.1095e-01, 0.9566e-02, 0.3412e-02, 0.5652e-03,-0.1382e-02,-0.1266e-02,-0.1566e-02, &
961  0.4563e-02, 0.1480e-01, 0.9566e-02, 0.3412e-02, 0.5443e-03,-0.1423e-02,-0.1199e-02,-0.1679e-02, 0.2261e-02, 0.1865e-01, &
962  0.9608e-02, 0.3454e-02, 0.4815e-03,-0.1423e-02,-0.1296e-02,-0.1555e-02, 0.2198e-02, 0.2250e-01, 0.9671e-02, 0.3412e-02, &
963  0.4187e-03,-0.1426e-02,-0.1472e-02,-0.1800e-02, 0.2072e-02, 0.2600e-01, 0.9734e-02, 0.3433e-02, 0.3977e-03,-0.1428e-02, &
964  -0.1541e-02,-0.1591e-02, 0.1987e-01, 0.8645e-02, 0.6280e-02, 0.1298e-02,-0.1151e-02,-0.1509e-02,-0.1662e-02,-0.1570e-02, &
965  0.4668e-02, 0.8373e-02, 0.3956e-02,-0.4187e-04,-0.1968e-02,-0.1624e-02,-0.1700e-02,-0.1947e-02, 0.9231e-02, 0.5694e-02, &
966  0.1444e-02,-0.2512e-03,-0.1827e-02,-0.1662e-02,-0.1576e-02,-0.1633e-02, 0.8666e-02, 0.3077e-02,-0.1737e-02,-0.1277e-02, &
967  -0.1507e-02,-0.1757e-02,-0.1612e-02,-0.1612e-02, 0.8164e-03,-0.4375e-02,-0.1884e-02,-0.1277e-02,-0.1564e-02,-0.1853e-02, &
968  -0.1591e-02,-0.1486e-02,-0.1486e-02,-0.2596e-02,-0.1633e-02,-0.1539e-02,-0.1662e-02,-0.1846e-02,-0.1423e-02,-0.1277e-02, &
969  -0.1423e-02,-0.2617e-02,-0.1005e-02,-0.1379e-02,-0.1687e-02,-0.1905e-02,-0.1528e-02,-0.1298e-02,-0.1675e-03,-0.1947e-02, &
970  -0.5024e-03,-0.1325e-02,-0.1696e-02,-0.1698e-02,-0.1486e-02,-0.1277e-02, 0.1047e-03,-0.1109e-02,-0.5861e-03,-0.1363e-02, &
971  -0.1620e-02,-0.1666e-02,-0.1507e-02,-0.9210e-03, 0.1047e-03,-0.1047e-02,-0.8394e-03,-0.1342e-02,-0.1591e-02,-0.1323e-02, &
972  -0.1340e-02,-0.9420e-03 /
973 
974  data((c18h2o(t_3,j,i),i=1,kg(18)),j=1,numpir) / &
975  -0.1085e-03, 0.2283e-05,-0.4719e-04,-0.3807e-06,-0.1522e-05,-0.3425e-05,-0.7612e-06, 0.1751e-05,-0.1766e-03, 0.1523e-05, &
976  -0.4719e-04,-0.7609e-06,-0.3807e-06,-0.3045e-05, 0.1599e-05, 0.8723e-05,-0.2443e-03, 0.1941e-04,-0.4757e-04,-0.1522e-05, &
977  -0.3806e-06,-0.1903e-05,-0.2778e-05, 0.1294e-04,-0.1838e-03, 0.8563e-04,-0.4757e-04,-0.1903e-05, 0.1142e-05,-0.2664e-05, &
978  -0.6090e-06, 0.1321e-04,-0.1161e-03, 0.1526e-03,-0.4757e-04,-0.2664e-05,-0.3805e-06,-0.3806e-05,-0.2093e-05, 0.2253e-04, &
979  -0.4795e-04, 0.9248e-04,-0.4757e-04,-0.1903e-05, 0.0000e+00,-0.3045e-05,-0.7992e-06, 0.1393e-04,-0.9134e-05, 0.2246e-04, &
980  -0.4834e-04,-0.2664e-05, 0.3804e-06,-0.5328e-05,-0.1510e-05, 0.1465e-04,-0.1028e-04,-0.4757e-04,-0.4948e-04,-0.1142e-05, &
981  0.7614e-06,-0.4910e-05,-0.5709e-06, 0.1477e-04,-0.1256e-04,-0.1066e-03,-0.4910e-04,-0.1523e-05,-0.3805e-06,-0.3121e-05, &
982  -0.2512e-05, 0.1142e-04,-0.7878e-04,-0.2664e-05,-0.8373e-05,-0.7612e-06, 0.1104e-04,-0.3311e-05,-0.1979e-05, 0.5709e-05, &
983  -0.2626e-04,-0.4872e-04,-0.3808e-06,-0.2283e-05, 0.2284e-05,-0.3349e-05,-0.4034e-05, 0.7231e-05,-0.4910e-04, 0.1599e-04, &
984  0.1256e-04,-0.7612e-05, 0.1180e-05,-0.1815e-05,-0.7193e-05, 0.3045e-05, 0.1576e-09, 0.6470e-05,-0.1408e-04,-0.1903e-05, &
985  0.1522e-05,-0.4746e-05,-0.4948e-05, 0.3806e-06, 0.9020e-04, 0.5214e-04, 0.6090e-05,-0.1104e-04, 0.1180e-05,-0.2778e-05, &
986  -0.6090e-05,-0.2664e-05,-0.6737e-04,-0.1218e-04,-0.3806e-05,-0.5214e-05,-0.1066e-05,-0.1294e-05,-0.3045e-05,-0.2664e-05, &
987  -0.4643e-04, 0.1713e-04,-0.1218e-04,-0.6204e-05,-0.2360e-05,-0.1979e-05,-0.1903e-05,-0.3806e-05,-0.3045e-04,-0.1256e-04, &
988  -0.9134e-05,-0.6508e-05,-0.1027e-05,-0.7993e-06,-0.1142e-05,-0.7992e-05,-0.3616e-04,-0.1028e-04,-0.1066e-04,-0.6051e-05, &
989  0.1066e-05,-0.1751e-05,-0.2284e-05,-0.2284e-05,-0.3920e-04,-0.9895e-05,-0.1321e-04,-0.3844e-05,-0.2055e-05,-0.2512e-05, &
990  -0.3806e-05,-0.3425e-05 /
991 
992 
993 contains
994 
995 subroutine gases ( ncol , nlm , ib , ig &
996  , pp , dp , tt , rmix &
997  , o3mix , umco2 , umch4 , umn2o &
998  , hk , tg , pkd , ip1 &
999  , ip2 &
1000  )
1001 
1002 !-----------------------------------------------------------------------
1003 ! Calculates the non-gray spectral absorption due to H2O, O3, CO2,
1004 ! CHA4, and N2O at long and short wavelengths.
1005 
1006  ! INPUT ARGUMENTS:
1007  integer (kind=int_kind), intent(in):: &
1008  ncol &! Length of sub-domain.
1009  ,nlm &! Number of layers.
1010  ,ib &! Index of spectral interval.
1011  ,ig ! Index of k-distribution.
1012 
1013  integer (kind=int_kind), intent(in), dimension(ncol,nlm):: &
1014  ip1 &!Used in conjunction with pressure weighting.
1015  ,ip2 !Used in conjunction with pressure weighting.
1016 
1017  real (kind=dbl_kind), intent(in):: &
1018  umco2 &!Concentration of CO2 (ppm).
1019  ,umch4 &!Concentration of CH4 (???).
1020  ,umn2o !Concentration of N2O (???).
1021 
1022  real (kind=dbl_kind), intent(in), dimension(ncol,nlm):: &
1023  dp &!Layer Thickness (hPa).
1024  ,tt &!Layer Temperature (K).
1025  ,rmix &!Layer Water Vapor Mixing Ratio (kg/kg).
1026  ,o3mix &!Layer Ozone Mixing Ratio (kg/kg).
1027  ,pkd !Used in conjunction with pressure weighting (-).
1028 
1029  real (kind=dbl_kind), intent(in), dimension(ncol,nlm+1):: &
1030  pp !Level Pressure (hPa).
1031 
1032  ! OUTPUT ARGUMENTS:
1033  real (kind=dbl_kind), intent(out):: &
1034  hk !Weight (-).
1035 
1036  real (kind=dbl_kind), intent(out), dimension(ncol,nlm):: &
1037  tg !Transmission (-).
1038 
1039  ! LOCAL VARIABLES:
1040 
1041  real (kind=dbl_kind), dimension(ncol,nlm):: &
1042  fkg , fkga , fkgb , pq &
1043  ,tg1 , tg2 , tg3
1044 
1045 !-----------------------------------------------------------------------
1046  select case (ib)
1047  case (1)
1048  ! (50000-14500 cm^-1 ): nongray absorption by O3.
1049  ! 619.618 is the solar energy contained in that band (Wm^-2).
1050  call qopo3s(fk1o3(ig),dp,o3mix,tg)
1051  hk = 619.618*hk1(ig)
1052 
1053  case (2)
1054  ! (14500-7700 cm^-1): nongray absorption by H2O.
1055  ! 484.295 is the solar energy contained in the band (Wm^-2).
1056  call qk(ncol,nlm,c2h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1057  call qoph2o(fkg,dp,rmix,tg)
1058  hk = 484.295*hk2(ig)
1059 
1060  case (3)
1061  ! (7700-5250 cm^-1): nongray absorption by H2O.
1062  ! 149.845 is the solar energy contained in the band (Wm^-2).
1063  call qk(ncol,nlm,c3h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1064  call qoph2o(fkg,dp,rmix,tg)
1065  hk = 149.845*hk3(ig)
1066 
1067  case (4)
1068  ! (5250-4000 cm^-1): nongray absorption by H2O.
1069  ! 48.7302 is the solar energy contained in the band (Wm^-2).
1070  call qk(ncol,nlm,c4h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1071  call qoph2o(fkg,dp,rmix,tg)
1072  hk = 48.7302*hk4(ig)
1073 
1074  case (5)
1075  ! (4000-2850 cm^-1): nongray absorption by H2O.
1076  ! 31.6576 is the solar energy contained in the band (Wm^-2).
1077  call qk(ncol,nlm,c5h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1078  call qoph2o(fkg,dp,rmix,tg)
1079  hk = 31.6576*hk5(ig)
1080 
1081  case (6)
1082  ! (2850-2500 cm^-1): nongray absorption by H2O.
1083  ! 5.79927 is the solar energy contained in the band (Wm^-2).
1084  call qk(ncol,nlm,c6h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1085  call qoph2o(fkg,dp,rmix,tg)
1086  hk = 5.79927*hk6(ig)
1087 
1088  case (7)
1089  ! (2200-1900 cm^-1): nongray absorption by H2O.
1090  call qk(ncol,nlm,c7h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1091  call qoph2o(fkg,dp,rmix,tg)
1092  hk = hk7(ig)
1093 
1094  case (8)
1095  ! (1900-1700 cm^-1): nongray absorption by H2O.
1096  call qk(ncol,nlm,c8h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1097  call qoph2o(fkg,dp,rmix,tg)
1098  hk = hk8(ig)
1099 
1100  case (9)
1101  ! (1700-1400 cm^-1): nongray absorption by H2O.
1102  call qk(ncol,nlm,c9h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1103  call qoph2o(fkg,dp,rmix,tg)
1104  hk = hk9(ig)
1105 
1106  case (10)
1107  ! (1400-1250 cm^-1): overlapping absorption by H2O, CH4, and N2O
1108  ! using approach one of Fu(1991).
1109  call qk(ncol,nlm,c10h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1110  call qoph2o(fkg,dp,rmix,tg1)
1111  call qk(ncol,nlm,c10ch4,tt,fkg,pkd,ip1,ip2)
1112  call qopch4(fkg,dp,tg2)
1113  call qk(ncol,nlm,c10n2o,tt,fkg,pkd,ip1,ip2)
1114  call qopn2o(fkg,dp,tg3)
1115 
1116  tg = tg1 + tg2/1.6*umch4 + tg3/0.28*umn2o
1117  hk = hk10(ig)
1118 
1119  case (11)
1120  ! (1250-1100 cm^-1): overlapping absorption by H2O, CH4, and N2O
1121  ! using approach one of Fu(1991).
1122  call qk(ncol,nlm,c11h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1123  call qoph2o(fkg,dp,rmix,tg1)
1124  call qk(ncol,nlm,c11ch4,tt,fkg,pkd,ip1,ip2)
1125  call qopch4(fkg,dp,tg2)
1126  call qk(ncol,nlm,c11n2o,tt,fkg,pkd,ip1,ip2)
1127  call qopn2o(fkg,dp,tg3)
1128 
1129  tg = tg1 + tg2/1.6*umch4 + tg3/0.28*umn2o
1130  hk = hk11(ig)
1131 
1132  case(12)
1133  ! (1100-980 cm^-1): overlapping absorption by H2O and O3
1134  ! using approach one of Fu(1991).
1135  call qkio3(ncol,nlm,c12o3(:,:,ig),tt,fkg,pkd,ip1,ip2)
1136  call qopo3i(fkg,dp,o3mix,tg1)
1137  call qk(ncol,nlm,c12h2o,tt,fkg,pkd,ip1,ip2)
1138  call qoph2o(fkg,dp,rmix,tg2)
1139 
1140  tg = tg1 + tg2
1141  hk = hk12(ig)
1142 
1143  case (13)
1144  ! (980-800 cm^-1 ): nongray absorption by H2O.
1145  call qk(ncol,nlm,c13h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1146  call qoph2o(fkg,dp,rmix,tg)
1147  hk = hk13(ig)
1148 
1149  case (14)
1150  ! (800-670 cm^-1): overlapping absorption by H2O and CO2
1151  ! using approach two of Fu(1991).
1152 
1153  where (pp .ge. 63.1)
1154  pq = rmix
1155  elsewhere
1156  pq = 0._dbl_kind
1157  end where
1158 
1159  call qk(ncol,nlm,c14hca(:,:,ig),tt,fkga,pkd,ip1,ip2)
1160  call qk(ncol,nlm,c14hcb(:,:,ig),tt,fkgb,pkd,ip1,ip2)
1161 
1162  fkg = fkga/330.0*umco2+pq*fkgb
1163  call qophc(fkg,dp,tg)
1164  hk = hk14(ig)
1165 
1166  case (15)
1167  ! (670-540 cm^-1): overlapping absorption by H2O and CO2
1168  ! using approach two of Fu(1991).
1169 
1170  where (pp .ge. 63.1)
1171  pq = rmix
1172  elsewhere
1173  pq = 0._dbl_kind
1174  end where
1175  call qk(ncol,nlm,c15hca(:,:,ig),tt,fkga,pkd,ip1,ip2)
1176  call qk(ncol,nlm,c15hcb(:,:,ig),tt,fkgb,pkd,ip1,ip2)
1177 
1178  fkg = fkga/330.0*umco2+pq*fkgb
1179  call qophc(fkg,dp,tg)
1180  hk = hk15(ig)
1181 
1182  case (16)
1183  ! (540-400 cm^-1 ): nongray absorption by H2O.
1184  call qk(ncol,nlm,c16h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1185  call qoph2o(fkg,dp,rmix,tg)
1186  hk = hk16(ig)
1187 
1188  case (17)
1189  ! (400-280 cm^-1 ): nongray absorption by H2O.
1190  call qk(ncol,nlm,c17h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1191  call qoph2o(fkg,dp,rmix,tg)
1192  hk = hk17(ig)
1193 
1194  case (18)
1195  ! (280-000 cm^-1 ): nongray absorption by H2O.
1196  call qk(ncol,nlm,c18h2o(:,:,ig),tt,fkg,pkd,ip1,ip2)
1197  call qoph2o(fkg,dp,rmix,tg)
1198  hk = hk18(ig)
1199 
1200  case default
1201  stop
1202  end select
1203 
1204  return
1205 end subroutine gases
1206 
1207 
1208 subroutine qopch4 (fkg, dp, tg)
1209 !-----------------------------------------------------------------------
1210 ! Absorption by CH4 in the 1400-1250cm^-1 and 1250-1100cm^-1
1211 ! spectral intervals.
1212 
1213  use bugsrad_physconst, only: molar_volume,gravity,mw_dry_air
1214  implicit none
1215 
1216 ! ARGUMENT LIST VARIABLES:
1217 ! INPUT ARGUMENTS:
1218 ! ----------------
1219  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1220  dp &!Layer Thickness (hPa).
1221  ,fkg !Gaseous absorption coefficients (cm-atm)^-1.
1222 
1223 ! OUTPUT ARGUMENTS:
1224 ! -----------------
1225  real (kind=dbl_kind), intent(out), dimension(:,:):: &
1226  tg !Transmission (-).
1227 
1228  real (kind=dbl_kind) :: &
1229  ch4_conc !ppv
1230 
1231  ch4_conc = 1.6e-6_dbl_kind
1232  tg = fkg*dp*molar_volume*10.0/gravity*ch4_conc/mw_dry_air
1233 
1234  return
1235 end subroutine qopch4
1236 
1237 
1238 subroutine qopn2o (fkg , dp, tg)
1239 !-----------------------------------------------------------------------
1240 ! Absorption by N2O in the 1400-1250cm^-1 and 1250-1100cm^-1
1241 ! spectral intervals.
1242 
1243  use bugsrad_physconst, only: molar_volume, gravity, mw_dry_air
1244  implicit none
1245 
1246  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1247  dp & !Layer Thickness (hPa).
1248  ,fkg !Gaseous absorption coefficients (cm-atm)^-1.
1249 
1250  ! OUTPUT ARGUMENTS:
1251  real (kind=dbl_kind), intent(out), dimension(:,:):: &
1252  tg !Transmission (-).
1253 
1254  real (kind = dbl_kind):: &
1255  n2o_conc !ppv
1256 
1257  n2o_conc = 0.28e-6_dbl_kind
1258  tg = fkg*dp*molar_volume*10.0/gravity*n2o_conc/mw_dry_air
1259 
1260  return
1261 end subroutine qopn2o
1262 
1263 
1264 subroutine qoph2o (fkg , dp , rmix , tg)
1265 !-----------------------------------------------------------------------
1266 ! Absorption by H2O in each SW and LW spectral interval.
1267 
1268  use bugsrad_physconst, only: molar_volume, gravity, mw_h2o
1269  implicit none
1270 
1271  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1272  dp & !Layer Thickness (hPa).
1273  ,rmix & !Water Vapor Mixing Ratio (kg/kg).
1274  ,fkg !Gaseous absorption coefficients (cm-atm)^-1.
1275 
1276  ! OUTPUT ARGUMENTS:
1277  real (kind=dbl_kind), intent(out), dimension(:,:):: &
1278  tg !Transmisssion (-).
1279 
1280  tg = fkg*rmix*dp*molar_volume/mw_h2o*10.0/gravity
1281 
1282  return
1283 end subroutine qoph2o
1284 
1285 
1286 subroutine qopo3i (fkg , dp , o3mix , tg )
1287 !-----------------------------------------------------------------------
1288 ! Absorption by O3 in the infrared band (1100-980 cm^-1).
1289 
1290  use bugsrad_physconst, only: molar_volume, gravity, mw_o3
1291  implicit none
1292 
1293  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1294  dp & !Layer Thickness (hPa).
1295  ,o3mix & !Ozone Mixing Ratio (kg/kg).
1296  ,fkg !Gaseous absorption coefficients (cm-atm)^-1.
1297 
1298  ! OUTPUT ARGUMENTS:
1299  real (kind=dbl_kind), intent(out), dimension(:,:):: &
1300  tg !Transmission (-).
1301 
1302  tg = fkg*o3mix*dp*molar_volume/mw_o3*10.0/gravity
1303 
1304  return
1305 end subroutine qopo3i
1306 
1307 
1308 subroutine qopo3s (fk , dp , o3mix , tg )
1309 !-----------------------------------------------------------------------
1310 ! Computes non-gray absorption of O3.
1311 
1312  use bugsrad_physconst, only: molar_volume, gravity, mw_o3
1313  implicit none
1314 
1315  real (kind=dbl_kind), intent(in):: &
1316  fk !
1317 
1318  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1319  dp & !Layer Thickness (hPa).
1320  ,o3mix !Ozone Mixing Ratio (kg/kg).
1321 
1322  ! OUTPUT ARGUMENTS:
1323  real (kind=dbl_kind), intent(out), dimension(:,:):: &
1324  tg !Transmission (-).
1325 
1326  tg = o3mix*dp*fk*molar_volume/mw_o3*10.0/gravity
1327 
1328  return
1329 end subroutine qopo3s
1330 
1331 
1332 subroutine qophc(fkg,dp,tg)
1333 !-----------------------------------------------------------------------
1334 ! Overlapping of H2O and CO2 bands.
1335 
1336  implicit none
1337 
1338  ! INPUT ARGUMENTS
1339  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1340  dp & !Layer Thickness (hPa).
1341  ,fkg !Gaseous absorption coefficients (cm-atm)^-1.
1342 
1343  ! OUTPUT ARGUMENTS:
1344  real (kind=dbl_kind), intent(out), dimension(:,:):: &
1345  tg !Transmission (-).
1346 
1347  tg = fkg*dp
1348  return
1349 end subroutine qophc
1350 
1351 subroutine qk(ncol,nlm,coefk,tt,fkg,pkd,ip1,ip2)
1352 !-----------------------------------------------------------------------
1353 ! Interpolates the correlated k coeffients to the temperatures tt
1354 ! using the form
1355 ! ln k = a + b * ( t - 245 ) + c * ( t - 245 ) ** 2
1356 ! and linearly interpolates in pressure if needed.
1357 
1358  implicit none
1359 
1360  ! INPUT ARGUMENTS:
1361  integer (kind=int_kind), intent(in):: &
1362  ncol & !Length of sub-domain.
1363  ,nlm !Number of layers.
1364 
1365  integer (kind=int_kind), intent(in), dimension(:,:):: &
1366  ip1 & !Used in conjunction with pressure weigthing.
1367  ,ip2 !Used in conjunction with pressure weigthing.
1368 
1369  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1370  coefk !Pre-computed coefficients.
1371 
1372  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1373  tt & !Layer temperature (K).
1374  ,pkd !Scaled pressure (hPa).
1375 
1376  ! OUTPUT ARGUMENTS:
1377  real (kind=dbl_kind), intent(out), dimension(:,:):: &
1378  fkg !Gaseous absorption coefficients (cm-atm)^-1.
1379 
1380  ! LOCAL VARIABLES:
1381  integer (kind=int_kind):: &
1382  i, l, i1, i2
1383 
1384  real (kind=dbl_kind):: &
1385  x1,x2,y1
1386 
1387 
1388  do l = 1, nlm
1389  do i = 1, ncol
1390  i1=ip1(i,l)
1391  i2=ip2(i,l)
1392  y1 = tt(i,l) - 245.0
1393  if (y1.lt.-65.) y1 = -65.
1394  if (y1.gt.75.) y1 = 75.
1395  x1 = exp(coefk(1,i1)+coefk(2,i1)*y1+coefk(3,i1)*y1**2)
1396  if(i1.eq.i2) then
1397  fkg(i,l) = x1 * pkd(i,l)
1398  else
1399  x2 = exp(coefk(1,i2)+coefk(2,i2)*y1+coefk(3,i2)*y1**2)
1400  fkg(i,l) = x1 + ( x2 - x1 ) * pkd(i,l)
1401  endif
1402  enddo
1403  enddo
1404 
1405  return
1406 end subroutine qk
1407 
1408 subroutine qkio3 (ncol , nlm , coefk , tt , fkg , pkd , ip1 , ip2)
1409 ! Computes the gaseous absorption coefficients in units of (cm-atm)^-1
1410 ! for a given cumulative probability in nlm layers. coefk are
1411 ! the coefficients used to calculate the absorption coefficient at the
1412 ! temperature t for the 19 pressures by
1413 ! ln k = a + b * ( t - 250 ) + c * ( t - 250 ) ** 2
1414 ! and the absorption coefficient at conditions other than those nineteen
1415 ! pressures is interpolated linearly with pressure (Fu, 1991).
1416 
1417  implicit none
1418 
1419  ! INPUT ARGUMENTS:
1420  integer (kind=int_kind), intent(in):: &
1421  ncol & !Length of sub-domain.
1422  ,nlm !Number of layers.
1423 
1424  integer (kind=int_kind), intent(in), dimension(:,:):: &
1425  ip1 & !Used in conjunction with pressure weigthing.
1426  ,ip2 !Used in conjunction with pressure weigthing.
1427 
1428  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1429  coefk !Pre-computed coefficients.
1430 
1431  real (kind=dbl_kind), intent(in), dimension(:,:):: &
1432  tt & !Layer temperature (K).
1433  ,pkd !Scaled pressure (hPa).
1434 
1435 
1436  ! OUTPUT ARGUMENTS:
1437  real (kind=dbl_kind), intent(out), dimension(:,:):: &
1438  fkg !Gaseous absorption coefficients (cm-atm)^-1.
1439 
1440 
1441  ! LOCAL VARIABLES:
1442  integer (kind=int_kind):: &
1443  i, l, i1, i2
1444 
1445  real (kind=dbl_kind):: &
1446  x1,x2,y1
1447 
1448  do l = 1, nlm
1449  do i = 1, ncol
1450  i1 = ip1(i,l)
1451  i2 =ip2(i,l)
1452  y1 = tt(i,l) - 250.0
1453  if (y1.lt.-70.) y1 = -70.
1454  if (y1.gt.75.) y1 = 70.
1455  x1 = exp(coefk(1,i1)+coefk(2,i1)*y1+coefk(3,i1)*y1**2)
1456  if(i1.eq.i2) then
1457  fkg(i,l) = x1 * pkd(i,l)
1458  else
1459  x2 = exp(coefk(1,i2)+coefk(2,i2)*y1+coefk(3,i2)*y1**2)
1460  fkg(i,l) = x1 + ( x2 - x1 ) * pkd(i,l)
1461  endif
1462  enddo
1463  enddo
1464 
1465  return
1466 end subroutine qkio3
1467 
1468 subroutine pscale (ncol , nlm , ppl, stanp, pkd , ip1 , ip2 )
1469 !-----------------------------------------------------------------------
1470 ! pscale is used to scale pressure for use in the k-distribu
1471 ! tion routine.
1472 
1473  implicit none
1474 
1475  ! INPUT ARGUMENTS:
1476  integer (kind=int_kind), intent(in):: &
1477  ncol & !Length of sub-domain.
1478  ,nlm !Number of layers.
1479 
1480  real (kind=dbl_kind), dimension(:,:):: &
1481  ppl !Pressure (hPa).
1482 
1483  real (kind=dbl_kind), dimension(:):: &
1484  stanp
1485 
1486  ! OUTPUT ARGUMENTS:
1487  integer (kind=int_kind), dimension(:,:):: &
1488  ip1 & !pointer.
1489  ,ip2 !pointer.
1490 
1491  real (kind=dbl_kind), dimension(:,:):: &
1492  pkd !Weighted pressure (hPa).
1493 
1494  ! LOCAL VARIABLES:
1495  integer (kind=int_kind):: &
1496  n_stanp
1497 
1498  integer (kind=int_kind):: &
1499  i1, icol, l
1500 
1501  n_stanp = size(stanp)
1502 
1503  do icol = 1, ncol
1504  i1 = 1
1505  do l = 1, nlm
1506  if( ppl(icol,l) .lt. stanp(1) ) then
1507  pkd(icol,l) = ppl(icol,l) / stanp(1)
1508  ip1(icol,l) = 1
1509  ip2(icol,l) = 1
1510  elseif ( ppl(icol,l) .ge. stanp(n_stanp) ) then
1511  pkd(icol,l) = (ppl(icol,l) - stanp(n_stanp-1)) &
1512  / (stanp(n_stanp) - stanp(n_stanp-1))
1513  ip1(icol,l) = n_stanp - 1
1514  ip2(icol,l) = n_stanp
1515  else
1516 30 continue
1517  if ( ppl(icol,l) .ge. stanp(i1) ) goto 20
1518  pkd(icol,l) = (ppl(icol,l)-stanp(i1-1)) &
1519  / (stanp(i1)-stanp(i1-1))
1520  ip1(icol,l) = i1-1
1521  ip2(icol,l) = i1
1522  goto 5
1523 20 i1 = i1 + 1
1524  goto 30
1525  endif
1526 5 enddo
1527  enddo
1528 
1529  return
1530 end subroutine pscale
1531 
1532 
1533 !-----------------------------------------------------------------------
1534 end module gases_ckd
Generated on Tue Aug 13 2013 10:38:25 for BUGSrad by   doxygen 1.8.3.1