source: LMDZ6/branches/LMDZ_cdrag_LSCE/libf/phylmd/ecrad/rrtm_taumol9.F90 @ 5428

Last change on this file since 5428 was 3908, checked in by idelkadi, 4 years ago

Online implementation of the radiative transfer code ECRAD in the LMDZ model.

  • Inclusion of the ecrad directory containing the sources of the ECRAD code
    • interface routine : radiation_scheme.F90
  • Adaptation of compilation scripts :
    • compilation under CPP key CPP_ECRAD
    • compilation with option "-rad ecard" or "-ecard true"
    • The "-rad old/rtm/ecran" build option will need to replace the "-rrtm true" and "-ecrad true" options in the future.
  • Runing LMDZ simulations with ecrad, you need :
    • logical key iflag_rrtm = 2 in physiq.def
    • namelist_ecrad (DefLists?)
    • the directory "data" containing the configuration files is temporarily placed in ../libfphylmd/ecrad/
  • Compilation and execution are tested in the 1D case. The repository under svn would allow to continue the implementation work: tests, verification of the results, ...
File size: 14.1 KB
Line 
1!----------------------------------------------------------------------------
2SUBROUTINE RRTM_TAUMOL9 (KIDIA,KFDIA,KLEV,P_TAU,&
3 & P_TAUAERL,P_FAC00,P_FAC01,P_FAC10,P_FAC11,P_FORFAC,P_FORFRAC,K_INDFOR,K_JP,K_JT,K_JT1,P_ONEMINUS,&
4 & P_COLH2O,P_COLN2O,P_COLCH4,P_COLDRY,K_LAYTROP,K_LAYSWTCH,K_LAYLOW,P_SELFFAC,P_SELFFRAC,K_INDSELF,PFRAC, &
5 & PRAT_H2OCH4,PRAT_H2OCH4_1,PMINORFRAC,KINDMINOR) 
6
7!     BAND 9:  1180-1390 cm-1 (low - H2O,CH4; high - CH4)
8
9!     AUTHOR.
10!     -------
11!      JJMorcrette, ECMWF
12
13!     MODIFICATIONS.
14!     --------------
15!      M.Hamrud      01-Oct-2003 CY28 Cleaning
16!      NEC           25-Oct-2007 Optimisations
17!      JJMorcrette 20110613 flexible number of g-points
18!      201306 ABozzo updated to rrtmg v4.85
19!     band 9:  1180-1390 cm-1 (low key - h2o,ch4; low minor - n2o)
20!                             (high key - ch4; high minor - n2o)
21! ---------------------------------------------------------------------------
22
23USE PARKIND1  ,ONLY : JPIM     ,JPRB
24USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
25
26USE PARRRTM  , ONLY : JPBAND
27USE YOERRTM  , ONLY : JPGPT  ,NG9   ,NGS8
28USE YOERRTWN , ONLY :      NSPA   ,NSPB
29USE YOERRTA9 , ONLY : ABSA   ,ABSB   ,FRACREFA, FRACREFB,SELFREF,FORREF,KA_MN2O, KB_MN2O   
30USE YOERRTRF, ONLY : CHI_MLS
31
32IMPLICIT NONE
33
34INTEGER(KIND=JPIM),INTENT(IN)    :: KIDIA
35INTEGER(KIND=JPIM),INTENT(IN)    :: KFDIA
36INTEGER(KIND=JPIM),INTENT(IN)    :: KLEV
37REAL(KIND=JPRB)   ,INTENT(OUT)   :: P_TAU(KIDIA:KFDIA,JPGPT,KLEV)
38REAL(KIND=JPRB)   ,INTENT(IN)    :: P_TAUAERL(KIDIA:KFDIA,KLEV,JPBAND)
39REAL(KIND=JPRB)   ,INTENT(IN)    :: P_FAC00(KIDIA:KFDIA,KLEV)
40REAL(KIND=JPRB)   ,INTENT(IN)    :: P_FAC01(KIDIA:KFDIA,KLEV)
41REAL(KIND=JPRB)   ,INTENT(IN)    :: P_FAC10(KIDIA:KFDIA,KLEV)
42REAL(KIND=JPRB)   ,INTENT(IN)    :: P_FAC11(KIDIA:KFDIA,KLEV)
43INTEGER(KIND=JPIM),INTENT(IN)    :: K_JP(KIDIA:KFDIA,KLEV)
44INTEGER(KIND=JPIM),INTENT(IN)    :: K_JT(KIDIA:KFDIA,KLEV)
45INTEGER(KIND=JPIM),INTENT(IN)    :: K_JT1(KIDIA:KFDIA,KLEV)
46REAL(KIND=JPRB)   ,INTENT(IN)    :: P_ONEMINUS
47REAL(KIND=JPRB)   ,INTENT(IN)    :: P_COLH2O(KIDIA:KFDIA,KLEV)
48REAL(KIND=JPRB)   ,INTENT(IN)    :: P_COLN2O(KIDIA:KFDIA,KLEV)
49REAL(KIND=JPRB)   ,INTENT(IN)    :: P_COLCH4(KIDIA:KFDIA,KLEV)
50REAL(KIND=JPRB)   ,INTENT(IN)    :: P_COLDRY(KIDIA:KFDIA,KLEV)
51INTEGER(KIND=JPIM),INTENT(IN)    :: K_LAYTROP(KIDIA:KFDIA)
52INTEGER(KIND=JPIM),INTENT(IN)    :: K_LAYSWTCH(KIDIA:KFDIA)
53INTEGER(KIND=JPIM),INTENT(IN)    :: K_LAYLOW(KIDIA:KFDIA)
54REAL(KIND=JPRB)   ,INTENT(IN)    :: P_SELFFAC(KIDIA:KFDIA,KLEV)
55REAL(KIND=JPRB)   ,INTENT(IN)    :: P_SELFFRAC(KIDIA:KFDIA,KLEV)
56INTEGER(KIND=JPIM),INTENT(IN)    :: K_INDSELF(KIDIA:KFDIA,KLEV)
57REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFRAC(KIDIA:KFDIA,JPGPT,KLEV)
58
59REAL(KIND=JPRB)   ,INTENT(IN)   :: PRAT_H2OCH4(KIDIA:KFDIA,KLEV)
60REAL(KIND=JPRB)   ,INTENT(IN)   :: PRAT_H2OCH4_1(KIDIA:KFDIA,KLEV)
61INTEGER(KIND=JPIM),INTENT(IN)   :: K_INDFOR(KIDIA:KFDIA,KLEV)
62REAL(KIND=JPRB)   ,INTENT(IN)   :: P_FORFAC(KIDIA:KFDIA,KLEV)
63REAL(KIND=JPRB)   ,INTENT(IN)   :: P_FORFRAC(KIDIA:KFDIA,KLEV)
64REAL(KIND=JPRB)   ,INTENT(IN)   :: PMINORFRAC(KIDIA:KFDIA,KLEV)
65INTEGER(KIND=JPIM),INTENT(IN)   :: KINDMINOR(KIDIA:KFDIA,KLEV)
66
67
68! ---------------------------------------------------------------------------
69
70INTEGER(KIND=JPIM) :: IND0(KLEV),IND1(KLEV),INDS(KLEV),INDF(KLEV),INDM(KLEV)
71INTEGER(KIND=JPIM) :: IG, JS, JLAY, JS1, JMN2O, JPL
72INTEGER(KIND=JPIM) :: JLON
73
74REAL(KIND=JPRB) :: Z_SPECCOMB(KLEV),Z_SPECCOMB1(KLEV),Z_SPECCOMB_MN2O(KLEV), &
75& Z_SPECCOMB_PLANCK(KLEV)
76REAL(KIND=JPRB) :: ZREFRAT_PLANCK_A, ZREFRAT_M_A
77
78REAL(KIND=JPRB) :: Z_FS, Z_SPECMULT, Z_SPECPARM,  &
79 & Z_FS1, Z_SPECMULT1, Z_SPECPARM1,   &
80 & Z_FMN2O, Z_SPECMULT_MN2O, Z_SPECPARM_MN2O,   &
81 & Z_FPL, Z_SPECMULT_PLANCK, Z_SPECPARM_PLANCK
82
83REAL(KIND=JPRB) :: ZADJFAC,ZADJCOLN2O(KIDIA:KFDIA,KLEV),ZRATN2O,ZCHI_N2O
84REAL(KIND=JPRB) ::  Z_FAC000, Z_FAC100, Z_FAC200,&
85 & Z_FAC010, Z_FAC110, Z_FAC210, &
86 & Z_FAC001, Z_FAC101, Z_FAC201, &
87 & Z_FAC011, Z_FAC111, Z_FAC211
88REAL(KIND=JPRB) :: ZP, ZP4, ZFK0, ZFK1, ZFK2
89REAL(KIND=JPRB) :: ZTAUFOR,ZTAUSELF,ZN2OM1,ZN2OM2,ZABSN2O,ZTAU_MAJOR,ZTAU_MAJOR1
90
91REAL(KIND=JPRB) :: ZHOOK_HANDLE
92
93ASSOCIATE(NFLEVG=>KLEV)
94IF (LHOOK) CALL DR_HOOK('RRTM_TAUMOL9',0,ZHOOK_HANDLE)
95
96! Minor gas mapping level :
97!     lower - n2o, p = 706.272 mbar, t = 278.94 k
98!     upper - n2o, p = 95.58 mbar, t = 215.7 k
99
100! Calculate reference ratio to be used in calculation of Planck
101! fraction in lower/upper atmosphere.
102
103! P = 212 mb
104      ZREFRAT_PLANCK_A = CHI_MLS(1,9)/CHI_MLS(6,9)
105
106! P = 706.272 mb
107      ZREFRAT_M_A = CHI_MLS(1,3)/CHI_MLS(6,3)
108
109! Compute the optical depth by interpolating in ln(pressure),
110! temperature, and appropriate species.  Below laytrop, the water
111! vapor self-continuum and foreign continuum is interpolated
112! (in temperature) separately. 
113
114DO JLAY = 1, KLEV
115  DO JLON = KIDIA, KFDIA
116
117    IF (JLAY <= K_LAYTROP(JLON)) THEN
118      Z_SPECCOMB(JLAY) = P_COLH2O(JLON,JLAY) + PRAT_H2OCH4(JLON,JLAY)*P_COLCH4(JLON,JLAY)
119      Z_SPECPARM = P_COLH2O(JLON,JLAY)/Z_SPECCOMB(JLAY)
120      Z_SPECPARM=MIN(P_ONEMINUS,Z_SPECPARM)
121      Z_SPECMULT = 8._JPRB*(Z_SPECPARM)
122      JS = 1 + INT(Z_SPECMULT)
123      Z_FS = MOD(Z_SPECMULT,1.0_JPRB)
124
125      Z_SPECCOMB1(JLAY) = P_COLH2O(JLON,JLAY) + PRAT_H2OCH4_1(JLON,JLAY)*P_COLCH4(JLON,JLAY)
126      Z_SPECPARM1 = P_COLH2O(JLON,JLAY)/Z_SPECCOMB1(JLAY)
127      IF (Z_SPECPARM1 >= P_ONEMINUS) Z_SPECPARM1 = P_ONEMINUS
128      Z_SPECMULT1 = 8._JPRB*(Z_SPECPARM1)
129      JS1 = 1 + INT(Z_SPECMULT1)
130      Z_FS1 = MOD(Z_SPECMULT1,1.0_JPRB)
131
132      Z_SPECCOMB_MN2O(JLAY) = P_COLH2O(JLON,JLAY) + ZREFRAT_M_A*P_COLCH4(JLON,JLAY)
133      Z_SPECPARM_MN2O = P_COLH2O(JLON,JLAY)/Z_SPECCOMB_MN2O(JLAY)
134      IF (Z_SPECPARM_MN2O >= P_ONEMINUS) Z_SPECPARM_MN2O = P_ONEMINUS
135      Z_SPECMULT_MN2O = 8._JPRB*Z_SPECPARM_MN2O
136      JMN2O = 1 + INT(Z_SPECMULT_MN2O)
137      Z_FMN2O = MOD(Z_SPECMULT_MN2O,1.0_JPRB)
138!  In atmospheres where the amount of N2O is too great to be considered
139!  a minor species, adjust the column amount of N2O by an empirical factor
140!  to obtain the proper contribution.
141      ZCHI_N2O = P_COLN2O(JLON,JLAY)/P_COLDRY(JLON,JLAY)
142      ZRATN2O = 1.E20_JPRB*ZCHI_N2O/CHI_MLS(4,K_JP(JLON,JLAY)+1)
143      IF (ZRATN2O > 1.5_JPRB) THEN
144         ZADJFAC = 0.5_JPRB+(ZRATN2O-0.5_JPRB)**0.65_JPRB
145         ZADJCOLN2O(JLON,JLAY) = ZADJFAC*CHI_MLS(4,K_JP(JLON,JLAY)+1)*P_COLDRY(JLON,JLAY)*1.E-20_JPRB
146      ELSE
147         ZADJCOLN2O(JLON,JLAY) = P_COLN2O(JLON,JLAY)
148      ENDIF
149
150      Z_SPECCOMB_PLANCK(JLAY) = P_COLH2O(JLON,JLAY)+ZREFRAT_PLANCK_A*P_COLCH4(JLON,JLAY)
151      Z_SPECPARM_PLANCK = P_COLH2O(JLON,JLAY)/Z_SPECCOMB_PLANCK(JLAY)
152      IF (Z_SPECPARM_PLANCK >= P_ONEMINUS) Z_SPECPARM_PLANCK=P_ONEMINUS
153      Z_SPECMULT_PLANCK = 8._JPRB*Z_SPECPARM_PLANCK
154      JPL= 1 + INT(Z_SPECMULT_PLANCK)
155      Z_FPL = MOD(Z_SPECMULT_PLANCK,1.0_JPRB)
156
157
158      IND0(JLAY) = ((K_JP(JLON,JLAY)-1)*5+(K_JT(JLON,JLAY)-1))*NSPA(9) + JS
159      IND1(JLAY) = (K_JP(JLON,JLAY)*5+(K_JT1(JLON,JLAY)-1))*NSPA(9) + JS1
160      INDS(JLAY) = K_INDSELF(JLON,JLAY)
161      INDF(JLAY) = K_INDFOR(JLON,JLAY)
162      INDM(JLAY) = KINDMINOR(JLON,JLAY)
163
164         IF (Z_SPECPARM < 0.125_JPRB) THEN
165            ZP = Z_FS - 1
166            ZP4 = ZP**4
167            ZFK0 = ZP4
168            ZFK1 = 1 - ZP - 2.0_JPRB*ZP4
169            ZFK2 = ZP + ZP4
170            Z_FAC000 = ZFK0*P_FAC00(JLON,JLAY)
171            Z_FAC100 = ZFK1*P_FAC00(JLON,JLAY)
172            Z_FAC200 = ZFK2*P_FAC00(JLON,JLAY)
173            Z_FAC010 = ZFK0*P_FAC10(JLON,JLAY)
174            Z_FAC110 = ZFK1*P_FAC10(JLON,JLAY)
175            Z_FAC210 = ZFK2*P_FAC10(JLON,JLAY)
176         ELSEIF (Z_SPECPARM > 0.875_JPRB) THEN
177            ZP = -Z_FS
178            ZP4 = ZP**4
179            ZFK0 = ZP4
180            ZFK1 = 1 - ZP - 2.0_JPRB*ZP4
181            ZFK2 = ZP + ZP4
182            Z_FAC000 = ZFK0*P_FAC00(JLON,JLAY)
183            Z_FAC100 = ZFK1*P_FAC00(JLON,JLAY)
184            Z_FAC200 = ZFK2*P_FAC00(JLON,JLAY)
185            Z_FAC010 = ZFK0*P_FAC10(JLON,JLAY)
186            Z_FAC110 = ZFK1*P_FAC10(JLON,JLAY)
187            Z_FAC210 = ZFK2*P_FAC10(JLON,JLAY)
188         ELSE
189            Z_FAC000 = (1._JPRB - Z_FS) * P_FAC00(JLON,JLAY)
190            Z_FAC010 = (1._JPRB - Z_FS) * P_FAC10(JLON,JLAY)
191            Z_FAC100 = Z_FS * P_FAC00(JLON,JLAY)
192            Z_FAC110 = Z_FS * P_FAC10(JLON,JLAY)
193         ENDIF
194         IF (Z_SPECPARM1 < 0.125_JPRB) THEN
195            ZP = Z_FS1 - 1
196            ZP4 = ZP**4
197            ZFK0 = ZP4
198            ZFK1 = 1 - ZP - 2.0_JPRB*ZP4
199            ZFK2 = ZP + ZP4
200            Z_FAC001 = ZFK0*P_FAC01(JLON,JLAY)
201            Z_FAC101 = ZFK1*P_FAC01(JLON,JLAY)
202            Z_FAC201 = ZFK2*P_FAC01(JLON,JLAY)
203            Z_FAC011 = ZFK0*P_FAC11(JLON,JLAY)
204            Z_FAC111 = ZFK1*P_FAC11(JLON,JLAY)
205            Z_FAC211 = ZFK2*P_FAC11(JLON,JLAY)
206         ELSEIF (Z_SPECPARM1 > 0.875_JPRB) THEN
207            ZP = -Z_FS1
208            ZP4 = ZP**4
209            ZFK0 = ZP4
210            ZFK1 = 1 - ZP - 2.0_JPRB*ZP4
211            ZFK2 = ZP + ZP4
212            Z_FAC001 = ZFK0*P_FAC01(JLON,JLAY)
213            Z_FAC101 = ZFK1*P_FAC01(JLON,JLAY)
214            Z_FAC201 = ZFK2*P_FAC01(JLON,JLAY)
215            Z_FAC011 = ZFK0*P_FAC11(JLON,JLAY)
216            Z_FAC111 = ZFK1*P_FAC11(JLON,JLAY)
217            Z_FAC211 = ZFK2*P_FAC11(JLON,JLAY)
218         ELSE
219            Z_FAC001 = (1._JPRB - Z_FS1) * P_FAC01(JLON,JLAY)
220            Z_FAC011 = (1._JPRB - Z_FS1) * P_FAC11(JLON,JLAY)
221            Z_FAC101 = Z_FS1 * P_FAC01(JLON,JLAY)
222            Z_FAC111 = Z_FS1 * P_FAC11(JLON,JLAY)
223         ENDIF
224
225
226
227!-- DS_000515
228!CDIR UNROLL=NG9
229      DO IG = 1, NG9
230!-- DS_000515
231        ZTAUSELF = P_SELFFAC(JLON,JLAY)* (SELFREF(INDS(JLAY),IG) + P_SELFFRAC(JLON,JLAY) * &
232          &      (SELFREF(INDS(JLAY)+1,IG) - SELFREF(INDS(JLAY),IG)))
233        ZTAUFOR = P_FORFAC(JLON,JLAY) * (FORREF(INDF(JLAY),IG) + P_FORFRAC(JLON,JLAY) * &
234          &      (FORREF(INDF(JLAY)+1,IG) - FORREF(INDF(JLAY),IG)))
235        ZN2OM1 = KA_MN2O(JMN2O,INDM(JLAY),IG) + Z_FMN2O * &
236          &      (KA_MN2O(JMN2O+1,INDM(JLAY),IG) - KA_MN2O(JMN2O,INDM(JLAY),IG))
237        ZN2OM2 = KA_MN2O(JMN2O,INDM(JLAY)+1,IG) + Z_FMN2O * &
238          &      (KA_MN2O(JMN2O+1,INDM(JLAY)+1,IG) - KA_MN2O(JMN2O,INDM(JLAY)+1,IG))
239        ZABSN2O = ZN2OM1 + PMINORFRAC(JLON,JLAY) * (ZN2OM2 - ZN2OM1)
240
241   IF (Z_SPECPARM < 0.125_JPRB) THEN
242               ZTAU_MAJOR = Z_SPECCOMB(JLAY) * &
243                 &  (Z_FAC000 * ABSA(IND0(JLAY),IG) + &
244                 &  Z_FAC100 * ABSA(IND0(JLAY)+1,IG) + &
245                 &  Z_FAC200 * ABSA(IND0(JLAY)+2,IG) + &
246                 &  Z_FAC010 * ABSA(IND0(JLAY)+9,IG) + &
247                 &  Z_FAC110 * ABSA(IND0(JLAY)+10,IG) + &
248                 &  Z_FAC210 * ABSA(IND0(JLAY)+11,IG))
249            ELSEIF (Z_SPECPARM > 0.875_JPRB) THEN
250               ZTAU_MAJOR = Z_SPECCOMB(JLAY) * &
251                 &  (Z_FAC200 * ABSA(IND0(JLAY)-1,IG) + &
252                 &  Z_FAC100 * ABSA(IND0(JLAY),IG) + &
253                 &  Z_FAC000 * ABSA(IND0(JLAY)+1,IG) + &
254                 &  Z_FAC210 * ABSA(IND0(JLAY)+8,IG) + &
255                 &  Z_FAC110 * ABSA(IND0(JLAY)+9,IG) + &
256                 &  Z_FAC010 * ABSA(IND0(JLAY)+10,IG))
257            ELSE
258               ZTAU_MAJOR = Z_SPECCOMB(JLAY) * &
259                 &  (Z_FAC000 * ABSA(IND0(JLAY),IG) + &
260                 &  Z_FAC100 * ABSA(IND0(JLAY)+1,IG) + &
261                 &  Z_FAC010 * ABSA(IND0(JLAY)+9,IG) + &
262                 &  Z_FAC110 * ABSA(IND0(JLAY)+10,IG))
263            ENDIF
264
265            IF (Z_SPECPARM1 < 0.125_JPRB) THEN
266               ZTAU_MAJOR1 = Z_SPECCOMB1(JLAY) * &
267                 &  (Z_FAC001 * ABSA(IND1(JLAY),IG) + &
268                 &  Z_FAC101 * ABSA(IND1(JLAY)+1,IG) + &
269                 &  Z_FAC201 * ABSA(IND1(JLAY)+2,IG) + &
270                 &  Z_FAC011 * ABSA(IND1(JLAY)+9,IG) + &
271                 &  Z_FAC111 * ABSA(IND1(JLAY)+10,IG) + &
272                 &  Z_FAC211 * ABSA(IND1(JLAY)+11,IG))
273            ELSEIF (Z_SPECPARM1 > 0.875_JPRB) THEN
274               ZTAU_MAJOR1 = Z_SPECCOMB1(JLAY) * &
275                 &  (Z_FAC201 * ABSA(IND1(JLAY)-1,IG) + &
276                 &  Z_FAC101 * ABSA(IND1(JLAY),IG) + &
277                 &  Z_FAC001 * ABSA(IND1(JLAY)+1,IG) + &
278                 &  Z_FAC211 * ABSA(IND1(JLAY)+8,IG) + &
279                 &  Z_FAC111 * ABSA(IND1(JLAY)+9,IG) + &
280                 &  Z_FAC011 * ABSA(IND1(JLAY)+10,IG))
281            ELSE
282               ZTAU_MAJOR1 = Z_SPECCOMB1(JLAY) * &
283                 &  (Z_FAC001 * ABSA(IND1(JLAY),IG) +  &
284                 &  Z_FAC101 * ABSA(IND1(JLAY)+1,IG) + &
285                 &  Z_FAC011 * ABSA(IND1(JLAY)+9,IG) + &
286                 &  Z_FAC111 * ABSA(IND1(JLAY)+10,IG))
287            ENDIF
288
289
290        P_TAU(JLON,NGS8+IG,JLAY) = ZTAU_MAJOR + ZTAU_MAJOR1 &
291         &       + ZTAUSELF + ZTAUFOR &
292         &       + ZADJCOLN2O(JLON,JLAY)*ZABSN2O &
293         & + P_TAUAERL(JLON,JLAY,9) 
294        PFRAC(JLON,NGS8+IG,JLAY) = FRACREFA(IG,JPL) + Z_FPL *&
295         & (FRACREFA(IG,JPL+1) - FRACREFA(IG,JPL))
296      ENDDO
297    ENDIF
298
299    IF (JLAY > K_LAYTROP(JLON)) THEN
300
301!  In atmospheres where the amount of N2O is too great to be considered
302!  a minor species, adjust the column amount of N2O by an empirical factor
303!  to obtain the proper contribution.
304      ZCHI_N2O = P_COLN2O(JLON,JLAY)/P_COLDRY(JLON,JLAY)
305      ZRATN2O = 1.E20_JPRB*ZCHI_N2O/CHI_MLS(4,K_JP(JLON,JLAY)+1)
306      IF (ZRATN2O > 1.5_JPRB) THEN
307         ZADJFAC = 0.5_JPRB+(ZRATN2O-0.5_JPRB)**0.65_JPRB
308         ZADJCOLN2O(JLON,JLAY) = ZADJFAC*CHI_MLS(4,K_JP(JLON,JLAY)+1)*P_COLDRY(JLON,JLAY)*1.E-20_JPRB
309      ELSE
310         ZADJCOLN2O(JLON,JLAY) = P_COLN2O(JLON,JLAY)
311      ENDIF
312
313      IND0(JLAY) = ((K_JP(JLON,JLAY)-13)*5+(K_JT(JLON,JLAY)-1))*NSPB(9) + 1
314      IND1(JLAY) = ((K_JP(JLON,JLAY)-12)*5+(K_JT1(JLON,JLAY)-1))*NSPB(9) + 1
315      INDM(JLAY) = KINDMINOR(JLON,JLAY)
316
317
318!-- JJM_000517
319!CDIR UNROLL=NG9
320      DO IG = 1, NG9
321        ZABSN2O = KB_MN2O(INDM(JLAY),IG) + PMINORFRAC(JLON,JLAY) * &
322         &       (KB_MN2O(INDM(JLAY)+1,IG) - KB_MN2O(INDM(JLAY),IG))
323
324!-- JJM_000517
325        P_TAU(JLON,NGS8+IG,JLAY) = P_COLCH4(JLON,JLAY) *&
326         & (P_FAC00(JLON,JLAY) * ABSB(IND0(JLAY)  ,IG) +&
327         & P_FAC10(JLON,JLAY) * ABSB(IND0(JLAY)+1,IG) +&
328         & P_FAC01(JLON,JLAY) * ABSB(IND1(JLAY)  ,IG) +&
329         & P_FAC11(JLON,JLAY) * ABSB(IND1(JLAY)+1,IG))&
330         & + ZADJCOLN2O(JLON,JLAY)*ZABSN2O &
331         & + P_TAUAERL(JLON,JLAY,9) 
332        PFRAC(JLON,NGS8+IG,JLAY) = FRACREFB(IG)
333      ENDDO
334    ENDIF
335  ENDDO
336ENDDO
337
338IF (LHOOK) CALL DR_HOOK('RRTM_TAUMOL9',1,ZHOOK_HANDLE)
339
340END ASSOCIATE
341END SUBROUTINE RRTM_TAUMOL9
Note: See TracBrowser for help on using the repository browser.