3 &tauaerl,fac00,fac01,fac10,fac11,jp,jt,jt1,oneminus,&
4 &colh2o,coln2o,laytrop,selffac,selffrac,indself,pfrac)
37 real_b :: fac00(
jplay)
38 real_b :: fac01(
jplay)
39 real_b :: fac10(
jplay)
40 real_b :: fac11(
jplay)
43 integer_m :: jp(
jplay)
44 integer_m :: jt(
jplay)
45 integer_m :: jt1(
jplay)
51 real_b :: colh2o(
jplay)
52 real_b :: coln2o(
jplay)
56 real_b :: selffac(
jplay)
57 real_b :: selffrac(
jplay)
58 integer_m :: indself(
jplay)
63 integer_m :: ijs(
jplay)
68 integer_m :: ig, js, lay
71 real_b :: fac000, fac001, fac010, fac011, fac100, fac101,&
72 &fac110, fac111, fs, specmult, specparm
82 speccomb(lay) = colh2o(lay) +
strrat*coln2o(lay)
83 specparm = colh2o(lay)/speccomb(lay)
84 specparm=min(oneminus,specparm)
85 specmult = 8._jprb*(specparm)
86 js = 1 + int(specmult)
87 fs = mod(specmult,_one_)
88 ind0(lay) = ((jp(lay)-1)*5+(jt(lay)-1))*
nspa(13) + js
89 ind1(lay) = (jp(lay)*5+(jt1(lay)-1))*
nspa(13) + js
90 inds(lay) = indself(lay)
114 tau(
ngs12+ig,lay) = speccomb(lay) *&
116 &( (1. - fs) *( fac00(lay) *
absa(ind0(lay) ,ig) + &
117 & fac10(lay) *
absa(ind0(lay)+ 9,ig) + &
118 & fac01(lay) *
absa(ind1(lay) ,ig) + &
119 & fac11(lay) *
absa(ind1(lay)+ 9,ig))+ &
120 & fs *( fac00(lay) *
absa(ind0(lay)+ 1,ig) + &
121 & fac10(lay) *
absa(ind0(lay)+10,ig) + &
122 & fac01(lay) *
absa(ind1(lay)+ 1,ig) + &
123 & fac11(lay) *
absa(ind1(lay)+10,ig))) + &
134 &selffac(lay) * (
selfref(inds(lay),ig) + &
145 DO lay = laytrop+1,
klev
147 tau(
ngs12+ig,lay) = tauaerl(lay,13)
148 pfrac(
ngs12+ig,lay) = _zero_
real(kind=jprb), dimension(9, 5, 13, ng13) ka
subroutine rrtm_taumol13(KLEV, P_TAU, P_TAUAERL, P_FAC00, P_FAC01, P_FAC10, P_FAC11, K_JP, K_JT, K_JT1, P_ONEMINUS, P_COLH2O, P_COLN2O, K_LAYTROP, P_SELFFAC, P_SELFFRAC, K_INDSELF, PFRAC)
real(kind=jprb), dimension(ng13, 9) fracrefa
integer(kind=jpim), parameter jpgpt
integer(kind=jpim), parameter jpband
real(kind=jprb), dimension(10, ng13) selfref
integer(kind=jpim), parameter ngs12
integer(kind=jpim), dimension(16) nspb
integer(kind=jpim), dimension(16) ng
integer(kind=jpim), parameter ng13
integer(kind=jpim), dimension(16) nspa
real(kind=jprb), dimension(585, ng13) absa
integer(kind=jpim), parameter jplay
integer(kind=jpim), parameter jpxsec