[2027] | 1 | ! |
---|
| 2 | ! $Id: rrtm_ecrt_140gp.F90 2626 2016-09-15 14:20:56Z adurocher $ |
---|
| 3 | ! |
---|
[1989] | 4 | !****************** SUBROUTINE RRTM_ECRT_140GP ************************** |
---|
| 5 | |
---|
| 6 | SUBROUTINE RRTM_ECRT_140GP & |
---|
| 7 | & ( K_IPLON, klon , klev, kcld,& |
---|
| 8 | & paer , paph , pap,& |
---|
| 9 | & pts , pth , pt,& |
---|
| 10 | & P_ZEMIS, P_ZEMIW,& |
---|
| 11 | & pq , pcco2, pozn, pcldf, ptaucld, ptclear,& |
---|
[2146] | 12 | & P_CLDFRAC,P_TAUCLD,& |
---|
| 13 | & PTAU_LW,& |
---|
| 14 | & P_COLDRY,P_WKL,P_WX,& |
---|
| 15 | & P_TAUAERL,PAVEL,P_TAVEL,PZ,P_TZ,P_TBOUND,K_NLAYERS,P_SEMISS,K_IREFLECT ) |
---|
[1989] | 16 | |
---|
| 17 | ! Reformatted for F90 by JJMorcrette, ECMWF, 980714 |
---|
| 18 | |
---|
| 19 | ! Read in atmospheric profile from ECMWF radiation code, and prepare it |
---|
| 20 | ! for use in RRTM. Set other RRTM input parameters. Values are passed |
---|
| 21 | ! back through existing RRTM arrays and commons. |
---|
| 22 | |
---|
| 23 | !- Modifications |
---|
| 24 | |
---|
| 25 | ! 2000-05-15 Deborah Salmond Speed-up |
---|
| 26 | |
---|
| 27 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
---|
| 28 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
---|
| 29 | |
---|
| 30 | USE PARRRTM , ONLY : JPBAND ,JPXSEC ,JPLAY ,& |
---|
| 31 | & JPINPX |
---|
[2146] | 32 | USE YOERAD , ONLY : NLW ,NOVLP |
---|
[2626] | 33 | !MPL/IM 20160915 on prend GES de phylmd USE YOERDI , ONLY : RCH4 ,RN2O ,RCFC11 ,RCFC12 |
---|
[1989] | 34 | USE YOESW , ONLY : RAER |
---|
| 35 | |
---|
| 36 | !------------------------------Arguments-------------------------------- |
---|
| 37 | |
---|
[2146] | 38 | IMPLICIT NONE |
---|
[2027] | 39 | |
---|
| 40 | |
---|
[1989] | 41 | INTEGER(KIND=JPIM),INTENT(IN) :: KLON! Number of atmospheres (longitudes) |
---|
| 42 | INTEGER(KIND=JPIM),INTENT(IN) :: KLEV! Number of atmospheric layers |
---|
| 43 | INTEGER(KIND=JPIM),INTENT(IN) :: K_IPLON |
---|
| 44 | INTEGER(KIND=JPIM),INTENT(OUT) :: KCLD |
---|
| 45 | REAL(KIND=JPRB) ,INTENT(IN) :: PAER(KLON,6,KLEV) ! Aerosol optical thickness |
---|
| 46 | REAL(KIND=JPRB) ,INTENT(IN) :: PAPH(KLON,KLEV+1) ! Interface pressures (Pa) |
---|
| 47 | REAL(KIND=JPRB) ,INTENT(IN) :: PAP(KLON,KLEV) ! Layer pressures (Pa) |
---|
| 48 | REAL(KIND=JPRB) ,INTENT(IN) :: PTS(KLON) ! Surface temperature (K) |
---|
| 49 | REAL(KIND=JPRB) ,INTENT(IN) :: PTH(KLON,KLEV+1) ! Interface temperatures (K) |
---|
| 50 | REAL(KIND=JPRB) ,INTENT(IN) :: PT(KLON,KLEV) ! Layer temperature (K) |
---|
| 51 | REAL(KIND=JPRB) ,INTENT(IN) :: P_ZEMIS(KLON) ! Non-window surface emissivity |
---|
| 52 | REAL(KIND=JPRB) ,INTENT(IN) :: P_ZEMIW(KLON) ! Window surface emissivity |
---|
| 53 | REAL(KIND=JPRB) ,INTENT(IN) :: PQ(KLON,KLEV) ! H2O specific humidity (mmr) |
---|
| 54 | REAL(KIND=JPRB) ,INTENT(IN) :: PCCO2 ! CO2 mass mixing ratio |
---|
| 55 | REAL(KIND=JPRB) ,INTENT(IN) :: POZN(KLON,KLEV) ! O3 mass mixing ratio |
---|
| 56 | REAL(KIND=JPRB) ,INTENT(IN) :: PCLDF(KLON,KLEV) ! Cloud fraction |
---|
| 57 | REAL(KIND=JPRB) ,INTENT(IN) :: PTAUCLD(KLON,KLEV,JPBAND) ! Cloud optical depth |
---|
[2146] | 58 | !--C.Kleinschmitt |
---|
| 59 | REAL(KIND=JPRB) ,INTENT(IN) :: PTAU_LW(KLON,KLEV,NLW) ! LW Optical depth of aerosols |
---|
| 60 | !--end |
---|
[1989] | 61 | REAL(KIND=JPRB) ,INTENT(OUT) :: PTCLEAR |
---|
| 62 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_CLDFRAC(JPLAY) ! Cloud fraction |
---|
| 63 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_TAUCLD(JPLAY,JPBAND) ! Spectral optical thickness |
---|
| 64 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_COLDRY(JPLAY) |
---|
| 65 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_WKL(JPINPX,JPLAY) |
---|
| 66 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_WX(JPXSEC,JPLAY) ! Amount of trace gases |
---|
| 67 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_TAUAERL(JPLAY,JPBAND) |
---|
| 68 | REAL(KIND=JPRB) ,INTENT(OUT) :: PAVEL(JPLAY) |
---|
| 69 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_TAVEL(JPLAY) |
---|
| 70 | REAL(KIND=JPRB) ,INTENT(OUT) :: PZ(0:JPLAY) |
---|
| 71 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_TZ(0:JPLAY) |
---|
| 72 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_TBOUND |
---|
| 73 | INTEGER(KIND=JPIM),INTENT(OUT) :: K_NLAYERS |
---|
| 74 | REAL(KIND=JPRB) ,INTENT(OUT) :: P_SEMISS(JPBAND) |
---|
| 75 | INTEGER(KIND=JPIM),INTENT(OUT) :: K_IREFLECT |
---|
| 76 | ! real rch4 ! CH4 mass mixing ratio |
---|
| 77 | ! real rn2o ! N2O mass mixing ratio |
---|
| 78 | ! real rcfc11 ! CFC11 mass mixing ratio |
---|
| 79 | ! real rcfc12 ! CFC12 mass mixing ratio |
---|
| 80 | !- from AER |
---|
| 81 | !- from PROFILE |
---|
| 82 | !- from SURFACE |
---|
| 83 | REAL(KIND=JPRB) :: ztauaer(5) |
---|
| 84 | REAL(KIND=JPRB) :: zc1j(0:klev) ! total cloud from top and level k |
---|
| 85 | REAL(KIND=JPRB) :: Z_AMD ! Effective molecular weight of dry air (g/mol) |
---|
| 86 | REAL(KIND=JPRB) :: Z_AMW ! Molecular weight of water vapor (g/mol) |
---|
| 87 | REAL(KIND=JPRB) :: Z_AMCO2 ! Molecular weight of carbon dioxide (g/mol) |
---|
| 88 | REAL(KIND=JPRB) :: Z_AMO ! Molecular weight of ozone (g/mol) |
---|
| 89 | REAL(KIND=JPRB) :: Z_AMCH4 ! Molecular weight of methane (g/mol) |
---|
| 90 | REAL(KIND=JPRB) :: Z_AMN2O ! Molecular weight of nitrous oxide (g/mol) |
---|
| 91 | REAL(KIND=JPRB) :: Z_AMC11 ! Molecular weight of CFC11 (g/mol) - CFCL3 |
---|
| 92 | REAL(KIND=JPRB) :: Z_AMC12 ! Molecular weight of CFC12 (g/mol) - CF2CL2 |
---|
| 93 | REAL(KIND=JPRB) :: Z_AVGDRO ! Avogadro's number (molecules/mole) |
---|
| 94 | REAL(KIND=JPRB) :: Z_GRAVIT ! Gravitational acceleration (cm/sec2) |
---|
| 95 | |
---|
| 96 | ! Atomic weights for conversion from mass to volume mixing ratios; these |
---|
| 97 | ! are the same values used in ECRT to assure accurate conversion to vmr |
---|
| 98 | data Z_AMD / 28.970_JPRB / |
---|
| 99 | data Z_AMW / 18.0154_JPRB / |
---|
| 100 | data Z_AMCO2 / 44.011_JPRB / |
---|
| 101 | data Z_AMO / 47.9982_JPRB / |
---|
| 102 | data Z_AMCH4 / 16.043_JPRB / |
---|
| 103 | data Z_AMN2O / 44.013_JPRB / |
---|
| 104 | data Z_AMC11 / 137.3686_JPRB / |
---|
| 105 | data Z_AMC12 / 120.9140_JPRB / |
---|
| 106 | data Z_AVGDRO/ 6.02214E23_JPRB / |
---|
| 107 | data Z_GRAVIT/ 9.80665E02_JPRB / |
---|
| 108 | |
---|
| 109 | INTEGER(KIND=JPIM) :: IATM, IMOL, IXMAX, J1, J2, JAE, JB, JK, JL, I_L |
---|
| 110 | INTEGER(KIND=JPIM) :: I_NMOL, I_NXMOL |
---|
| 111 | |
---|
| 112 | REAL(KIND=JPRB) :: Z_AMM, ZCLDLY, ZCLEAR, ZCLOUD, ZEPSEC |
---|
| 113 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
---|
| 114 | |
---|
[2626] | 115 | !MPL/IM 20160915 on prend GES de phylmd |
---|
| 116 | #include "clesphys.h" |
---|
[1989] | 117 | ! *** |
---|
| 118 | |
---|
| 119 | ! *** mji |
---|
| 120 | ! Initialize all molecular amounts and aerosol optical depths to zero here, |
---|
| 121 | ! then pass ECRT amounts into RRTM arrays below. |
---|
| 122 | |
---|
| 123 | ! DATA ZWKL /MAXPRDW*0.0/ |
---|
| 124 | ! DATA ZWX /MAXPROD*0.0/ |
---|
| 125 | ! DATA KREFLECT /0/ |
---|
| 126 | |
---|
| 127 | ! Activate cross section molecules: |
---|
| 128 | ! NXMOL - number of cross-sections input by user |
---|
| 129 | ! IXINDX(I) - index of cross-section molecule corresponding to Ith |
---|
| 130 | ! cross-section specified by user |
---|
| 131 | ! = 0 -- not allowed in RRTM |
---|
| 132 | ! = 1 -- CCL4 |
---|
| 133 | ! = 2 -- CFC11 |
---|
| 134 | ! = 3 -- CFC12 |
---|
| 135 | ! = 4 -- CFC22 |
---|
| 136 | ! DATA KXMOL /2/ |
---|
| 137 | ! DATA KXINDX /0,2,3,0,31*0/ |
---|
| 138 | |
---|
| 139 | ! IREFLECT=KREFLECT |
---|
| 140 | ! NXMOL=KXMOL |
---|
| 141 | |
---|
| 142 | IF (LHOOK) CALL DR_HOOK('RRTM_ECRT_140GP',0,ZHOOK_HANDLE) |
---|
| 143 | K_IREFLECT=0 |
---|
| 144 | I_NXMOL=2 |
---|
| 145 | |
---|
| 146 | DO J1=1,35 |
---|
| 147 | ! IXINDX(J1)=0 |
---|
| 148 | DO J2=1,KLEV |
---|
| 149 | P_WKL(J1,J2)=0.0_JPRB |
---|
| 150 | ENDDO |
---|
| 151 | ENDDO |
---|
| 152 | !IXINDX(2)=2 |
---|
| 153 | !IXINDX(3)=3 |
---|
| 154 | |
---|
| 155 | ! Set parameters needed for RRTM execution: |
---|
| 156 | IATM = 0 |
---|
| 157 | ! IXSECT = 1 |
---|
| 158 | ! NUMANGS = 0 |
---|
| 159 | ! IOUT = -1 |
---|
| 160 | IXMAX = 4 |
---|
| 161 | |
---|
| 162 | ! Bands 6,7,8 are considered the 'window' and allowed to have a |
---|
| 163 | ! different surface emissivity (as in ECMWF). Eli wrote this part.... |
---|
| 164 | P_SEMISS(1) = P_ZEMIS(K_IPLON) |
---|
| 165 | P_SEMISS(2) = P_ZEMIS(K_IPLON) |
---|
| 166 | P_SEMISS(3) = P_ZEMIS(K_IPLON) |
---|
| 167 | P_SEMISS(4) = P_ZEMIS(K_IPLON) |
---|
| 168 | P_SEMISS(5) = P_ZEMIS(K_IPLON) |
---|
| 169 | P_SEMISS(6) = P_ZEMIW(K_IPLON) |
---|
| 170 | P_SEMISS(7) = P_ZEMIW(K_IPLON) |
---|
| 171 | P_SEMISS(8) = P_ZEMIW(K_IPLON) |
---|
| 172 | P_SEMISS(9) = P_ZEMIS(K_IPLON) |
---|
| 173 | P_SEMISS(10) = P_ZEMIS(K_IPLON) |
---|
| 174 | P_SEMISS(11) = P_ZEMIS(K_IPLON) |
---|
| 175 | P_SEMISS(12) = P_ZEMIS(K_IPLON) |
---|
| 176 | P_SEMISS(13) = P_ZEMIS(K_IPLON) |
---|
| 177 | P_SEMISS(14) = P_ZEMIS(K_IPLON) |
---|
| 178 | P_SEMISS(15) = P_ZEMIS(K_IPLON) |
---|
| 179 | P_SEMISS(16) = P_ZEMIS(K_IPLON) |
---|
| 180 | |
---|
| 181 | ! Set surface temperature. |
---|
| 182 | |
---|
| 183 | P_TBOUND = pts(K_IPLON) |
---|
| 184 | |
---|
| 185 | ! Install ECRT arrays into RRTM arrays for pressure, temperature, |
---|
| 186 | ! and molecular amounts. Pressures are converted from Pascals |
---|
| 187 | ! (ECRT) to mb (RRTM). H2O, CO2, O3 and trace gas amounts are |
---|
| 188 | ! converted from mass mixing ratio to volume mixing ratio. CO2 |
---|
| 189 | ! converted with same dry air and CO2 molecular weights used in |
---|
| 190 | ! ECRT to assure correct conversion back to the proper CO2 vmr. |
---|
| 191 | ! The dry air column COLDRY (in molec/cm2) is calculated from |
---|
| 192 | ! the level pressures PZ (in mb) based on the hydrostatic equation |
---|
| 193 | ! and includes a correction to account for H2O in the layer. The |
---|
| 194 | ! molecular weight of moist air (amm) is calculated for each layer. |
---|
| 195 | ! Note: RRTM levels count from bottom to top, while the ECRT input |
---|
| 196 | ! variables count from the top down and must be reversed here. |
---|
| 197 | |
---|
| 198 | K_NLAYERS = klev |
---|
| 199 | I_NMOL = 6 |
---|
| 200 | PZ(0) = paph(K_IPLON,klev+1)/100._JPRB |
---|
| 201 | P_TZ(0) = pth(K_IPLON,klev+1) |
---|
| 202 | DO I_L = 1, KLEV |
---|
| 203 | PAVEL(I_L) = pap(K_IPLON,KLEV-I_L+1)/100._JPRB |
---|
| 204 | P_TAVEL(I_L) = pt(K_IPLON,KLEV-I_L+1) |
---|
| 205 | PZ(I_L) = paph(K_IPLON,KLEV-I_L+1)/100._JPRB |
---|
| 206 | P_TZ(I_L) = pth(K_IPLON,KLEV-I_L+1) |
---|
| 207 | P_WKL(1,I_L) = pq(K_IPLON,KLEV-I_L+1)*Z_AMD/Z_AMW |
---|
| 208 | P_WKL(2,I_L) = pcco2*Z_AMD/Z_AMCO2 |
---|
| 209 | P_WKL(3,I_L) = pozn(K_IPLON,KLEV-I_L+1)*Z_AMD/Z_AMO |
---|
| 210 | P_WKL(4,I_L) = rn2o*Z_AMD/Z_AMN2O |
---|
| 211 | P_WKL(6,I_L) = rch4*Z_AMD/Z_AMCH4 |
---|
| 212 | Z_AMM = (1-P_WKL(1,I_L))*Z_AMD + P_WKL(1,I_L)*Z_AMW |
---|
| 213 | P_COLDRY(I_L) = (PZ(I_L-1)-PZ(I_L))*1.E3_JPRB*Z_AVGDRO/(Z_GRAVIT*Z_AMM*(1+P_WKL(1,I_L))) |
---|
| 214 | ENDDO |
---|
| 215 | |
---|
| 216 | !- Fill RRTM aerosol arrays with operational ECMWF aerosols, |
---|
| 217 | ! do the mixing and distribute over the 16 spectral intervals |
---|
| 218 | |
---|
| 219 | DO I_L=1,KLEV |
---|
| 220 | JK=KLEV-I_L+1 |
---|
| 221 | ! DO JAE=1,5 |
---|
| 222 | JAE=1 |
---|
| 223 | ZTAUAER(JAE) =& |
---|
| 224 | & RAER(JAE,1)*PAER(K_IPLON,1,JK)+RAER(JAE,2)*PAER(K_IPLON,2,JK)& |
---|
| 225 | & +RAER(JAE,3)*PAER(K_IPLON,3,JK)+RAER(JAE,4)*PAER(K_IPLON,4,JK)& |
---|
| 226 | & +RAER(JAE,5)*PAER(K_IPLON,5,JK)+RAER(JAE,6)*PAER(K_IPLON,6,JK) |
---|
| 227 | P_TAUAERL(I_L, 1)=ZTAUAER(1) |
---|
| 228 | P_TAUAERL(I_L, 2)=ZTAUAER(1) |
---|
| 229 | JAE=2 |
---|
| 230 | ZTAUAER(JAE) =& |
---|
| 231 | & RAER(JAE,1)*PAER(K_IPLON,1,JK)+RAER(JAE,2)*PAER(K_IPLON,2,JK)& |
---|
| 232 | & +RAER(JAE,3)*PAER(K_IPLON,3,JK)+RAER(JAE,4)*PAER(K_IPLON,4,JK)& |
---|
| 233 | & +RAER(JAE,5)*PAER(K_IPLON,5,JK)+RAER(JAE,6)*PAER(K_IPLON,6,JK) |
---|
| 234 | P_TAUAERL(I_L, 3)=ZTAUAER(2) |
---|
| 235 | P_TAUAERL(I_L, 4)=ZTAUAER(2) |
---|
| 236 | P_TAUAERL(I_L, 5)=ZTAUAER(2) |
---|
| 237 | JAE=3 |
---|
| 238 | ZTAUAER(JAE) =& |
---|
| 239 | & RAER(JAE,1)*PAER(K_IPLON,1,JK)+RAER(JAE,2)*PAER(K_IPLON,2,JK)& |
---|
| 240 | & +RAER(JAE,3)*PAER(K_IPLON,3,JK)+RAER(JAE,4)*PAER(K_IPLON,4,JK)& |
---|
| 241 | & +RAER(JAE,5)*PAER(K_IPLON,5,JK)+RAER(JAE,6)*PAER(K_IPLON,6,JK) |
---|
| 242 | P_TAUAERL(I_L, 6)=ZTAUAER(3) |
---|
| 243 | P_TAUAERL(I_L, 8)=ZTAUAER(3) |
---|
| 244 | P_TAUAERL(I_L, 9)=ZTAUAER(3) |
---|
| 245 | JAE=4 |
---|
| 246 | ZTAUAER(JAE) =& |
---|
| 247 | & RAER(JAE,1)*PAER(K_IPLON,1,JK)+RAER(JAE,2)*PAER(K_IPLON,2,JK)& |
---|
| 248 | & +RAER(JAE,3)*PAER(K_IPLON,3,JK)+RAER(JAE,4)*PAER(K_IPLON,4,JK)& |
---|
| 249 | & +RAER(JAE,5)*PAER(K_IPLON,5,JK)+RAER(JAE,6)*PAER(K_IPLON,6,JK) |
---|
| 250 | P_TAUAERL(I_L, 7)=ZTAUAER(4) |
---|
| 251 | JAE=5 |
---|
| 252 | ZTAUAER(JAE) =& |
---|
| 253 | & RAER(JAE,1)*PAER(K_IPLON,1,JK)+RAER(JAE,2)*PAER(K_IPLON,2,JK)& |
---|
| 254 | & +RAER(JAE,3)*PAER(K_IPLON,3,JK)+RAER(JAE,4)*PAER(K_IPLON,4,JK)& |
---|
| 255 | & +RAER(JAE,5)*PAER(K_IPLON,5,JK)+RAER(JAE,6)*PAER(K_IPLON,6,JK) |
---|
| 256 | ! END DO |
---|
| 257 | P_TAUAERL(I_L,10)=ZTAUAER(5) |
---|
| 258 | P_TAUAERL(I_L,11)=ZTAUAER(5) |
---|
| 259 | P_TAUAERL(I_L,12)=ZTAUAER(5) |
---|
| 260 | P_TAUAERL(I_L,13)=ZTAUAER(5) |
---|
| 261 | P_TAUAERL(I_L,14)=ZTAUAER(5) |
---|
| 262 | P_TAUAERL(I_L,15)=ZTAUAER(5) |
---|
| 263 | P_TAUAERL(I_L,16)=ZTAUAER(5) |
---|
| 264 | ENDDO |
---|
[2146] | 265 | !--Use LW AOD from own Mie calculations (C. Kleinschmitt) |
---|
| 266 | DO I_L=1,KLEV |
---|
| 267 | JK=KLEV-I_L+1 |
---|
[2152] | 268 | DO JAE=1, NLW |
---|
[2146] | 269 | P_TAUAERL(I_L,JAE) = MAX( PTAU_LW(K_IPLON, JK, JAE), 1e-30 ) |
---|
| 270 | ENDDO |
---|
| 271 | ENDDO |
---|
| 272 | !--end C. Kleinschmitt |
---|
[1989] | 273 | |
---|
| 274 | DO J2=1,KLEV |
---|
| 275 | DO J1=1,JPXSEC |
---|
| 276 | P_WX(J1,J2)=0.0_JPRB |
---|
| 277 | ENDDO |
---|
| 278 | ENDDO |
---|
| 279 | |
---|
| 280 | DO I_L = 1, KLEV |
---|
| 281 | !- Set cross section molecule amounts from ECRT; convert to vmr |
---|
| 282 | P_WX(2,I_L) = rcfc11*Z_AMD/Z_AMC11 |
---|
| 283 | P_WX(3,I_L) = rcfc12*Z_AMD/Z_AMC12 |
---|
| 284 | P_WX(2,I_L) = P_COLDRY(I_L) * P_WX(2,I_L) * 1.E-20_JPRB |
---|
| 285 | P_WX(3,I_L) = P_COLDRY(I_L) * P_WX(3,I_L) * 1.E-20_JPRB |
---|
| 286 | |
---|
| 287 | !- Here, all molecules in WKL and WX are in volume mixing ratio; convert to |
---|
| 288 | ! molec/cm2 based on COLDRY for use in RRTM |
---|
| 289 | |
---|
| 290 | DO IMOL = 1, I_NMOL |
---|
| 291 | P_WKL(IMOL,I_L) = P_COLDRY(I_L) * P_WKL(IMOL,I_L) |
---|
| 292 | ENDDO |
---|
| 293 | |
---|
| 294 | ! DO IX = 1,JPXSEC |
---|
| 295 | ! IF (IXINDX(IX) /= 0) THEN |
---|
| 296 | ! WX(IXINDX(IX),L) = COLDRY(L) * WX(IX,L) * 1.E-20_JPRB |
---|
| 297 | ! ENDIF |
---|
| 298 | ! END DO |
---|
| 299 | |
---|
| 300 | ENDDO |
---|
| 301 | |
---|
| 302 | !- Approximate treatment for various cloud overlaps |
---|
| 303 | ZCLEAR=1.0_JPRB |
---|
| 304 | ZCLOUD=0.0_JPRB |
---|
| 305 | ZC1J(0)=0.0_JPRB |
---|
| 306 | ZEPSEC=1.E-03_JPRB |
---|
| 307 | JL=K_IPLON |
---|
| 308 | |
---|
| 309 | !++MODIFCODE |
---|
| 310 | IF ((NOVLP == 1).OR.(NOVLP ==6).OR.(NOVLP ==8)) THEN |
---|
| 311 | !--MODIFCODE |
---|
| 312 | |
---|
| 313 | DO JK=1,KLEV |
---|
| 314 | IF (pcldf(JL,JK) > ZEPSEC) THEN |
---|
| 315 | ZCLDLY=pcldf(JL,JK) |
---|
| 316 | ZCLEAR=ZCLEAR & |
---|
| 317 | & *(1.0_JPRB-MAX( ZCLDLY , ZCLOUD ))& |
---|
| 318 | & /(1.0_JPRB-MIN( ZCLOUD , 1.0_JPRB-ZEPSEC )) |
---|
| 319 | ZCLOUD = ZCLDLY |
---|
| 320 | ZC1J(JK)= 1.0_JPRB - ZCLEAR |
---|
| 321 | ELSE |
---|
| 322 | ZCLDLY=0.0_JPRB |
---|
| 323 | ZCLEAR=ZCLEAR & |
---|
| 324 | & *(1.0_JPRB-MAX( ZCLDLY , ZCLOUD ))& |
---|
| 325 | & /(1.0_JPRB-MIN( ZCLOUD , 1.0_JPRB-ZEPSEC )) |
---|
| 326 | ZCLOUD = ZCLDLY |
---|
| 327 | ZC1J(JK)= 1.0_JPRB - ZCLEAR |
---|
| 328 | ENDIF |
---|
| 329 | ENDDO |
---|
| 330 | |
---|
| 331 | !++MODIFCODE |
---|
| 332 | ELSEIF ((NOVLP == 2).OR.(NOVLP ==7)) THEN |
---|
| 333 | !--MODIFCODE |
---|
| 334 | |
---|
| 335 | DO JK=1,KLEV |
---|
| 336 | IF (pcldf(JL,JK) > ZEPSEC) THEN |
---|
| 337 | ZCLDLY=pcldf(JL,JK) |
---|
| 338 | ZCLOUD = MAX( ZCLDLY , ZCLOUD ) |
---|
| 339 | ZC1J(JK) = ZCLOUD |
---|
| 340 | ELSE |
---|
| 341 | ZCLDLY=0.0_JPRB |
---|
| 342 | ZCLOUD = MAX( ZCLDLY , ZCLOUD ) |
---|
| 343 | ZC1J(JK) = ZCLOUD |
---|
| 344 | ENDIF |
---|
| 345 | ENDDO |
---|
| 346 | |
---|
| 347 | !++MODIFCODE |
---|
| 348 | ELSEIF ((NOVLP == 3).OR.(NOVLP ==5)) THEN |
---|
| 349 | !--MODIFCODE |
---|
| 350 | |
---|
| 351 | DO JK=1,KLEV |
---|
| 352 | IF (pcldf(JL,JK) > ZEPSEC) THEN |
---|
| 353 | ZCLDLY=pcldf(JL,JK) |
---|
| 354 | ZCLEAR = ZCLEAR * (1.0_JPRB-ZCLDLY) |
---|
| 355 | ZCLOUD = 1.0_JPRB - ZCLEAR |
---|
| 356 | ZC1J(JK) = ZCLOUD |
---|
| 357 | ELSE |
---|
| 358 | ZCLDLY=0.0_JPRB |
---|
| 359 | ZCLEAR = ZCLEAR * (1.0_JPRB-ZCLDLY) |
---|
| 360 | ZCLOUD = 1.0_JPRB - ZCLEAR |
---|
| 361 | ZC1J(JK) = ZCLOUD |
---|
| 362 | ENDIF |
---|
| 363 | ENDDO |
---|
| 364 | |
---|
| 365 | ELSEIF (NOVLP == 4) THEN |
---|
| 366 | |
---|
| 367 | ENDIF |
---|
| 368 | PTCLEAR=1.0_JPRB-ZC1J(KLEV) |
---|
| 369 | |
---|
| 370 | ! Transfer cloud fraction and cloud optical depth to RRTM arrays; |
---|
| 371 | ! invert array index for pcldf to go from bottom to top for RRTM |
---|
| 372 | |
---|
| 373 | !- clear-sky column |
---|
| 374 | IF (PTCLEAR > 1.0_JPRB-ZEPSEC) THEN |
---|
| 375 | KCLD=0 |
---|
| 376 | DO I_L = 1, KLEV |
---|
| 377 | P_CLDFRAC(I_L) = 0.0_JPRB |
---|
| 378 | ENDDO |
---|
| 379 | DO JB=1,JPBAND |
---|
| 380 | DO I_L=1,KLEV |
---|
| 381 | P_TAUCLD(I_L,JB) = 0.0_JPRB |
---|
| 382 | ENDDO |
---|
| 383 | ENDDO |
---|
| 384 | |
---|
| 385 | ELSE |
---|
| 386 | |
---|
| 387 | !- cloudy column |
---|
| 388 | ! The diffusivity factor (Savijarvi, 1997) on the cloud optical |
---|
| 389 | ! thickness TAUCLD has already been applied in RADLSW |
---|
| 390 | |
---|
| 391 | KCLD=1 |
---|
| 392 | DO I_L=1,KLEV |
---|
| 393 | P_CLDFRAC(I_L) = pcldf(K_IPLON,I_L) |
---|
| 394 | ENDDO |
---|
| 395 | DO JB=1,JPBAND |
---|
| 396 | DO I_L=1,KLEV |
---|
| 397 | P_TAUCLD(I_L,JB) = ptaucld(K_IPLON,I_L,JB) |
---|
| 398 | ENDDO |
---|
| 399 | ENDDO |
---|
| 400 | |
---|
| 401 | ENDIF |
---|
| 402 | |
---|
| 403 | ! ------------------------------------------------------------------ |
---|
[2626] | 404 | |
---|
[1989] | 405 | IF (LHOOK) CALL DR_HOOK('RRTM_ECRT_140GP',1,ZHOOK_HANDLE) |
---|
| 406 | END SUBROUTINE RRTM_ECRT_140GP |
---|