1 | ! path: $Source: /storm/rc1/cvsroot/rc/rrtmg_lw/src/rrtmg_lw_init.f90,v $ |
---|
2 | ! author: $Author: miacono $ |
---|
3 | ! revision: $Revision: 1.6 $ |
---|
4 | ! created: $Date: 2011/04/08 20:25:01 $ |
---|
5 | ! |
---|
6 | module rrtmg_lw_init |
---|
7 | |
---|
8 | ! -------------------------------------------------------------------------- |
---|
9 | ! | | |
---|
10 | ! | Copyright 2002-2009, Atmospheric & Environmental Research, Inc. (AER). | |
---|
11 | ! | This software may be used, copied, or redistributed as long as it is | |
---|
12 | ! | not sold and this copyright notice is reproduced on each copy made. | |
---|
13 | ! | This model is provided as is without any express or implied warranties. | |
---|
14 | ! | (http://www.rtweb.aer.com/) | |
---|
15 | ! | | |
---|
16 | ! -------------------------------------------------------------------------- |
---|
17 | |
---|
18 | ! ------- Modules ------- |
---|
19 | use parkind, only : im => kind_im, rb => kind_rb |
---|
20 | use rrlw_wvn |
---|
21 | use rrtmg_lw_setcoef, only: lwatmref, lwavplank, lwavplankderiv |
---|
22 | |
---|
23 | implicit none |
---|
24 | |
---|
25 | contains |
---|
26 | |
---|
27 | ! ************************************************************************** |
---|
28 | subroutine rrtmg_lw_ini(cpdair) |
---|
29 | ! ************************************************************************** |
---|
30 | ! |
---|
31 | ! Original version: Michael J. Iacono; July, 1998 |
---|
32 | ! First revision for GCMs: September, 1998 |
---|
33 | ! Second revision for RRTM_V3.0: September, 2002 |
---|
34 | ! |
---|
35 | ! This subroutine performs calculations necessary for the initialization |
---|
36 | ! of the longwave model. Lookup tables are computed for use in the LW |
---|
37 | ! radiative transfer, and input absorption coefficient data for each |
---|
38 | ! spectral band are reduced from 256 g-point intervals to 140. |
---|
39 | ! ************************************************************************** |
---|
40 | |
---|
41 | use parrrtm, only : mg, nbndlw, ngptlw |
---|
42 | use rrlw_tbl, only: ntbl, tblint, pade, bpade, tau_tbl, exp_tbl, tfn_tbl |
---|
43 | use rrlw_vsn, only: hvrini, hnamini |
---|
44 | |
---|
45 | real(kind=rb), intent(in) :: cpdair ! Specific heat capacity of dry air |
---|
46 | ! at constant pressure at 273 K |
---|
47 | ! (J kg-1 K-1) |
---|
48 | |
---|
49 | ! ------- Local ------- |
---|
50 | |
---|
51 | integer(kind=im) :: itr, ibnd, igc, ig, ind, ipr |
---|
52 | integer(kind=im) :: igcsm, iprsm |
---|
53 | |
---|
54 | real(kind=rb) :: wtsum, wtsm(mg) ! |
---|
55 | real(kind=rb) :: tfn ! |
---|
56 | |
---|
57 | real(kind=rb), parameter :: expeps = 1.e-20_rb ! Smallest value for exponential table |
---|
58 | |
---|
59 | ! ------- Definitions ------- |
---|
60 | ! Arrays for 10000-point look-up tables: |
---|
61 | ! TAU_TBL Clear-sky optical depth (used in cloudy radiative transfer) |
---|
62 | ! EXP_TBL Exponential lookup table for ransmittance |
---|
63 | ! TFN_TBL Tau transition function; i.e. the transition of the Planck |
---|
64 | ! function from that for the mean layer temperature to that for |
---|
65 | ! the layer boundary temperature as a function of optical depth. |
---|
66 | ! The "linear in tau" method is used to make the table. |
---|
67 | ! PADE Pade approximation constant (= 0.278) |
---|
68 | ! BPADE Inverse of the Pade approximation constant |
---|
69 | ! |
---|
70 | |
---|
71 | hvrini = '$Revision: 1.6 $' |
---|
72 | |
---|
73 | ! Initialize model data |
---|
74 | call lwdatinit(cpdair) |
---|
75 | call lwcmbdat ! g-point interval reduction data |
---|
76 | call lwcldpr ! cloud optical properties |
---|
77 | call lwatmref ! reference MLS profile |
---|
78 | call lwavplank ! Planck function |
---|
79 | call lwavplankderiv ! Planck function derivative wrt temp |
---|
80 | call lw_kgb01 ! molecular absorption coefficients |
---|
81 | call lw_kgb02 |
---|
82 | call lw_kgb03 |
---|
83 | call lw_kgb04 |
---|
84 | call lw_kgb05 |
---|
85 | call lw_kgb06 |
---|
86 | call lw_kgb07 |
---|
87 | call lw_kgb08 |
---|
88 | call lw_kgb09 |
---|
89 | call lw_kgb10 |
---|
90 | call lw_kgb11 |
---|
91 | call lw_kgb12 |
---|
92 | call lw_kgb13 |
---|
93 | call lw_kgb14 |
---|
94 | call lw_kgb15 |
---|
95 | call lw_kgb16 |
---|
96 | |
---|
97 | ! Compute lookup tables for transmittance, tau transition function, |
---|
98 | ! and clear sky tau (for the cloudy sky radiative transfer). Tau is |
---|
99 | ! computed as a function of the tau transition function, transmittance |
---|
100 | ! is calculated as a function of tau, and the tau transition function |
---|
101 | ! is calculated using the linear in tau formulation at values of tau |
---|
102 | ! above 0.01. TF is approximated as tau/6 for tau < 0.01. All tables |
---|
103 | ! are computed at intervals of 0.001. The inverse of the constant used |
---|
104 | ! in the Pade approximation to the tau transition function is set to b. |
---|
105 | |
---|
106 | tau_tbl(0) = 0.0_rb |
---|
107 | tau_tbl(ntbl) = 1.e10_rb |
---|
108 | exp_tbl(0) = 1.0_rb |
---|
109 | exp_tbl(ntbl) = expeps |
---|
110 | tfn_tbl(0) = 0.0_rb |
---|
111 | tfn_tbl(ntbl) = 1.0_rb |
---|
112 | bpade = 1.0_rb / pade |
---|
113 | do itr = 1, ntbl-1 |
---|
114 | tfn = real(itr) / real(ntbl) |
---|
115 | tau_tbl(itr) = bpade * tfn / (1._rb - tfn) |
---|
116 | exp_tbl(itr) = exp(-tau_tbl(itr)) |
---|
117 | if (exp_tbl(itr) .le. expeps) exp_tbl(itr) = expeps |
---|
118 | if (tau_tbl(itr) .lt. 0.06_rb) then |
---|
119 | tfn_tbl(itr) = tau_tbl(itr)/6._rb |
---|
120 | else |
---|
121 | tfn_tbl(itr) = 1._rb-2._rb*((1._rb/tau_tbl(itr))-(exp_tbl(itr)/(1.-exp_tbl(itr)))) |
---|
122 | endif |
---|
123 | enddo |
---|
124 | |
---|
125 | ! Perform g-point reduction from 16 per band (256 total points) to |
---|
126 | ! a band dependant number (140 total points) for all absorption |
---|
127 | ! coefficient input data and Planck fraction input data. |
---|
128 | ! Compute relative weighting for new g-point combinations. |
---|
129 | |
---|
130 | igcsm = 0 |
---|
131 | do ibnd = 1,nbndlw |
---|
132 | iprsm = 0 |
---|
133 | if (ngc(ibnd).lt.mg) then |
---|
134 | do igc = 1,ngc(ibnd) |
---|
135 | igcsm = igcsm + 1 |
---|
136 | wtsum = 0._rb |
---|
137 | do ipr = 1, ngn(igcsm) |
---|
138 | iprsm = iprsm + 1 |
---|
139 | wtsum = wtsum + wt(iprsm) |
---|
140 | enddo |
---|
141 | wtsm(igc) = wtsum |
---|
142 | enddo |
---|
143 | do ig = 1, ng(ibnd) |
---|
144 | ind = (ibnd-1)*mg + ig |
---|
145 | rwgt(ind) = wt(ig)/wtsm(ngm(ind)) |
---|
146 | enddo |
---|
147 | else |
---|
148 | do ig = 1, ng(ibnd) |
---|
149 | igcsm = igcsm + 1 |
---|
150 | ind = (ibnd-1)*mg + ig |
---|
151 | rwgt(ind) = 1.0_rb |
---|
152 | enddo |
---|
153 | endif |
---|
154 | enddo |
---|
155 | |
---|
156 | ! Reduce g-points for absorption coefficient data in each LW spectral band. |
---|
157 | |
---|
158 | call cmbgb1 |
---|
159 | call cmbgb2 |
---|
160 | call cmbgb3 |
---|
161 | call cmbgb4 |
---|
162 | call cmbgb5 |
---|
163 | call cmbgb6 |
---|
164 | call cmbgb7 |
---|
165 | call cmbgb8 |
---|
166 | call cmbgb9 |
---|
167 | call cmbgb10 |
---|
168 | call cmbgb11 |
---|
169 | call cmbgb12 |
---|
170 | call cmbgb13 |
---|
171 | call cmbgb14 |
---|
172 | call cmbgb15 |
---|
173 | call cmbgb16 |
---|
174 | |
---|
175 | end subroutine rrtmg_lw_ini |
---|
176 | |
---|
177 | !*************************************************************************** |
---|
178 | subroutine lwdatinit(cpdair) |
---|
179 | !*************************************************************************** |
---|
180 | |
---|
181 | ! --------- Modules ---------- |
---|
182 | |
---|
183 | use parrrtm, only : maxxsec, maxinpx |
---|
184 | use rrlw_con, only: heatfac, grav, planck, boltz, & |
---|
185 | clight, avogad, alosmt, gascon, radcn1, radcn2, & |
---|
186 | sbcnst, secdy |
---|
187 | use rrlw_vsn |
---|
188 | |
---|
189 | save |
---|
190 | |
---|
191 | real(kind=rb), intent(in) :: cpdair ! Specific heat capacity of dry air |
---|
192 | ! at constant pressure at 273 K |
---|
193 | ! (J kg-1 K-1) |
---|
194 | |
---|
195 | ! Longwave spectral band limits (wavenumbers) |
---|
196 | wavenum1(:) = (/ 10._rb, 350._rb, 500._rb, 630._rb, 700._rb, 820._rb, & |
---|
197 | 980._rb,1080._rb,1180._rb,1390._rb,1480._rb,1800._rb, & |
---|
198 | 2080._rb,2250._rb,2380._rb,2600._rb/) |
---|
199 | wavenum2(:) = (/350._rb, 500._rb, 630._rb, 700._rb, 820._rb, 980._rb, & |
---|
200 | 1080._rb,1180._rb,1390._rb,1480._rb,1800._rb,2080._rb, & |
---|
201 | 2250._rb,2380._rb,2600._rb,3250._rb/) |
---|
202 | delwave(:) = (/340._rb, 150._rb, 130._rb, 70._rb, 120._rb, 160._rb, & |
---|
203 | 100._rb, 100._rb, 210._rb, 90._rb, 320._rb, 280._rb, & |
---|
204 | 170._rb, 130._rb, 220._rb, 650._rb/) |
---|
205 | |
---|
206 | ! Spectral band information |
---|
207 | ng(:) = (/16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16/) |
---|
208 | nspa(:) = (/1,1,9,9,9,1,9,1,9,1,1,9,9,1,9,9/) |
---|
209 | nspb(:) = (/1,1,5,5,5,0,1,1,1,1,1,0,0,1,0,0/) |
---|
210 | |
---|
211 | ! nxmol - number of cross-sections input by user |
---|
212 | ! ixindx(i) - index of cross-section molecule corresponding to Ith |
---|
213 | ! cross-section specified by user |
---|
214 | ! = 0 -- not allowed in rrtm |
---|
215 | ! = 1 -- ccl4 |
---|
216 | ! = 2 -- cfc11 |
---|
217 | ! = 3 -- cfc12 |
---|
218 | ! = 4 -- cfc22 |
---|
219 | nxmol = 4 |
---|
220 | ixindx(1) = 1 |
---|
221 | ixindx(2) = 2 |
---|
222 | ixindx(3) = 3 |
---|
223 | ixindx(4) = 4 |
---|
224 | ixindx(5:maxinpx) = 0 |
---|
225 | |
---|
226 | ! Fundamental physical constants from NIST 2002 |
---|
227 | |
---|
228 | grav = 9.8066_rb ! Acceleration of gravity |
---|
229 | ! (m s-2) |
---|
230 | planck = 6.62606876e-27_rb ! Planck constant |
---|
231 | ! (ergs s; g cm2 s-1) |
---|
232 | boltz = 1.3806503e-16_rb ! Boltzmann constant |
---|
233 | ! (ergs K-1; g cm2 s-2 K-1) |
---|
234 | clight = 2.99792458e+10_rb ! Speed of light in a vacuum |
---|
235 | ! (cm s-1) |
---|
236 | avogad = 6.02214199e+23_rb ! Avogadro constant |
---|
237 | ! (mol-1) |
---|
238 | alosmt = 2.6867775e+19_rb ! Loschmidt constant |
---|
239 | ! (cm-3) |
---|
240 | gascon = 8.31447200e+07_rb ! Molar gas constant |
---|
241 | ! (ergs mol-1 K-1) |
---|
242 | radcn1 = 1.191042722e-12_rb ! First radiation constant |
---|
243 | ! (W cm2 sr-1) |
---|
244 | radcn2 = 1.4387752_rb ! Second radiation constant |
---|
245 | ! (cm K) |
---|
246 | sbcnst = 5.670400e-04_rb ! Stefan-Boltzmann constant |
---|
247 | ! (W cm-2 K-4) |
---|
248 | secdy = 8.6400e4_rb ! Number of seconds per day |
---|
249 | ! (s d-1) |
---|
250 | ! |
---|
251 | ! units are generally cgs |
---|
252 | ! |
---|
253 | ! The first and second radiation constants are taken from NIST. |
---|
254 | ! They were previously obtained from the relations: |
---|
255 | ! radcn1 = 2.*planck*clight*clight*1.e-07 |
---|
256 | ! radcn2 = planck*clight/boltz |
---|
257 | |
---|
258 | ! Heatfac is the factor by which delta-flux / delta-pressure is |
---|
259 | ! multiplied, with flux in W/m-2 and pressure in mbar, to get |
---|
260 | ! the heating rate in units of degrees/day. It is equal to: |
---|
261 | ! Original value: |
---|
262 | ! (g)x(#sec/day)x(1e-5)/(specific heat of air at const. p) |
---|
263 | ! Here, cpdair (1.004) is in units of J g-1 K-1, and the |
---|
264 | ! constant (1.e-5) converts mb to Pa and g-1 to kg-1. |
---|
265 | ! = (9.8066)(86400)(1e-5)/(1.004) |
---|
266 | ! heatfac = 8.4391_rb |
---|
267 | ! |
---|
268 | ! Modified value for consistency with CAM3: |
---|
269 | ! (g)x(#sec/day)x(1e-5)/(specific heat of air at const. p) |
---|
270 | ! Here, cpdair (1.00464) is in units of J g-1 K-1, and the |
---|
271 | ! constant (1.e-5) converts mb to Pa and g-1 to kg-1. |
---|
272 | ! = (9.80616)(86400)(1e-5)/(1.00464) |
---|
273 | ! heatfac = 8.43339130434_rb |
---|
274 | ! |
---|
275 | ! Calculated value: |
---|
276 | ! (grav) x (#sec/day) / (specific heat of dry air at const. p x 1.e2) |
---|
277 | ! Here, cpdair is in units of J kg-1 K-1, and the constant (1.e2) |
---|
278 | ! converts mb to Pa when heatfac is multiplied by W m-2 mb-1. |
---|
279 | heatfac = grav * secdy / (cpdair * 1.e2_rb) |
---|
280 | |
---|
281 | end subroutine lwdatinit |
---|
282 | |
---|
283 | !*************************************************************************** |
---|
284 | subroutine lwcmbdat |
---|
285 | !*************************************************************************** |
---|
286 | |
---|
287 | save |
---|
288 | |
---|
289 | ! ------- Definitions ------- |
---|
290 | ! Arrays for the g-point reduction from 256 to 140 for the 16 LW bands: |
---|
291 | ! This mapping from 256 to 140 points has been carefully selected to |
---|
292 | ! minimize the effect on the resulting fluxes and cooling rates, and |
---|
293 | ! caution should be used if the mapping is modified. The full 256 |
---|
294 | ! g-point set can be restored with ngptlw=256, ngc=16*16, ngn=256*1., etc. |
---|
295 | ! ngptlw The total number of new g-points |
---|
296 | ! ngc The number of new g-points in each band |
---|
297 | ! ngs The cumulative sum of new g-points for each band |
---|
298 | ! ngm The index of each new g-point relative to the original |
---|
299 | ! 16 g-points for each band. |
---|
300 | ! ngn The number of original g-points that are combined to make |
---|
301 | ! each new g-point in each band. |
---|
302 | ! ngb The band index for each new g-point. |
---|
303 | ! wt RRTM weights for 16 g-points. |
---|
304 | |
---|
305 | ! ------- Data statements ------- |
---|
306 | ngc(:) = (/10,12,16,14,16,8,12,8,12,6,8,8,4,2,2,2/) |
---|
307 | ngs(:) = (/10,22,38,52,68,76,88,96,108,114,122,130,134,136,138,140/) |
---|
308 | ngm(:) = (/1,2,3,3,4,4,5,5,6,6,7,7,8,8,9,10, & ! band 1 |
---|
309 | 1,2,3,4,5,6,7,8,9,9,10,10,11,11,12,12, & ! band 2 |
---|
310 | 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! band 3 |
---|
311 | 1,2,3,4,5,6,7,8,9,10,11,12,13,14,14,14, & ! band 4 |
---|
312 | 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, & ! band 5 |
---|
313 | 1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8, & ! band 6 |
---|
314 | 1,1,2,2,3,4,5,6,7,8,9,10,11,11,12,12, & ! band 7 |
---|
315 | 1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8, & ! band 8 |
---|
316 | 1,2,3,4,5,6,7,8,9,9,10,10,11,11,12,12, & ! band 9 |
---|
317 | 1,1,2,2,3,3,4,4,5,5,5,5,6,6,6,6, & ! band 10 |
---|
318 | 1,2,3,3,4,4,5,5,6,6,7,7,7,8,8,8, & ! band 11 |
---|
319 | 1,2,3,4,5,5,6,6,7,7,7,7,8,8,8,8, & ! band 12 |
---|
320 | 1,1,1,2,2,2,3,3,3,3,4,4,4,4,4,4, & ! band 13 |
---|
321 | 1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2, & ! band 14 |
---|
322 | 1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2, & ! band 15 |
---|
323 | 1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2/) ! band 16 |
---|
324 | ngn(:) = (/1,1,2,2,2,2,2,2,1,1, & ! band 1 |
---|
325 | 1,1,1,1,1,1,1,1,2,2,2,2, & ! band 2 |
---|
326 | 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! band 3 |
---|
327 | 1,1,1,1,1,1,1,1,1,1,1,1,1,3, & ! band 4 |
---|
328 | 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, & ! band 5 |
---|
329 | 2,2,2,2,2,2,2,2, & ! band 6 |
---|
330 | 2,2,1,1,1,1,1,1,1,1,2,2, & ! band 7 |
---|
331 | 2,2,2,2,2,2,2,2, & ! band 8 |
---|
332 | 1,1,1,1,1,1,1,1,2,2,2,2, & ! band 9 |
---|
333 | 2,2,2,2,4,4, & ! band 10 |
---|
334 | 1,1,2,2,2,2,3,3, & ! band 11 |
---|
335 | 1,1,1,1,2,2,4,4, & ! band 12 |
---|
336 | 3,3,4,6, & ! band 13 |
---|
337 | 8,8, & ! band 14 |
---|
338 | 8,8, & ! band 15 |
---|
339 | 4,12/) ! band 16 |
---|
340 | ngb(:) = (/1,1,1,1,1,1,1,1,1,1, & ! band 1 |
---|
341 | 2,2,2,2,2,2,2,2,2,2,2,2, & ! band 2 |
---|
342 | 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3, & ! band 3 |
---|
343 | 4,4,4,4,4,4,4,4,4,4,4,4,4,4, & ! band 4 |
---|
344 | 5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5, & ! band 5 |
---|
345 | 6,6,6,6,6,6,6,6, & ! band 6 |
---|
346 | 7,7,7,7,7,7,7,7,7,7,7,7, & ! band 7 |
---|
347 | 8,8,8,8,8,8,8,8, & ! band 8 |
---|
348 | 9,9,9,9,9,9,9,9,9,9,9,9, & ! band 9 |
---|
349 | 10,10,10,10,10,10, & ! band 10 |
---|
350 | 11,11,11,11,11,11,11,11, & ! band 11 |
---|
351 | 12,12,12,12,12,12,12,12, & ! band 12 |
---|
352 | 13,13,13,13, & ! band 13 |
---|
353 | 14,14, & ! band 14 |
---|
354 | 15,15, & ! band 15 |
---|
355 | 16,16/) ! band 16 |
---|
356 | wt(:) = (/ 0.1527534276_rb, 0.1491729617_rb, 0.1420961469_rb, & |
---|
357 | 0.1316886544_rb, 0.1181945205_rb, 0.1019300893_rb, & |
---|
358 | 0.0832767040_rb, 0.0626720116_rb, 0.0424925000_rb, & |
---|
359 | 0.0046269894_rb, 0.0038279891_rb, 0.0030260086_rb, & |
---|
360 | 0.0022199750_rb, 0.0014140010_rb, 0.0005330000_rb, & |
---|
361 | 0.0000750000_rb/) |
---|
362 | |
---|
363 | end subroutine lwcmbdat |
---|
364 | |
---|
365 | !*************************************************************************** |
---|
366 | subroutine cmbgb1 |
---|
367 | !*************************************************************************** |
---|
368 | ! |
---|
369 | ! Original version: MJIacono; July 1998 |
---|
370 | ! Revision for GCMs: MJIacono; September 1998 |
---|
371 | ! Revision for RRTMG: MJIacono, September 2002 |
---|
372 | ! Revision for F90 reformatting: MJIacono, June 2006 |
---|
373 | ! |
---|
374 | ! The subroutines CMBGB1->CMBGB16 input the absorption coefficient |
---|
375 | ! data for each band, which are defined for 16 g-points and 16 spectral |
---|
376 | ! bands. The data are combined with appropriate weighting following the |
---|
377 | ! g-point mapping arrays specified in RRTMINIT. Plank fraction data |
---|
378 | ! in arrays FRACREFA and FRACREFB are combined without weighting. All |
---|
379 | ! g-point reduced data are put into new arrays for use in RRTM. |
---|
380 | ! |
---|
381 | ! band 1: 10-350 cm-1 (low key - h2o; low minor - n2) |
---|
382 | ! (high key - h2o; high minor - n2) |
---|
383 | ! note: previous versions of rrtm band 1: |
---|
384 | ! 10-250 cm-1 (low - h2o; high - h2o) |
---|
385 | !*************************************************************************** |
---|
386 | |
---|
387 | use parrrtm, only : mg, nbndlw, ngptlw, ng1 |
---|
388 | use rrlw_kg01, only: fracrefao, fracrefbo, kao, kbo, kao_mn2, kbo_mn2, & |
---|
389 | selfrefo, forrefo, & |
---|
390 | fracrefa, fracrefb, absa, ka, absb, kb, ka_mn2, kb_mn2, & |
---|
391 | selfref, forref |
---|
392 | |
---|
393 | ! ------- Local ------- |
---|
394 | integer(kind=im) :: jt, jp, igc, ipr, iprsm |
---|
395 | real(kind=rb) :: sumk, sumk1, sumk2, sumf1, sumf2 |
---|
396 | |
---|
397 | |
---|
398 | do jt = 1,5 |
---|
399 | do jp = 1,13 |
---|
400 | iprsm = 0 |
---|
401 | do igc = 1,ngc(1) |
---|
402 | sumk = 0. |
---|
403 | do ipr = 1, ngn(igc) |
---|
404 | iprsm = iprsm + 1 |
---|
405 | sumk = sumk + kao(jt,jp,iprsm)*rwgt(iprsm) |
---|
406 | enddo |
---|
407 | ka(jt,jp,igc) = sumk |
---|
408 | enddo |
---|
409 | enddo |
---|
410 | do jp = 13,59 |
---|
411 | iprsm = 0 |
---|
412 | do igc = 1,ngc(1) |
---|
413 | sumk = 0. |
---|
414 | do ipr = 1, ngn(igc) |
---|
415 | iprsm = iprsm + 1 |
---|
416 | sumk = sumk + kbo(jt,jp,iprsm)*rwgt(iprsm) |
---|
417 | enddo |
---|
418 | kb(jt,jp,igc) = sumk |
---|
419 | enddo |
---|
420 | enddo |
---|
421 | enddo |
---|
422 | |
---|
423 | do jt = 1,10 |
---|
424 | iprsm = 0 |
---|
425 | do igc = 1,ngc(1) |
---|
426 | sumk = 0. |
---|
427 | do ipr = 1, ngn(igc) |
---|
428 | iprsm = iprsm + 1 |
---|
429 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm) |
---|
430 | enddo |
---|
431 | selfref(jt,igc) = sumk |
---|
432 | enddo |
---|
433 | enddo |
---|
434 | |
---|
435 | do jt = 1,4 |
---|
436 | iprsm = 0 |
---|
437 | do igc = 1,ngc(1) |
---|
438 | sumk = 0. |
---|
439 | do ipr = 1, ngn(igc) |
---|
440 | iprsm = iprsm + 1 |
---|
441 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm) |
---|
442 | enddo |
---|
443 | forref(jt,igc) = sumk |
---|
444 | enddo |
---|
445 | enddo |
---|
446 | |
---|
447 | do jt = 1,19 |
---|
448 | iprsm = 0 |
---|
449 | do igc = 1,ngc(1) |
---|
450 | sumk1 = 0. |
---|
451 | sumk2 = 0. |
---|
452 | do ipr = 1, ngn(igc) |
---|
453 | iprsm = iprsm + 1 |
---|
454 | sumk1 = sumk1 + kao_mn2(jt,iprsm)*rwgt(iprsm) |
---|
455 | sumk2 = sumk2 + kbo_mn2(jt,iprsm)*rwgt(iprsm) |
---|
456 | enddo |
---|
457 | ka_mn2(jt,igc) = sumk1 |
---|
458 | kb_mn2(jt,igc) = sumk2 |
---|
459 | enddo |
---|
460 | enddo |
---|
461 | |
---|
462 | iprsm = 0 |
---|
463 | do igc = 1,ngc(1) |
---|
464 | sumf1 = 0. |
---|
465 | sumf2 = 0. |
---|
466 | do ipr = 1, ngn(igc) |
---|
467 | iprsm = iprsm + 1 |
---|
468 | sumf1= sumf1+ fracrefao(iprsm) |
---|
469 | sumf2= sumf2+ fracrefbo(iprsm) |
---|
470 | enddo |
---|
471 | fracrefa(igc) = sumf1 |
---|
472 | fracrefb(igc) = sumf2 |
---|
473 | enddo |
---|
474 | |
---|
475 | end subroutine cmbgb1 |
---|
476 | |
---|
477 | !*************************************************************************** |
---|
478 | subroutine cmbgb2 |
---|
479 | !*************************************************************************** |
---|
480 | ! |
---|
481 | ! band 2: 350-500 cm-1 (low key - h2o; high key - h2o) |
---|
482 | ! |
---|
483 | ! note: previous version of rrtm band 2: |
---|
484 | ! 250 - 500 cm-1 (low - h2o; high - h2o) |
---|
485 | !*************************************************************************** |
---|
486 | |
---|
487 | use parrrtm, only : mg, nbndlw, ngptlw, ng2 |
---|
488 | use rrlw_kg02, only: fracrefao, fracrefbo, kao, kbo, selfrefo, forrefo, & |
---|
489 | fracrefa, fracrefb, absa, ka, absb, kb, selfref, forref |
---|
490 | |
---|
491 | ! ------- Local ------- |
---|
492 | integer(kind=im) :: jt, jp, igc, ipr, iprsm |
---|
493 | real(kind=rb) :: sumk, sumf1, sumf2 |
---|
494 | |
---|
495 | |
---|
496 | do jt = 1,5 |
---|
497 | do jp = 1,13 |
---|
498 | iprsm = 0 |
---|
499 | do igc = 1,ngc(2) |
---|
500 | sumk = 0. |
---|
501 | do ipr = 1, ngn(ngs(1)+igc) |
---|
502 | iprsm = iprsm + 1 |
---|
503 | sumk = sumk + kao(jt,jp,iprsm)*rwgt(iprsm+16) |
---|
504 | enddo |
---|
505 | ka(jt,jp,igc) = sumk |
---|
506 | enddo |
---|
507 | enddo |
---|
508 | do jp = 13,59 |
---|
509 | iprsm = 0 |
---|
510 | do igc = 1,ngc(2) |
---|
511 | sumk = 0. |
---|
512 | do ipr = 1, ngn(ngs(1)+igc) |
---|
513 | iprsm = iprsm + 1 |
---|
514 | sumk = sumk + kbo(jt,jp,iprsm)*rwgt(iprsm+16) |
---|
515 | enddo |
---|
516 | kb(jt,jp,igc) = sumk |
---|
517 | enddo |
---|
518 | enddo |
---|
519 | enddo |
---|
520 | |
---|
521 | do jt = 1,10 |
---|
522 | iprsm = 0 |
---|
523 | do igc = 1,ngc(2) |
---|
524 | sumk = 0. |
---|
525 | do ipr = 1, ngn(ngs(1)+igc) |
---|
526 | iprsm = iprsm + 1 |
---|
527 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+16) |
---|
528 | enddo |
---|
529 | selfref(jt,igc) = sumk |
---|
530 | enddo |
---|
531 | enddo |
---|
532 | |
---|
533 | do jt = 1,4 |
---|
534 | iprsm = 0 |
---|
535 | do igc = 1,ngc(2) |
---|
536 | sumk = 0. |
---|
537 | do ipr = 1, ngn(ngs(1)+igc) |
---|
538 | iprsm = iprsm + 1 |
---|
539 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+16) |
---|
540 | enddo |
---|
541 | forref(jt,igc) = sumk |
---|
542 | enddo |
---|
543 | enddo |
---|
544 | |
---|
545 | iprsm = 0 |
---|
546 | do igc = 1,ngc(2) |
---|
547 | sumf1 = 0. |
---|
548 | sumf2 = 0. |
---|
549 | do ipr = 1, ngn(ngs(1)+igc) |
---|
550 | iprsm = iprsm + 1 |
---|
551 | sumf1= sumf1+ fracrefao(iprsm) |
---|
552 | sumf2= sumf2+ fracrefbo(iprsm) |
---|
553 | enddo |
---|
554 | fracrefa(igc) = sumf1 |
---|
555 | fracrefb(igc) = sumf2 |
---|
556 | enddo |
---|
557 | |
---|
558 | end subroutine cmbgb2 |
---|
559 | |
---|
560 | !*************************************************************************** |
---|
561 | subroutine cmbgb3 |
---|
562 | !*************************************************************************** |
---|
563 | ! |
---|
564 | ! band 3: 500-630 cm-1 (low key - h2o,co2; low minor - n2o) |
---|
565 | ! (high key - h2o,co2; high minor - n2o) |
---|
566 | ! |
---|
567 | ! old band 3: 500-630 cm-1 (low - h2o,co2; high - h2o,co2) |
---|
568 | !*************************************************************************** |
---|
569 | |
---|
570 | use parrrtm, only : mg, nbndlw, ngptlw, ng3 |
---|
571 | use rrlw_kg03, only: fracrefao, fracrefbo, kao, kbo, kao_mn2o, kbo_mn2o, & |
---|
572 | selfrefo, forrefo, & |
---|
573 | fracrefa, fracrefb, absa, ka, absb, kb, ka_mn2o, kb_mn2o, & |
---|
574 | selfref, forref |
---|
575 | |
---|
576 | ! ------- Local ------- |
---|
577 | integer(kind=im) :: jn, jt, jp, igc, ipr, iprsm |
---|
578 | real(kind=rb) :: sumk, sumf |
---|
579 | |
---|
580 | |
---|
581 | do jn = 1,9 |
---|
582 | do jt = 1,5 |
---|
583 | do jp = 1,13 |
---|
584 | iprsm = 0 |
---|
585 | do igc = 1,ngc(3) |
---|
586 | sumk = 0. |
---|
587 | do ipr = 1, ngn(ngs(2)+igc) |
---|
588 | iprsm = iprsm + 1 |
---|
589 | sumk = sumk + kao(jn,jt,jp,iprsm)*rwgt(iprsm+32) |
---|
590 | enddo |
---|
591 | ka(jn,jt,jp,igc) = sumk |
---|
592 | enddo |
---|
593 | enddo |
---|
594 | enddo |
---|
595 | enddo |
---|
596 | do jn = 1,5 |
---|
597 | do jt = 1,5 |
---|
598 | do jp = 13,59 |
---|
599 | iprsm = 0 |
---|
600 | do igc = 1,ngc(3) |
---|
601 | sumk = 0. |
---|
602 | do ipr = 1, ngn(ngs(2)+igc) |
---|
603 | iprsm = iprsm + 1 |
---|
604 | sumk = sumk + kbo(jn,jt,jp,iprsm)*rwgt(iprsm+32) |
---|
605 | enddo |
---|
606 | kb(jn,jt,jp,igc) = sumk |
---|
607 | enddo |
---|
608 | enddo |
---|
609 | enddo |
---|
610 | enddo |
---|
611 | |
---|
612 | do jn = 1,9 |
---|
613 | do jt = 1,19 |
---|
614 | iprsm = 0 |
---|
615 | do igc = 1,ngc(3) |
---|
616 | sumk = 0. |
---|
617 | do ipr = 1, ngn(ngs(2)+igc) |
---|
618 | iprsm = iprsm + 1 |
---|
619 | sumk = sumk + kao_mn2o(jn,jt,iprsm)*rwgt(iprsm+32) |
---|
620 | enddo |
---|
621 | ka_mn2o(jn,jt,igc) = sumk |
---|
622 | enddo |
---|
623 | enddo |
---|
624 | enddo |
---|
625 | |
---|
626 | do jn = 1,5 |
---|
627 | do jt = 1,19 |
---|
628 | iprsm = 0 |
---|
629 | do igc = 1,ngc(3) |
---|
630 | sumk = 0. |
---|
631 | do ipr = 1, ngn(ngs(2)+igc) |
---|
632 | iprsm = iprsm + 1 |
---|
633 | sumk = sumk + kbo_mn2o(jn,jt,iprsm)*rwgt(iprsm+32) |
---|
634 | enddo |
---|
635 | kb_mn2o(jn,jt,igc) = sumk |
---|
636 | enddo |
---|
637 | enddo |
---|
638 | enddo |
---|
639 | |
---|
640 | do jt = 1,10 |
---|
641 | iprsm = 0 |
---|
642 | do igc = 1,ngc(3) |
---|
643 | sumk = 0. |
---|
644 | do ipr = 1, ngn(ngs(2)+igc) |
---|
645 | iprsm = iprsm + 1 |
---|
646 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+32) |
---|
647 | enddo |
---|
648 | selfref(jt,igc) = sumk |
---|
649 | enddo |
---|
650 | enddo |
---|
651 | |
---|
652 | do jt = 1,4 |
---|
653 | iprsm = 0 |
---|
654 | do igc = 1,ngc(3) |
---|
655 | sumk = 0. |
---|
656 | do ipr = 1, ngn(ngs(2)+igc) |
---|
657 | iprsm = iprsm + 1 |
---|
658 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+32) |
---|
659 | enddo |
---|
660 | forref(jt,igc) = sumk |
---|
661 | enddo |
---|
662 | enddo |
---|
663 | |
---|
664 | do jp = 1,9 |
---|
665 | iprsm = 0 |
---|
666 | do igc = 1,ngc(3) |
---|
667 | sumf = 0. |
---|
668 | do ipr = 1, ngn(ngs(2)+igc) |
---|
669 | iprsm = iprsm + 1 |
---|
670 | sumf = sumf + fracrefao(iprsm,jp) |
---|
671 | enddo |
---|
672 | fracrefa(igc,jp) = sumf |
---|
673 | enddo |
---|
674 | enddo |
---|
675 | |
---|
676 | do jp = 1,5 |
---|
677 | iprsm = 0 |
---|
678 | do igc = 1,ngc(3) |
---|
679 | sumf = 0. |
---|
680 | do ipr = 1, ngn(ngs(2)+igc) |
---|
681 | iprsm = iprsm + 1 |
---|
682 | sumf = sumf + fracrefbo(iprsm,jp) |
---|
683 | enddo |
---|
684 | fracrefb(igc,jp) = sumf |
---|
685 | enddo |
---|
686 | enddo |
---|
687 | |
---|
688 | end subroutine cmbgb3 |
---|
689 | |
---|
690 | !*************************************************************************** |
---|
691 | subroutine cmbgb4 |
---|
692 | !*************************************************************************** |
---|
693 | ! |
---|
694 | ! band 4: 630-700 cm-1 (low key - h2o,co2; high key - o3,co2) |
---|
695 | ! |
---|
696 | ! old band 4: 630-700 cm-1 (low - h2o,co2; high - o3,co2) |
---|
697 | !*************************************************************************** |
---|
698 | |
---|
699 | use parrrtm, only : mg, nbndlw, ngptlw, ng4 |
---|
700 | use rrlw_kg04, only: fracrefao, fracrefbo, kao, kbo, selfrefo, forrefo, & |
---|
701 | fracrefa, fracrefb, absa, ka, absb, kb, selfref, forref |
---|
702 | |
---|
703 | ! ------- Local ------- |
---|
704 | integer(kind=im) :: jn, jt, jp, igc, ipr, iprsm |
---|
705 | real(kind=rb) :: sumk, sumf |
---|
706 | |
---|
707 | |
---|
708 | do jn = 1,9 |
---|
709 | do jt = 1,5 |
---|
710 | do jp = 1,13 |
---|
711 | iprsm = 0 |
---|
712 | do igc = 1,ngc(4) |
---|
713 | sumk = 0. |
---|
714 | do ipr = 1, ngn(ngs(3)+igc) |
---|
715 | iprsm = iprsm + 1 |
---|
716 | sumk = sumk + kao(jn,jt,jp,iprsm)*rwgt(iprsm+48) |
---|
717 | enddo |
---|
718 | ka(jn,jt,jp,igc) = sumk |
---|
719 | enddo |
---|
720 | enddo |
---|
721 | enddo |
---|
722 | enddo |
---|
723 | do jn = 1,5 |
---|
724 | do jt = 1,5 |
---|
725 | do jp = 13,59 |
---|
726 | iprsm = 0 |
---|
727 | do igc = 1,ngc(4) |
---|
728 | sumk = 0. |
---|
729 | do ipr = 1, ngn(ngs(3)+igc) |
---|
730 | iprsm = iprsm + 1 |
---|
731 | sumk = sumk + kbo(jn,jt,jp,iprsm)*rwgt(iprsm+48) |
---|
732 | enddo |
---|
733 | kb(jn,jt,jp,igc) = sumk |
---|
734 | enddo |
---|
735 | enddo |
---|
736 | enddo |
---|
737 | enddo |
---|
738 | |
---|
739 | do jt = 1,10 |
---|
740 | iprsm = 0 |
---|
741 | do igc = 1,ngc(4) |
---|
742 | sumk = 0. |
---|
743 | do ipr = 1, ngn(ngs(3)+igc) |
---|
744 | iprsm = iprsm + 1 |
---|
745 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+48) |
---|
746 | enddo |
---|
747 | selfref(jt,igc) = sumk |
---|
748 | enddo |
---|
749 | enddo |
---|
750 | |
---|
751 | do jt = 1,4 |
---|
752 | iprsm = 0 |
---|
753 | do igc = 1,ngc(4) |
---|
754 | sumk = 0. |
---|
755 | do ipr = 1, ngn(ngs(3)+igc) |
---|
756 | iprsm = iprsm + 1 |
---|
757 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+48) |
---|
758 | enddo |
---|
759 | forref(jt,igc) = sumk |
---|
760 | enddo |
---|
761 | enddo |
---|
762 | |
---|
763 | do jp = 1,9 |
---|
764 | iprsm = 0 |
---|
765 | do igc = 1,ngc(4) |
---|
766 | sumf = 0. |
---|
767 | do ipr = 1, ngn(ngs(3)+igc) |
---|
768 | iprsm = iprsm + 1 |
---|
769 | sumf = sumf + fracrefao(iprsm,jp) |
---|
770 | enddo |
---|
771 | fracrefa(igc,jp) = sumf |
---|
772 | enddo |
---|
773 | enddo |
---|
774 | |
---|
775 | do jp = 1,5 |
---|
776 | iprsm = 0 |
---|
777 | do igc = 1,ngc(4) |
---|
778 | sumf = 0. |
---|
779 | do ipr = 1, ngn(ngs(3)+igc) |
---|
780 | iprsm = iprsm + 1 |
---|
781 | sumf = sumf + fracrefbo(iprsm,jp) |
---|
782 | enddo |
---|
783 | fracrefb(igc,jp) = sumf |
---|
784 | enddo |
---|
785 | enddo |
---|
786 | |
---|
787 | end subroutine cmbgb4 |
---|
788 | |
---|
789 | !*************************************************************************** |
---|
790 | subroutine cmbgb5 |
---|
791 | !*************************************************************************** |
---|
792 | ! |
---|
793 | ! band 5: 700-820 cm-1 (low key - h2o,co2; low minor - o3, ccl4) |
---|
794 | ! (high key - o3,co2) |
---|
795 | ! |
---|
796 | ! old band 5: 700-820 cm-1 (low - h2o,co2; high - o3,co2) |
---|
797 | !*************************************************************************** |
---|
798 | |
---|
799 | use parrrtm, only : mg, nbndlw, ngptlw, ng5 |
---|
800 | use rrlw_kg05, only: fracrefao, fracrefbo, kao, kbo, kao_mo3, ccl4o, & |
---|
801 | selfrefo, forrefo, & |
---|
802 | fracrefa, fracrefb, absa, ka, absb, kb, ka_mo3, ccl4, & |
---|
803 | selfref, forref |
---|
804 | |
---|
805 | ! ------- Local ------- |
---|
806 | integer(kind=im) :: jn, jt, jp, igc, ipr, iprsm |
---|
807 | real(kind=rb) :: sumk, sumf |
---|
808 | |
---|
809 | |
---|
810 | do jn = 1,9 |
---|
811 | do jt = 1,5 |
---|
812 | do jp = 1,13 |
---|
813 | iprsm = 0 |
---|
814 | do igc = 1,ngc(5) |
---|
815 | sumk = 0. |
---|
816 | do ipr = 1, ngn(ngs(4)+igc) |
---|
817 | iprsm = iprsm + 1 |
---|
818 | sumk = sumk + kao(jn,jt,jp,iprsm)*rwgt(iprsm+64) |
---|
819 | enddo |
---|
820 | ka(jn,jt,jp,igc) = sumk |
---|
821 | enddo |
---|
822 | enddo |
---|
823 | enddo |
---|
824 | enddo |
---|
825 | do jn = 1,5 |
---|
826 | do jt = 1,5 |
---|
827 | do jp = 13,59 |
---|
828 | iprsm = 0 |
---|
829 | do igc = 1,ngc(5) |
---|
830 | sumk = 0. |
---|
831 | do ipr = 1, ngn(ngs(4)+igc) |
---|
832 | iprsm = iprsm + 1 |
---|
833 | sumk = sumk + kbo(jn,jt,jp,iprsm)*rwgt(iprsm+64) |
---|
834 | enddo |
---|
835 | kb(jn,jt,jp,igc) = sumk |
---|
836 | enddo |
---|
837 | enddo |
---|
838 | enddo |
---|
839 | enddo |
---|
840 | |
---|
841 | do jn = 1,9 |
---|
842 | do jt = 1,19 |
---|
843 | iprsm = 0 |
---|
844 | do igc = 1,ngc(5) |
---|
845 | sumk = 0. |
---|
846 | do ipr = 1, ngn(ngs(4)+igc) |
---|
847 | iprsm = iprsm + 1 |
---|
848 | sumk = sumk + kao_mo3(jn,jt,iprsm)*rwgt(iprsm+64) |
---|
849 | enddo |
---|
850 | ka_mo3(jn,jt,igc) = sumk |
---|
851 | enddo |
---|
852 | enddo |
---|
853 | enddo |
---|
854 | |
---|
855 | do jt = 1,10 |
---|
856 | iprsm = 0 |
---|
857 | do igc = 1,ngc(5) |
---|
858 | sumk = 0. |
---|
859 | do ipr = 1, ngn(ngs(4)+igc) |
---|
860 | iprsm = iprsm + 1 |
---|
861 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+64) |
---|
862 | enddo |
---|
863 | selfref(jt,igc) = sumk |
---|
864 | enddo |
---|
865 | enddo |
---|
866 | |
---|
867 | do jt = 1,4 |
---|
868 | iprsm = 0 |
---|
869 | do igc = 1,ngc(5) |
---|
870 | sumk = 0. |
---|
871 | do ipr = 1, ngn(ngs(4)+igc) |
---|
872 | iprsm = iprsm + 1 |
---|
873 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+64) |
---|
874 | enddo |
---|
875 | forref(jt,igc) = sumk |
---|
876 | enddo |
---|
877 | enddo |
---|
878 | |
---|
879 | do jp = 1,9 |
---|
880 | iprsm = 0 |
---|
881 | do igc = 1,ngc(5) |
---|
882 | sumf = 0. |
---|
883 | do ipr = 1, ngn(ngs(4)+igc) |
---|
884 | iprsm = iprsm + 1 |
---|
885 | sumf = sumf + fracrefao(iprsm,jp) |
---|
886 | enddo |
---|
887 | fracrefa(igc,jp) = sumf |
---|
888 | enddo |
---|
889 | enddo |
---|
890 | |
---|
891 | do jp = 1,5 |
---|
892 | iprsm = 0 |
---|
893 | do igc = 1,ngc(5) |
---|
894 | sumf = 0. |
---|
895 | do ipr = 1, ngn(ngs(4)+igc) |
---|
896 | iprsm = iprsm + 1 |
---|
897 | sumf = sumf + fracrefbo(iprsm,jp) |
---|
898 | enddo |
---|
899 | fracrefb(igc,jp) = sumf |
---|
900 | enddo |
---|
901 | enddo |
---|
902 | |
---|
903 | iprsm = 0 |
---|
904 | do igc = 1,ngc(5) |
---|
905 | sumk = 0. |
---|
906 | do ipr = 1, ngn(ngs(4)+igc) |
---|
907 | iprsm = iprsm + 1 |
---|
908 | sumk = sumk + ccl4o(iprsm)*rwgt(iprsm+64) |
---|
909 | enddo |
---|
910 | ccl4(igc) = sumk |
---|
911 | enddo |
---|
912 | |
---|
913 | end subroutine cmbgb5 |
---|
914 | |
---|
915 | !*************************************************************************** |
---|
916 | subroutine cmbgb6 |
---|
917 | !*************************************************************************** |
---|
918 | ! |
---|
919 | ! band 6: 820-980 cm-1 (low key - h2o; low minor - co2) |
---|
920 | ! (high key - nothing; high minor - cfc11, cfc12) |
---|
921 | ! |
---|
922 | ! old band 6: 820-980 cm-1 (low - h2o; high - nothing) |
---|
923 | !*************************************************************************** |
---|
924 | |
---|
925 | use parrrtm, only : mg, nbndlw, ngptlw, ng6 |
---|
926 | use rrlw_kg06, only: fracrefao, kao, kao_mco2, cfc11adjo, cfc12o, & |
---|
927 | selfrefo, forrefo, & |
---|
928 | fracrefa, absa, ka, ka_mco2, cfc11adj, cfc12, & |
---|
929 | selfref, forref |
---|
930 | |
---|
931 | ! ------- Local ------- |
---|
932 | integer(kind=im) :: jt, jp, igc, ipr, iprsm |
---|
933 | real(kind=rb) :: sumk, sumf, sumk1, sumk2 |
---|
934 | |
---|
935 | |
---|
936 | do jt = 1,5 |
---|
937 | do jp = 1,13 |
---|
938 | iprsm = 0 |
---|
939 | do igc = 1,ngc(6) |
---|
940 | sumk = 0. |
---|
941 | do ipr = 1, ngn(ngs(5)+igc) |
---|
942 | iprsm = iprsm + 1 |
---|
943 | sumk = sumk + kao(jt,jp,iprsm)*rwgt(iprsm+80) |
---|
944 | enddo |
---|
945 | ka(jt,jp,igc) = sumk |
---|
946 | enddo |
---|
947 | enddo |
---|
948 | enddo |
---|
949 | |
---|
950 | do jt = 1,19 |
---|
951 | iprsm = 0 |
---|
952 | do igc = 1,ngc(6) |
---|
953 | sumk = 0. |
---|
954 | do ipr = 1, ngn(ngs(5)+igc) |
---|
955 | iprsm = iprsm + 1 |
---|
956 | sumk = sumk + kao_mco2(jt,iprsm)*rwgt(iprsm+80) |
---|
957 | enddo |
---|
958 | ka_mco2(jt,igc) = sumk |
---|
959 | enddo |
---|
960 | enddo |
---|
961 | |
---|
962 | do jt = 1,10 |
---|
963 | iprsm = 0 |
---|
964 | do igc = 1,ngc(6) |
---|
965 | sumk = 0. |
---|
966 | do ipr = 1, ngn(ngs(5)+igc) |
---|
967 | iprsm = iprsm + 1 |
---|
968 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+80) |
---|
969 | enddo |
---|
970 | selfref(jt,igc) = sumk |
---|
971 | enddo |
---|
972 | enddo |
---|
973 | |
---|
974 | do jt = 1,4 |
---|
975 | iprsm = 0 |
---|
976 | do igc = 1,ngc(6) |
---|
977 | sumk = 0. |
---|
978 | do ipr = 1, ngn(ngs(5)+igc) |
---|
979 | iprsm = iprsm + 1 |
---|
980 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+80) |
---|
981 | enddo |
---|
982 | forref(jt,igc) = sumk |
---|
983 | enddo |
---|
984 | enddo |
---|
985 | |
---|
986 | iprsm = 0 |
---|
987 | do igc = 1,ngc(6) |
---|
988 | sumf = 0. |
---|
989 | sumk1= 0. |
---|
990 | sumk2= 0. |
---|
991 | do ipr = 1, ngn(ngs(5)+igc) |
---|
992 | iprsm = iprsm + 1 |
---|
993 | sumf = sumf + fracrefao(iprsm) |
---|
994 | sumk1= sumk1+ cfc11adjo(iprsm)*rwgt(iprsm+80) |
---|
995 | sumk2= sumk2+ cfc12o(iprsm)*rwgt(iprsm+80) |
---|
996 | enddo |
---|
997 | fracrefa(igc) = sumf |
---|
998 | cfc11adj(igc) = sumk1 |
---|
999 | cfc12(igc) = sumk2 |
---|
1000 | enddo |
---|
1001 | |
---|
1002 | end subroutine cmbgb6 |
---|
1003 | |
---|
1004 | !*************************************************************************** |
---|
1005 | subroutine cmbgb7 |
---|
1006 | !*************************************************************************** |
---|
1007 | ! |
---|
1008 | ! band 7: 980-1080 cm-1 (low key - h2o,o3; low minor - co2) |
---|
1009 | ! (high key - o3; high minor - co2) |
---|
1010 | ! |
---|
1011 | ! old band 7: 980-1080 cm-1 (low - h2o,o3; high - o3) |
---|
1012 | !*************************************************************************** |
---|
1013 | |
---|
1014 | use parrrtm, only : mg, nbndlw, ngptlw, ng7 |
---|
1015 | use rrlw_kg07, only: fracrefao, fracrefbo, kao, kbo, kao_mco2, kbo_mco2, & |
---|
1016 | selfrefo, forrefo, & |
---|
1017 | fracrefa, fracrefb, absa, ka, absb, kb, ka_mco2, kb_mco2, & |
---|
1018 | selfref, forref |
---|
1019 | |
---|
1020 | ! ------- Local ------- |
---|
1021 | integer(kind=im) :: jn, jt, jp, igc, ipr, iprsm |
---|
1022 | real(kind=rb) :: sumk, sumf |
---|
1023 | |
---|
1024 | |
---|
1025 | do jn = 1,9 |
---|
1026 | do jt = 1,5 |
---|
1027 | do jp = 1,13 |
---|
1028 | iprsm = 0 |
---|
1029 | do igc = 1,ngc(7) |
---|
1030 | sumk = 0. |
---|
1031 | do ipr = 1, ngn(ngs(6)+igc) |
---|
1032 | iprsm = iprsm + 1 |
---|
1033 | sumk = sumk + kao(jn,jt,jp,iprsm)*rwgt(iprsm+96) |
---|
1034 | enddo |
---|
1035 | ka(jn,jt,jp,igc) = sumk |
---|
1036 | enddo |
---|
1037 | enddo |
---|
1038 | enddo |
---|
1039 | enddo |
---|
1040 | do jt = 1,5 |
---|
1041 | do jp = 13,59 |
---|
1042 | iprsm = 0 |
---|
1043 | do igc = 1,ngc(7) |
---|
1044 | sumk = 0. |
---|
1045 | do ipr = 1, ngn(ngs(6)+igc) |
---|
1046 | iprsm = iprsm + 1 |
---|
1047 | sumk = sumk + kbo(jt,jp,iprsm)*rwgt(iprsm+96) |
---|
1048 | enddo |
---|
1049 | kb(jt,jp,igc) = sumk |
---|
1050 | enddo |
---|
1051 | enddo |
---|
1052 | enddo |
---|
1053 | |
---|
1054 | do jn = 1,9 |
---|
1055 | do jt = 1,19 |
---|
1056 | iprsm = 0 |
---|
1057 | do igc = 1,ngc(7) |
---|
1058 | sumk = 0. |
---|
1059 | do ipr = 1, ngn(ngs(6)+igc) |
---|
1060 | iprsm = iprsm + 1 |
---|
1061 | sumk = sumk + kao_mco2(jn,jt,iprsm)*rwgt(iprsm+96) |
---|
1062 | enddo |
---|
1063 | ka_mco2(jn,jt,igc) = sumk |
---|
1064 | enddo |
---|
1065 | enddo |
---|
1066 | enddo |
---|
1067 | |
---|
1068 | do jt = 1,19 |
---|
1069 | iprsm = 0 |
---|
1070 | do igc = 1,ngc(7) |
---|
1071 | sumk = 0. |
---|
1072 | do ipr = 1, ngn(ngs(6)+igc) |
---|
1073 | iprsm = iprsm + 1 |
---|
1074 | sumk = sumk + kbo_mco2(jt,iprsm)*rwgt(iprsm+96) |
---|
1075 | enddo |
---|
1076 | kb_mco2(jt,igc) = sumk |
---|
1077 | enddo |
---|
1078 | enddo |
---|
1079 | |
---|
1080 | do jt = 1,10 |
---|
1081 | iprsm = 0 |
---|
1082 | do igc = 1,ngc(7) |
---|
1083 | sumk = 0. |
---|
1084 | do ipr = 1, ngn(ngs(6)+igc) |
---|
1085 | iprsm = iprsm + 1 |
---|
1086 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+96) |
---|
1087 | enddo |
---|
1088 | selfref(jt,igc) = sumk |
---|
1089 | enddo |
---|
1090 | enddo |
---|
1091 | |
---|
1092 | do jt = 1,4 |
---|
1093 | iprsm = 0 |
---|
1094 | do igc = 1,ngc(7) |
---|
1095 | sumk = 0. |
---|
1096 | do ipr = 1, ngn(ngs(6)+igc) |
---|
1097 | iprsm = iprsm + 1 |
---|
1098 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+96) |
---|
1099 | enddo |
---|
1100 | forref(jt,igc) = sumk |
---|
1101 | enddo |
---|
1102 | enddo |
---|
1103 | |
---|
1104 | do jp = 1,9 |
---|
1105 | iprsm = 0 |
---|
1106 | do igc = 1,ngc(7) |
---|
1107 | sumf = 0. |
---|
1108 | do ipr = 1, ngn(ngs(6)+igc) |
---|
1109 | iprsm = iprsm + 1 |
---|
1110 | sumf = sumf + fracrefao(iprsm,jp) |
---|
1111 | enddo |
---|
1112 | fracrefa(igc,jp) = sumf |
---|
1113 | enddo |
---|
1114 | enddo |
---|
1115 | |
---|
1116 | iprsm = 0 |
---|
1117 | do igc = 1,ngc(7) |
---|
1118 | sumf = 0. |
---|
1119 | do ipr = 1, ngn(ngs(6)+igc) |
---|
1120 | iprsm = iprsm + 1 |
---|
1121 | sumf = sumf + fracrefbo(iprsm) |
---|
1122 | enddo |
---|
1123 | fracrefb(igc) = sumf |
---|
1124 | enddo |
---|
1125 | |
---|
1126 | end subroutine cmbgb7 |
---|
1127 | |
---|
1128 | !*************************************************************************** |
---|
1129 | subroutine cmbgb8 |
---|
1130 | !*************************************************************************** |
---|
1131 | ! |
---|
1132 | ! band 8: 1080-1180 cm-1 (low key - h2o; low minor - co2,o3,n2o) |
---|
1133 | ! (high key - o3; high minor - co2, n2o) |
---|
1134 | ! |
---|
1135 | ! old band 8: 1080-1180 cm-1 (low (i.e.>~300mb) - h2o; high - o3) |
---|
1136 | !*************************************************************************** |
---|
1137 | |
---|
1138 | use parrrtm, only : mg, nbndlw, ngptlw, ng8 |
---|
1139 | use rrlw_kg08, only: fracrefao, fracrefbo, kao, kao_mco2, kao_mn2o, & |
---|
1140 | kao_mo3, kbo, kbo_mco2, kbo_mn2o, selfrefo, forrefo, & |
---|
1141 | cfc12o, cfc22adjo, & |
---|
1142 | fracrefa, fracrefb, absa, ka, ka_mco2, ka_mn2o, & |
---|
1143 | ka_mo3, absb, kb, kb_mco2, kb_mn2o, selfref, forref, & |
---|
1144 | cfc12, cfc22adj |
---|
1145 | |
---|
1146 | ! ------- Local ------- |
---|
1147 | integer(kind=im) :: jt, jp, igc, ipr, iprsm |
---|
1148 | real(kind=rb) :: sumk, sumk1, sumk2, sumk3, sumk4, sumk5, sumf1, sumf2 |
---|
1149 | |
---|
1150 | |
---|
1151 | do jt = 1,5 |
---|
1152 | do jp = 1,13 |
---|
1153 | iprsm = 0 |
---|
1154 | do igc = 1,ngc(8) |
---|
1155 | sumk = 0. |
---|
1156 | do ipr = 1, ngn(ngs(7)+igc) |
---|
1157 | iprsm = iprsm + 1 |
---|
1158 | sumk = sumk + kao(jt,jp,iprsm)*rwgt(iprsm+112) |
---|
1159 | enddo |
---|
1160 | ka(jt,jp,igc) = sumk |
---|
1161 | enddo |
---|
1162 | enddo |
---|
1163 | enddo |
---|
1164 | do jt = 1,5 |
---|
1165 | do jp = 13,59 |
---|
1166 | iprsm = 0 |
---|
1167 | do igc = 1,ngc(8) |
---|
1168 | sumk = 0. |
---|
1169 | do ipr = 1, ngn(ngs(7)+igc) |
---|
1170 | iprsm = iprsm + 1 |
---|
1171 | sumk = sumk + kbo(jt,jp,iprsm)*rwgt(iprsm+112) |
---|
1172 | enddo |
---|
1173 | kb(jt,jp,igc) = sumk |
---|
1174 | enddo |
---|
1175 | enddo |
---|
1176 | enddo |
---|
1177 | |
---|
1178 | do jt = 1,10 |
---|
1179 | iprsm = 0 |
---|
1180 | do igc = 1,ngc(8) |
---|
1181 | sumk = 0. |
---|
1182 | do ipr = 1, ngn(ngs(7)+igc) |
---|
1183 | iprsm = iprsm + 1 |
---|
1184 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+112) |
---|
1185 | enddo |
---|
1186 | selfref(jt,igc) = sumk |
---|
1187 | enddo |
---|
1188 | enddo |
---|
1189 | |
---|
1190 | do jt = 1,4 |
---|
1191 | iprsm = 0 |
---|
1192 | do igc = 1,ngc(8) |
---|
1193 | sumk = 0. |
---|
1194 | do ipr = 1, ngn(ngs(7)+igc) |
---|
1195 | iprsm = iprsm + 1 |
---|
1196 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+112) |
---|
1197 | enddo |
---|
1198 | forref(jt,igc) = sumk |
---|
1199 | enddo |
---|
1200 | enddo |
---|
1201 | |
---|
1202 | do jt = 1,19 |
---|
1203 | iprsm = 0 |
---|
1204 | do igc = 1,ngc(8) |
---|
1205 | sumk1 = 0. |
---|
1206 | sumk2 = 0. |
---|
1207 | sumk3 = 0. |
---|
1208 | sumk4 = 0. |
---|
1209 | sumk5 = 0. |
---|
1210 | do ipr = 1, ngn(ngs(7)+igc) |
---|
1211 | iprsm = iprsm + 1 |
---|
1212 | sumk1 = sumk1 + kao_mco2(jt,iprsm)*rwgt(iprsm+112) |
---|
1213 | sumk2 = sumk2 + kbo_mco2(jt,iprsm)*rwgt(iprsm+112) |
---|
1214 | sumk3 = sumk3 + kao_mo3(jt,iprsm)*rwgt(iprsm+112) |
---|
1215 | sumk4 = sumk4 + kao_mn2o(jt,iprsm)*rwgt(iprsm+112) |
---|
1216 | sumk5 = sumk5 + kbo_mn2o(jt,iprsm)*rwgt(iprsm+112) |
---|
1217 | enddo |
---|
1218 | ka_mco2(jt,igc) = sumk1 |
---|
1219 | kb_mco2(jt,igc) = sumk2 |
---|
1220 | ka_mo3(jt,igc) = sumk3 |
---|
1221 | ka_mn2o(jt,igc) = sumk4 |
---|
1222 | kb_mn2o(jt,igc) = sumk5 |
---|
1223 | enddo |
---|
1224 | enddo |
---|
1225 | |
---|
1226 | iprsm = 0 |
---|
1227 | do igc = 1,ngc(8) |
---|
1228 | sumf1= 0. |
---|
1229 | sumf2= 0. |
---|
1230 | sumk1= 0. |
---|
1231 | sumk2= 0. |
---|
1232 | do ipr = 1, ngn(ngs(7)+igc) |
---|
1233 | iprsm = iprsm + 1 |
---|
1234 | sumf1= sumf1+ fracrefao(iprsm) |
---|
1235 | sumf2= sumf2+ fracrefbo(iprsm) |
---|
1236 | sumk1= sumk1+ cfc12o(iprsm)*rwgt(iprsm+112) |
---|
1237 | sumk2= sumk2+ cfc22adjo(iprsm)*rwgt(iprsm+112) |
---|
1238 | enddo |
---|
1239 | fracrefa(igc) = sumf1 |
---|
1240 | fracrefb(igc) = sumf2 |
---|
1241 | cfc12(igc) = sumk1 |
---|
1242 | cfc22adj(igc) = sumk2 |
---|
1243 | enddo |
---|
1244 | |
---|
1245 | end subroutine cmbgb8 |
---|
1246 | |
---|
1247 | !*************************************************************************** |
---|
1248 | subroutine cmbgb9 |
---|
1249 | !*************************************************************************** |
---|
1250 | ! |
---|
1251 | ! band 9: 1180-1390 cm-1 (low key - h2o,ch4; low minor - n2o) |
---|
1252 | ! (high key - ch4; high minor - n2o)! |
---|
1253 | |
---|
1254 | ! old band 9: 1180-1390 cm-1 (low - h2o,ch4; high - ch4) |
---|
1255 | !*************************************************************************** |
---|
1256 | |
---|
1257 | use parrrtm, only : mg, nbndlw, ngptlw, ng9 |
---|
1258 | use rrlw_kg09, only: fracrefao, fracrefbo, kao, kao_mn2o, & |
---|
1259 | kbo, kbo_mn2o, selfrefo, forrefo, & |
---|
1260 | fracrefa, fracrefb, absa, ka, ka_mn2o, & |
---|
1261 | absb, kb, kb_mn2o, selfref, forref |
---|
1262 | |
---|
1263 | ! ------- Local ------- |
---|
1264 | integer(kind=im) :: jn, jt, jp, igc, ipr, iprsm |
---|
1265 | real(kind=rb) :: sumk, sumf |
---|
1266 | |
---|
1267 | |
---|
1268 | do jn = 1,9 |
---|
1269 | do jt = 1,5 |
---|
1270 | do jp = 1,13 |
---|
1271 | iprsm = 0 |
---|
1272 | do igc = 1,ngc(9) |
---|
1273 | sumk = 0. |
---|
1274 | do ipr = 1, ngn(ngs(8)+igc) |
---|
1275 | iprsm = iprsm + 1 |
---|
1276 | sumk = sumk + kao(jn,jt,jp,iprsm)*rwgt(iprsm+128) |
---|
1277 | enddo |
---|
1278 | ka(jn,jt,jp,igc) = sumk |
---|
1279 | enddo |
---|
1280 | enddo |
---|
1281 | enddo |
---|
1282 | enddo |
---|
1283 | |
---|
1284 | do jt = 1,5 |
---|
1285 | do jp = 13,59 |
---|
1286 | iprsm = 0 |
---|
1287 | do igc = 1,ngc(9) |
---|
1288 | sumk = 0. |
---|
1289 | do ipr = 1, ngn(ngs(8)+igc) |
---|
1290 | iprsm = iprsm + 1 |
---|
1291 | sumk = sumk + kbo(jt,jp,iprsm)*rwgt(iprsm+128) |
---|
1292 | enddo |
---|
1293 | kb(jt,jp,igc) = sumk |
---|
1294 | enddo |
---|
1295 | enddo |
---|
1296 | enddo |
---|
1297 | |
---|
1298 | do jn = 1,9 |
---|
1299 | do jt = 1,19 |
---|
1300 | iprsm = 0 |
---|
1301 | do igc = 1,ngc(9) |
---|
1302 | sumk = 0. |
---|
1303 | do ipr = 1, ngn(ngs(8)+igc) |
---|
1304 | iprsm = iprsm + 1 |
---|
1305 | sumk = sumk + kao_mn2o(jn,jt,iprsm)*rwgt(iprsm+128) |
---|
1306 | enddo |
---|
1307 | ka_mn2o(jn,jt,igc) = sumk |
---|
1308 | enddo |
---|
1309 | enddo |
---|
1310 | enddo |
---|
1311 | |
---|
1312 | do jt = 1,19 |
---|
1313 | iprsm = 0 |
---|
1314 | do igc = 1,ngc(9) |
---|
1315 | sumk = 0. |
---|
1316 | do ipr = 1, ngn(ngs(8)+igc) |
---|
1317 | iprsm = iprsm + 1 |
---|
1318 | sumk = sumk + kbo_mn2o(jt,iprsm)*rwgt(iprsm+128) |
---|
1319 | enddo |
---|
1320 | kb_mn2o(jt,igc) = sumk |
---|
1321 | enddo |
---|
1322 | enddo |
---|
1323 | |
---|
1324 | do jt = 1,10 |
---|
1325 | iprsm = 0 |
---|
1326 | do igc = 1,ngc(9) |
---|
1327 | sumk = 0. |
---|
1328 | do ipr = 1, ngn(ngs(8)+igc) |
---|
1329 | iprsm = iprsm + 1 |
---|
1330 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+128) |
---|
1331 | enddo |
---|
1332 | selfref(jt,igc) = sumk |
---|
1333 | enddo |
---|
1334 | enddo |
---|
1335 | |
---|
1336 | do jt = 1,4 |
---|
1337 | iprsm = 0 |
---|
1338 | do igc = 1,ngc(9) |
---|
1339 | sumk = 0. |
---|
1340 | do ipr = 1, ngn(ngs(8)+igc) |
---|
1341 | iprsm = iprsm + 1 |
---|
1342 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+128) |
---|
1343 | enddo |
---|
1344 | forref(jt,igc) = sumk |
---|
1345 | enddo |
---|
1346 | enddo |
---|
1347 | |
---|
1348 | do jp = 1,9 |
---|
1349 | iprsm = 0 |
---|
1350 | do igc = 1,ngc(9) |
---|
1351 | sumf = 0. |
---|
1352 | do ipr = 1, ngn(ngs(8)+igc) |
---|
1353 | iprsm = iprsm + 1 |
---|
1354 | sumf = sumf + fracrefao(iprsm,jp) |
---|
1355 | enddo |
---|
1356 | fracrefa(igc,jp) = sumf |
---|
1357 | enddo |
---|
1358 | enddo |
---|
1359 | |
---|
1360 | iprsm = 0 |
---|
1361 | do igc = 1,ngc(9) |
---|
1362 | sumf = 0. |
---|
1363 | do ipr = 1, ngn(ngs(8)+igc) |
---|
1364 | iprsm = iprsm + 1 |
---|
1365 | sumf = sumf + fracrefbo(iprsm) |
---|
1366 | enddo |
---|
1367 | fracrefb(igc) = sumf |
---|
1368 | enddo |
---|
1369 | |
---|
1370 | end subroutine cmbgb9 |
---|
1371 | |
---|
1372 | !*************************************************************************** |
---|
1373 | subroutine cmbgb10 |
---|
1374 | !*************************************************************************** |
---|
1375 | ! |
---|
1376 | ! band 10: 1390-1480 cm-1 (low key - h2o; high key - h2o) |
---|
1377 | ! |
---|
1378 | ! old band 10: 1390-1480 cm-1 (low - h2o; high - h2o) |
---|
1379 | !*************************************************************************** |
---|
1380 | |
---|
1381 | use parrrtm, only : mg, nbndlw, ngptlw, ng10 |
---|
1382 | use rrlw_kg10, only: fracrefao, fracrefbo, kao, kbo, & |
---|
1383 | selfrefo, forrefo, & |
---|
1384 | fracrefa, fracrefb, absa, ka, absb, kb, & |
---|
1385 | selfref, forref |
---|
1386 | |
---|
1387 | ! ------- Local ------- |
---|
1388 | integer(kind=im) :: jt, jp, igc, ipr, iprsm |
---|
1389 | real(kind=rb) :: sumk, sumf1, sumf2 |
---|
1390 | |
---|
1391 | |
---|
1392 | do jt = 1,5 |
---|
1393 | do jp = 1,13 |
---|
1394 | iprsm = 0 |
---|
1395 | do igc = 1,ngc(10) |
---|
1396 | sumk = 0. |
---|
1397 | do ipr = 1, ngn(ngs(9)+igc) |
---|
1398 | iprsm = iprsm + 1 |
---|
1399 | sumk = sumk + kao(jt,jp,iprsm)*rwgt(iprsm+144) |
---|
1400 | enddo |
---|
1401 | ka(jt,jp,igc) = sumk |
---|
1402 | enddo |
---|
1403 | enddo |
---|
1404 | enddo |
---|
1405 | |
---|
1406 | do jt = 1,5 |
---|
1407 | do jp = 13,59 |
---|
1408 | iprsm = 0 |
---|
1409 | do igc = 1,ngc(10) |
---|
1410 | sumk = 0. |
---|
1411 | do ipr = 1, ngn(ngs(9)+igc) |
---|
1412 | iprsm = iprsm + 1 |
---|
1413 | sumk = sumk + kbo(jt,jp,iprsm)*rwgt(iprsm+144) |
---|
1414 | enddo |
---|
1415 | kb(jt,jp,igc) = sumk |
---|
1416 | enddo |
---|
1417 | enddo |
---|
1418 | enddo |
---|
1419 | |
---|
1420 | do jt = 1,10 |
---|
1421 | iprsm = 0 |
---|
1422 | do igc = 1,ngc(10) |
---|
1423 | sumk = 0. |
---|
1424 | do ipr = 1, ngn(ngs(9)+igc) |
---|
1425 | iprsm = iprsm + 1 |
---|
1426 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+144) |
---|
1427 | enddo |
---|
1428 | selfref(jt,igc) = sumk |
---|
1429 | enddo |
---|
1430 | enddo |
---|
1431 | |
---|
1432 | do jt = 1,4 |
---|
1433 | iprsm = 0 |
---|
1434 | do igc = 1,ngc(10) |
---|
1435 | sumk = 0. |
---|
1436 | do ipr = 1, ngn(ngs(9)+igc) |
---|
1437 | iprsm = iprsm + 1 |
---|
1438 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+144) |
---|
1439 | enddo |
---|
1440 | forref(jt,igc) = sumk |
---|
1441 | enddo |
---|
1442 | enddo |
---|
1443 | |
---|
1444 | iprsm = 0 |
---|
1445 | do igc = 1,ngc(10) |
---|
1446 | sumf1= 0. |
---|
1447 | sumf2= 0. |
---|
1448 | do ipr = 1, ngn(ngs(9)+igc) |
---|
1449 | iprsm = iprsm + 1 |
---|
1450 | sumf1= sumf1+ fracrefao(iprsm) |
---|
1451 | sumf2= sumf2+ fracrefbo(iprsm) |
---|
1452 | enddo |
---|
1453 | fracrefa(igc) = sumf1 |
---|
1454 | fracrefb(igc) = sumf2 |
---|
1455 | enddo |
---|
1456 | |
---|
1457 | end subroutine cmbgb10 |
---|
1458 | |
---|
1459 | !*************************************************************************** |
---|
1460 | subroutine cmbgb11 |
---|
1461 | !*************************************************************************** |
---|
1462 | ! |
---|
1463 | ! band 11: 1480-1800 cm-1 (low - h2o; low minor - o2) |
---|
1464 | ! (high key - h2o; high minor - o2) |
---|
1465 | ! |
---|
1466 | ! old band 11: 1480-1800 cm-1 (low - h2o; low minor - o2) |
---|
1467 | ! (high key - h2o; high minor - o2) |
---|
1468 | !*************************************************************************** |
---|
1469 | |
---|
1470 | use parrrtm, only : mg, nbndlw, ngptlw, ng11 |
---|
1471 | use rrlw_kg11, only: fracrefao, fracrefbo, kao, kao_mo2, & |
---|
1472 | kbo, kbo_mo2, selfrefo, forrefo, & |
---|
1473 | fracrefa, fracrefb, absa, ka, ka_mo2, & |
---|
1474 | absb, kb, kb_mo2, selfref, forref |
---|
1475 | |
---|
1476 | ! ------- Local ------- |
---|
1477 | integer(kind=im) :: jt, jp, igc, ipr, iprsm |
---|
1478 | real(kind=rb) :: sumk, sumk1, sumk2, sumf1, sumf2 |
---|
1479 | |
---|
1480 | |
---|
1481 | do jt = 1,5 |
---|
1482 | do jp = 1,13 |
---|
1483 | iprsm = 0 |
---|
1484 | do igc = 1,ngc(11) |
---|
1485 | sumk = 0. |
---|
1486 | do ipr = 1, ngn(ngs(10)+igc) |
---|
1487 | iprsm = iprsm + 1 |
---|
1488 | sumk = sumk + kao(jt,jp,iprsm)*rwgt(iprsm+160) |
---|
1489 | enddo |
---|
1490 | ka(jt,jp,igc) = sumk |
---|
1491 | enddo |
---|
1492 | enddo |
---|
1493 | enddo |
---|
1494 | do jt = 1,5 |
---|
1495 | do jp = 13,59 |
---|
1496 | iprsm = 0 |
---|
1497 | do igc = 1,ngc(11) |
---|
1498 | sumk = 0. |
---|
1499 | do ipr = 1, ngn(ngs(10)+igc) |
---|
1500 | iprsm = iprsm + 1 |
---|
1501 | sumk = sumk + kbo(jt,jp,iprsm)*rwgt(iprsm+160) |
---|
1502 | enddo |
---|
1503 | kb(jt,jp,igc) = sumk |
---|
1504 | enddo |
---|
1505 | enddo |
---|
1506 | enddo |
---|
1507 | |
---|
1508 | do jt = 1,19 |
---|
1509 | iprsm = 0 |
---|
1510 | do igc = 1,ngc(11) |
---|
1511 | sumk1 = 0. |
---|
1512 | sumk2 = 0. |
---|
1513 | do ipr = 1, ngn(ngs(10)+igc) |
---|
1514 | iprsm = iprsm + 1 |
---|
1515 | sumk1 = sumk1 + kao_mo2(jt,iprsm)*rwgt(iprsm+160) |
---|
1516 | sumk2 = sumk2 + kbo_mo2(jt,iprsm)*rwgt(iprsm+160) |
---|
1517 | enddo |
---|
1518 | ka_mo2(jt,igc) = sumk1 |
---|
1519 | kb_mo2(jt,igc) = sumk2 |
---|
1520 | enddo |
---|
1521 | enddo |
---|
1522 | |
---|
1523 | do jt = 1,10 |
---|
1524 | iprsm = 0 |
---|
1525 | do igc = 1,ngc(11) |
---|
1526 | sumk = 0. |
---|
1527 | do ipr = 1, ngn(ngs(10)+igc) |
---|
1528 | iprsm = iprsm + 1 |
---|
1529 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+160) |
---|
1530 | enddo |
---|
1531 | selfref(jt,igc) = sumk |
---|
1532 | enddo |
---|
1533 | enddo |
---|
1534 | |
---|
1535 | do jt = 1,4 |
---|
1536 | iprsm = 0 |
---|
1537 | do igc = 1,ngc(11) |
---|
1538 | sumk = 0. |
---|
1539 | do ipr = 1, ngn(ngs(10)+igc) |
---|
1540 | iprsm = iprsm + 1 |
---|
1541 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+160) |
---|
1542 | enddo |
---|
1543 | forref(jt,igc) = sumk |
---|
1544 | enddo |
---|
1545 | enddo |
---|
1546 | |
---|
1547 | iprsm = 0 |
---|
1548 | do igc = 1,ngc(11) |
---|
1549 | sumf1= 0. |
---|
1550 | sumf2= 0. |
---|
1551 | do ipr = 1, ngn(ngs(10)+igc) |
---|
1552 | iprsm = iprsm + 1 |
---|
1553 | sumf1= sumf1+ fracrefao(iprsm) |
---|
1554 | sumf2= sumf2+ fracrefbo(iprsm) |
---|
1555 | enddo |
---|
1556 | fracrefa(igc) = sumf1 |
---|
1557 | fracrefb(igc) = sumf2 |
---|
1558 | enddo |
---|
1559 | |
---|
1560 | end subroutine cmbgb11 |
---|
1561 | |
---|
1562 | !*************************************************************************** |
---|
1563 | subroutine cmbgb12 |
---|
1564 | !*************************************************************************** |
---|
1565 | ! |
---|
1566 | ! band 12: 1800-2080 cm-1 (low - h2o,co2; high - nothing) |
---|
1567 | ! |
---|
1568 | ! old band 12: 1800-2080 cm-1 (low - h2o,co2; high - nothing) |
---|
1569 | !*************************************************************************** |
---|
1570 | |
---|
1571 | use parrrtm, only : mg, nbndlw, ngptlw, ng12 |
---|
1572 | use rrlw_kg12, only: fracrefao, kao, selfrefo, forrefo, & |
---|
1573 | fracrefa, absa, ka, selfref, forref |
---|
1574 | |
---|
1575 | ! ------- Local ------- |
---|
1576 | integer(kind=im) :: jn, jt, jp, igc, ipr, iprsm |
---|
1577 | real(kind=rb) :: sumk, sumf |
---|
1578 | |
---|
1579 | |
---|
1580 | do jn = 1,9 |
---|
1581 | do jt = 1,5 |
---|
1582 | do jp = 1,13 |
---|
1583 | iprsm = 0 |
---|
1584 | do igc = 1,ngc(12) |
---|
1585 | sumk = 0. |
---|
1586 | do ipr = 1, ngn(ngs(11)+igc) |
---|
1587 | iprsm = iprsm + 1 |
---|
1588 | sumk = sumk + kao(jn,jt,jp,iprsm)*rwgt(iprsm+176) |
---|
1589 | enddo |
---|
1590 | ka(jn,jt,jp,igc) = sumk |
---|
1591 | enddo |
---|
1592 | enddo |
---|
1593 | enddo |
---|
1594 | enddo |
---|
1595 | |
---|
1596 | do jt = 1,10 |
---|
1597 | iprsm = 0 |
---|
1598 | do igc = 1,ngc(12) |
---|
1599 | sumk = 0. |
---|
1600 | do ipr = 1, ngn(ngs(11)+igc) |
---|
1601 | iprsm = iprsm + 1 |
---|
1602 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+176) |
---|
1603 | enddo |
---|
1604 | selfref(jt,igc) = sumk |
---|
1605 | enddo |
---|
1606 | enddo |
---|
1607 | |
---|
1608 | do jt = 1,4 |
---|
1609 | iprsm = 0 |
---|
1610 | do igc = 1,ngc(12) |
---|
1611 | sumk = 0. |
---|
1612 | do ipr = 1, ngn(ngs(11)+igc) |
---|
1613 | iprsm = iprsm + 1 |
---|
1614 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+176) |
---|
1615 | enddo |
---|
1616 | forref(jt,igc) = sumk |
---|
1617 | enddo |
---|
1618 | enddo |
---|
1619 | |
---|
1620 | do jp = 1,9 |
---|
1621 | iprsm = 0 |
---|
1622 | do igc = 1,ngc(12) |
---|
1623 | sumf = 0. |
---|
1624 | do ipr = 1, ngn(ngs(11)+igc) |
---|
1625 | iprsm = iprsm + 1 |
---|
1626 | sumf = sumf + fracrefao(iprsm,jp) |
---|
1627 | enddo |
---|
1628 | fracrefa(igc,jp) = sumf |
---|
1629 | enddo |
---|
1630 | enddo |
---|
1631 | |
---|
1632 | end subroutine cmbgb12 |
---|
1633 | |
---|
1634 | !*************************************************************************** |
---|
1635 | subroutine cmbgb13 |
---|
1636 | !*************************************************************************** |
---|
1637 | ! |
---|
1638 | ! band 13: 2080-2250 cm-1 (low key - h2o,n2o; high minor - o3 minor) |
---|
1639 | ! |
---|
1640 | ! old band 13: 2080-2250 cm-1 (low - h2o,n2o; high - nothing) |
---|
1641 | !*************************************************************************** |
---|
1642 | |
---|
1643 | use parrrtm, only : mg, nbndlw, ngptlw, ng13 |
---|
1644 | use rrlw_kg13, only: fracrefao, fracrefbo, kao, kao_mco2, kao_mco, & |
---|
1645 | kbo_mo3, selfrefo, forrefo, & |
---|
1646 | fracrefa, fracrefb, absa, ka, ka_mco2, ka_mco, & |
---|
1647 | kb_mo3, selfref, forref |
---|
1648 | |
---|
1649 | ! ------- Local ------- |
---|
1650 | integer(kind=im) :: jn, jt, jp, igc, ipr, iprsm |
---|
1651 | real(kind=rb) :: sumk, sumk1, sumk2, sumf |
---|
1652 | |
---|
1653 | |
---|
1654 | do jn = 1,9 |
---|
1655 | do jt = 1,5 |
---|
1656 | do jp = 1,13 |
---|
1657 | iprsm = 0 |
---|
1658 | do igc = 1,ngc(13) |
---|
1659 | sumk = 0. |
---|
1660 | do ipr = 1, ngn(ngs(12)+igc) |
---|
1661 | iprsm = iprsm + 1 |
---|
1662 | sumk = sumk + kao(jn,jt,jp,iprsm)*rwgt(iprsm+192) |
---|
1663 | enddo |
---|
1664 | ka(jn,jt,jp,igc) = sumk |
---|
1665 | enddo |
---|
1666 | enddo |
---|
1667 | enddo |
---|
1668 | enddo |
---|
1669 | |
---|
1670 | do jn = 1,9 |
---|
1671 | do jt = 1,19 |
---|
1672 | iprsm = 0 |
---|
1673 | do igc = 1,ngc(13) |
---|
1674 | sumk1 = 0. |
---|
1675 | sumk2 = 0. |
---|
1676 | do ipr = 1, ngn(ngs(12)+igc) |
---|
1677 | iprsm = iprsm + 1 |
---|
1678 | sumk1 = sumk1 + kao_mco2(jn,jt,iprsm)*rwgt(iprsm+192) |
---|
1679 | sumk2 = sumk2 + kao_mco(jn,jt,iprsm)*rwgt(iprsm+192) |
---|
1680 | enddo |
---|
1681 | ka_mco2(jn,jt,igc) = sumk1 |
---|
1682 | ka_mco(jn,jt,igc) = sumk2 |
---|
1683 | enddo |
---|
1684 | enddo |
---|
1685 | enddo |
---|
1686 | |
---|
1687 | do jt = 1,19 |
---|
1688 | iprsm = 0 |
---|
1689 | do igc = 1,ngc(13) |
---|
1690 | sumk = 0. |
---|
1691 | do ipr = 1, ngn(ngs(12)+igc) |
---|
1692 | iprsm = iprsm + 1 |
---|
1693 | sumk = sumk + kbo_mo3(jt,iprsm)*rwgt(iprsm+192) |
---|
1694 | enddo |
---|
1695 | kb_mo3(jt,igc) = sumk |
---|
1696 | enddo |
---|
1697 | enddo |
---|
1698 | |
---|
1699 | do jt = 1,10 |
---|
1700 | iprsm = 0 |
---|
1701 | do igc = 1,ngc(13) |
---|
1702 | sumk = 0. |
---|
1703 | do ipr = 1, ngn(ngs(12)+igc) |
---|
1704 | iprsm = iprsm + 1 |
---|
1705 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+192) |
---|
1706 | enddo |
---|
1707 | selfref(jt,igc) = sumk |
---|
1708 | enddo |
---|
1709 | enddo |
---|
1710 | |
---|
1711 | do jt = 1,4 |
---|
1712 | iprsm = 0 |
---|
1713 | do igc = 1,ngc(13) |
---|
1714 | sumk = 0. |
---|
1715 | do ipr = 1, ngn(ngs(12)+igc) |
---|
1716 | iprsm = iprsm + 1 |
---|
1717 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+192) |
---|
1718 | enddo |
---|
1719 | forref(jt,igc) = sumk |
---|
1720 | enddo |
---|
1721 | enddo |
---|
1722 | |
---|
1723 | iprsm = 0 |
---|
1724 | do igc = 1,ngc(13) |
---|
1725 | sumf = 0. |
---|
1726 | do ipr = 1, ngn(ngs(12)+igc) |
---|
1727 | iprsm = iprsm + 1 |
---|
1728 | sumf = sumf + fracrefbo(iprsm) |
---|
1729 | enddo |
---|
1730 | fracrefb(igc) = sumf |
---|
1731 | enddo |
---|
1732 | |
---|
1733 | do jp = 1,9 |
---|
1734 | iprsm = 0 |
---|
1735 | do igc = 1,ngc(13) |
---|
1736 | sumf = 0. |
---|
1737 | do ipr = 1, ngn(ngs(12)+igc) |
---|
1738 | iprsm = iprsm + 1 |
---|
1739 | sumf = sumf + fracrefao(iprsm,jp) |
---|
1740 | enddo |
---|
1741 | fracrefa(igc,jp) = sumf |
---|
1742 | enddo |
---|
1743 | enddo |
---|
1744 | |
---|
1745 | end subroutine cmbgb13 |
---|
1746 | |
---|
1747 | !*************************************************************************** |
---|
1748 | subroutine cmbgb14 |
---|
1749 | !*************************************************************************** |
---|
1750 | ! |
---|
1751 | ! band 14: 2250-2380 cm-1 (low - co2; high - co2) |
---|
1752 | ! |
---|
1753 | ! old band 14: 2250-2380 cm-1 (low - co2; high - co2) |
---|
1754 | !*************************************************************************** |
---|
1755 | |
---|
1756 | use parrrtm, only : mg, nbndlw, ngptlw, ng14 |
---|
1757 | use rrlw_kg14, only: fracrefao, fracrefbo, kao, kbo, & |
---|
1758 | selfrefo, forrefo, & |
---|
1759 | fracrefa, fracrefb, absa, ka, absb, kb, & |
---|
1760 | selfref, forref |
---|
1761 | |
---|
1762 | ! ------- Local ------- |
---|
1763 | integer(kind=im) :: jt, jp, igc, ipr, iprsm |
---|
1764 | real(kind=rb) :: sumk, sumf1, sumf2 |
---|
1765 | |
---|
1766 | |
---|
1767 | do jt = 1,5 |
---|
1768 | do jp = 1,13 |
---|
1769 | iprsm = 0 |
---|
1770 | do igc = 1,ngc(14) |
---|
1771 | sumk = 0. |
---|
1772 | do ipr = 1, ngn(ngs(13)+igc) |
---|
1773 | iprsm = iprsm + 1 |
---|
1774 | sumk = sumk + kao(jt,jp,iprsm)*rwgt(iprsm+208) |
---|
1775 | enddo |
---|
1776 | ka(jt,jp,igc) = sumk |
---|
1777 | enddo |
---|
1778 | enddo |
---|
1779 | enddo |
---|
1780 | |
---|
1781 | do jt = 1,5 |
---|
1782 | do jp = 13,59 |
---|
1783 | iprsm = 0 |
---|
1784 | do igc = 1,ngc(14) |
---|
1785 | sumk = 0. |
---|
1786 | do ipr = 1, ngn(ngs(13)+igc) |
---|
1787 | iprsm = iprsm + 1 |
---|
1788 | sumk = sumk + kbo(jt,jp,iprsm)*rwgt(iprsm+208) |
---|
1789 | enddo |
---|
1790 | kb(jt,jp,igc) = sumk |
---|
1791 | enddo |
---|
1792 | enddo |
---|
1793 | enddo |
---|
1794 | |
---|
1795 | do jt = 1,10 |
---|
1796 | iprsm = 0 |
---|
1797 | do igc = 1,ngc(14) |
---|
1798 | sumk = 0. |
---|
1799 | do ipr = 1, ngn(ngs(13)+igc) |
---|
1800 | iprsm = iprsm + 1 |
---|
1801 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+208) |
---|
1802 | enddo |
---|
1803 | selfref(jt,igc) = sumk |
---|
1804 | enddo |
---|
1805 | enddo |
---|
1806 | |
---|
1807 | do jt = 1,4 |
---|
1808 | iprsm = 0 |
---|
1809 | do igc = 1,ngc(14) |
---|
1810 | sumk = 0. |
---|
1811 | do ipr = 1, ngn(ngs(13)+igc) |
---|
1812 | iprsm = iprsm + 1 |
---|
1813 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+208) |
---|
1814 | enddo |
---|
1815 | forref(jt,igc) = sumk |
---|
1816 | enddo |
---|
1817 | enddo |
---|
1818 | |
---|
1819 | iprsm = 0 |
---|
1820 | do igc = 1,ngc(14) |
---|
1821 | sumf1= 0. |
---|
1822 | sumf2= 0. |
---|
1823 | do ipr = 1, ngn(ngs(13)+igc) |
---|
1824 | iprsm = iprsm + 1 |
---|
1825 | sumf1= sumf1+ fracrefao(iprsm) |
---|
1826 | sumf2= sumf2+ fracrefbo(iprsm) |
---|
1827 | enddo |
---|
1828 | fracrefa(igc) = sumf1 |
---|
1829 | fracrefb(igc) = sumf2 |
---|
1830 | enddo |
---|
1831 | |
---|
1832 | end subroutine cmbgb14 |
---|
1833 | |
---|
1834 | !*************************************************************************** |
---|
1835 | subroutine cmbgb15 |
---|
1836 | !*************************************************************************** |
---|
1837 | ! |
---|
1838 | ! band 15: 2380-2600 cm-1 (low - n2o,co2; low minor - n2) |
---|
1839 | ! (high - nothing) |
---|
1840 | ! |
---|
1841 | ! old band 15: 2380-2600 cm-1 (low - n2o,co2; high - nothing) |
---|
1842 | !*************************************************************************** |
---|
1843 | |
---|
1844 | use parrrtm, only : mg, nbndlw, ngptlw, ng15 |
---|
1845 | use rrlw_kg15, only: fracrefao, kao, kao_mn2, selfrefo, forrefo, & |
---|
1846 | fracrefa, absa, ka, ka_mn2, selfref, forref |
---|
1847 | |
---|
1848 | ! ------- Local ------- |
---|
1849 | integer(kind=im) :: jn, jt, jp, igc, ipr, iprsm |
---|
1850 | real(kind=rb) :: sumk, sumf |
---|
1851 | |
---|
1852 | |
---|
1853 | do jn = 1,9 |
---|
1854 | do jt = 1,5 |
---|
1855 | do jp = 1,13 |
---|
1856 | iprsm = 0 |
---|
1857 | do igc = 1,ngc(15) |
---|
1858 | sumk = 0. |
---|
1859 | do ipr = 1, ngn(ngs(14)+igc) |
---|
1860 | iprsm = iprsm + 1 |
---|
1861 | sumk = sumk + kao(jn,jt,jp,iprsm)*rwgt(iprsm+224) |
---|
1862 | enddo |
---|
1863 | ka(jn,jt,jp,igc) = sumk |
---|
1864 | enddo |
---|
1865 | enddo |
---|
1866 | enddo |
---|
1867 | enddo |
---|
1868 | |
---|
1869 | do jn = 1,9 |
---|
1870 | do jt = 1,19 |
---|
1871 | iprsm = 0 |
---|
1872 | do igc = 1,ngc(15) |
---|
1873 | sumk = 0. |
---|
1874 | do ipr = 1, ngn(ngs(14)+igc) |
---|
1875 | iprsm = iprsm + 1 |
---|
1876 | sumk = sumk + kao_mn2(jn,jt,iprsm)*rwgt(iprsm+224) |
---|
1877 | enddo |
---|
1878 | ka_mn2(jn,jt,igc) = sumk |
---|
1879 | enddo |
---|
1880 | enddo |
---|
1881 | enddo |
---|
1882 | |
---|
1883 | do jt = 1,10 |
---|
1884 | iprsm = 0 |
---|
1885 | do igc = 1,ngc(15) |
---|
1886 | sumk = 0. |
---|
1887 | do ipr = 1, ngn(ngs(14)+igc) |
---|
1888 | iprsm = iprsm + 1 |
---|
1889 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+224) |
---|
1890 | enddo |
---|
1891 | selfref(jt,igc) = sumk |
---|
1892 | enddo |
---|
1893 | enddo |
---|
1894 | |
---|
1895 | do jt = 1,4 |
---|
1896 | iprsm = 0 |
---|
1897 | do igc = 1,ngc(15) |
---|
1898 | sumk = 0. |
---|
1899 | do ipr = 1, ngn(ngs(14)+igc) |
---|
1900 | iprsm = iprsm + 1 |
---|
1901 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+224) |
---|
1902 | enddo |
---|
1903 | forref(jt,igc) = sumk |
---|
1904 | enddo |
---|
1905 | enddo |
---|
1906 | |
---|
1907 | do jp = 1,9 |
---|
1908 | iprsm = 0 |
---|
1909 | do igc = 1,ngc(15) |
---|
1910 | sumf = 0. |
---|
1911 | do ipr = 1, ngn(ngs(14)+igc) |
---|
1912 | iprsm = iprsm + 1 |
---|
1913 | sumf = sumf + fracrefao(iprsm,jp) |
---|
1914 | enddo |
---|
1915 | fracrefa(igc,jp) = sumf |
---|
1916 | enddo |
---|
1917 | enddo |
---|
1918 | |
---|
1919 | end subroutine cmbgb15 |
---|
1920 | |
---|
1921 | !*************************************************************************** |
---|
1922 | subroutine cmbgb16 |
---|
1923 | !*************************************************************************** |
---|
1924 | ! |
---|
1925 | ! band 16: 2600-3250 cm-1 (low key- h2o,ch4; high key - ch4) |
---|
1926 | ! |
---|
1927 | ! old band 16: 2600-3000 cm-1 (low - h2o,ch4; high - nothing) |
---|
1928 | !*************************************************************************** |
---|
1929 | |
---|
1930 | use parrrtm, only : mg, nbndlw, ngptlw, ng16 |
---|
1931 | use rrlw_kg16, only: fracrefao, fracrefbo, kao, kbo, selfrefo, forrefo, & |
---|
1932 | fracrefa, fracrefb, absa, ka, absb, kb, selfref, forref |
---|
1933 | |
---|
1934 | ! ------- Local ------- |
---|
1935 | integer(kind=im) :: jn, jt, jp, igc, ipr, iprsm |
---|
1936 | real(kind=rb) :: sumk, sumf |
---|
1937 | |
---|
1938 | |
---|
1939 | do jn = 1,9 |
---|
1940 | do jt = 1,5 |
---|
1941 | do jp = 1,13 |
---|
1942 | iprsm = 0 |
---|
1943 | do igc = 1,ngc(16) |
---|
1944 | sumk = 0. |
---|
1945 | do ipr = 1, ngn(ngs(15)+igc) |
---|
1946 | iprsm = iprsm + 1 |
---|
1947 | sumk = sumk + kao(jn,jt,jp,iprsm)*rwgt(iprsm+240) |
---|
1948 | enddo |
---|
1949 | ka(jn,jt,jp,igc) = sumk |
---|
1950 | enddo |
---|
1951 | enddo |
---|
1952 | enddo |
---|
1953 | enddo |
---|
1954 | |
---|
1955 | do jt = 1,5 |
---|
1956 | do jp = 13,59 |
---|
1957 | iprsm = 0 |
---|
1958 | do igc = 1,ngc(16) |
---|
1959 | sumk = 0. |
---|
1960 | do ipr = 1, ngn(ngs(15)+igc) |
---|
1961 | iprsm = iprsm + 1 |
---|
1962 | sumk = sumk + kbo(jt,jp,iprsm)*rwgt(iprsm+240) |
---|
1963 | enddo |
---|
1964 | kb(jt,jp,igc) = sumk |
---|
1965 | enddo |
---|
1966 | enddo |
---|
1967 | enddo |
---|
1968 | |
---|
1969 | do jt = 1,10 |
---|
1970 | iprsm = 0 |
---|
1971 | do igc = 1,ngc(16) |
---|
1972 | sumk = 0. |
---|
1973 | do ipr = 1, ngn(ngs(15)+igc) |
---|
1974 | iprsm = iprsm + 1 |
---|
1975 | sumk = sumk + selfrefo(jt,iprsm)*rwgt(iprsm+240) |
---|
1976 | enddo |
---|
1977 | selfref(jt,igc) = sumk |
---|
1978 | enddo |
---|
1979 | enddo |
---|
1980 | |
---|
1981 | do jt = 1,4 |
---|
1982 | iprsm = 0 |
---|
1983 | do igc = 1,ngc(16) |
---|
1984 | sumk = 0. |
---|
1985 | do ipr = 1, ngn(ngs(15)+igc) |
---|
1986 | iprsm = iprsm + 1 |
---|
1987 | sumk = sumk + forrefo(jt,iprsm)*rwgt(iprsm+240) |
---|
1988 | enddo |
---|
1989 | forref(jt,igc) = sumk |
---|
1990 | enddo |
---|
1991 | enddo |
---|
1992 | |
---|
1993 | iprsm = 0 |
---|
1994 | do igc = 1,ngc(16) |
---|
1995 | sumf = 0. |
---|
1996 | do ipr = 1, ngn(ngs(15)+igc) |
---|
1997 | iprsm = iprsm + 1 |
---|
1998 | sumf = sumf + fracrefbo(iprsm) |
---|
1999 | enddo |
---|
2000 | fracrefb(igc) = sumf |
---|
2001 | enddo |
---|
2002 | |
---|
2003 | do jp = 1,9 |
---|
2004 | iprsm = 0 |
---|
2005 | do igc = 1,ngc(16) |
---|
2006 | sumf = 0. |
---|
2007 | do ipr = 1, ngn(ngs(15)+igc) |
---|
2008 | iprsm = iprsm + 1 |
---|
2009 | sumf = sumf + fracrefao(iprsm,jp) |
---|
2010 | enddo |
---|
2011 | fracrefa(igc,jp) = sumf |
---|
2012 | enddo |
---|
2013 | enddo |
---|
2014 | |
---|
2015 | end subroutine cmbgb16 |
---|
2016 | |
---|
2017 | !*************************************************************************** |
---|
2018 | subroutine lwcldpr |
---|
2019 | !*************************************************************************** |
---|
2020 | |
---|
2021 | ! --------- Modules ---------- |
---|
2022 | |
---|
2023 | use rrlw_cld, only: abscld1, absliq0, absliq1, & |
---|
2024 | absice0, absice1, absice2, absice3 |
---|
2025 | |
---|
2026 | save |
---|
2027 | |
---|
2028 | ! ABSCLDn is the liquid water absorption coefficient (m2/g). |
---|
2029 | ! For INFLAG = 1. |
---|
2030 | abscld1 = 0.0602410_rb |
---|
2031 | ! |
---|
2032 | ! Everything below is for INFLAG = 2. |
---|
2033 | |
---|
2034 | ! ABSICEn(J,IB) are the parameters needed to compute the liquid water |
---|
2035 | ! absorption coefficient in spectral region IB for ICEFLAG=n. The units |
---|
2036 | ! of ABSICEn(1,IB) are m2/g and ABSICEn(2,IB) has units (microns (m2/g)). |
---|
2037 | ! For ICEFLAG = 0. |
---|
2038 | |
---|
2039 | absice0(:)= (/0.005_rb, 1.0_rb/) |
---|
2040 | |
---|
2041 | ! For ICEFLAG = 1. |
---|
2042 | absice1(1,:) = (/0.0036_rb, 0.0068_rb, 0.0003_rb, 0.0016_rb, 0.0020_rb/) |
---|
2043 | absice1(2,:) = (/1.136_rb , 0.600_rb , 1.338_rb , 1.166_rb , 1.118_rb /) |
---|
2044 | |
---|
2045 | ! For ICEFLAG = 2. In each band, the absorption |
---|
2046 | ! coefficients are listed for a range of effective radii from 5.0 |
---|
2047 | ! to 131.0 microns in increments of 3.0 microns. |
---|
2048 | ! Spherical Ice Particle Parameterization |
---|
2049 | ! absorption units (abs coef/iwc): [(m^-1)/(g m^-3)] |
---|
2050 | absice2(:,1) = (/ & |
---|
2051 | ! band 1 |
---|
2052 | 7.798999e-02_rb,6.340479e-02_rb,5.417973e-02_rb,4.766245e-02_rb,4.272663e-02_rb, & |
---|
2053 | 3.880939e-02_rb,3.559544e-02_rb,3.289241e-02_rb,3.057511e-02_rb,2.855800e-02_rb, & |
---|
2054 | 2.678022e-02_rb,2.519712e-02_rb,2.377505e-02_rb,2.248806e-02_rb,2.131578e-02_rb, & |
---|
2055 | 2.024194e-02_rb,1.925337e-02_rb,1.833926e-02_rb,1.749067e-02_rb,1.670007e-02_rb, & |
---|
2056 | 1.596113e-02_rb,1.526845e-02_rb,1.461739e-02_rb,1.400394e-02_rb,1.342462e-02_rb, & |
---|
2057 | 1.287639e-02_rb,1.235656e-02_rb,1.186279e-02_rb,1.139297e-02_rb,1.094524e-02_rb, & |
---|
2058 | 1.051794e-02_rb,1.010956e-02_rb,9.718755e-03_rb,9.344316e-03_rb,8.985139e-03_rb, & |
---|
2059 | 8.640223e-03_rb,8.308656e-03_rb,7.989606e-03_rb,7.682312e-03_rb,7.386076e-03_rb, & |
---|
2060 | 7.100255e-03_rb,6.824258e-03_rb,6.557540e-03_rb/) |
---|
2061 | absice2(:,2) = (/ & |
---|
2062 | ! band 2 |
---|
2063 | 2.784879e-02_rb,2.709863e-02_rb,2.619165e-02_rb,2.529230e-02_rb,2.443225e-02_rb, & |
---|
2064 | 2.361575e-02_rb,2.284021e-02_rb,2.210150e-02_rb,2.139548e-02_rb,2.071840e-02_rb, & |
---|
2065 | 2.006702e-02_rb,1.943856e-02_rb,1.883064e-02_rb,1.824120e-02_rb,1.766849e-02_rb, & |
---|
2066 | 1.711099e-02_rb,1.656737e-02_rb,1.603647e-02_rb,1.551727e-02_rb,1.500886e-02_rb, & |
---|
2067 | 1.451045e-02_rb,1.402132e-02_rb,1.354084e-02_rb,1.306842e-02_rb,1.260355e-02_rb, & |
---|
2068 | 1.214575e-02_rb,1.169460e-02_rb,1.124971e-02_rb,1.081072e-02_rb,1.037731e-02_rb, & |
---|
2069 | 9.949167e-03_rb,9.526021e-03_rb,9.107615e-03_rb,8.693714e-03_rb,8.284096e-03_rb, & |
---|
2070 | 7.878558e-03_rb,7.476910e-03_rb,7.078974e-03_rb,6.684586e-03_rb,6.293589e-03_rb, & |
---|
2071 | 5.905839e-03_rb,5.521200e-03_rb,5.139543e-03_rb/) |
---|
2072 | absice2(:,3) = (/ & |
---|
2073 | ! band 3 |
---|
2074 | 1.065397e-01_rb,8.005726e-02_rb,6.546428e-02_rb,5.589131e-02_rb,4.898681e-02_rb, & |
---|
2075 | 4.369932e-02_rb,3.947901e-02_rb,3.600676e-02_rb,3.308299e-02_rb,3.057561e-02_rb, & |
---|
2076 | 2.839325e-02_rb,2.647040e-02_rb,2.475872e-02_rb,2.322164e-02_rb,2.183091e-02_rb, & |
---|
2077 | 2.056430e-02_rb,1.940407e-02_rb,1.833586e-02_rb,1.734787e-02_rb,1.643034e-02_rb, & |
---|
2078 | 1.557512e-02_rb,1.477530e-02_rb,1.402501e-02_rb,1.331924e-02_rb,1.265364e-02_rb, & |
---|
2079 | 1.202445e-02_rb,1.142838e-02_rb,1.086257e-02_rb,1.032445e-02_rb,9.811791e-03_rb, & |
---|
2080 | 9.322587e-03_rb,8.855053e-03_rb,8.407591e-03_rb,7.978763e-03_rb,7.567273e-03_rb, & |
---|
2081 | 7.171949e-03_rb,6.791728e-03_rb,6.425642e-03_rb,6.072809e-03_rb,5.732424e-03_rb, & |
---|
2082 | 5.403748e-03_rb,5.086103e-03_rb,4.778865e-03_rb/) |
---|
2083 | absice2(:,4) = (/ & |
---|
2084 | ! band 4 |
---|
2085 | 1.804566e-01_rb,1.168987e-01_rb,8.680442e-02_rb,6.910060e-02_rb,5.738174e-02_rb, & |
---|
2086 | 4.902332e-02_rb,4.274585e-02_rb,3.784923e-02_rb,3.391734e-02_rb,3.068690e-02_rb, & |
---|
2087 | 2.798301e-02_rb,2.568480e-02_rb,2.370600e-02_rb,2.198337e-02_rb,2.046940e-02_rb, & |
---|
2088 | 1.912777e-02_rb,1.793016e-02_rb,1.685420e-02_rb,1.588193e-02_rb,1.499882e-02_rb, & |
---|
2089 | 1.419293e-02_rb,1.345440e-02_rb,1.277496e-02_rb,1.214769e-02_rb,1.156669e-02_rb, & |
---|
2090 | 1.102694e-02_rb,1.052412e-02_rb,1.005451e-02_rb,9.614854e-03_rb,9.202335e-03_rb, & |
---|
2091 | 8.814470e-03_rb,8.449077e-03_rb,8.104223e-03_rb,7.778195e-03_rb,7.469466e-03_rb, & |
---|
2092 | 7.176671e-03_rb,6.898588e-03_rb,6.634117e-03_rb,6.382264e-03_rb,6.142134e-03_rb, & |
---|
2093 | 5.912913e-03_rb,5.693862e-03_rb,5.484308e-03_rb/) |
---|
2094 | absice2(:,5) = (/ & |
---|
2095 | ! band 5 |
---|
2096 | 2.131806e-01_rb,1.311372e-01_rb,9.407171e-02_rb,7.299442e-02_rb,5.941273e-02_rb, & |
---|
2097 | 4.994043e-02_rb,4.296242e-02_rb,3.761113e-02_rb,3.337910e-02_rb,2.994978e-02_rb, & |
---|
2098 | 2.711556e-02_rb,2.473461e-02_rb,2.270681e-02_rb,2.095943e-02_rb,1.943839e-02_rb, & |
---|
2099 | 1.810267e-02_rb,1.692057e-02_rb,1.586719e-02_rb,1.492275e-02_rb,1.407132e-02_rb, & |
---|
2100 | 1.329989e-02_rb,1.259780e-02_rb,1.195618e-02_rb,1.136761e-02_rb,1.082583e-02_rb, & |
---|
2101 | 1.032552e-02_rb,9.862158e-03_rb,9.431827e-03_rb,9.031157e-03_rb,8.657217e-03_rb, & |
---|
2102 | 8.307449e-03_rb,7.979609e-03_rb,7.671724e-03_rb,7.382048e-03_rb,7.109032e-03_rb, & |
---|
2103 | 6.851298e-03_rb,6.607615e-03_rb,6.376881e-03_rb,6.158105e-03_rb,5.950394e-03_rb, & |
---|
2104 | 5.752942e-03_rb,5.565019e-03_rb,5.385963e-03_rb/) |
---|
2105 | absice2(:,6) = (/ & |
---|
2106 | ! band 6 |
---|
2107 | 1.546177e-01_rb,1.039251e-01_rb,7.910347e-02_rb,6.412429e-02_rb,5.399997e-02_rb, & |
---|
2108 | 4.664937e-02_rb,4.104237e-02_rb,3.660781e-02_rb,3.300218e-02_rb,3.000586e-02_rb, & |
---|
2109 | 2.747148e-02_rb,2.529633e-02_rb,2.340647e-02_rb,2.174723e-02_rb,2.027731e-02_rb, & |
---|
2110 | 1.896487e-02_rb,1.778492e-02_rb,1.671761e-02_rb,1.574692e-02_rb,1.485978e-02_rb, & |
---|
2111 | 1.404543e-02_rb,1.329489e-02_rb,1.260066e-02_rb,1.195636e-02_rb,1.135657e-02_rb, & |
---|
2112 | 1.079664e-02_rb,1.027257e-02_rb,9.780871e-03_rb,9.318505e-03_rb,8.882815e-03_rb, & |
---|
2113 | 8.471458e-03_rb,8.082364e-03_rb,7.713696e-03_rb,7.363817e-03_rb,7.031264e-03_rb, & |
---|
2114 | 6.714725e-03_rb,6.413021e-03_rb,6.125086e-03_rb,5.849958e-03_rb,5.586764e-03_rb, & |
---|
2115 | 5.334707e-03_rb,5.093066e-03_rb,4.861179e-03_rb/) |
---|
2116 | absice2(:,7) = (/ & |
---|
2117 | ! band 7 |
---|
2118 | 7.583404e-02_rb,6.181558e-02_rb,5.312027e-02_rb,4.696039e-02_rb,4.225986e-02_rb, & |
---|
2119 | 3.849735e-02_rb,3.538340e-02_rb,3.274182e-02_rb,3.045798e-02_rb,2.845343e-02_rb, & |
---|
2120 | 2.667231e-02_rb,2.507353e-02_rb,2.362606e-02_rb,2.230595e-02_rb,2.109435e-02_rb, & |
---|
2121 | 1.997617e-02_rb,1.893916e-02_rb,1.797328e-02_rb,1.707016e-02_rb,1.622279e-02_rb, & |
---|
2122 | 1.542523e-02_rb,1.467241e-02_rb,1.395997e-02_rb,1.328414e-02_rb,1.264164e-02_rb, & |
---|
2123 | 1.202958e-02_rb,1.144544e-02_rb,1.088697e-02_rb,1.035218e-02_rb,9.839297e-03_rb, & |
---|
2124 | 9.346733e-03_rb,8.873057e-03_rb,8.416980e-03_rb,7.977335e-03_rb,7.553066e-03_rb, & |
---|
2125 | 7.143210e-03_rb,6.746888e-03_rb,6.363297e-03_rb,5.991700e-03_rb,5.631422e-03_rb, & |
---|
2126 | 5.281840e-03_rb,4.942378e-03_rb,4.612505e-03_rb/) |
---|
2127 | absice2(:,8) = (/ & |
---|
2128 | ! band 8 |
---|
2129 | 9.022185e-02_rb,6.922700e-02_rb,5.710674e-02_rb,4.898377e-02_rb,4.305946e-02_rb, & |
---|
2130 | 3.849553e-02_rb,3.484183e-02_rb,3.183220e-02_rb,2.929794e-02_rb,2.712627e-02_rb, & |
---|
2131 | 2.523856e-02_rb,2.357810e-02_rb,2.210286e-02_rb,2.078089e-02_rb,1.958747e-02_rb, & |
---|
2132 | 1.850310e-02_rb,1.751218e-02_rb,1.660205e-02_rb,1.576232e-02_rb,1.498440e-02_rb, & |
---|
2133 | 1.426107e-02_rb,1.358624e-02_rb,1.295474e-02_rb,1.236212e-02_rb,1.180456e-02_rb, & |
---|
2134 | 1.127874e-02_rb,1.078175e-02_rb,1.031106e-02_rb,9.864433e-03_rb,9.439878e-03_rb, & |
---|
2135 | 9.035637e-03_rb,8.650140e-03_rb,8.281981e-03_rb,7.929895e-03_rb,7.592746e-03_rb, & |
---|
2136 | 7.269505e-03_rb,6.959238e-03_rb,6.661100e-03_rb,6.374317e-03_rb,6.098185e-03_rb, & |
---|
2137 | 5.832059e-03_rb,5.575347e-03_rb,5.327504e-03_rb/) |
---|
2138 | absice2(:,9) = (/ & |
---|
2139 | ! band 9 |
---|
2140 | 1.294087e-01_rb,8.788217e-02_rb,6.728288e-02_rb,5.479720e-02_rb,4.635049e-02_rb, & |
---|
2141 | 4.022253e-02_rb,3.555576e-02_rb,3.187259e-02_rb,2.888498e-02_rb,2.640843e-02_rb, & |
---|
2142 | 2.431904e-02_rb,2.253038e-02_rb,2.098024e-02_rb,1.962267e-02_rb,1.842293e-02_rb, & |
---|
2143 | 1.735426e-02_rb,1.639571e-02_rb,1.553060e-02_rb,1.474552e-02_rb,1.402953e-02_rb, & |
---|
2144 | 1.337363e-02_rb,1.277033e-02_rb,1.221336e-02_rb,1.169741e-02_rb,1.121797e-02_rb, & |
---|
2145 | 1.077117e-02_rb,1.035369e-02_rb,9.962643e-03_rb,9.595509e-03_rb,9.250088e-03_rb, & |
---|
2146 | 8.924447e-03_rb,8.616876e-03_rb,8.325862e-03_rb,8.050057e-03_rb,7.788258e-03_rb, & |
---|
2147 | 7.539388e-03_rb,7.302478e-03_rb,7.076656e-03_rb,6.861134e-03_rb,6.655197e-03_rb, & |
---|
2148 | 6.458197e-03_rb,6.269543e-03_rb,6.088697e-03_rb/) |
---|
2149 | absice2(:,10) = (/ & |
---|
2150 | ! band 10 |
---|
2151 | 1.593628e-01_rb,1.014552e-01_rb,7.458955e-02_rb,5.903571e-02_rb,4.887582e-02_rb, & |
---|
2152 | 4.171159e-02_rb,3.638480e-02_rb,3.226692e-02_rb,2.898717e-02_rb,2.631256e-02_rb, & |
---|
2153 | 2.408925e-02_rb,2.221156e-02_rb,2.060448e-02_rb,1.921325e-02_rb,1.799699e-02_rb, & |
---|
2154 | 1.692456e-02_rb,1.597177e-02_rb,1.511961e-02_rb,1.435289e-02_rb,1.365933e-02_rb, & |
---|
2155 | 1.302890e-02_rb,1.245334e-02_rb,1.192576e-02_rb,1.144037e-02_rb,1.099230e-02_rb, & |
---|
2156 | 1.057739e-02_rb,1.019208e-02_rb,9.833302e-03_rb,9.498395e-03_rb,9.185047e-03_rb, & |
---|
2157 | 8.891237e-03_rb,8.615185e-03_rb,8.355325e-03_rb,8.110267e-03_rb,7.878778e-03_rb, & |
---|
2158 | 7.659759e-03_rb,7.452224e-03_rb,7.255291e-03_rb,7.068166e-03_rb,6.890130e-03_rb, & |
---|
2159 | 6.720536e-03_rb,6.558794e-03_rb,6.404371e-03_rb/) |
---|
2160 | absice2(:,11) = (/ & |
---|
2161 | ! band 11 |
---|
2162 | 1.656227e-01_rb,1.032129e-01_rb,7.487359e-02_rb,5.871431e-02_rb,4.828355e-02_rb, & |
---|
2163 | 4.099989e-02_rb,3.562924e-02_rb,3.150755e-02_rb,2.824593e-02_rb,2.560156e-02_rb, & |
---|
2164 | 2.341503e-02_rb,2.157740e-02_rb,2.001169e-02_rb,1.866199e-02_rb,1.748669e-02_rb, & |
---|
2165 | 1.645421e-02_rb,1.554015e-02_rb,1.472535e-02_rb,1.399457e-02_rb,1.333553e-02_rb, & |
---|
2166 | 1.273821e-02_rb,1.219440e-02_rb,1.169725e-02_rb,1.124104e-02_rb,1.082096e-02_rb, & |
---|
2167 | 1.043290e-02_rb,1.007336e-02_rb,9.739338e-03_rb,9.428223e-03_rb,9.137756e-03_rb, & |
---|
2168 | 8.865964e-03_rb,8.611115e-03_rb,8.371686e-03_rb,8.146330e-03_rb,7.933852e-03_rb, & |
---|
2169 | 7.733187e-03_rb,7.543386e-03_rb,7.363597e-03_rb,7.193056e-03_rb,7.031072e-03_rb, & |
---|
2170 | 6.877024e-03_rb,6.730348e-03_rb,6.590531e-03_rb/) |
---|
2171 | absice2(:,12) = (/ & |
---|
2172 | ! band 12 |
---|
2173 | 9.194591e-02_rb,6.446867e-02_rb,4.962034e-02_rb,4.042061e-02_rb,3.418456e-02_rb, & |
---|
2174 | 2.968856e-02_rb,2.629900e-02_rb,2.365572e-02_rb,2.153915e-02_rb,1.980791e-02_rb, & |
---|
2175 | 1.836689e-02_rb,1.714979e-02_rb,1.610900e-02_rb,1.520946e-02_rb,1.442476e-02_rb, & |
---|
2176 | 1.373468e-02_rb,1.312345e-02_rb,1.257858e-02_rb,1.209010e-02_rb,1.164990e-02_rb, & |
---|
2177 | 1.125136e-02_rb,1.088901e-02_rb,1.055827e-02_rb,1.025531e-02_rb,9.976896e-03_rb, & |
---|
2178 | 9.720255e-03_rb,9.483022e-03_rb,9.263160e-03_rb,9.058902e-03_rb,8.868710e-03_rb, & |
---|
2179 | 8.691240e-03_rb,8.525312e-03_rb,8.369886e-03_rb,8.224042e-03_rb,8.086961e-03_rb, & |
---|
2180 | 7.957917e-03_rb,7.836258e-03_rb,7.721400e-03_rb,7.612821e-03_rb,7.510045e-03_rb, & |
---|
2181 | 7.412648e-03_rb,7.320242e-03_rb,7.232476e-03_rb/) |
---|
2182 | absice2(:,13) = (/ & |
---|
2183 | ! band 13 |
---|
2184 | 1.437021e-01_rb,8.872535e-02_rb,6.392420e-02_rb,4.991833e-02_rb,4.096790e-02_rb, & |
---|
2185 | 3.477881e-02_rb,3.025782e-02_rb,2.681909e-02_rb,2.412102e-02_rb,2.195132e-02_rb, & |
---|
2186 | 2.017124e-02_rb,1.868641e-02_rb,1.743044e-02_rb,1.635529e-02_rb,1.542540e-02_rb, & |
---|
2187 | 1.461388e-02_rb,1.390003e-02_rb,1.326766e-02_rb,1.270395e-02_rb,1.219860e-02_rb, & |
---|
2188 | 1.174326e-02_rb,1.133107e-02_rb,1.095637e-02_rb,1.061442e-02_rb,1.030126e-02_rb, & |
---|
2189 | 1.001352e-02_rb,9.748340e-03_rb,9.503256e-03_rb,9.276155e-03_rb,9.065205e-03_rb, & |
---|
2190 | 8.868808e-03_rb,8.685571e-03_rb,8.514268e-03_rb,8.353820e-03_rb,8.203272e-03_rb, & |
---|
2191 | 8.061776e-03_rb,7.928578e-03_rb,7.803001e-03_rb,7.684443e-03_rb,7.572358e-03_rb, & |
---|
2192 | 7.466258e-03_rb,7.365701e-03_rb,7.270286e-03_rb/) |
---|
2193 | absice2(:,14) = (/ & |
---|
2194 | ! band 14 |
---|
2195 | 1.288870e-01_rb,8.160295e-02_rb,5.964745e-02_rb,4.703790e-02_rb,3.888637e-02_rb, & |
---|
2196 | 3.320115e-02_rb,2.902017e-02_rb,2.582259e-02_rb,2.330224e-02_rb,2.126754e-02_rb, & |
---|
2197 | 1.959258e-02_rb,1.819130e-02_rb,1.700289e-02_rb,1.598320e-02_rb,1.509942e-02_rb, & |
---|
2198 | 1.432666e-02_rb,1.364572e-02_rb,1.304156e-02_rb,1.250220e-02_rb,1.201803e-02_rb, & |
---|
2199 | 1.158123e-02_rb,1.118537e-02_rb,1.082513e-02_rb,1.049605e-02_rb,1.019440e-02_rb, & |
---|
2200 | 9.916989e-03_rb,9.661116e-03_rb,9.424457e-03_rb,9.205005e-03_rb,9.001022e-03_rb, & |
---|
2201 | 8.810992e-03_rb,8.633588e-03_rb,8.467646e-03_rb,8.312137e-03_rb,8.166151e-03_rb, & |
---|
2202 | 8.028878e-03_rb,7.899597e-03_rb,7.777663e-03_rb,7.662498e-03_rb,7.553581e-03_rb, & |
---|
2203 | 7.450444e-03_rb,7.352662e-03_rb,7.259851e-03_rb/) |
---|
2204 | absice2(:,15) = (/ & |
---|
2205 | ! band 15 |
---|
2206 | 8.254229e-02_rb,5.808787e-02_rb,4.492166e-02_rb,3.675028e-02_rb,3.119623e-02_rb, & |
---|
2207 | 2.718045e-02_rb,2.414450e-02_rb,2.177073e-02_rb,1.986526e-02_rb,1.830306e-02_rb, & |
---|
2208 | 1.699991e-02_rb,1.589698e-02_rb,1.495199e-02_rb,1.413374e-02_rb,1.341870e-02_rb, & |
---|
2209 | 1.278883e-02_rb,1.223002e-02_rb,1.173114e-02_rb,1.128322e-02_rb,1.087900e-02_rb, & |
---|
2210 | 1.051254e-02_rb,1.017890e-02_rb,9.873991e-03_rb,9.594347e-03_rb,9.337044e-03_rb, & |
---|
2211 | 9.099589e-03_rb,8.879842e-03_rb,8.675960e-03_rb,8.486341e-03_rb,8.309594e-03_rb, & |
---|
2212 | 8.144500e-03_rb,7.989986e-03_rb,7.845109e-03_rb,7.709031e-03_rb,7.581007e-03_rb, & |
---|
2213 | 7.460376e-03_rb,7.346544e-03_rb,7.238978e-03_rb,7.137201e-03_rb,7.040780e-03_rb, & |
---|
2214 | 6.949325e-03_rb,6.862483e-03_rb,6.779931e-03_rb/) |
---|
2215 | absice2(:,16) = (/ & |
---|
2216 | ! band 16 |
---|
2217 | 1.382062e-01_rb,8.643227e-02_rb,6.282935e-02_rb,4.934783e-02_rb,4.063891e-02_rb, & |
---|
2218 | 3.455591e-02_rb,3.007059e-02_rb,2.662897e-02_rb,2.390631e-02_rb,2.169972e-02_rb, & |
---|
2219 | 1.987596e-02_rb,1.834393e-02_rb,1.703924e-02_rb,1.591513e-02_rb,1.493679e-02_rb, & |
---|
2220 | 1.407780e-02_rb,1.331775e-02_rb,1.264061e-02_rb,1.203364e-02_rb,1.148655e-02_rb, & |
---|
2221 | 1.099099e-02_rb,1.054006e-02_rb,1.012807e-02_rb,9.750215e-03_rb,9.402477e-03_rb, & |
---|
2222 | 9.081428e-03_rb,8.784143e-03_rb,8.508107e-03_rb,8.251146e-03_rb,8.011373e-03_rb, & |
---|
2223 | 7.787140e-03_rb,7.577002e-03_rb,7.379687e-03_rb,7.194071e-03_rb,7.019158e-03_rb, & |
---|
2224 | 6.854061e-03_rb,6.697986e-03_rb,6.550224e-03_rb,6.410138e-03_rb,6.277153e-03_rb, & |
---|
2225 | 6.150751e-03_rb,6.030462e-03_rb,5.915860e-03_rb/) |
---|
2226 | |
---|
2227 | ! ICEFLAG = 3; Fu parameterization. Particle size 5 - 140 micron in |
---|
2228 | ! increments of 3 microns. |
---|
2229 | ! units = m2/g |
---|
2230 | ! Hexagonal Ice Particle Parameterization |
---|
2231 | ! absorption units (abs coef/iwc): [(m^-1)/(g m^-3)] |
---|
2232 | absice3(:,1) = (/ & |
---|
2233 | ! band 1 |
---|
2234 | 3.110649e-03_rb,4.666352e-02_rb,6.606447e-02_rb,6.531678e-02_rb,6.012598e-02_rb, & |
---|
2235 | 5.437494e-02_rb,4.906411e-02_rb,4.441146e-02_rb,4.040585e-02_rb,3.697334e-02_rb, & |
---|
2236 | 3.403027e-02_rb,3.149979e-02_rb,2.931596e-02_rb,2.742365e-02_rb,2.577721e-02_rb, & |
---|
2237 | 2.433888e-02_rb,2.307732e-02_rb,2.196644e-02_rb,2.098437e-02_rb,2.011264e-02_rb, & |
---|
2238 | 1.933561e-02_rb,1.863992e-02_rb,1.801407e-02_rb,1.744812e-02_rb,1.693346e-02_rb, & |
---|
2239 | 1.646252e-02_rb,1.602866e-02_rb,1.562600e-02_rb,1.524933e-02_rb,1.489399e-02_rb, & |
---|
2240 | 1.455580e-02_rb,1.423098e-02_rb,1.391612e-02_rb,1.360812e-02_rb,1.330413e-02_rb, & |
---|
2241 | 1.300156e-02_rb,1.269801e-02_rb,1.239127e-02_rb,1.207928e-02_rb,1.176014e-02_rb, & |
---|
2242 | 1.143204e-02_rb,1.109334e-02_rb,1.074243e-02_rb,1.037786e-02_rb,9.998198e-03_rb, & |
---|
2243 | 9.602126e-03_rb/) |
---|
2244 | absice3(:,2) = (/ & |
---|
2245 | ! band 2 |
---|
2246 | 3.984966e-04_rb,1.681097e-02_rb,2.627680e-02_rb,2.767465e-02_rb,2.700722e-02_rb, & |
---|
2247 | 2.579180e-02_rb,2.448677e-02_rb,2.323890e-02_rb,2.209096e-02_rb,2.104882e-02_rb, & |
---|
2248 | 2.010547e-02_rb,1.925003e-02_rb,1.847128e-02_rb,1.775883e-02_rb,1.710358e-02_rb, & |
---|
2249 | 1.649769e-02_rb,1.593449e-02_rb,1.540829e-02_rb,1.491429e-02_rb,1.444837e-02_rb, & |
---|
2250 | 1.400704e-02_rb,1.358729e-02_rb,1.318654e-02_rb,1.280258e-02_rb,1.243346e-02_rb, & |
---|
2251 | 1.207750e-02_rb,1.173325e-02_rb,1.139941e-02_rb,1.107487e-02_rb,1.075861e-02_rb, & |
---|
2252 | 1.044975e-02_rb,1.014753e-02_rb,9.851229e-03_rb,9.560240e-03_rb,9.274003e-03_rb, & |
---|
2253 | 8.992020e-03_rb,8.713845e-03_rb,8.439074e-03_rb,8.167346e-03_rb,7.898331e-03_rb, & |
---|
2254 | 7.631734e-03_rb,7.367286e-03_rb,7.104742e-03_rb,6.843882e-03_rb,6.584504e-03_rb, & |
---|
2255 | 6.326424e-03_rb/) |
---|
2256 | absice3(:,3) = (/ & |
---|
2257 | ! band 3 |
---|
2258 | 6.933163e-02_rb,8.540475e-02_rb,7.701816e-02_rb,6.771158e-02_rb,5.986953e-02_rb, & |
---|
2259 | 5.348120e-02_rb,4.824962e-02_rb,4.390563e-02_rb,4.024411e-02_rb,3.711404e-02_rb, & |
---|
2260 | 3.440426e-02_rb,3.203200e-02_rb,2.993478e-02_rb,2.806474e-02_rb,2.638464e-02_rb, & |
---|
2261 | 2.486516e-02_rb,2.348288e-02_rb,2.221890e-02_rb,2.105780e-02_rb,1.998687e-02_rb, & |
---|
2262 | 1.899552e-02_rb,1.807490e-02_rb,1.721750e-02_rb,1.641693e-02_rb,1.566773e-02_rb, & |
---|
2263 | 1.496515e-02_rb,1.430509e-02_rb,1.368398e-02_rb,1.309865e-02_rb,1.254634e-02_rb, & |
---|
2264 | 1.202456e-02_rb,1.153114e-02_rb,1.106409e-02_rb,1.062166e-02_rb,1.020224e-02_rb, & |
---|
2265 | 9.804381e-03_rb,9.426771e-03_rb,9.068205e-03_rb,8.727578e-03_rb,8.403876e-03_rb, & |
---|
2266 | 8.096160e-03_rb,7.803564e-03_rb,7.525281e-03_rb,7.260560e-03_rb,7.008697e-03_rb, & |
---|
2267 | 6.769036e-03_rb/) |
---|
2268 | absice3(:,4) = (/ & |
---|
2269 | ! band 4 |
---|
2270 | 1.765735e-01_rb,1.382700e-01_rb,1.095129e-01_rb,8.987475e-02_rb,7.591185e-02_rb, & |
---|
2271 | 6.554169e-02_rb,5.755500e-02_rb,5.122083e-02_rb,4.607610e-02_rb,4.181475e-02_rb, & |
---|
2272 | 3.822697e-02_rb,3.516432e-02_rb,3.251897e-02_rb,3.021073e-02_rb,2.817876e-02_rb, & |
---|
2273 | 2.637607e-02_rb,2.476582e-02_rb,2.331871e-02_rb,2.201113e-02_rb,2.082388e-02_rb, & |
---|
2274 | 1.974115e-02_rb,1.874983e-02_rb,1.783894e-02_rb,1.699922e-02_rb,1.622280e-02_rb, & |
---|
2275 | 1.550296e-02_rb,1.483390e-02_rb,1.421064e-02_rb,1.362880e-02_rb,1.308460e-02_rb, & |
---|
2276 | 1.257468e-02_rb,1.209611e-02_rb,1.164628e-02_rb,1.122287e-02_rb,1.082381e-02_rb, & |
---|
2277 | 1.044725e-02_rb,1.009154e-02_rb,9.755166e-03_rb,9.436783e-03_rb,9.135163e-03_rb, & |
---|
2278 | 8.849193e-03_rb,8.577856e-03_rb,8.320225e-03_rb,8.075451e-03_rb,7.842755e-03_rb, & |
---|
2279 | 7.621418e-03_rb/) |
---|
2280 | absice3(:,5) = (/ & |
---|
2281 | ! band 5 |
---|
2282 | 2.339673e-01_rb,1.692124e-01_rb,1.291656e-01_rb,1.033837e-01_rb,8.562949e-02_rb, & |
---|
2283 | 7.273526e-02_rb,6.298262e-02_rb,5.537015e-02_rb,4.927787e-02_rb,4.430246e-02_rb, & |
---|
2284 | 4.017061e-02_rb,3.669072e-02_rb,3.372455e-02_rb,3.116995e-02_rb,2.894977e-02_rb, & |
---|
2285 | 2.700471e-02_rb,2.528842e-02_rb,2.376420e-02_rb,2.240256e-02_rb,2.117959e-02_rb, & |
---|
2286 | 2.007567e-02_rb,1.907456e-02_rb,1.816271e-02_rb,1.732874e-02_rb,1.656300e-02_rb, & |
---|
2287 | 1.585725e-02_rb,1.520445e-02_rb,1.459852e-02_rb,1.403419e-02_rb,1.350689e-02_rb, & |
---|
2288 | 1.301260e-02_rb,1.254781e-02_rb,1.210941e-02_rb,1.169468e-02_rb,1.130118e-02_rb, & |
---|
2289 | 1.092675e-02_rb,1.056945e-02_rb,1.022757e-02_rb,9.899560e-03_rb,9.584021e-03_rb, & |
---|
2290 | 9.279705e-03_rb,8.985479e-03_rb,8.700322e-03_rb,8.423306e-03_rb,8.153590e-03_rb, & |
---|
2291 | 7.890412e-03_rb/) |
---|
2292 | absice3(:,6) = (/ & |
---|
2293 | ! band 6 |
---|
2294 | 1.145369e-01_rb,1.174566e-01_rb,9.917866e-02_rb,8.332990e-02_rb,7.104263e-02_rb, & |
---|
2295 | 6.153370e-02_rb,5.405472e-02_rb,4.806281e-02_rb,4.317918e-02_rb,3.913795e-02_rb, & |
---|
2296 | 3.574916e-02_rb,3.287437e-02_rb,3.041067e-02_rb,2.828017e-02_rb,2.642292e-02_rb, & |
---|
2297 | 2.479206e-02_rb,2.335051e-02_rb,2.206851e-02_rb,2.092195e-02_rb,1.989108e-02_rb, & |
---|
2298 | 1.895958e-02_rb,1.811385e-02_rb,1.734245e-02_rb,1.663573e-02_rb,1.598545e-02_rb, & |
---|
2299 | 1.538456e-02_rb,1.482700e-02_rb,1.430750e-02_rb,1.382150e-02_rb,1.336499e-02_rb, & |
---|
2300 | 1.293447e-02_rb,1.252685e-02_rb,1.213939e-02_rb,1.176968e-02_rb,1.141555e-02_rb, & |
---|
2301 | 1.107508e-02_rb,1.074655e-02_rb,1.042839e-02_rb,1.011923e-02_rb,9.817799e-03_rb, & |
---|
2302 | 9.522962e-03_rb,9.233688e-03_rb,8.949041e-03_rb,8.668171e-03_rb,8.390301e-03_rb, & |
---|
2303 | 8.114723e-03_rb/) |
---|
2304 | absice3(:,7) = (/ & |
---|
2305 | ! band 7 |
---|
2306 | 1.222345e-02_rb,5.344230e-02_rb,5.523465e-02_rb,5.128759e-02_rb,4.676925e-02_rb, & |
---|
2307 | 4.266150e-02_rb,3.910561e-02_rb,3.605479e-02_rb,3.342843e-02_rb,3.115052e-02_rb, & |
---|
2308 | 2.915776e-02_rb,2.739935e-02_rb,2.583499e-02_rb,2.443266e-02_rb,2.316681e-02_rb, & |
---|
2309 | 2.201687e-02_rb,2.096619e-02_rb,2.000112e-02_rb,1.911044e-02_rb,1.828481e-02_rb, & |
---|
2310 | 1.751641e-02_rb,1.679866e-02_rb,1.612598e-02_rb,1.549360e-02_rb,1.489742e-02_rb, & |
---|
2311 | 1.433392e-02_rb,1.380002e-02_rb,1.329305e-02_rb,1.281068e-02_rb,1.235084e-02_rb, & |
---|
2312 | 1.191172e-02_rb,1.149171e-02_rb,1.108936e-02_rb,1.070341e-02_rb,1.033271e-02_rb, & |
---|
2313 | 9.976220e-03_rb,9.633021e-03_rb,9.302273e-03_rb,8.983216e-03_rb,8.675161e-03_rb, & |
---|
2314 | 8.377478e-03_rb,8.089595e-03_rb,7.810986e-03_rb,7.541170e-03_rb,7.279706e-03_rb, & |
---|
2315 | 7.026186e-03_rb/) |
---|
2316 | absice3(:,8) = (/ & |
---|
2317 | ! band 8 |
---|
2318 | 6.711058e-02_rb,6.918198e-02_rb,6.127484e-02_rb,5.411944e-02_rb,4.836902e-02_rb, & |
---|
2319 | 4.375293e-02_rb,3.998077e-02_rb,3.683587e-02_rb,3.416508e-02_rb,3.186003e-02_rb, & |
---|
2320 | 2.984290e-02_rb,2.805671e-02_rb,2.645895e-02_rb,2.501733e-02_rb,2.370689e-02_rb, & |
---|
2321 | 2.250808e-02_rb,2.140532e-02_rb,2.038609e-02_rb,1.944018e-02_rb,1.855918e-02_rb, & |
---|
2322 | 1.773609e-02_rb,1.696504e-02_rb,1.624106e-02_rb,1.555990e-02_rb,1.491793e-02_rb, & |
---|
2323 | 1.431197e-02_rb,1.373928e-02_rb,1.319743e-02_rb,1.268430e-02_rb,1.219799e-02_rb, & |
---|
2324 | 1.173682e-02_rb,1.129925e-02_rb,1.088393e-02_rb,1.048961e-02_rb,1.011516e-02_rb, & |
---|
2325 | 9.759543e-03_rb,9.421813e-03_rb,9.101089e-03_rb,8.796559e-03_rb,8.507464e-03_rb, & |
---|
2326 | 8.233098e-03_rb,7.972798e-03_rb,7.725942e-03_rb,7.491940e-03_rb,7.270238e-03_rb, & |
---|
2327 | 7.060305e-03_rb/) |
---|
2328 | absice3(:,9) = (/ & |
---|
2329 | ! band 9 |
---|
2330 | 1.236780e-01_rb,9.222386e-02_rb,7.383997e-02_rb,6.204072e-02_rb,5.381029e-02_rb, & |
---|
2331 | 4.770678e-02_rb,4.296928e-02_rb,3.916131e-02_rb,3.601540e-02_rb,3.335878e-02_rb, & |
---|
2332 | 3.107493e-02_rb,2.908247e-02_rb,2.732282e-02_rb,2.575276e-02_rb,2.433968e-02_rb, & |
---|
2333 | 2.305852e-02_rb,2.188966e-02_rb,2.081757e-02_rb,1.982974e-02_rb,1.891599e-02_rb, & |
---|
2334 | 1.806794e-02_rb,1.727865e-02_rb,1.654227e-02_rb,1.585387e-02_rb,1.520924e-02_rb, & |
---|
2335 | 1.460476e-02_rb,1.403730e-02_rb,1.350416e-02_rb,1.300293e-02_rb,1.253153e-02_rb, & |
---|
2336 | 1.208808e-02_rb,1.167094e-02_rb,1.127862e-02_rb,1.090979e-02_rb,1.056323e-02_rb, & |
---|
2337 | 1.023786e-02_rb,9.932665e-03_rb,9.646744e-03_rb,9.379250e-03_rb,9.129409e-03_rb, & |
---|
2338 | 8.896500e-03_rb,8.679856e-03_rb,8.478852e-03_rb,8.292904e-03_rb,8.121463e-03_rb, & |
---|
2339 | 7.964013e-03_rb/) |
---|
2340 | absice3(:,10) = (/ & |
---|
2341 | ! band 10 |
---|
2342 | 1.655966e-01_rb,1.134205e-01_rb,8.714344e-02_rb,7.129241e-02_rb,6.063739e-02_rb, & |
---|
2343 | 5.294203e-02_rb,4.709309e-02_rb,4.247476e-02_rb,3.871892e-02_rb,3.559206e-02_rb, & |
---|
2344 | 3.293893e-02_rb,3.065226e-02_rb,2.865558e-02_rb,2.689288e-02_rb,2.532221e-02_rb, & |
---|
2345 | 2.391150e-02_rb,2.263582e-02_rb,2.147549e-02_rb,2.041476e-02_rb,1.944089e-02_rb, & |
---|
2346 | 1.854342e-02_rb,1.771371e-02_rb,1.694456e-02_rb,1.622989e-02_rb,1.556456e-02_rb, & |
---|
2347 | 1.494415e-02_rb,1.436491e-02_rb,1.382354e-02_rb,1.331719e-02_rb,1.284339e-02_rb, & |
---|
2348 | 1.239992e-02_rb,1.198486e-02_rb,1.159647e-02_rb,1.123323e-02_rb,1.089375e-02_rb, & |
---|
2349 | 1.057679e-02_rb,1.028124e-02_rb,1.000607e-02_rb,9.750376e-03_rb,9.513303e-03_rb, & |
---|
2350 | 9.294082e-03_rb,9.092003e-03_rb,8.906412e-03_rb,8.736702e-03_rb,8.582314e-03_rb, & |
---|
2351 | 8.442725e-03_rb/) |
---|
2352 | absice3(:,11) = (/ & |
---|
2353 | ! band 11 |
---|
2354 | 1.775615e-01_rb,1.180046e-01_rb,8.929607e-02_rb,7.233500e-02_rb,6.108333e-02_rb, & |
---|
2355 | 5.303642e-02_rb,4.696927e-02_rb,4.221206e-02_rb,3.836768e-02_rb,3.518576e-02_rb, & |
---|
2356 | 3.250063e-02_rb,3.019825e-02_rb,2.819758e-02_rb,2.643943e-02_rb,2.487953e-02_rb, & |
---|
2357 | 2.348414e-02_rb,2.222705e-02_rb,2.108762e-02_rb,2.004936e-02_rb,1.909892e-02_rb, & |
---|
2358 | 1.822539e-02_rb,1.741975e-02_rb,1.667449e-02_rb,1.598330e-02_rb,1.534084e-02_rb, & |
---|
2359 | 1.474253e-02_rb,1.418446e-02_rb,1.366325e-02_rb,1.317597e-02_rb,1.272004e-02_rb, & |
---|
2360 | 1.229321e-02_rb,1.189350e-02_rb,1.151915e-02_rb,1.116859e-02_rb,1.084042e-02_rb, & |
---|
2361 | 1.053338e-02_rb,1.024636e-02_rb,9.978326e-03_rb,9.728357e-03_rb,9.495613e-03_rb, & |
---|
2362 | 9.279327e-03_rb,9.078798e-03_rb,8.893383e-03_rb,8.722488e-03_rb,8.565568e-03_rb, & |
---|
2363 | 8.422115e-03_rb/) |
---|
2364 | absice3(:,12) = (/ & |
---|
2365 | ! band 12 |
---|
2366 | 9.465447e-02_rb,6.432047e-02_rb,5.060973e-02_rb,4.267283e-02_rb,3.741843e-02_rb, & |
---|
2367 | 3.363096e-02_rb,3.073531e-02_rb,2.842405e-02_rb,2.651789e-02_rb,2.490518e-02_rb, & |
---|
2368 | 2.351273e-02_rb,2.229056e-02_rb,2.120335e-02_rb,2.022541e-02_rb,1.933763e-02_rb, & |
---|
2369 | 1.852546e-02_rb,1.777763e-02_rb,1.708528e-02_rb,1.644134e-02_rb,1.584009e-02_rb, & |
---|
2370 | 1.527684e-02_rb,1.474774e-02_rb,1.424955e-02_rb,1.377957e-02_rb,1.333549e-02_rb, & |
---|
2371 | 1.291534e-02_rb,1.251743e-02_rb,1.214029e-02_rb,1.178265e-02_rb,1.144337e-02_rb, & |
---|
2372 | 1.112148e-02_rb,1.081609e-02_rb,1.052642e-02_rb,1.025178e-02_rb,9.991540e-03_rb, & |
---|
2373 | 9.745130e-03_rb,9.512038e-03_rb,9.291797e-03_rb,9.083980e-03_rb,8.888195e-03_rb, & |
---|
2374 | 8.704081e-03_rb,8.531306e-03_rb,8.369560e-03_rb,8.218558e-03_rb,8.078032e-03_rb, & |
---|
2375 | 7.947730e-03_rb/) |
---|
2376 | absice3(:,13) = (/ & |
---|
2377 | ! band 13 |
---|
2378 | 1.560311e-01_rb,9.961097e-02_rb,7.502949e-02_rb,6.115022e-02_rb,5.214952e-02_rb, & |
---|
2379 | 4.578149e-02_rb,4.099731e-02_rb,3.724174e-02_rb,3.419343e-02_rb,3.165356e-02_rb, & |
---|
2380 | 2.949251e-02_rb,2.762222e-02_rb,2.598073e-02_rb,2.452322e-02_rb,2.321642e-02_rb, & |
---|
2381 | 2.203516e-02_rb,2.096002e-02_rb,1.997579e-02_rb,1.907036e-02_rb,1.823401e-02_rb, & |
---|
2382 | 1.745879e-02_rb,1.673819e-02_rb,1.606678e-02_rb,1.544003e-02_rb,1.485411e-02_rb, & |
---|
2383 | 1.430574e-02_rb,1.379215e-02_rb,1.331092e-02_rb,1.285996e-02_rb,1.243746e-02_rb, & |
---|
2384 | 1.204183e-02_rb,1.167164e-02_rb,1.132567e-02_rb,1.100281e-02_rb,1.070207e-02_rb, & |
---|
2385 | 1.042258e-02_rb,1.016352e-02_rb,9.924197e-03_rb,9.703953e-03_rb,9.502199e-03_rb, & |
---|
2386 | 9.318400e-03_rb,9.152066e-03_rb,9.002749e-03_rb,8.870038e-03_rb,8.753555e-03_rb, & |
---|
2387 | 8.652951e-03_rb/) |
---|
2388 | absice3(:,14) = (/ & |
---|
2389 | ! band 14 |
---|
2390 | 1.559547e-01_rb,9.896700e-02_rb,7.441231e-02_rb,6.061469e-02_rb,5.168730e-02_rb, & |
---|
2391 | 4.537821e-02_rb,4.064106e-02_rb,3.692367e-02_rb,3.390714e-02_rb,3.139438e-02_rb, & |
---|
2392 | 2.925702e-02_rb,2.740783e-02_rb,2.578547e-02_rb,2.434552e-02_rb,2.305506e-02_rb, & |
---|
2393 | 2.188910e-02_rb,2.082842e-02_rb,1.985789e-02_rb,1.896553e-02_rb,1.814165e-02_rb, & |
---|
2394 | 1.737839e-02_rb,1.666927e-02_rb,1.600891e-02_rb,1.539279e-02_rb,1.481712e-02_rb, & |
---|
2395 | 1.427865e-02_rb,1.377463e-02_rb,1.330266e-02_rb,1.286068e-02_rb,1.244689e-02_rb, & |
---|
2396 | 1.205973e-02_rb,1.169780e-02_rb,1.135989e-02_rb,1.104492e-02_rb,1.075192e-02_rb, & |
---|
2397 | 1.048004e-02_rb,1.022850e-02_rb,9.996611e-03_rb,9.783753e-03_rb,9.589361e-03_rb, & |
---|
2398 | 9.412924e-03_rb,9.253977e-03_rb,9.112098e-03_rb,8.986903e-03_rb,8.878039e-03_rb, & |
---|
2399 | 8.785184e-03_rb/) |
---|
2400 | absice3(:,15) = (/ & |
---|
2401 | ! band 15 |
---|
2402 | 1.102926e-01_rb,7.176622e-02_rb,5.530316e-02_rb,4.606056e-02_rb,4.006116e-02_rb, & |
---|
2403 | 3.579628e-02_rb,3.256909e-02_rb,3.001360e-02_rb,2.791920e-02_rb,2.615617e-02_rb, & |
---|
2404 | 2.464023e-02_rb,2.331426e-02_rb,2.213817e-02_rb,2.108301e-02_rb,2.012733e-02_rb, & |
---|
2405 | 1.925493e-02_rb,1.845331e-02_rb,1.771269e-02_rb,1.702531e-02_rb,1.638493e-02_rb, & |
---|
2406 | 1.578648e-02_rb,1.522579e-02_rb,1.469940e-02_rb,1.420442e-02_rb,1.373841e-02_rb, & |
---|
2407 | 1.329931e-02_rb,1.288535e-02_rb,1.249502e-02_rb,1.212700e-02_rb,1.178015e-02_rb, & |
---|
2408 | 1.145348e-02_rb,1.114612e-02_rb,1.085730e-02_rb,1.058633e-02_rb,1.033263e-02_rb, & |
---|
2409 | 1.009564e-02_rb,9.874895e-03_rb,9.669960e-03_rb,9.480449e-03_rb,9.306014e-03_rb, & |
---|
2410 | 9.146339e-03_rb,9.001138e-03_rb,8.870154e-03_rb,8.753148e-03_rb,8.649907e-03_rb, & |
---|
2411 | 8.560232e-03_rb/) |
---|
2412 | absice3(:,16) = (/ & |
---|
2413 | ! band 16 |
---|
2414 | 1.688344e-01_rb,1.077072e-01_rb,7.994467e-02_rb,6.403862e-02_rb,5.369850e-02_rb, & |
---|
2415 | 4.641582e-02_rb,4.099331e-02_rb,3.678724e-02_rb,3.342069e-02_rb,3.065831e-02_rb, & |
---|
2416 | 2.834557e-02_rb,2.637680e-02_rb,2.467733e-02_rb,2.319286e-02_rb,2.188299e-02_rb, & |
---|
2417 | 2.071701e-02_rb,1.967121e-02_rb,1.872692e-02_rb,1.786931e-02_rb,1.708641e-02_rb, & |
---|
2418 | 1.636846e-02_rb,1.570743e-02_rb,1.509665e-02_rb,1.453052e-02_rb,1.400433e-02_rb, & |
---|
2419 | 1.351407e-02_rb,1.305631e-02_rb,1.262810e-02_rb,1.222688e-02_rb,1.185044e-02_rb, & |
---|
2420 | 1.149683e-02_rb,1.116436e-02_rb,1.085153e-02_rb,1.055701e-02_rb,1.027961e-02_rb, & |
---|
2421 | 1.001831e-02_rb,9.772141e-03_rb,9.540280e-03_rb,9.321966e-03_rb,9.116517e-03_rb, & |
---|
2422 | 8.923315e-03_rb,8.741803e-03_rb,8.571472e-03_rb,8.411860e-03_rb,8.262543e-03_rb, & |
---|
2423 | 8.123136e-03_rb/) |
---|
2424 | |
---|
2425 | ! For LIQFLAG = 0. |
---|
2426 | absliq0 = 0.0903614_rb |
---|
2427 | |
---|
2428 | ! For LIQFLAG = 1. In each band, the absorption |
---|
2429 | ! coefficients are listed for a range of effective radii from 2.5 |
---|
2430 | ! to 59.5 microns in increments of 1.0 micron. |
---|
2431 | absliq1(:, 1) = (/ & |
---|
2432 | ! band 1 |
---|
2433 | 1.64047e-03_rb, 6.90533e-02_rb, 7.72017e-02_rb, 7.78054e-02_rb, 7.69523e-02_rb, & |
---|
2434 | 7.58058e-02_rb, 7.46400e-02_rb, 7.35123e-02_rb, 7.24162e-02_rb, 7.13225e-02_rb, & |
---|
2435 | 6.99145e-02_rb, 6.66409e-02_rb, 6.36582e-02_rb, 6.09425e-02_rb, 5.84593e-02_rb, & |
---|
2436 | 5.61743e-02_rb, 5.40571e-02_rb, 5.20812e-02_rb, 5.02245e-02_rb, 4.84680e-02_rb, & |
---|
2437 | 4.67959e-02_rb, 4.51944e-02_rb, 4.36516e-02_rb, 4.21570e-02_rb, 4.07015e-02_rb, & |
---|
2438 | 3.92766e-02_rb, 3.78747e-02_rb, 3.64886e-02_rb, 3.53632e-02_rb, 3.41992e-02_rb, & |
---|
2439 | 3.31016e-02_rb, 3.20643e-02_rb, 3.10817e-02_rb, 3.01490e-02_rb, 2.92620e-02_rb, & |
---|
2440 | 2.84171e-02_rb, 2.76108e-02_rb, 2.68404e-02_rb, 2.61031e-02_rb, 2.53966e-02_rb, & |
---|
2441 | 2.47189e-02_rb, 2.40678e-02_rb, 2.34418e-02_rb, 2.28392e-02_rb, 2.22586e-02_rb, & |
---|
2442 | 2.16986e-02_rb, 2.11580e-02_rb, 2.06356e-02_rb, 2.01305e-02_rb, 1.96417e-02_rb, & |
---|
2443 | 1.91682e-02_rb, 1.87094e-02_rb, 1.82643e-02_rb, 1.78324e-02_rb, 1.74129e-02_rb, & |
---|
2444 | 1.70052e-02_rb, 1.66088e-02_rb, 1.62231e-02_rb/) |
---|
2445 | absliq1(:, 2) = (/ & |
---|
2446 | ! band 2 |
---|
2447 | 2.19486e-01_rb, 1.80687e-01_rb, 1.59150e-01_rb, 1.44731e-01_rb, 1.33703e-01_rb, & |
---|
2448 | 1.24355e-01_rb, 1.15756e-01_rb, 1.07318e-01_rb, 9.86119e-02_rb, 8.92739e-02_rb, & |
---|
2449 | 8.34911e-02_rb, 7.70773e-02_rb, 7.15240e-02_rb, 6.66615e-02_rb, 6.23641e-02_rb, & |
---|
2450 | 5.85359e-02_rb, 5.51020e-02_rb, 5.20032e-02_rb, 4.91916e-02_rb, 4.66283e-02_rb, & |
---|
2451 | 4.42813e-02_rb, 4.21236e-02_rb, 4.01330e-02_rb, 3.82905e-02_rb, 3.65797e-02_rb, & |
---|
2452 | 3.49869e-02_rb, 3.35002e-02_rb, 3.21090e-02_rb, 3.08957e-02_rb, 2.97601e-02_rb, & |
---|
2453 | 2.86966e-02_rb, 2.76984e-02_rb, 2.67599e-02_rb, 2.58758e-02_rb, 2.50416e-02_rb, & |
---|
2454 | 2.42532e-02_rb, 2.35070e-02_rb, 2.27997e-02_rb, 2.21284e-02_rb, 2.14904e-02_rb, & |
---|
2455 | 2.08834e-02_rb, 2.03051e-02_rb, 1.97536e-02_rb, 1.92271e-02_rb, 1.87239e-02_rb, & |
---|
2456 | 1.82425e-02_rb, 1.77816e-02_rb, 1.73399e-02_rb, 1.69162e-02_rb, 1.65094e-02_rb, & |
---|
2457 | 1.61187e-02_rb, 1.57430e-02_rb, 1.53815e-02_rb, 1.50334e-02_rb, 1.46981e-02_rb, & |
---|
2458 | 1.43748e-02_rb, 1.40628e-02_rb, 1.37617e-02_rb/) |
---|
2459 | absliq1(:, 3) = (/ & |
---|
2460 | ! band 3 |
---|
2461 | 2.95174e-01_rb, 2.34765e-01_rb, 1.98038e-01_rb, 1.72114e-01_rb, 1.52083e-01_rb, & |
---|
2462 | 1.35654e-01_rb, 1.21613e-01_rb, 1.09252e-01_rb, 9.81263e-02_rb, 8.79448e-02_rb, & |
---|
2463 | 8.12566e-02_rb, 7.44563e-02_rb, 6.86374e-02_rb, 6.36042e-02_rb, 5.92094e-02_rb, & |
---|
2464 | 5.53402e-02_rb, 5.19087e-02_rb, 4.88455e-02_rb, 4.60951e-02_rb, 4.36124e-02_rb, & |
---|
2465 | 4.13607e-02_rb, 3.93096e-02_rb, 3.74338e-02_rb, 3.57119e-02_rb, 3.41261e-02_rb, & |
---|
2466 | 3.26610e-02_rb, 3.13036e-02_rb, 3.00425e-02_rb, 2.88497e-02_rb, 2.78077e-02_rb, & |
---|
2467 | 2.68317e-02_rb, 2.59158e-02_rb, 2.50545e-02_rb, 2.42430e-02_rb, 2.34772e-02_rb, & |
---|
2468 | 2.27533e-02_rb, 2.20679e-02_rb, 2.14181e-02_rb, 2.08011e-02_rb, 2.02145e-02_rb, & |
---|
2469 | 1.96561e-02_rb, 1.91239e-02_rb, 1.86161e-02_rb, 1.81311e-02_rb, 1.76673e-02_rb, & |
---|
2470 | 1.72234e-02_rb, 1.67981e-02_rb, 1.63903e-02_rb, 1.59989e-02_rb, 1.56230e-02_rb, & |
---|
2471 | 1.52615e-02_rb, 1.49138e-02_rb, 1.45791e-02_rb, 1.42565e-02_rb, 1.39455e-02_rb, & |
---|
2472 | 1.36455e-02_rb, 1.33559e-02_rb, 1.30761e-02_rb/) |
---|
2473 | absliq1(:, 4) = (/ & |
---|
2474 | ! band 4 |
---|
2475 | 3.00925e-01_rb, 2.36949e-01_rb, 1.96947e-01_rb, 1.68692e-01_rb, 1.47190e-01_rb, & |
---|
2476 | 1.29986e-01_rb, 1.15719e-01_rb, 1.03568e-01_rb, 9.30028e-02_rb, 8.36658e-02_rb, & |
---|
2477 | 7.71075e-02_rb, 7.07002e-02_rb, 6.52284e-02_rb, 6.05024e-02_rb, 5.63801e-02_rb, & |
---|
2478 | 5.27534e-02_rb, 4.95384e-02_rb, 4.66690e-02_rb, 4.40925e-02_rb, 4.17664e-02_rb, & |
---|
2479 | 3.96559e-02_rb, 3.77326e-02_rb, 3.59727e-02_rb, 3.43561e-02_rb, 3.28662e-02_rb, & |
---|
2480 | 3.14885e-02_rb, 3.02110e-02_rb, 2.90231e-02_rb, 2.78948e-02_rb, 2.69109e-02_rb, & |
---|
2481 | 2.59884e-02_rb, 2.51217e-02_rb, 2.43058e-02_rb, 2.35364e-02_rb, 2.28096e-02_rb, & |
---|
2482 | 2.21218e-02_rb, 2.14700e-02_rb, 2.08515e-02_rb, 2.02636e-02_rb, 1.97041e-02_rb, & |
---|
2483 | 1.91711e-02_rb, 1.86625e-02_rb, 1.81769e-02_rb, 1.77126e-02_rb, 1.72683e-02_rb, & |
---|
2484 | 1.68426e-02_rb, 1.64344e-02_rb, 1.60427e-02_rb, 1.56664e-02_rb, 1.53046e-02_rb, & |
---|
2485 | 1.49565e-02_rb, 1.46214e-02_rb, 1.42985e-02_rb, 1.39871e-02_rb, 1.36866e-02_rb, & |
---|
2486 | 1.33965e-02_rb, 1.31162e-02_rb, 1.28453e-02_rb/) |
---|
2487 | absliq1(:, 5) = (/ & |
---|
2488 | ! band 5 |
---|
2489 | 2.64691e-01_rb, 2.12018e-01_rb, 1.78009e-01_rb, 1.53539e-01_rb, 1.34721e-01_rb, & |
---|
2490 | 1.19580e-01_rb, 1.06996e-01_rb, 9.62772e-02_rb, 8.69710e-02_rb, 7.87670e-02_rb, & |
---|
2491 | 7.29272e-02_rb, 6.70920e-02_rb, 6.20977e-02_rb, 5.77732e-02_rb, 5.39910e-02_rb, & |
---|
2492 | 5.06538e-02_rb, 4.76866e-02_rb, 4.50301e-02_rb, 4.26374e-02_rb, 4.04704e-02_rb, & |
---|
2493 | 3.84981e-02_rb, 3.66948e-02_rb, 3.50394e-02_rb, 3.35141e-02_rb, 3.21038e-02_rb, & |
---|
2494 | 3.07957e-02_rb, 2.95788e-02_rb, 2.84438e-02_rb, 2.73790e-02_rb, 2.64390e-02_rb, & |
---|
2495 | 2.55565e-02_rb, 2.47263e-02_rb, 2.39437e-02_rb, 2.32047e-02_rb, 2.25056e-02_rb, & |
---|
2496 | 2.18433e-02_rb, 2.12149e-02_rb, 2.06177e-02_rb, 2.00495e-02_rb, 1.95081e-02_rb, & |
---|
2497 | 1.89917e-02_rb, 1.84984e-02_rb, 1.80269e-02_rb, 1.75755e-02_rb, 1.71431e-02_rb, & |
---|
2498 | 1.67283e-02_rb, 1.63303e-02_rb, 1.59478e-02_rb, 1.55801e-02_rb, 1.52262e-02_rb, & |
---|
2499 | 1.48853e-02_rb, 1.45568e-02_rb, 1.42400e-02_rb, 1.39342e-02_rb, 1.36388e-02_rb, & |
---|
2500 | 1.33533e-02_rb, 1.30773e-02_rb, 1.28102e-02_rb/) |
---|
2501 | absliq1(:, 6) = (/ & |
---|
2502 | ! band 6 |
---|
2503 | 8.81182e-02_rb, 1.06745e-01_rb, 9.79753e-02_rb, 8.99625e-02_rb, 8.35200e-02_rb, & |
---|
2504 | 7.81899e-02_rb, 7.35939e-02_rb, 6.94696e-02_rb, 6.56266e-02_rb, 6.19148e-02_rb, & |
---|
2505 | 5.83355e-02_rb, 5.49306e-02_rb, 5.19642e-02_rb, 4.93325e-02_rb, 4.69659e-02_rb, & |
---|
2506 | 4.48148e-02_rb, 4.28431e-02_rb, 4.10231e-02_rb, 3.93332e-02_rb, 3.77563e-02_rb, & |
---|
2507 | 3.62785e-02_rb, 3.48882e-02_rb, 3.35758e-02_rb, 3.23333e-02_rb, 3.11536e-02_rb, & |
---|
2508 | 3.00310e-02_rb, 2.89601e-02_rb, 2.79365e-02_rb, 2.70502e-02_rb, 2.62618e-02_rb, & |
---|
2509 | 2.55025e-02_rb, 2.47728e-02_rb, 2.40726e-02_rb, 2.34013e-02_rb, 2.27583e-02_rb, & |
---|
2510 | 2.21422e-02_rb, 2.15522e-02_rb, 2.09869e-02_rb, 2.04453e-02_rb, 1.99260e-02_rb, & |
---|
2511 | 1.94280e-02_rb, 1.89501e-02_rb, 1.84913e-02_rb, 1.80506e-02_rb, 1.76270e-02_rb, & |
---|
2512 | 1.72196e-02_rb, 1.68276e-02_rb, 1.64500e-02_rb, 1.60863e-02_rb, 1.57357e-02_rb, & |
---|
2513 | 1.53975e-02_rb, 1.50710e-02_rb, 1.47558e-02_rb, 1.44511e-02_rb, 1.41566e-02_rb, & |
---|
2514 | 1.38717e-02_rb, 1.35960e-02_rb, 1.33290e-02_rb/) |
---|
2515 | absliq1(:, 7) = (/ & |
---|
2516 | ! band 7 |
---|
2517 | 4.32174e-02_rb, 7.36078e-02_rb, 6.98340e-02_rb, 6.65231e-02_rb, 6.41948e-02_rb, & |
---|
2518 | 6.23551e-02_rb, 6.06638e-02_rb, 5.88680e-02_rb, 5.67124e-02_rb, 5.38629e-02_rb, & |
---|
2519 | 4.99579e-02_rb, 4.86289e-02_rb, 4.70120e-02_rb, 4.52854e-02_rb, 4.35466e-02_rb, & |
---|
2520 | 4.18480e-02_rb, 4.02169e-02_rb, 3.86658e-02_rb, 3.71992e-02_rb, 3.58168e-02_rb, & |
---|
2521 | 3.45155e-02_rb, 3.32912e-02_rb, 3.21390e-02_rb, 3.10538e-02_rb, 3.00307e-02_rb, & |
---|
2522 | 2.90651e-02_rb, 2.81524e-02_rb, 2.72885e-02_rb, 2.62821e-02_rb, 2.55744e-02_rb, & |
---|
2523 | 2.48799e-02_rb, 2.42029e-02_rb, 2.35460e-02_rb, 2.29108e-02_rb, 2.22981e-02_rb, & |
---|
2524 | 2.17079e-02_rb, 2.11402e-02_rb, 2.05945e-02_rb, 2.00701e-02_rb, 1.95663e-02_rb, & |
---|
2525 | 1.90824e-02_rb, 1.86174e-02_rb, 1.81706e-02_rb, 1.77411e-02_rb, 1.73281e-02_rb, & |
---|
2526 | 1.69307e-02_rb, 1.65483e-02_rb, 1.61801e-02_rb, 1.58254e-02_rb, 1.54835e-02_rb, & |
---|
2527 | 1.51538e-02_rb, 1.48358e-02_rb, 1.45288e-02_rb, 1.42322e-02_rb, 1.39457e-02_rb, & |
---|
2528 | 1.36687e-02_rb, 1.34008e-02_rb, 1.31416e-02_rb/) |
---|
2529 | absliq1(:, 8) = (/ & |
---|
2530 | ! band 8 |
---|
2531 | 1.41881e-01_rb, 7.15419e-02_rb, 6.30335e-02_rb, 6.11132e-02_rb, 6.01931e-02_rb, & |
---|
2532 | 5.92420e-02_rb, 5.78968e-02_rb, 5.58876e-02_rb, 5.28923e-02_rb, 4.84462e-02_rb, & |
---|
2533 | 4.60839e-02_rb, 4.56013e-02_rb, 4.45410e-02_rb, 4.31866e-02_rb, 4.17026e-02_rb, & |
---|
2534 | 4.01850e-02_rb, 3.86892e-02_rb, 3.72461e-02_rb, 3.58722e-02_rb, 3.45749e-02_rb, & |
---|
2535 | 3.33564e-02_rb, 3.22155e-02_rb, 3.11494e-02_rb, 3.01541e-02_rb, 2.92253e-02_rb, & |
---|
2536 | 2.83584e-02_rb, 2.75488e-02_rb, 2.67925e-02_rb, 2.57692e-02_rb, 2.50704e-02_rb, & |
---|
2537 | 2.43918e-02_rb, 2.37350e-02_rb, 2.31005e-02_rb, 2.24888e-02_rb, 2.18996e-02_rb, & |
---|
2538 | 2.13325e-02_rb, 2.07870e-02_rb, 2.02623e-02_rb, 1.97577e-02_rb, 1.92724e-02_rb, & |
---|
2539 | 1.88056e-02_rb, 1.83564e-02_rb, 1.79241e-02_rb, 1.75079e-02_rb, 1.71070e-02_rb, & |
---|
2540 | 1.67207e-02_rb, 1.63482e-02_rb, 1.59890e-02_rb, 1.56424e-02_rb, 1.53077e-02_rb, & |
---|
2541 | 1.49845e-02_rb, 1.46722e-02_rb, 1.43702e-02_rb, 1.40782e-02_rb, 1.37955e-02_rb, & |
---|
2542 | 1.35219e-02_rb, 1.32569e-02_rb, 1.30000e-02_rb/) |
---|
2543 | absliq1(:, 9) = (/ & |
---|
2544 | ! band 9 |
---|
2545 | 6.72726e-02_rb, 6.61013e-02_rb, 6.47866e-02_rb, 6.33780e-02_rb, 6.18985e-02_rb, & |
---|
2546 | 6.03335e-02_rb, 5.86136e-02_rb, 5.65876e-02_rb, 5.39839e-02_rb, 5.03536e-02_rb, & |
---|
2547 | 4.71608e-02_rb, 4.63630e-02_rb, 4.50313e-02_rb, 4.34526e-02_rb, 4.17876e-02_rb, & |
---|
2548 | 4.01261e-02_rb, 3.85171e-02_rb, 3.69860e-02_rb, 3.55442e-02_rb, 3.41954e-02_rb, & |
---|
2549 | 3.29384e-02_rb, 3.17693e-02_rb, 3.06832e-02_rb, 2.96745e-02_rb, 2.87374e-02_rb, & |
---|
2550 | 2.78662e-02_rb, 2.70557e-02_rb, 2.63008e-02_rb, 2.52450e-02_rb, 2.45424e-02_rb, & |
---|
2551 | 2.38656e-02_rb, 2.32144e-02_rb, 2.25885e-02_rb, 2.19873e-02_rb, 2.14099e-02_rb, & |
---|
2552 | 2.08554e-02_rb, 2.03230e-02_rb, 1.98116e-02_rb, 1.93203e-02_rb, 1.88482e-02_rb, & |
---|
2553 | 1.83944e-02_rb, 1.79578e-02_rb, 1.75378e-02_rb, 1.71335e-02_rb, 1.67440e-02_rb, & |
---|
2554 | 1.63687e-02_rb, 1.60069e-02_rb, 1.56579e-02_rb, 1.53210e-02_rb, 1.49958e-02_rb, & |
---|
2555 | 1.46815e-02_rb, 1.43778e-02_rb, 1.40841e-02_rb, 1.37999e-02_rb, 1.35249e-02_rb, & |
---|
2556 | 1.32585e-02_rb, 1.30004e-02_rb, 1.27502e-02_rb/) |
---|
2557 | absliq1(:,10) = (/ & |
---|
2558 | ! band 10 |
---|
2559 | 7.97040e-02_rb, 7.63844e-02_rb, 7.36499e-02_rb, 7.13525e-02_rb, 6.93043e-02_rb, & |
---|
2560 | 6.72807e-02_rb, 6.50227e-02_rb, 6.22395e-02_rb, 5.86093e-02_rb, 5.37815e-02_rb, & |
---|
2561 | 5.14682e-02_rb, 4.97214e-02_rb, 4.77392e-02_rb, 4.56961e-02_rb, 4.36858e-02_rb, & |
---|
2562 | 4.17569e-02_rb, 3.99328e-02_rb, 3.82224e-02_rb, 3.66265e-02_rb, 3.51416e-02_rb, & |
---|
2563 | 3.37617e-02_rb, 3.24798e-02_rb, 3.12887e-02_rb, 3.01812e-02_rb, 2.91505e-02_rb, & |
---|
2564 | 2.81900e-02_rb, 2.72939e-02_rb, 2.64568e-02_rb, 2.54165e-02_rb, 2.46832e-02_rb, & |
---|
2565 | 2.39783e-02_rb, 2.33017e-02_rb, 2.26531e-02_rb, 2.20314e-02_rb, 2.14359e-02_rb, & |
---|
2566 | 2.08653e-02_rb, 2.03187e-02_rb, 1.97947e-02_rb, 1.92924e-02_rb, 1.88106e-02_rb, & |
---|
2567 | 1.83483e-02_rb, 1.79043e-02_rb, 1.74778e-02_rb, 1.70678e-02_rb, 1.66735e-02_rb, & |
---|
2568 | 1.62941e-02_rb, 1.59286e-02_rb, 1.55766e-02_rb, 1.52371e-02_rb, 1.49097e-02_rb, & |
---|
2569 | 1.45937e-02_rb, 1.42885e-02_rb, 1.39936e-02_rb, 1.37085e-02_rb, 1.34327e-02_rb, & |
---|
2570 | 1.31659e-02_rb, 1.29075e-02_rb, 1.26571e-02_rb/) |
---|
2571 | absliq1(:,11) = (/ & |
---|
2572 | ! band 11 |
---|
2573 | 1.49438e-01_rb, 1.33535e-01_rb, 1.21542e-01_rb, 1.11743e-01_rb, 1.03263e-01_rb, & |
---|
2574 | 9.55774e-02_rb, 8.83382e-02_rb, 8.12943e-02_rb, 7.42533e-02_rb, 6.70609e-02_rb, & |
---|
2575 | 6.38761e-02_rb, 5.97788e-02_rb, 5.59841e-02_rb, 5.25318e-02_rb, 4.94132e-02_rb, & |
---|
2576 | 4.66014e-02_rb, 4.40644e-02_rb, 4.17706e-02_rb, 3.96910e-02_rb, 3.77998e-02_rb, & |
---|
2577 | 3.60742e-02_rb, 3.44947e-02_rb, 3.30442e-02_rb, 3.17079e-02_rb, 3.04730e-02_rb, & |
---|
2578 | 2.93283e-02_rb, 2.82642e-02_rb, 2.72720e-02_rb, 2.61789e-02_rb, 2.53277e-02_rb, & |
---|
2579 | 2.45237e-02_rb, 2.37635e-02_rb, 2.30438e-02_rb, 2.23615e-02_rb, 2.17140e-02_rb, & |
---|
2580 | 2.10987e-02_rb, 2.05133e-02_rb, 1.99557e-02_rb, 1.94241e-02_rb, 1.89166e-02_rb, & |
---|
2581 | 1.84317e-02_rb, 1.79679e-02_rb, 1.75238e-02_rb, 1.70983e-02_rb, 1.66901e-02_rb, & |
---|
2582 | 1.62983e-02_rb, 1.59219e-02_rb, 1.55599e-02_rb, 1.52115e-02_rb, 1.48761e-02_rb, & |
---|
2583 | 1.45528e-02_rb, 1.42411e-02_rb, 1.39402e-02_rb, 1.36497e-02_rb, 1.33690e-02_rb, & |
---|
2584 | 1.30976e-02_rb, 1.28351e-02_rb, 1.25810e-02_rb/) |
---|
2585 | absliq1(:,12) = (/ & |
---|
2586 | ! band 12 |
---|
2587 | 3.71985e-02_rb, 3.88586e-02_rb, 3.99070e-02_rb, 4.04351e-02_rb, 4.04610e-02_rb, & |
---|
2588 | 3.99834e-02_rb, 3.89953e-02_rb, 3.74886e-02_rb, 3.54551e-02_rb, 3.28870e-02_rb, & |
---|
2589 | 3.32576e-02_rb, 3.22444e-02_rb, 3.12384e-02_rb, 3.02584e-02_rb, 2.93146e-02_rb, & |
---|
2590 | 2.84120e-02_rb, 2.75525e-02_rb, 2.67361e-02_rb, 2.59618e-02_rb, 2.52280e-02_rb, & |
---|
2591 | 2.45327e-02_rb, 2.38736e-02_rb, 2.32487e-02_rb, 2.26558e-02_rb, 2.20929e-02_rb, & |
---|
2592 | 2.15579e-02_rb, 2.10491e-02_rb, 2.05648e-02_rb, 1.99749e-02_rb, 1.95704e-02_rb, & |
---|
2593 | 1.91731e-02_rb, 1.87839e-02_rb, 1.84032e-02_rb, 1.80315e-02_rb, 1.76689e-02_rb, & |
---|
2594 | 1.73155e-02_rb, 1.69712e-02_rb, 1.66362e-02_rb, 1.63101e-02_rb, 1.59928e-02_rb, & |
---|
2595 | 1.56842e-02_rb, 1.53840e-02_rb, 1.50920e-02_rb, 1.48080e-02_rb, 1.45318e-02_rb, & |
---|
2596 | 1.42631e-02_rb, 1.40016e-02_rb, 1.37472e-02_rb, 1.34996e-02_rb, 1.32586e-02_rb, & |
---|
2597 | 1.30239e-02_rb, 1.27954e-02_rb, 1.25728e-02_rb, 1.23559e-02_rb, 1.21445e-02_rb, & |
---|
2598 | 1.19385e-02_rb, 1.17376e-02_rb, 1.15417e-02_rb/) |
---|
2599 | absliq1(:,13) = (/ & |
---|
2600 | ! band 13 |
---|
2601 | 3.11868e-02_rb, 4.48357e-02_rb, 4.90224e-02_rb, 4.96406e-02_rb, 4.86806e-02_rb, & |
---|
2602 | 4.69610e-02_rb, 4.48630e-02_rb, 4.25795e-02_rb, 4.02138e-02_rb, 3.78236e-02_rb, & |
---|
2603 | 3.74266e-02_rb, 3.60384e-02_rb, 3.47074e-02_rb, 3.34434e-02_rb, 3.22499e-02_rb, & |
---|
2604 | 3.11264e-02_rb, 3.00704e-02_rb, 2.90784e-02_rb, 2.81463e-02_rb, 2.72702e-02_rb, & |
---|
2605 | 2.64460e-02_rb, 2.56698e-02_rb, 2.49381e-02_rb, 2.42475e-02_rb, 2.35948e-02_rb, & |
---|
2606 | 2.29774e-02_rb, 2.23925e-02_rb, 2.18379e-02_rb, 2.11793e-02_rb, 2.07076e-02_rb, & |
---|
2607 | 2.02470e-02_rb, 1.97981e-02_rb, 1.93613e-02_rb, 1.89367e-02_rb, 1.85243e-02_rb, & |
---|
2608 | 1.81240e-02_rb, 1.77356e-02_rb, 1.73588e-02_rb, 1.69935e-02_rb, 1.66392e-02_rb, & |
---|
2609 | 1.62956e-02_rb, 1.59624e-02_rb, 1.56393e-02_rb, 1.53259e-02_rb, 1.50219e-02_rb, & |
---|
2610 | 1.47268e-02_rb, 1.44404e-02_rb, 1.41624e-02_rb, 1.38925e-02_rb, 1.36302e-02_rb, & |
---|
2611 | 1.33755e-02_rb, 1.31278e-02_rb, 1.28871e-02_rb, 1.26530e-02_rb, 1.24253e-02_rb, & |
---|
2612 | 1.22038e-02_rb, 1.19881e-02_rb, 1.17782e-02_rb/) |
---|
2613 | absliq1(:,14) = (/ & |
---|
2614 | ! band 14 |
---|
2615 | 1.58988e-02_rb, 3.50652e-02_rb, 4.00851e-02_rb, 4.07270e-02_rb, 3.98101e-02_rb, & |
---|
2616 | 3.83306e-02_rb, 3.66829e-02_rb, 3.50327e-02_rb, 3.34497e-02_rb, 3.19609e-02_rb, & |
---|
2617 | 3.13712e-02_rb, 3.03348e-02_rb, 2.93415e-02_rb, 2.83973e-02_rb, 2.75037e-02_rb, & |
---|
2618 | 2.66604e-02_rb, 2.58654e-02_rb, 2.51161e-02_rb, 2.44100e-02_rb, 2.37440e-02_rb, & |
---|
2619 | 2.31154e-02_rb, 2.25215e-02_rb, 2.19599e-02_rb, 2.14282e-02_rb, 2.09242e-02_rb, & |
---|
2620 | 2.04459e-02_rb, 1.99915e-02_rb, 1.95594e-02_rb, 1.90254e-02_rb, 1.86598e-02_rb, & |
---|
2621 | 1.82996e-02_rb, 1.79455e-02_rb, 1.75983e-02_rb, 1.72584e-02_rb, 1.69260e-02_rb, & |
---|
2622 | 1.66013e-02_rb, 1.62843e-02_rb, 1.59752e-02_rb, 1.56737e-02_rb, 1.53799e-02_rb, & |
---|
2623 | 1.50936e-02_rb, 1.48146e-02_rb, 1.45429e-02_rb, 1.42782e-02_rb, 1.40203e-02_rb, & |
---|
2624 | 1.37691e-02_rb, 1.35243e-02_rb, 1.32858e-02_rb, 1.30534e-02_rb, 1.28270e-02_rb, & |
---|
2625 | 1.26062e-02_rb, 1.23909e-02_rb, 1.21810e-02_rb, 1.19763e-02_rb, 1.17766e-02_rb, & |
---|
2626 | 1.15817e-02_rb, 1.13915e-02_rb, 1.12058e-02_rb/) |
---|
2627 | absliq1(:,15) = (/ & |
---|
2628 | ! band 15 |
---|
2629 | 5.02079e-03_rb, 2.17615e-02_rb, 2.55449e-02_rb, 2.59484e-02_rb, 2.53650e-02_rb, & |
---|
2630 | 2.45281e-02_rb, 2.36843e-02_rb, 2.29159e-02_rb, 2.22451e-02_rb, 2.16716e-02_rb, & |
---|
2631 | 2.11451e-02_rb, 2.05817e-02_rb, 2.00454e-02_rb, 1.95372e-02_rb, 1.90567e-02_rb, & |
---|
2632 | 1.86028e-02_rb, 1.81742e-02_rb, 1.77693e-02_rb, 1.73866e-02_rb, 1.70244e-02_rb, & |
---|
2633 | 1.66815e-02_rb, 1.63563e-02_rb, 1.60477e-02_rb, 1.57544e-02_rb, 1.54755e-02_rb, & |
---|
2634 | 1.52097e-02_rb, 1.49564e-02_rb, 1.47146e-02_rb, 1.43684e-02_rb, 1.41728e-02_rb, & |
---|
2635 | 1.39762e-02_rb, 1.37797e-02_rb, 1.35838e-02_rb, 1.33891e-02_rb, 1.31961e-02_rb, & |
---|
2636 | 1.30051e-02_rb, 1.28164e-02_rb, 1.26302e-02_rb, 1.24466e-02_rb, 1.22659e-02_rb, & |
---|
2637 | 1.20881e-02_rb, 1.19131e-02_rb, 1.17412e-02_rb, 1.15723e-02_rb, 1.14063e-02_rb, & |
---|
2638 | 1.12434e-02_rb, 1.10834e-02_rb, 1.09264e-02_rb, 1.07722e-02_rb, 1.06210e-02_rb, & |
---|
2639 | 1.04725e-02_rb, 1.03269e-02_rb, 1.01839e-02_rb, 1.00436e-02_rb, 9.90593e-03_rb, & |
---|
2640 | 9.77080e-03_rb, 9.63818e-03_rb, 9.50800e-03_rb/) |
---|
2641 | absliq1(:,16) = (/ & |
---|
2642 | ! band 16 |
---|
2643 | 5.64971e-02_rb, 9.04736e-02_rb, 8.11726e-02_rb, 7.05450e-02_rb, 6.20052e-02_rb, & |
---|
2644 | 5.54286e-02_rb, 5.03503e-02_rb, 4.63791e-02_rb, 4.32290e-02_rb, 4.06959e-02_rb, & |
---|
2645 | 3.74690e-02_rb, 3.52964e-02_rb, 3.33799e-02_rb, 3.16774e-02_rb, 3.01550e-02_rb, & |
---|
2646 | 2.87856e-02_rb, 2.75474e-02_rb, 2.64223e-02_rb, 2.53953e-02_rb, 2.44542e-02_rb, & |
---|
2647 | 2.35885e-02_rb, 2.27894e-02_rb, 2.20494e-02_rb, 2.13622e-02_rb, 2.07222e-02_rb, & |
---|
2648 | 2.01246e-02_rb, 1.95654e-02_rb, 1.90408e-02_rb, 1.84398e-02_rb, 1.80021e-02_rb, & |
---|
2649 | 1.75816e-02_rb, 1.71775e-02_rb, 1.67889e-02_rb, 1.64152e-02_rb, 1.60554e-02_rb, & |
---|
2650 | 1.57089e-02_rb, 1.53751e-02_rb, 1.50531e-02_rb, 1.47426e-02_rb, 1.44428e-02_rb, & |
---|
2651 | 1.41532e-02_rb, 1.38734e-02_rb, 1.36028e-02_rb, 1.33410e-02_rb, 1.30875e-02_rb, & |
---|
2652 | 1.28420e-02_rb, 1.26041e-02_rb, 1.23735e-02_rb, 1.21497e-02_rb, 1.19325e-02_rb, & |
---|
2653 | 1.17216e-02_rb, 1.15168e-02_rb, 1.13177e-02_rb, 1.11241e-02_rb, 1.09358e-02_rb, & |
---|
2654 | 1.07525e-02_rb, 1.05741e-02_rb, 1.04003e-02_rb/) |
---|
2655 | |
---|
2656 | end subroutine lwcldpr |
---|
2657 | |
---|
2658 | end module rrtmg_lw_init |
---|
2659 | |
---|