Home

Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

rrtmg_lw_init.f90 @ 4358

Last change on this file since 4358 was 1585, checked in by maronga, 10 years ago
Added support for RRTMG radiation code
File size: 107.7 KB

Rev	Line
[1585]	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.planckclightclight1.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=1616, ngn=2561., 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

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: palm/trunk/LIB/rrtmg/rrtmg_lw_init.f90 @ 4358

Download in other formats: