| 1 | SUBROUTINE RADLSW & |
|---|
| 2 | & ( KIDIA, KFDIA , KLON , KLEV , KMODE, KAER,& |
|---|
| 3 | & PRII0,& |
|---|
| 4 | & PAER , PALBD , PALBP, PAPH , PAP,& |
|---|
| 5 | & PCCNL, PCCNO,& |
|---|
| 6 | & PCCO2, PCLFR , PDP , PEMIS, PEMIW , PLSM , PMU0, POZON,& |
|---|
| 7 | & PQ , PQIWP , PQLWP, PQS , PQRAIN, PRAINT,& |
|---|
| 8 | & PTH , PT , PTS , PNBAS, PNTOP,& |
|---|
| 9 | & PREF_LIQ, PREF_ICE,& |
|---|
| 10 | & PEMIT, PFCT , PFLT , PFCS , PFLS,& |
|---|
| 11 | & PFRSOD,PSUDU , PUVDF, PPARF, PPARCF, PTINCF,& |
|---|
| 12 | & PSFSWDIR, PSFSWDIF,PFSDNN,PFSDNV ,& |
|---|
| 13 | & LRDUST,PPIZA_DST,PCGA_DST,PTAUREL_DST,& |
|---|
| 14 | & PTAU_LW,& |
|---|
| 15 | & PFLUX,PFLUC,PFSDN ,PFSUP , PFSCDN , PFSCUP) |
|---|
| 16 | |
|---|
| 17 | use write_field_phy |
|---|
| 18 | ! Temporary fix waiting for cleaner interface (or not) |
|---|
| 19 | USE clesphys_mod_h, ONLY: NSW |
|---|
| 20 | |
|---|
| 21 | !**** *RADLSW* - INTERFACE TO ECMWF LW AND SW RADIATION SCHEMES |
|---|
| 22 | |
|---|
| 23 | ! PURPOSE. |
|---|
| 24 | ! -------- |
|---|
| 25 | ! CONTROLS RADIATION COMPUTATIONS |
|---|
| 26 | |
|---|
| 27 | !** INTERFACE. |
|---|
| 28 | ! ---------- |
|---|
| 29 | |
|---|
| 30 | ! EXPLICIT ARGUMENTS : |
|---|
| 31 | ! -------------------- |
|---|
| 32 | ! PAER : (KLON,6,KLEV) ; OPTICAL THICKNESS OF THE AEROSOLS |
|---|
| 33 | ! PALBD : (KLON,NSW) ; SURF. SW ALBEDO FOR DIFFUSE RADIATION |
|---|
| 34 | ! PALBP : (KLON,NSW) ; SURF. SW ALBEDO FOR PARALLEL RADIATION |
|---|
| 35 | ! PAPH : (KLON,KLEV+1) ; HALF LEVEL PRESSURE |
|---|
| 36 | ! PAP : (KLON,KLEV) ; FULL LEVEL PRESSURE |
|---|
| 37 | ! PCCNL : (KLON) ; CCN CONCENTRATION OVER LAND |
|---|
| 38 | ! PCCNO : (KLON) ; CCN CONCENTRATION OVER OCEAN |
|---|
| 39 | ! PCCO2 : ; CONCENTRATION IN CO2 (KG/KG) |
|---|
| 40 | ! PCLFR : (KLON,KLEV) ; CLOUD FRACTIONAL COVER |
|---|
| 41 | ! PDP : (KLON,KLEV) ; LAYER PRESSURE THICKNESS |
|---|
| 42 | ! PEMIS : (KLON) ; SURFACE LW EMISSIVITY |
|---|
| 43 | ! PEMIW : (KLON) ; SURFACE LW WINDOW EMISSIVITY |
|---|
| 44 | ! PLSM : (KLON) ; LAND-SEA MASK |
|---|
| 45 | ! PMU0 : (KLON) ; SOLAR ANGLE |
|---|
| 46 | ! PNBAS : (KLON) ; INDEX OF BASE OF CONVECTIVE LAYER |
|---|
| 47 | ! PNTOP : (KLON) ; INDEX OF TOP OF CONVECTIVE LAYER |
|---|
| 48 | ! POZON : (KLON,KLEV) ; OZONE AMOUNT in LAYER (KG/KG*PA) |
|---|
| 49 | ! PQ : (KLON,KLEV) ; SPECIFIC HUMIDITY KG/KG |
|---|
| 50 | ! PQIWP : (KLON,KLEV) ; SOLID WATER KG/KG |
|---|
| 51 | ! PQLWP : (KLON,KLEV) ; LIQUID WATER KG/KG |
|---|
| 52 | ! PQS : (KLON,KLEV) ; SATURATION WATER VAPOR KG/KG |
|---|
| 53 | ! PQRAIN : (KLON,KLEV) ; RAIN WATER KG/KG |
|---|
| 54 | ! PRAINT : (KLON,KLEV) ; RAIN RATE (m/s) |
|---|
| 55 | ! PTH : (KLON,KLEV+1) ; HALF LEVEL TEMPERATURE |
|---|
| 56 | ! PT : (KLON,KLEV) ; FULL LEVEL TEMPERATURE |
|---|
| 57 | ! PTS : (KLON) ; SURFACE TEMPERATURE |
|---|
| 58 | ! LDDUST ; Dust properties switch |
|---|
| 59 | ! PPIZA_DST : (KPROMA,KLEV,NSW); Single scattering albedo of dust |
|---|
| 60 | ! PCGA_DST : (KPROMA,KLEV,NSW); Assymetry factor for dust |
|---|
| 61 | ! PTAUREL_DST: (KPROMA,KLEV,NSW); Optical depth of dust relative to at 550nm |
|---|
| 62 | ! PTAU_LW (KPROMA,KLEV,NLW); LW Optical depth of aerosols |
|---|
| 63 | ! PREF_LIQ (KPROMA,KLEV) ; Liquid droplet radius (um) |
|---|
| 64 | ! PREF_ICE (KPROMA,KLEV) ; Ice crystal radius (um) |
|---|
| 65 | ! ==== OUTPUTS === |
|---|
| 66 | ! PFCT : (KLON,KLEV+1) ; CLEAR-SKY LW NET FLUXES |
|---|
| 67 | ! PFLT : (KLON,KLEV+1) ; TOTAL LW NET FLUXES |
|---|
| 68 | ! PFCS : (KLON,KLEV+1) ; CLEAR-SKY SW NET FLUXES |
|---|
| 69 | ! PFLS : (KLON,KLEV+1) ; TOTAL SW NET FLUXES |
|---|
| 70 | ! PFRSOD : (KLON) ; TOTAL-SKY SURFACE SW DOWNWARD FLUX |
|---|
| 71 | ! PEMIT : (KLON) ; SURFACE TOTAL LONGWAVE EMISSIVITY |
|---|
| 72 | ! PSUDU : (KLON) ; SOLAR RADIANCE IN SUN'S DIRECTION |
|---|
| 73 | ! PPARF : (KLON) ; PHOTOSYNTHETICALLY ACTIVE RADIATION |
|---|
| 74 | ! PUVDF : (KLON) ; UV(-B) RADIATION |
|---|
| 75 | ! PPARCF : (KLON) ; CLEAR-SKY PHOTOSYNTHETICALLY ACTIVE RADIATION |
|---|
| 76 | ! PTINCF : (KLON) ; TOA INCIDENT SOLAR RADIATION |
|---|
| 77 | ! Ajout flux LW et SW montants et descendants, et ciel clair (MPL 19.12.08) |
|---|
| 78 | ! PFLUX : (KLON,2,KLEV+1) ; LW total sky flux (1=up, 2=down) |
|---|
| 79 | ! PFLUC : (KLON,2,KLEV+1) ; LW clear sky flux (1=up, 2=down) |
|---|
| 80 | ! PFSDN(KLON,KLEV+1) ; SW total sky flux down |
|---|
| 81 | ! PFSUP(KLON,KLEV+1) ; SW total sky flux up |
|---|
| 82 | ! PFSCDN(KLON,KLEV+1) ; SW clear sky flux down |
|---|
| 83 | ! PFSCUP(KLON,KLEV+1) ; SW clear sky flux up |
|---|
| 84 | |
|---|
| 85 | |
|---|
| 86 | |
|---|
| 87 | ! IMPLICIT ARGUMENTS : NONE |
|---|
| 88 | ! -------------------- |
|---|
| 89 | |
|---|
| 90 | ! METHOD. |
|---|
| 91 | ! ------- |
|---|
| 92 | ! SEE DOCUMENTATION |
|---|
| 93 | |
|---|
| 94 | ! EXTERNALS. |
|---|
| 95 | ! ---------- |
|---|
| 96 | |
|---|
| 97 | ! REFERENCE. |
|---|
| 98 | ! ---------- |
|---|
| 99 | ! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS |
|---|
| 100 | |
|---|
| 101 | ! AUTHORS. |
|---|
| 102 | ! -------- |
|---|
| 103 | ! J.-J. MORCRETTE *ECMWF* |
|---|
| 104 | |
|---|
| 105 | ! MODIFICATIONS. |
|---|
| 106 | ! -------------- |
|---|
| 107 | ! ORIGINAL : 88-02-04 |
|---|
| 108 | ! J.-J. MORCRETTE 94-11-15 DIRECT/DIFFUSE SURFACE ALBEDO |
|---|
| 109 | ! 08/96: J.-J. Morcrette/Ph. Dandin: tests of eff. radius param. |
|---|
| 110 | ! 9909 : JJMorcrette effect.radius + inhomogeneity factors |
|---|
| 111 | ! JJMorcrette 990128 : sunshine duration |
|---|
| 112 | ! JJMorcrette : 990831 RRTM-140gp |
|---|
| 113 | ! JJMorcrette : 010112 Sun-Rikus ice particle Diameter |
|---|
| 114 | ! JJMorcrette : 010301 cleaning liq/ice cloud optical properties |
|---|
| 115 | ! JJMorcrette : 011005 CCN --> Re liquid water clouds |
|---|
| 116 | ! JJMorcrette : 011108 Safety checks |
|---|
| 117 | ! JJMorcrette : 011108 Safety checks |
|---|
| 118 | ! DJSalmond : 020211 Check before R-To-R |
|---|
| 119 | ! JJMorcrette : 020901 PAR & UV |
|---|
| 120 | ! M.Hamrud 01-Oct-2003 CY28 Cleaning |
|---|
| 121 | ! JJMorcrette : 050402 New sets of optical properties (NB: inactive) |
|---|
| 122 | ! Y.Seity 04-11-18 : add 4 arguments for AROME externalized surface |
|---|
| 123 | ! Y.Seity 05-10-10 : add 3 optional arg. for dust SW properties |
|---|
| 124 | ! JJMorcrette 20060721 PP of clear-sky PAR and TOA incident solar radiation |
|---|
| 125 | !----------------------------------------------------------------------- |
|---|
| 126 | |
|---|
| 127 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
|---|
| 128 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
|---|
| 129 | |
|---|
| 130 | USE YOMCST , ONLY : RG ,RD ,RTT ,RPI |
|---|
| 131 | !USE YOERAD , ONLY : NSW ,LRRTM ,LCCNL ,LCCNO, LDIFFC, & |
|---|
| 132 | ! NSW mis dans .def MPL 20140211 |
|---|
| 133 | USE YOERAD , ONLY : NLW, LRRTM ,LCCNL ,LCCNO, LDIFFC, & |
|---|
| 134 | & NRADIP , NRADLP , NICEOPT, NLIQOPT, NINHOM ,NLAYINH ,& |
|---|
| 135 | & RCCNLND, RCCNSEA, RLWINHF, RSWINHF, RRe2De ,& |
|---|
| 136 | & LEDBUG |
|---|
| 137 | USE YOELW , ONLY : NSIL ,NTRA ,NUA ,TSTAND ,XP |
|---|
| 138 | USE YOESW , ONLY : RYFWCA ,RYFWCB ,RYFWCC ,RYFWCD ,& |
|---|
| 139 | & RYFWCE ,RYFWCF ,REBCUA ,REBCUB ,REBCUC ,& |
|---|
| 140 | & REBCUD ,REBCUE ,REBCUF ,REBCUI ,REBCUJ ,& |
|---|
| 141 | & REBCUG ,REBCUH ,RHSAVI ,RFULIO ,RFLAA0 ,& |
|---|
| 142 | & RFLAA1 ,RFLBB0 ,RFLBB1 ,RFLBB2 ,RFLBB3 ,& |
|---|
| 143 | & RFLCC0 ,RFLCC1 ,RFLCC2 ,RFLCC3 ,RFLDD0 ,& |
|---|
| 144 | & RFLDD1 ,RFLDD2 ,RFLDD3 ,RFUETA ,RFUETB ,RFUETC ,RASWCA ,& |
|---|
| 145 | & RASWCB ,RASWCC ,RASWCD ,RASWCE ,RASWCF ,& |
|---|
| 146 | & RFUAA0 ,RFUAA1 ,RFUBB0 ,RFUBB1 ,RFUBB2 ,& |
|---|
| 147 | & RFUBB3 ,RFUCC0 ,RFUCC1 ,RFUCC2 ,RFUCC3 ,& |
|---|
| 148 | & RLILIA ,RLILIB |
|---|
| 149 | USE YOERDU , ONLY : NUAER ,NTRAER ,REPLOG ,REPSC ,REPSCW ,DIFF |
|---|
| 150 | !USE YOETHF , ONLY : RTICE |
|---|
| 151 | USE YOEPHLI , ONLY : LPHYLIN |
|---|
| 152 | USE YOERRTWN , ONLY : DELWAVE ,TOTPLNK |
|---|
| 153 | |
|---|
| 154 | USE YOMLUN_IFSAUX , ONLY : NULOUT |
|---|
| 155 | USE YOMCT3 , ONLY : NSTEP |
|---|
| 156 | |
|---|
| 157 | IMPLICIT NONE |
|---|
| 158 | |
|---|
| 159 | !!include "clesphys.h" |
|---|
| 160 | !!include "clesrrtm.h" |
|---|
| 161 | include "YOETHF.h" |
|---|
| 162 | INTEGER(KIND=JPIM),INTENT(IN) :: KLON |
|---|
| 163 | INTEGER(KIND=JPIM),INTENT(IN) :: KLEV |
|---|
| 164 | INTEGER(KIND=JPIM),INTENT(IN) :: KIDIA |
|---|
| 165 | INTEGER(KIND=JPIM),INTENT(IN) :: KFDIA |
|---|
| 166 | INTEGER(KIND=JPIM),INTENT(IN) :: KMODE |
|---|
| 167 | INTEGER(KIND=JPIM),INTENT(IN) :: KAER |
|---|
| 168 | REAL(KIND=JPRB) ,INTENT(IN) :: PRII0 |
|---|
| 169 | REAL(KIND=JPRB) ,INTENT(IN) :: PAER(KLON,6,KLEV) |
|---|
| 170 | REAL(KIND=JPRB) ,INTENT(IN) :: PALBD(KLON,NSW) |
|---|
| 171 | REAL(KIND=JPRB) ,INTENT(IN) :: PALBP(KLON,NSW) |
|---|
| 172 | REAL(KIND=JPRB) ,INTENT(IN) :: PAPH(KLON,KLEV+1) |
|---|
| 173 | REAL(KIND=JPRB) ,INTENT(IN) :: PAP(KLON,KLEV) |
|---|
| 174 | REAL(KIND=JPRB) ,INTENT(IN) :: PCCNL(KLON) |
|---|
| 175 | REAL(KIND=JPRB) ,INTENT(IN) :: PCCNO(KLON) |
|---|
| 176 | REAL(KIND=JPRB) ,INTENT(IN) :: PCCO2 |
|---|
| 177 | REAL(KIND=JPRB) ,INTENT(IN) :: PCLFR(KLON,KLEV) |
|---|
| 178 | REAL(KIND=JPRB) ,INTENT(IN) :: PDP(KLON,KLEV) |
|---|
| 179 | REAL(KIND=JPRB) ,INTENT(IN) :: PEMIS(KLON) |
|---|
| 180 | REAL(KIND=JPRB) ,INTENT(IN) :: PEMIW(KLON) |
|---|
| 181 | REAL(KIND=JPRB) ,INTENT(IN) :: PLSM(KLON) |
|---|
| 182 | REAL(KIND=JPRB) ,INTENT(IN) :: PMU0(KLON) |
|---|
| 183 | REAL(KIND=JPRB) ,INTENT(IN) :: POZON(KLON,KLEV) |
|---|
| 184 | REAL(KIND=JPRB) ,INTENT(IN) :: PQ(KLON,KLEV) |
|---|
| 185 | REAL(KIND=JPRB) ,INTENT(IN) :: PQIWP(KLON,KLEV) |
|---|
| 186 | REAL(KIND=JPRB) ,INTENT(IN) :: PQLWP(KLON,KLEV) |
|---|
| 187 | REAL(KIND=JPRB) ,INTENT(IN) :: PQS(KLON,KLEV) |
|---|
| 188 | REAL(KIND=JPRB) :: PQRAIN(KLON,KLEV) ! Argument NOT used |
|---|
| 189 | REAL(KIND=JPRB) :: PRAINT(KLON,KLEV) ! Argument NOT used |
|---|
| 190 | REAL(KIND=JPRB) ,INTENT(IN) :: PTH(KLON,KLEV+1) |
|---|
| 191 | REAL(KIND=JPRB) ,INTENT(IN) :: PT(KLON,KLEV) |
|---|
| 192 | REAL(KIND=JPRB) ,INTENT(IN) :: PTS(KLON) |
|---|
| 193 | REAL(KIND=JPRB) ,INTENT(IN) :: PNBAS(KLON) |
|---|
| 194 | REAL(KIND=JPRB) ,INTENT(IN) :: PNTOP(KLON) |
|---|
| 195 | LOGICAL ,INTENT(IN) :: LRDUST |
|---|
| 196 | REAL(KIND=JPRB) ,INTENT(IN) :: PPIZA_DST(KLON,KLEV,NSW) |
|---|
| 197 | REAL(KIND=JPRB) ,INTENT(IN) :: PCGA_DST(KLON,KLEV,NSW) |
|---|
| 198 | REAL(KIND=JPRB) ,INTENT(IN) :: PTAUREL_DST(KLON,KLEV,NSW) |
|---|
| 199 | !--C.Kleinschmitt |
|---|
| 200 | REAL(KIND=JPRB) ,INTENT(IN) :: PTAU_LW(KLON,KLEV,NLW) |
|---|
| 201 | !--end |
|---|
| 202 | REAL(KIND=JPRB) ,INTENT(IN) :: PREF_LIQ(KLON,KLEV) |
|---|
| 203 | REAL(KIND=JPRB) ,INTENT(IN) :: PREF_ICE(KLON,KLEV) |
|---|
| 204 | REAL(KIND=JPRB) ,INTENT(OUT) :: PEMIT(KLON) |
|---|
| 205 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFCT(KLON,KLEV+1) |
|---|
| 206 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFLT(KLON,KLEV+1) |
|---|
| 207 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFCS(KLON,KLEV+1) |
|---|
| 208 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFLS(KLON,KLEV+1) |
|---|
| 209 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFRSOD(KLON) |
|---|
| 210 | REAL(KIND=JPRB) ,INTENT(OUT) :: PSUDU(KLON) |
|---|
| 211 | REAL(KIND=JPRB) ,INTENT(OUT) :: PUVDF(KLON) |
|---|
| 212 | REAL(KIND=JPRB) ,INTENT(OUT) :: PPARF(KLON) |
|---|
| 213 | REAL(KIND=JPRB) ,INTENT(OUT) :: PPARCF(KLON), PTINCF(KLON) |
|---|
| 214 | REAL(KIND=JPRB) ,INTENT(OUT) :: PSFSWDIR(KLON,NSW) |
|---|
| 215 | REAL(KIND=JPRB) ,INTENT(OUT) :: PSFSWDIF(KLON,NSW) |
|---|
| 216 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSDNN(KLON) |
|---|
| 217 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSDNV(KLON) |
|---|
| 218 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFLUX(KLON,2,KLEV+1) ! LW total sky flux (1=up, 2=down) |
|---|
| 219 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFLUC(KLON,2,KLEV+1) ! LW clear sky flux (1=up, 2=down) |
|---|
| 220 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSDN(KLON,KLEV+1) ! SW total sky flux down |
|---|
| 221 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSUP(KLON,KLEV+1) ! SW total sky flux up |
|---|
| 222 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSCDN(KLON,KLEV+1) ! SW clear sky flux down |
|---|
| 223 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSCUP(KLON,KLEV+1) ! SW clear sky flux up |
|---|
| 224 | |
|---|
| 225 | |
|---|
| 226 | ! ----------------------------------------------------------------- |
|---|
| 227 | |
|---|
| 228 | !* 0.1 ARGUMENTS. |
|---|
| 229 | ! ---------- |
|---|
| 230 | ! ==== COMPUTED IN RADLSW === |
|---|
| 231 | ! ----------------------------------------------------------------- |
|---|
| 232 | |
|---|
| 233 | !* 0.2 LOCAL ARRAYS. |
|---|
| 234 | ! ------------- |
|---|
| 235 | ! ----------------------------------------------------------------- |
|---|
| 236 | |
|---|
| 237 | !-- ARRAYS FOR LOCAL VARIABLES ----------------------------------------- |
|---|
| 238 | |
|---|
| 239 | INTEGER(KIND=JPIM) :: IBAS(KLON) , ITOP(KLON) |
|---|
| 240 | |
|---|
| 241 | REAL(KIND=JPRB) ::& |
|---|
| 242 | & ZALBD(KLON,NSW) , ZALBP(KLON,NSW)& |
|---|
| 243 | & , ZCG(KLON,NSW,KLEV) , ZOMEGA(KLON,NSW,KLEV)& |
|---|
| 244 | & , ZTAU (KLON,NSW,KLEV) & |
|---|
| 245 | & , ZTAUCLD(KLON,KLEV,16), ZTCLEAR(KLON) |
|---|
| 246 | REAL(KIND=JPRB) ::& |
|---|
| 247 | & ZCLDLD(KLON,KLEV) , ZCLDLU(KLON,KLEV)& |
|---|
| 248 | & , ZCLDSW(KLON,KLEV) , ZCLD0(KLON,KLEV)& |
|---|
| 249 | & , ZDT0(KLON) & |
|---|
| 250 | & , ZEMIS(KLON) , ZEMIW(KLON)& |
|---|
| 251 | & , ZFIWP(KLON) , ZFLWP(KLON) , ZFRWP(KLON)& |
|---|
| 252 | & , ZIWC(KLON) , ZLWC(KLON)& |
|---|
| 253 | !cc , ZRWC(KLON) |
|---|
| 254 | & , ZMU0(KLON) , ZOZ(KLON,KLEV) , ZOZN(KLON,KLEV)& |
|---|
| 255 | & , ZPMB(KLON,KLEV+1) , ZPSOL(KLON)& |
|---|
| 256 | & , ZTAVE (KLON,KLEV) , ZTL(KLON,KLEV+1)& |
|---|
| 257 | & , ZVIEW(KLON) |
|---|
| 258 | REAL(KIND=JPRB) ::& |
|---|
| 259 | & ZFCDWN(KLON,KLEV+1), ZFCUP(KLON,KLEV+1)& |
|---|
| 260 | & , ZFSDWN(KLON,KLEV+1), ZFSUP(KLON,KLEV+1)& |
|---|
| 261 | & , ZFSUPN(KLON) , ZFSUPV(KLON)& |
|---|
| 262 | & , ZFCUPN(KLON) , ZFCUPV(KLON)& |
|---|
| 263 | & , ZFSDNN(KLON) , ZFSDNV(KLON)& |
|---|
| 264 | & , ZFCDNN(KLON) , ZFCDNV(KLON)& |
|---|
| 265 | & , ZDIRFS(KLON,NSW) , ZDIFFS(KLON,NSW) |
|---|
| 266 | REAL(KIND=JPRB) ::& |
|---|
| 267 | & ZALFICE(KLON) , ZGAMICE(KLON) , ZBICE(KLON) , ZDESR(KLON)& |
|---|
| 268 | & , ZRADIP(KLON) , ZRADLP(KLON) & |
|---|
| 269 | !cc , ZRADRD(KLON) |
|---|
| 270 | & , ZRAINT(KLON) , ZRES(KLON)& |
|---|
| 271 | & , ZTICE(KLON) , ZEMIT(KLON), ZBICFU(KLON)& |
|---|
| 272 | & , ZKICFU(KLON) |
|---|
| 273 | REAL(KIND=JPRB) :: ZSUDU(KLON) , ZPARF(KLON) , ZUVDF(KLON), ZPARCF(KLON) |
|---|
| 274 | INTEGER(KIND=JPIM) :: IKL, JK, JKL, JKLP1, JKP1, JL, JNU, JRTM, JSW, INDLAY |
|---|
| 275 | |
|---|
| 276 | REAL(KIND=JPRB) :: ZASYMX, ZDIFFD, ZGI, ZGL, ZGR, ZIWGKG, ZLWGKG,& |
|---|
| 277 | & ZMSAID, ZMSAIU, ZMSALD, ZMSALU, ZRSAIA, ZRSAID, ZRSAIE, ZRSAIF, ZRSAIG, ZRSALD, & |
|---|
| 278 | & ZMULTI, ZMULTL, ZOI , ZOL, & |
|---|
| 279 | & ZOMGMX, ZOR, ZRMUZ, ZRWGKG, ZTAUD, ZTAUMX, ZTEMPC, & |
|---|
| 280 | & ZTOI, ZTOL, ZTOR, ZZFIWP, ZZFLWP, ZDPOG, ZPODT |
|---|
| 281 | |
|---|
| 282 | REAL(KIND=JPRB) :: ZALND, ZASEA, ZD, ZDEN, ZNTOT, ZNUM, ZRATIO, Z1RADI, & |
|---|
| 283 | & Z1RADL, ZBETAI, ZOMGI, ZOMGP, ZFDEL, ZTCELS, ZFSR, ZAIWC, & |
|---|
| 284 | & ZBIWC, ZTBLAY, ZADDPLK, ZPLANCK, ZEXTCF, Z1MOMG, & |
|---|
| 285 | & ZDefRe, ZRefDe, ZVI , ZMABSD |
|---|
| 286 | |
|---|
| 287 | !REAL(KIND=JPRB) :: ZAVDP(KLON), ZAVTO(KLON), ZSQTO(KLON) |
|---|
| 288 | REAL(KIND=JPRB) :: ZAVTO(KLON), ZSQTO(KLON) |
|---|
| 289 | REAL(KIND=JPRB) :: ZSQUAR(KLON,KLEV), ZVARIA(KLON,KLEV) |
|---|
| 290 | INTEGER(KIND=JPIM) :: IKI, JKI, JEXPLR, JXPLDN |
|---|
| 291 | LOGICAL :: LLDEBUG |
|---|
| 292 | |
|---|
| 293 | |
|---|
| 294 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
|---|
| 295 | |
|---|
| 296 | #include "lw.intfb.h" |
|---|
| 297 | #include "rrtm_rrtm_140gp.intfb.h" |
|---|
| 298 | #include "sw.intfb.h" |
|---|
| 299 | |
|---|
| 300 | ! ----------------------------------------------------------------- |
|---|
| 301 | |
|---|
| 302 | !* 1. SET-UP INPUT QUANTITIES FOR RADIATION |
|---|
| 303 | ! ------------------------------------- |
|---|
| 304 | |
|---|
| 305 | IF (LHOOK) CALL DR_HOOK('RADLSW',0,ZHOOK_HANDLE) |
|---|
| 306 | |
|---|
| 307 | LLDEBUG=.FALSE. |
|---|
| 308 | ZRefDe = RRe2De |
|---|
| 309 | ZDefRe = 1.0_JPRB / ZRefDe |
|---|
| 310 | |
|---|
| 311 | DO JL = KIDIA,KFDIA |
|---|
| 312 | ZFCUP(JL,KLEV+1) = 0.0_JPRB |
|---|
| 313 | ZFCDWN(JL,KLEV+1) = REPLOG |
|---|
| 314 | ZFSUP(JL,KLEV+1) = 0.0_JPRB |
|---|
| 315 | ZFSDWN(JL,KLEV+1) = REPLOG |
|---|
| 316 | PFLUX(JL,1,KLEV+1) = 0.0_JPRB |
|---|
| 317 | PFLUX(JL,2,KLEV+1) = 0.0_JPRB |
|---|
| 318 | PFLUC(JL,1,KLEV+1) = 0.0_JPRB |
|---|
| 319 | PFLUC(JL,2,KLEV+1) = 0.0_JPRB |
|---|
| 320 | ZFSDNN(JL) = 0.0_JPRB |
|---|
| 321 | ZFSDNV(JL) = 0.0_JPRB |
|---|
| 322 | ZFCDNN(JL) = 0.0_JPRB |
|---|
| 323 | ZFCDNV(JL) = 0.0_JPRB |
|---|
| 324 | ZFSUPN(JL) = 0.0_JPRB |
|---|
| 325 | ZFSUPV(JL) = 0.0_JPRB |
|---|
| 326 | ZFCUPN(JL) = 0.0_JPRB |
|---|
| 327 | ZFCUPV(JL) = 0.0_JPRB |
|---|
| 328 | ZPSOL(JL) = PAPH(JL,KLEV+1) |
|---|
| 329 | ZPMB(JL,1) = ZPSOL(JL) / 100.0_JPRB |
|---|
| 330 | ZDT0(JL) = PTS(JL) - PTH(JL,KLEV+1) |
|---|
| 331 | PSUDU(JL) = 0.0_JPRB |
|---|
| 332 | PPARF(JL) = 0.0_JPRB |
|---|
| 333 | PPARCF(JL)= 0.0_JPRB |
|---|
| 334 | PUVDF(JL) = 0.0_JPRB |
|---|
| 335 | PSFSWDIR(JL,:)=0.0_JPRB |
|---|
| 336 | PSFSWDIF(JL,:)=0.0_JPRB |
|---|
| 337 | IBAS(JL) = INT ( 0.01_JPRB + PNBAS(JL) ) |
|---|
| 338 | ITOP(JL) = INT ( 0.01_JPRB + PNTOP(JL) ) |
|---|
| 339 | ENDDO |
|---|
| 340 | |
|---|
| 341 | !* 1.1 INITIALIZE VARIOUS FIELDS |
|---|
| 342 | ! ------------------------- |
|---|
| 343 | |
|---|
| 344 | DO JSW=1,NSW |
|---|
| 345 | DO JL = KIDIA,KFDIA |
|---|
| 346 | ZALBD(JL,JSW)=PALBD(JL,JSW) |
|---|
| 347 | ZALBP(JL,JSW)=PALBP(JL,JSW) |
|---|
| 348 | ENDDO |
|---|
| 349 | ENDDO |
|---|
| 350 | DO JL = KIDIA,KFDIA |
|---|
| 351 | ZEMIS(JL) =PEMIS(JL) |
|---|
| 352 | ZEMIW(JL) =PEMIW(JL) |
|---|
| 353 | ZMU0(JL) =PMU0(JL) |
|---|
| 354 | ENDDO |
|---|
| 355 | |
|---|
| 356 | DO JK = 1 , KLEV |
|---|
| 357 | JKP1 = JK + 1 |
|---|
| 358 | JKL = KLEV+ 1 - JK |
|---|
| 359 | JKLP1 = JKL + 1 |
|---|
| 360 | DO JL = KIDIA,KFDIA |
|---|
| 361 | ZPMB(JL,JK+1)=PAPH(JL,JKL)/100.0_JPRB |
|---|
| 362 | |
|---|
| 363 | !-- ZOZ in cm.atm for SW scheme |
|---|
| 364 | ZOZ(JL,JK) = POZON(JL,JKL) * 46.6968_JPRB / RG |
|---|
| 365 | |
|---|
| 366 | ZCLD0(JL,JK) = 0.0_JPRB |
|---|
| 367 | ZFCUP(JL,JK) = 0.0_JPRB |
|---|
| 368 | ZFCDWN(JL,JK) = 0.0_JPRB |
|---|
| 369 | ZFSUP(JL,JK) = 0.0_JPRB |
|---|
| 370 | ZFSDWN(JL,JK) = 0.0_JPRB |
|---|
| 371 | PFLUX(JL,1,JK) = 0.0_JPRB |
|---|
| 372 | PFLUX(JL,2,JK) = 0.0_JPRB |
|---|
| 373 | PFLUC(JL,1,JK) = 0.0_JPRB |
|---|
| 374 | PFLUC(JL,2,JK) = 0.0_JPRB |
|---|
| 375 | ENDDO |
|---|
| 376 | ENDDO |
|---|
| 377 | |
|---|
| 378 | DO JK=1,KLEV |
|---|
| 379 | JKL=KLEV+1-JK |
|---|
| 380 | JKLP1=JKL+1 |
|---|
| 381 | DO JL=KIDIA,KFDIA |
|---|
| 382 | ZTL(JL,JK)=PTH(JL,JKLP1) |
|---|
| 383 | ZTAVE(JL,JK)=PT(JL,JKL) |
|---|
| 384 | ENDDO |
|---|
| 385 | ENDDO |
|---|
| 386 | DO JL=KIDIA,KFDIA |
|---|
| 387 | ZTL(JL,KLEV+1)= PTH(JL,1) |
|---|
| 388 | ZPMB(JL,KLEV+1) = PAPH(JL,1)/100.0_JPRB |
|---|
| 389 | ENDDO |
|---|
| 390 | !*** |
|---|
| 391 | |
|---|
| 392 | ! ------------------------------------------------------------------ |
|---|
| 393 | |
|---|
| 394 | !* 2. CLOUD AND AEROSOL PARAMETERS |
|---|
| 395 | ! ---------------------------- |
|---|
| 396 | |
|---|
| 397 | DO JK = 1 , KLEV |
|---|
| 398 | IKL = KLEV + 1 - JK |
|---|
| 399 | |
|---|
| 400 | ! 2.1 INITIALIZE OPTICAL PROPERTIES TO CLEAR SKY VALUES |
|---|
| 401 | ! ------------------------------------------------- |
|---|
| 402 | |
|---|
| 403 | DO JSW = 1,NSW |
|---|
| 404 | DO JL = KIDIA,KFDIA |
|---|
| 405 | ZTAU(JL,JSW,JK) = 0.0_JPRB |
|---|
| 406 | ZOMEGA(JL,JSW,JK)= 1.0_JPRB |
|---|
| 407 | ZCG(JL,JSW,JK) = 0.0_JPRB |
|---|
| 408 | ENDDO |
|---|
| 409 | ENDDO |
|---|
| 410 | DO JL = KIDIA,KFDIA |
|---|
| 411 | ZCLDSW(JL,JK) = 0.0_JPRB |
|---|
| 412 | ZCLDLD(JL,JK) = 0.0_JPRB |
|---|
| 413 | ZCLDLU(JL,JK) = 0.0_JPRB |
|---|
| 414 | ENDDO |
|---|
| 415 | |
|---|
| 416 | ! 2.2 CLOUD ICE AND LIQUID CONTENT AND PATH |
|---|
| 417 | ! ------------------------------------- |
|---|
| 418 | |
|---|
| 419 | DO JL = KIDIA,KFDIA |
|---|
| 420 | |
|---|
| 421 | ! --- LIQUID WATER CONTENT (g.m-3) AND LIQUID WATER PATH (g.m-2) |
|---|
| 422 | IF (PCLFR(JL,IKL) > REPSC ) THEN |
|---|
| 423 | ZLWGKG=MAX(PQLWP(JL,IKL)*1000.0_JPRB,0.0_JPRB) |
|---|
| 424 | ZIWGKG=MAX(PQIWP(JL,IKL)*1000.0_JPRB,0.0_JPRB) |
|---|
| 425 | ZLWGKG=ZLWGKG/PCLFR(JL,IKL) |
|---|
| 426 | ZIWGKG=ZIWGKG/PCLFR(JL,IKL) |
|---|
| 427 | ELSE |
|---|
| 428 | ZLWGKG=0.0_JPRB |
|---|
| 429 | ZIWGKG=0.0_JPRB |
|---|
| 430 | ENDIF |
|---|
| 431 | ZRWGKG=0.0_JPRB |
|---|
| 432 | ZRAINT(JL)=0.0_JPRB |
|---|
| 433 | |
|---|
| 434 | ! --- RAIN LIQUID WATER CONTENT (g.m-3) AND LIQUID WATER PATH (g.m-2) |
|---|
| 435 | ! IF (PRAINT(JL,IKL) >= REPSCW) THEN |
|---|
| 436 | ! ZRWGKG=MAX(PQRAIN(JL,IKL)*1000., 0.0) |
|---|
| 437 | ! ZRAINT(JL)=PRAINT(JL,IKL)*3600.*1000. |
|---|
| 438 | !- no radiative effect of rain (for the moment) |
|---|
| 439 | ! ZRWGKG=0. |
|---|
| 440 | ! ZRAINT(JL)=0. |
|---|
| 441 | ! =========================================================== |
|---|
| 442 | |
|---|
| 443 | ! Modifications Martin et al. |
|---|
| 444 | ! ELSE |
|---|
| 445 | ! ENDIF |
|---|
| 446 | ZDPOG=PDP(JL,IKL)/RG |
|---|
| 447 | ZFLWP(JL)= ZLWGKG*ZDPOG |
|---|
| 448 | ZFIWP(JL)= ZIWGKG*ZDPOG |
|---|
| 449 | ZFRWP(JL)= ZRWGKG*ZDPOG |
|---|
| 450 | ZPODT=PAP(JL,IKL)/(RD*PT(JL,IKL)) |
|---|
| 451 | ZLWC(JL)=ZLWGKG*ZPODT |
|---|
| 452 | ZIWC(JL)=ZIWGKG*ZPODT |
|---|
| 453 | ! ZRWC(JL)=ZRWGKG*ZPODT |
|---|
| 454 | |
|---|
| 455 | ENDDO |
|---|
| 456 | DO JL = KIDIA,KFDIA |
|---|
| 457 | ! --- EFFECTIVE RADIUS FOR WATER, ICE AND RAIN PARTICLES |
|---|
| 458 | |
|---|
| 459 | ! very old parametrization as f(pressure) |
|---|
| 460 | |
|---|
| 461 | IF (NRADLP == 0) THEN |
|---|
| 462 | !-- very old parametrization as f(pressure) ERA-15 |
|---|
| 463 | ZRADLP(JL)=10.0_JPRB + (100000.0_JPRB-PAP(JL,IKL))*3.5_JPRB |
|---|
| 464 | |
|---|
| 465 | ELSEIF (NRADLP == 1) THEN |
|---|
| 466 | ! simple distinction between land (10) and ocean (13) Zhang and Rossow |
|---|
| 467 | IF (PLSM(JL) < 0.5_JPRB) THEN |
|---|
| 468 | ZRADLP(JL)=13.0_JPRB |
|---|
| 469 | ELSE |
|---|
| 470 | ZRADLP(JL)=10.0_JPRB |
|---|
| 471 | ENDIF |
|---|
| 472 | |
|---|
| 473 | ELSEIF (NRADLP == 2) THEN |
|---|
| 474 | !-- based on Martin et al., 1994, JAS |
|---|
| 475 | IF (PLSM(JL) < 0.5_JPRB) THEN |
|---|
| 476 | IF (LCCNO) THEN |
|---|
| 477 | ! ZASEA=50.0_JPRB |
|---|
| 478 | ZASEA=PCCNO(JL) |
|---|
| 479 | ELSE |
|---|
| 480 | ZASEA=RCCNSEA |
|---|
| 481 | ENDIF |
|---|
| 482 | ZD=0.33_JPRB |
|---|
| 483 | ZNTOT=-1.15E-03_JPRB*ZASEA*ZASEA+0.963_JPRB*ZASEA+5.30_JPRB |
|---|
| 484 | ELSE |
|---|
| 485 | IF (LCCNL) THEN |
|---|
| 486 | ! ZALND=900.0_JPRB |
|---|
| 487 | ZALND=PCCNL(JL) |
|---|
| 488 | ELSE |
|---|
| 489 | ZALND=RCCNLND |
|---|
| 490 | ENDIF |
|---|
| 491 | ZD=0.43_JPRB |
|---|
| 492 | ZNTOT=-2.10E-04_JPRB*ZALND*ZALND+0.568_JPRB*ZALND-27.9_JPRB |
|---|
| 493 | ENDIF |
|---|
| 494 | ZNUM=3.0_JPRB*ZLWC(JL)*(1.0_JPRB+3.0_JPRB*ZD*ZD)**2 |
|---|
| 495 | ZDEN=4.0_JPRB*RPI*ZNTOT*(1.0_JPRB+ZD*ZD)**3 |
|---|
| 496 | IF((ZNUM/ZDEN) > REPLOG)THEN |
|---|
| 497 | ZRADLP(JL)=100.0_JPRB*EXP(0.333_JPRB*LOG(ZNUM/ZDEN)) |
|---|
| 498 | ZRADLP(JL)=MAX(ZRADLP(JL), 4.0_JPRB) |
|---|
| 499 | ZRADLP(JL)=MIN(ZRADLP(JL),16.0_JPRB) |
|---|
| 500 | ELSE |
|---|
| 501 | ZRADLP(JL)=4.0_JPRB |
|---|
| 502 | ENDIF |
|---|
| 503 | |
|---|
| 504 | ELSEIF (NRADLP == 3) THEN |
|---|
| 505 | ! one uses the cloud droplet radius from newmicro |
|---|
| 506 | ! IKL or JK ?? - I think IKL but needs to be verified > ref_liq_i |
|---|
| 507 | ! (inverted) is used in the call of RECMWF_AERO in radlwsw_m.F90, |
|---|
| 508 | ! so everything is fine - JBM 6/2019 |
|---|
| 509 | ZRADLP(JL)=PREF_LIQ(JL,IKL) |
|---|
| 510 | ENDIF |
|---|
| 511 | |
|---|
| 512 | ! =========================================================== |
|---|
| 513 | ! ___________________________________________________________ |
|---|
| 514 | |
|---|
| 515 | ! rain drop from : unused as ZRAINT is 0. |
|---|
| 516 | ! ZRADRD(JL)=500.0_JPRB*ZRAINT(JL)**0.22_JPRB |
|---|
| 517 | ! IF (ZFLWP(JL).GT.0.) THEN |
|---|
| 518 | ! ZRADRD(JL)=ZRADLP(JL)+ZRADRD(JL) |
|---|
| 519 | ! ENDIF |
|---|
| 520 | |
|---|
| 521 | ENDDO |
|---|
| 522 | DO JL = KIDIA,KFDIA |
|---|
| 523 | |
|---|
| 524 | ! diagnosing the ice particle effective radius/diameter |
|---|
| 525 | |
|---|
| 526 | !- ice particle effective radius =f(T) from Liou and Ou (1994) |
|---|
| 527 | |
|---|
| 528 | IF (PT(JL,IKL) < RTICE) THEN |
|---|
| 529 | ZTEMPC=PT(JL,IKL)-RTT |
|---|
| 530 | ELSE |
|---|
| 531 | ZTEMPC=RTICE-RTT |
|---|
| 532 | ENDIF |
|---|
| 533 | ZRADIP(JL)=326.3_JPRB+ZTEMPC*(12.42_JPRB + ZTEMPC*(0.197_JPRB + ZTEMPC*& |
|---|
| 534 | & 0.0012_JPRB)) |
|---|
| 535 | |
|---|
| 536 | IF (NRADIP == 0) THEN |
|---|
| 537 | !-- fixed 40 micron effective radius |
|---|
| 538 | ZRADIP(JL)= 40.0_JPRB |
|---|
| 539 | ZDESR(JL) = ZDefRe * ZRADIP(JL) |
|---|
| 540 | |
|---|
| 541 | ELSEIF (NRADIP == 1) THEN |
|---|
| 542 | |
|---|
| 543 | !-- old formulation based on Liou & Ou (1994) temperature (40-130microns) |
|---|
| 544 | ZRADIP(JL)=MAX(ZRADIP(JL),40.0_JPRB) |
|---|
| 545 | ZDESR(JL) = ZDefRe * ZRADIP(JL) |
|---|
| 546 | |
|---|
| 547 | ELSEIF (NRADIP == 2) THEN |
|---|
| 548 | !-- formulation following Jakob, Klein modifications to ice content |
|---|
| 549 | ZRADIP(JL)=MAX(ZRADIP(JL),30.0_JPRB) |
|---|
| 550 | ZRADIP(JL)=MIN(ZRADIP(JL),60.0_JPRB) |
|---|
| 551 | ZDESR(JL)= ZDefRe * ZRADIP(JL) |
|---|
| 552 | |
|---|
| 553 | ELSEIF (NRADIP == 3 ) THEN |
|---|
| 554 | |
|---|
| 555 | !- ice particle effective radius =f(T,IWC) from Sun and Rikus (1999) |
|---|
| 556 | ! revised by Sun (2001) |
|---|
| 557 | IF (ZIWC(JL) > 0.0_JPRB ) THEN |
|---|
| 558 | ZTEMPC = PT(JL,IKL)-83.15_JPRB |
|---|
| 559 | ZTCELS = PT(JL,IKL)-RTT |
|---|
| 560 | ZFSR = 1.2351_JPRB +0.0105_JPRB * ZTCELS |
|---|
| 561 | ! Sun, 2001 (corrected from Sun & Rikus, 1999) |
|---|
| 562 | ZAIWC = 45.8966_JPRB * ZIWC(JL)**0.2214_JPRB |
|---|
| 563 | ZBIWC = 0.7957_JPRB * ZIWC(JL)**0.2535_JPRB |
|---|
| 564 | ZDESR(JL) = ZFSR * (ZAIWC + ZBIWC*ZTEMPC) |
|---|
| 565 | !-new ZDESR(JL) = MIN ( MAX( ZDESR(JL), 30.0_JPRB), 155.0_JPRB) |
|---|
| 566 | ZDESR(JL) = MIN ( MAX( ZDESR(JL), 45.0_JPRB), 350.0_JPRB) |
|---|
| 567 | ZRADIP(JL)= ZRefDe * ZDESR(JL) |
|---|
| 568 | ELSE |
|---|
| 569 | ! ZDESR(JL) = 92.5_JPRB |
|---|
| 570 | ZDESR(JL) = 80.0_JPRB |
|---|
| 571 | ZRADIP(JL)= ZRefDe * ZDESR(JL) |
|---|
| 572 | ENDIF |
|---|
| 573 | |
|---|
| 574 | ELSEIF (NRADIP == 4 ) THEN |
|---|
| 575 | ! one uses the cloud droplet radius from newmicro |
|---|
| 576 | ! IKL or JK ?? - I think IKL but needs to be verified |
|---|
| 577 | ZRADIP(JL)=PREF_ICE(JL,IKL) |
|---|
| 578 | ENDIF |
|---|
| 579 | |
|---|
| 580 | ENDDO |
|---|
| 581 | |
|---|
| 582 | ! 2.3 CLOUD SHORTWAVE OPTICAL PROPERTIES |
|---|
| 583 | ! ---------------------------------- |
|---|
| 584 | |
|---|
| 585 | ! ------------------------- |
|---|
| 586 | ! --+ SW OPTICAL PARAMETERS + Water clouds after Fouquart (1987) |
|---|
| 587 | ! ------------------------- Ice clouds (Ebert, Curry, 1992) |
|---|
| 588 | |
|---|
| 589 | DO JSW=1,NSW |
|---|
| 590 | DO JL = KIDIA,KFDIA |
|---|
| 591 | ZTOL=0.0_JPRB |
|---|
| 592 | ZGL =0.0_JPRB |
|---|
| 593 | ZOL =0.0_JPRB |
|---|
| 594 | ZTOI=0.0_JPRB |
|---|
| 595 | ZGI =0.0_JPRB |
|---|
| 596 | ZOI =0.0_JPRB |
|---|
| 597 | ZTOR=0.0_JPRB |
|---|
| 598 | ZGR =0.0_JPRB |
|---|
| 599 | ZOR =0.0_JPRB |
|---|
| 600 | IF (ZFLWP(JL)+ZFIWP(JL)+ZFRWP(JL) > 2.0_JPRB * REPSCW ) THEN |
|---|
| 601 | IF (ZFLWP(JL) >= REPSCW ) THEN |
|---|
| 602 | IF (NLIQOPT /= 0 ) THEN |
|---|
| 603 | !-- SW: Slingo, 1989 |
|---|
| 604 | ZTOL = ZFLWP(JL)*(RASWCA(JSW)+RASWCB(JSW)/ZRADLP(JL)) |
|---|
| 605 | ZGL = RASWCE(JSW)+RASWCF(JSW)*ZRADLP(JL) |
|---|
| 606 | ZOL = 1. - RASWCC(JSW)-RASWCD(JSW)*ZRADLP(JL) |
|---|
| 607 | ELSE |
|---|
| 608 | !-- SW: Fouquart, 1991 |
|---|
| 609 | ZTOL = ZFLWP(JL)*(RYFWCA(JSW)+RYFWCB(JSW)/ZRADLP(JL)) |
|---|
| 610 | ZGL = RYFWCF(JSW) |
|---|
| 611 | ! ZOL = RYFWCC(JSW)-RYFWCD(JSW)*EXP(-RYFWCE(JSW)*ZTOL) |
|---|
| 612 | !-- NB: RSWINHF is there simply for making the CY29R2 branch bit compatible with |
|---|
| 613 | ! the previous. Should be cleaned when RRTM_SW becomes active |
|---|
| 614 | ZOL = RYFWCC(JSW)-RYFWCD(JSW)*EXP(-RYFWCE(JSW)*ZTOL*RSWINHF) |
|---|
| 615 | ENDIF |
|---|
| 616 | ENDIF |
|---|
| 617 | |
|---|
| 618 | IF (ZFIWP(JL) >= REPSCW ) THEN |
|---|
| 619 | IF (NICEOPT <= 1) THEN |
|---|
| 620 | !-- SW: Ebert-Curry |
|---|
| 621 | ZTOI = ZFIWP(JL)*(REBCUA(JSW)+REBCUB(JSW)/ZRADIP(JL)) |
|---|
| 622 | ZGI = REBCUE(JSW)+REBCUF(JSW)*ZRADIP(JL) |
|---|
| 623 | ZOI = 1.0_JPRB - REBCUC(JSW)-REBCUD(JSW)*ZRADIP(JL) |
|---|
| 624 | |
|---|
| 625 | ELSEIF (NICEOPT == 2) THEN |
|---|
| 626 | !-- SW: Fu-Liou 1993 |
|---|
| 627 | Z1RADI = 1.0_JPRB / ZDESR(JL) |
|---|
| 628 | ZBETAI = RFLAA0(JSW)+Z1RADI* RFLAA1(JSW) |
|---|
| 629 | ZTOI = ZFIWP(JL) * ZBETAI |
|---|
| 630 | ZOMGI= RFLBB0(JSW)+ZRADIP(JL)*(RFLBB1(JSW) + ZRADIP(JL) & |
|---|
| 631 | & *(RFLBB2(JSW)+ZRADIP(JL)* RFLBB3(JSW) )) |
|---|
| 632 | ZOI = 1.0_JPRB - ZOMGI |
|---|
| 633 | ZOMGP= RFLCC0(JSW)+ZRADIP(JL)*(RFLCC1(JSW) + ZRADIP(JL) & |
|---|
| 634 | & *(RFLCC2(JSW)+ZRADIP(JL)* RFLCC3(JSW) )) |
|---|
| 635 | ZFDEL= RFLDD0(JSW)+ZRADIP(JL)*(RFLDD1(JSW) + ZRADIP(JL) & |
|---|
| 636 | & *(RFLDD2(JSW)+ZRADIP(JL)* RFLDD3(JSW) )) |
|---|
| 637 | ZGI = ((1.0_JPRB -ZFDEL)*ZOMGP + ZFDEL*3.0_JPRB) / 3.0_JPRB |
|---|
| 638 | |
|---|
| 639 | ELSEIF (NICEOPT == 3) THEN |
|---|
| 640 | !-- SW: Fu 1996 |
|---|
| 641 | Z1RADI = 1.0_JPRB / ZDESR(JL) |
|---|
| 642 | ZBETAI = RFUAA0(JSW)+Z1RADI* RFUAA1(JSW) |
|---|
| 643 | ZTOI = ZFIWP(JL) * ZBETAI |
|---|
| 644 | ZOMGI= RFUBB0(JSW)+ZDESR(JL)*(RFUBB1(JSW) + ZDESR(JL) & |
|---|
| 645 | & *(RFUBB2(JSW)+ZDESR(JL)* RFUBB3(JSW) )) |
|---|
| 646 | ZOI = 1.0_JPRB - ZOMGI |
|---|
| 647 | ZGI = RFUCC0(JSW)+ZDESR(JL)*(RFUCC1(JSW) + ZDESR(JL) & |
|---|
| 648 | & *(RFUCC2(JSW)+ZDESR(JL)* RFUCC3(JSW) )) |
|---|
| 649 | ZGI = MIN(1.0_JPRB, ZGI) |
|---|
| 650 | |
|---|
| 651 | ENDIF |
|---|
| 652 | ENDIF |
|---|
| 653 | |
|---|
| 654 | ! IF (ZFRWP(JL) >= REPSCW ) THEN |
|---|
| 655 | ! ZTOR= ZFRWP(JL)*0.003_JPRB * ZRAINT(JL)**(-0.22_JPRB) |
|---|
| 656 | ! ZOR = 1.0_JPRB - RROMA(JSW)*ZRAINT(JL)**RROMB(JSW) |
|---|
| 657 | ! ZGR = RRASY(JSW) |
|---|
| 658 | ! ENDIF |
|---|
| 659 | |
|---|
| 660 | ! - MIX of WATER and ICE CLOUDS |
|---|
| 661 | ZTAUMX= ZTOL + ZTOI + ZTOR |
|---|
| 662 | ZOMGMX= ZTOL*ZOL + ZTOI*ZOI + ZTOR*ZOR |
|---|
| 663 | ZASYMX= ZTOL*ZOL*ZGL + ZTOI*ZOI*ZGI + ZTOR*ZOR*ZGR |
|---|
| 664 | |
|---|
| 665 | ZASYMX= ZASYMX/ZOMGMX |
|---|
| 666 | ZOMGMX= ZOMGMX/ZTAUMX |
|---|
| 667 | |
|---|
| 668 | ! --- SW FINAL CLOUD OPTICAL PARAMETERS |
|---|
| 669 | |
|---|
| 670 | ZCLDSW(JL,JK) = PCLFR(JL,IKL) |
|---|
| 671 | ZTAU(JL,JSW,JK) = ZTAUMX |
|---|
| 672 | ZOMEGA(JL,JSW,JK)= ZOMGMX |
|---|
| 673 | ZCG(JL,JSW,JK) = ZASYMX |
|---|
| 674 | ENDIF |
|---|
| 675 | ENDDO |
|---|
| 676 | ENDDO |
|---|
| 677 | |
|---|
| 678 | IF(LLDEBUG) THEN |
|---|
| 679 | call writefield_phy("radlsw_ztau",ztau(:,1,:),klev) |
|---|
| 680 | ENDIF |
|---|
| 681 | |
|---|
| 682 | ! 2.4 CLOUD LONGWAVE OPTICAL PROPERTIES FOR EC-OPE |
|---|
| 683 | ! -------------------------------------------- |
|---|
| 684 | |
|---|
| 685 | ! ------------------------- |
|---|
| 686 | ! --+ LW OPTICAL PARAMETERS + Water (and Ice) from Smith and Shi (1992) |
|---|
| 687 | ! ------------------------- Ice clouds (Ebert, Curry, 1992) |
|---|
| 688 | |
|---|
| 689 | IF (.NOT.LRRTM) THEN |
|---|
| 690 | |
|---|
| 691 | DO JL = KIDIA,KFDIA |
|---|
| 692 | ZALFICE(JL)=0.0_JPRB |
|---|
| 693 | ZGAMICE(JL)=0.0_JPRB |
|---|
| 694 | ZBICE(JL)=0.0_JPRB |
|---|
| 695 | ZTICE(JL)=(PT(JL,IKL)-TSTAND)/TSTAND |
|---|
| 696 | IF (NICEOPT == 1) THEN |
|---|
| 697 | ZBICFU(JL)=1.0_JPRB |
|---|
| 698 | ELSE |
|---|
| 699 | ZBICFU(JL)=0.0_JPRB |
|---|
| 700 | ENDIF |
|---|
| 701 | ZKICFU(JL)=0.0_JPRB |
|---|
| 702 | ENDDO |
|---|
| 703 | |
|---|
| 704 | DO JNU= 1,NSIL |
|---|
| 705 | DO JL = KIDIA,KFDIA |
|---|
| 706 | ZRES(JL) = XP(1,JNU)+ZTICE(JL)*(XP(2,JNU)+ZTICE(JL)*(XP(3,& |
|---|
| 707 | & JNU)& |
|---|
| 708 | & +ZTICE(JL)*(XP(4,JNU)+ZTICE(JL)*(XP(5,JNU)+ZTICE(JL)*(XP(6,& |
|---|
| 709 | & JNU)& |
|---|
| 710 | & ))))) |
|---|
| 711 | ZBICE(JL) = ZBICE(JL) + ZRES(JL) |
|---|
| 712 | ZGAMICE(JL) = ZGAMICE(JL) + REBCUI(JNU)*ZRES(JL) |
|---|
| 713 | ZALFICE(JL) = ZALFICE(JL) + REBCUJ(JNU)*ZRES(JL) |
|---|
| 714 | ENDDO |
|---|
| 715 | ENDDO |
|---|
| 716 | |
|---|
| 717 | !-- Fu et al. (1998) with M'91 LW scheme |
|---|
| 718 | IF (NICEOPT == 2 .OR. NICEOPT == 3) THEN |
|---|
| 719 | DO JRTM=1,16 |
|---|
| 720 | DO JL=KIDIA,KFDIA |
|---|
| 721 | IF (PT(JL,IKL) < 160.0_JPRB) THEN |
|---|
| 722 | INDLAY=1 |
|---|
| 723 | ZTBLAY =PT(JL,IKL)-160.0_JPRB |
|---|
| 724 | ELSEIF (PT(JL,IKL) < 339.0_JPRB ) THEN |
|---|
| 725 | INDLAY=PT(JL,IKL)-159.0_JPRB |
|---|
| 726 | INDLAY=MAX(INDLAY,1) |
|---|
| 727 | ZTBLAY =PT(JL,IKL)-INT(PT(JL,IKL)) |
|---|
| 728 | ELSE |
|---|
| 729 | INDLAY=180 |
|---|
| 730 | ZTBLAY =PT(JL,IKL)-339.0_JPRB |
|---|
| 731 | ENDIF |
|---|
| 732 | ZADDPLK = TOTPLNK(INDLAY+1,JRTM)-TOTPLNK(INDLAY,JRTM) |
|---|
| 733 | ZPLANCK = DELWAVE(JRTM) * (TOTPLNK(INDLAY,JRTM) + ZTBLAY*ZADDPLK) |
|---|
| 734 | ZBICFU(JL) = ZBICFU(JL) + ZPLANCK |
|---|
| 735 | |
|---|
| 736 | IF (ZIWC(JL) > 0.0_JPRB ) THEN |
|---|
| 737 | ZRATIO = 1.0_JPRB / ZDESR(JL) |
|---|
| 738 | IF (NICEOPT == 2) THEN |
|---|
| 739 | ! ice cloud spectral emissivity a la Fu & Liou (1993) |
|---|
| 740 | ZMABSD = RFULIO(JRTM,1) + ZRATIO & |
|---|
| 741 | & *(RFULIO(JRTM,2) + ZRATIO*RFULIO(JRTM,3)) |
|---|
| 742 | |
|---|
| 743 | ! ice cloud spectral emissivity a la Fu et al (1998) |
|---|
| 744 | ELSEIF (NICEOPT == 3) THEN |
|---|
| 745 | ZMABSD = RFUETA(JRTM,1) + ZRATIO & |
|---|
| 746 | & *(RFUETA(JRTM,2) + ZRATIO*RFUETA(JRTM,3)) |
|---|
| 747 | ENDIF |
|---|
| 748 | ZKICFU(JL) = ZKICFU(JL)+ ZMABSD*ZPLANCK |
|---|
| 749 | ENDIF |
|---|
| 750 | ENDDO |
|---|
| 751 | ENDDO |
|---|
| 752 | ENDIF |
|---|
| 753 | |
|---|
| 754 | DO JL = KIDIA,KFDIA |
|---|
| 755 | ZGAMICE(JL) = ZGAMICE(JL) / ZBICE(JL) |
|---|
| 756 | ZALFICE(JL) = ZALFICE(JL) / ZBICE(JL) |
|---|
| 757 | ZKICFU(JL) = ZKICFU(JL) / ZBICFU(JL) |
|---|
| 758 | |
|---|
| 759 | IF (ZFLWP(JL)+ZFIWP(JL) > REPSCW) THEN |
|---|
| 760 | |
|---|
| 761 | IF (NLIQOPT == 0) THEN |
|---|
| 762 | ! water cloud emissivity a la Smith & Shi (1992) |
|---|
| 763 | ZMULTL=1.2_JPRB-0.006_JPRB*ZRADLP(JL) |
|---|
| 764 | ZMSALD= 0.158_JPRB*ZMULTL |
|---|
| 765 | ZMSALU= 0.130_JPRB*ZMULTL |
|---|
| 766 | |
|---|
| 767 | ELSE |
|---|
| 768 | ! water cloud emissivity a la Savijarvi (1997) |
|---|
| 769 | ZMSALU= 0.2441_JPRB-0.0105_JPRB*ZRADLP(JL) |
|---|
| 770 | ZMSALD= 1.2154_JPRB*ZMSALU |
|---|
| 771 | |
|---|
| 772 | ENDIF |
|---|
| 773 | |
|---|
| 774 | IF (NICEOPT == 0) THEN |
|---|
| 775 | ! ice cloud emissivity a la Smith & Shi (1992) |
|---|
| 776 | ZMULTI=1.2_JPRB-0.006_JPRB*ZRADIP(JL) |
|---|
| 777 | ZMSAID= 0.113_JPRB*ZMULTI |
|---|
| 778 | ZMSAIU= 0.093_JPRB*ZMULTI |
|---|
| 779 | |
|---|
| 780 | ELSEIF (NICEOPT == 1) THEN |
|---|
| 781 | ! ice cloud emissivity a la Ebert & Curry (1992) |
|---|
| 782 | ZMSAID= 1.66_JPRB*(ZALFICE(JL)+ZGAMICE(JL)/ZRADIP(JL)) |
|---|
| 783 | ZMSAIU= ZMSAID |
|---|
| 784 | |
|---|
| 785 | ELSEIF (NICEOPT == 2 .OR. NICEOPT == 3) THEN |
|---|
| 786 | ! ice cloud emissivity a la Fu & Liou (1993) or Fu et al. (1998) |
|---|
| 787 | ZMSAID= 1.66_JPRB*ZKICFU(JL) |
|---|
| 788 | ZMSAIU= ZMSAID |
|---|
| 789 | ENDIF |
|---|
| 790 | |
|---|
| 791 | IF (NINHOM == 1) THEN |
|---|
| 792 | ZZFLWP= ZFLWP(JL) * RLWINHF |
|---|
| 793 | ZZFIWP= ZFIWP(JL) * RLWINHF |
|---|
| 794 | ELSE |
|---|
| 795 | ZZFLWP= ZFLWP(JL) |
|---|
| 796 | ZZFIWP= ZFIWP(JL) |
|---|
| 797 | ENDIF |
|---|
| 798 | |
|---|
| 799 | ! effective cloudiness accounting for condensed water |
|---|
| 800 | ZCLDLD(JL,JK) = PCLFR(JL,IKL)*(1.0_JPRB-EXP(-ZMSALD*ZZFLWP-ZMSAID* & |
|---|
| 801 | & ZZFIWP)) |
|---|
| 802 | ZCLDLU(JL,JK) = PCLFR(JL,IKL)*(1.0_JPRB-EXP(-ZMSALU*ZZFLWP-ZMSAIU* & |
|---|
| 803 | & ZZFIWP)) |
|---|
| 804 | ENDIF |
|---|
| 805 | ENDDO |
|---|
| 806 | |
|---|
| 807 | ELSE |
|---|
| 808 | |
|---|
| 809 | ! 2.5 CLOUD LONGWAVE OPTICAL PROPERTIES FOR RRTM |
|---|
| 810 | ! ------------------------------------------ |
|---|
| 811 | |
|---|
| 812 | ! ------------------------- |
|---|
| 813 | ! --+ LW OPTICAL PARAMETERS + Water (and Ice) from Savijarvi (1998) |
|---|
| 814 | ! ------------------------- Ice clouds (Ebert, Curry, 1992) |
|---|
| 815 | |
|---|
| 816 | ! No need for a fixed diffusivity factor, accounted for spectrally below |
|---|
| 817 | ! The detailed spectral structure does not require defining upward and |
|---|
| 818 | ! downward effective optical properties |
|---|
| 819 | |
|---|
| 820 | DO JRTM=1,16 |
|---|
| 821 | DO JL = KIDIA,KFDIA |
|---|
| 822 | ZTAUCLD(JL,JK,JRTM) = 0.0_JPRB |
|---|
| 823 | ZMSALD = 0.0_JPRB |
|---|
| 824 | ZMSAID = 0.0_JPRB |
|---|
| 825 | |
|---|
| 826 | IF (ZFLWP(JL)+ZFIWP(JL) > REPSCW) THEN |
|---|
| 827 | |
|---|
| 828 | IF (NLIQOPT == 0 .OR. NLIQOPT >= 3 ) THEN |
|---|
| 829 | ! water cloud total emissivity a la Smith and Shi (1992) |
|---|
| 830 | ZMULTL=1.2_JPRB-0.006_JPRB*ZRADLP(JL) |
|---|
| 831 | ZRSALD= 0.144_JPRB*ZMULTL / 1.66_JPRB |
|---|
| 832 | |
|---|
| 833 | ELSEIF (NLIQOPT == 1) THEN |
|---|
| 834 | ! water cloud spectral emissivity a la Savijarvi (1997) |
|---|
| 835 | ZRSALD= RHSAVI(JRTM,1) + ZRADLP(JL)& |
|---|
| 836 | & *(RHSAVI(JRTM,2) + ZRADLP(JL)*RHSAVI(JRTM,3)) |
|---|
| 837 | |
|---|
| 838 | ELSEIF (NLIQOPT == 2) THEN |
|---|
| 839 | ! water cloud spectral emissivity a la Lindner and Li (2000) |
|---|
| 840 | Z1RADL = 1.0_JPRB / ZRADLP(JL) |
|---|
| 841 | ZEXTCF = RLILIA(JRTM,1)+ZRADLP(JL)*RLILIA(JRTM,2)+ Z1RADL*& |
|---|
| 842 | & (RLILIA(JRTM,3) + Z1RADL*(RLILIA(JRTM,4) + Z1RADL*& |
|---|
| 843 | & RLILIA(JRTM,5) )) |
|---|
| 844 | Z1MOMG = RLILIB(JRTM,1) + Z1RADL*RLILIB(JRTM,2) & |
|---|
| 845 | & + ZRADLP(JL) *(RLILIB(JRTM,3) + ZRADLP(JL)*RLILIB(JRTM,4) ) |
|---|
| 846 | ZRSALD = Z1MOMG * ZEXTCF |
|---|
| 847 | ENDIF |
|---|
| 848 | |
|---|
| 849 | IF (NICEOPT == 0) THEN |
|---|
| 850 | ! ice cloud spectral emissivity a la Smith & Shi (1992) |
|---|
| 851 | ZMULTI=1.2_JPRB-0.006_JPRB*ZRADIP(JL) |
|---|
| 852 | ZRSAID= 0.103_JPRB*ZMULTI / 1.66_JPRB |
|---|
| 853 | |
|---|
| 854 | ELSEIF (NICEOPT == 1) THEN |
|---|
| 855 | ! ice cloud spectral emissivity a la Ebert-Curry (1992) |
|---|
| 856 | ZRSAID= REBCUH(JRTM)+REBCUG(JRTM)/ZRADIP(JL) |
|---|
| 857 | |
|---|
| 858 | ELSEIF (NICEOPT == 2) THEN |
|---|
| 859 | ! ice cloud spectral emissivity a la Fu & Liou (1993) |
|---|
| 860 | Z1RADI = 1.0_JPRB / ZDESR(JL) |
|---|
| 861 | ZRSAID = RFULIO(JRTM,1) + Z1RADI & |
|---|
| 862 | & *(RFULIO(JRTM,2) + Z1RADI * RFULIO(JRTM,3)) |
|---|
| 863 | |
|---|
| 864 | ELSEIF (NICEOPT == 3) THEN |
|---|
| 865 | ! ice cloud spectral emissivity a la Fu et al (1998) including |
|---|
| 866 | ! parametrisation for LW scattering effect |
|---|
| 867 | Z1RADI = 1.0_JPRB / ZDESR(JL) |
|---|
| 868 | ZRSAIE = RFUETA(JRTM,1) + Z1RADI & |
|---|
| 869 | &*(RFUETA(JRTM,2) + Z1RADI * RFUETA(JRTM,3)) |
|---|
| 870 | ZRSAIA = Z1RADI*(RFUETB(JRTM,1) +ZDESR(JL)*( RFUETB(JRTM,2) +ZDESR(JL)*( RFUETB(JRTM,3) +ZDESR(JL)* RFUETB(JRTM,4)))) |
|---|
| 871 | ZRSAIG = RFUETC(JRTM,1) +ZDESR(JL)*( RFUETC(JRTM,2) +ZDESR(JL)*( RFUETC(JRTM,3) +ZDESR(JL)* RFUETC(JRTM,4))) |
|---|
| 872 | ZRSAIF = 0.5_JPRB + ZRSAIG*( 0.3738_JPRB + ZRSAIG*( 0.0076_JPRB + ZRSAIG*0.1186_JPRB ) ) |
|---|
| 873 | ZRSAID = (1.0_JPRB - ZRSAIA/ZRSAIE * ZRSAIF) * ZRSAIE |
|---|
| 874 | ENDIF |
|---|
| 875 | |
|---|
| 876 | ZTAUD = ZRSALD*ZFLWP(JL)+ZRSAID*ZFIWP(JL) |
|---|
| 877 | |
|---|
| 878 | ! Diffusivity correction within clouds a la Savijarvi |
|---|
| 879 | IF (LDIFFC) THEN |
|---|
| 880 | ZDIFFD=MIN(MAX(1.517_JPRB-0.156_JPRB*LOG(ZTAUD) , 1.0_JPRB), & |
|---|
| 881 | & 2.0_JPRB) |
|---|
| 882 | ELSE |
|---|
| 883 | ZDIFFD=1.66_JPRB |
|---|
| 884 | ENDIF |
|---|
| 885 | |
|---|
| 886 | ZTAUCLD(JL,JK,JRTM) = ZTAUD*ZDIFFD |
|---|
| 887 | ENDIF |
|---|
| 888 | |
|---|
| 889 | ENDDO |
|---|
| 890 | ENDDO |
|---|
| 891 | |
|---|
| 892 | ENDIF |
|---|
| 893 | |
|---|
| 894 | ENDDO |
|---|
| 895 | |
|---|
| 896 | NUAER = NUA |
|---|
| 897 | NTRAER = NTRA |
|---|
| 898 | |
|---|
| 899 | ! ------------------------------------------------------------------ |
|---|
| 900 | ! |
|---|
| 901 | ! 2.6 SCALING OF OPTICAL THICKNESS |
|---|
| 902 | ! SPECTRALLY, ACCOUNTING FOR VERTICAL VARIABILITY |
|---|
| 903 | |
|---|
| 904 | JEXPLR=NLAYINH |
|---|
| 905 | JXPLDN=2*JEXPLR+1 |
|---|
| 906 | |
|---|
| 907 | IF (NINHOM == 1) THEN |
|---|
| 908 | !-- simple scaling a la Tiedtke (1996) with RSWINHF in SW and RLWINHF in LW |
|---|
| 909 | DO JSW=1,NSW |
|---|
| 910 | DO JK=1,KLEV |
|---|
| 911 | DO JL=KIDIA,KFDIA |
|---|
| 912 | ZTAU(JL,JSW,JK)=ZTAU(JL,JSW,JK) * RSWINHF |
|---|
| 913 | ENDDO |
|---|
| 914 | ENDDO |
|---|
| 915 | ENDDO |
|---|
| 916 | |
|---|
| 917 | DO JRTM=1,16 |
|---|
| 918 | DO JK=1,KLEV |
|---|
| 919 | DO JL=KIDIA,KFDIA |
|---|
| 920 | ZTAUCLD(JL,JK,JRTM)=ZTAUCLD(JL,JK,JRTM) * RLWINHF |
|---|
| 921 | ENDDO |
|---|
| 922 | ENDDO |
|---|
| 923 | ENDDO |
|---|
| 924 | |
|---|
| 925 | ELSEIF (JEXPLR /= 0) THEN |
|---|
| 926 | DO JSW=1,NSW |
|---|
| 927 | DO JK=1,KLEV |
|---|
| 928 | DO JL=KIDIA,KFDIA |
|---|
| 929 | ZSQUAR(JL,JK)=0.0_JPRB |
|---|
| 930 | ZVARIA(JL,JK)=1.0_JPRB |
|---|
| 931 | ENDDO |
|---|
| 932 | ENDDO |
|---|
| 933 | !-- range should be defined from Hogan & Illingworth |
|---|
| 934 | DO JK=1+JEXPLR,KLEV-JEXPLR |
|---|
| 935 | DO JL=KIDIA,KFDIA |
|---|
| 936 | ! ZAVDP(JL)=0.0_JPRB |
|---|
| 937 | ZAVTO(JL)=0.0_JPRB |
|---|
| 938 | ZSQTO(JL)=0.0_JPRB |
|---|
| 939 | ENDDO |
|---|
| 940 | DO JKI=JK-JEXPLR,JK+JEXPLR |
|---|
| 941 | IKI=KLEV+1-JKI |
|---|
| 942 | DO JL=KIDIA,KFDIA |
|---|
| 943 | ! ZAVDP(JL)=ZAVDP(JL)+PDP(JL,IKI)/RG |
|---|
| 944 | ZAVTO(JL)=ZAVTO(JL)+ZTAU(JL,JSW,JKI) |
|---|
| 945 | ENDDO |
|---|
| 946 | ENDDO |
|---|
| 947 | DO JL=KIDIA,KFDIA |
|---|
| 948 | ! ZAVTO(JL)=ZAVTO(JL)/ZAVDP(JL) |
|---|
| 949 | ZAVTO(JL)=ZAVTO(JL)/JXPLDN |
|---|
| 950 | ENDDO |
|---|
| 951 | DO JKI=JK-JEXPLR,JK+JEXPLR |
|---|
| 952 | IKI=KLEV+1-JKI |
|---|
| 953 | DO JL=KIDIA,KFDIA |
|---|
| 954 | ! ZSQTO(JL)=ZSQTO(JL)+(ZTAU(JL,JSW,JKI)/PDP(JL,IKI)-ZAVTO(JL))**2 |
|---|
| 955 | ZSQTO(JL)=ZSQTO(JL)+(ZTAU(JL,JSW,JKI)-ZAVTO(JL))**2 |
|---|
| 956 | ENDDO |
|---|
| 957 | ENDDO |
|---|
| 958 | DO JL=KIDIA,KFDIA |
|---|
| 959 | ZSQTO(JL)=SQRT(ZSQTO(JL)/(JXPLDN*(JXPLDN-1))) |
|---|
| 960 | IF (ZAVTO(JL) > 0.0_JPRB) THEN |
|---|
| 961 | ZVARIA(JL,JK)=(ZSQTO(JL)/ZAVTO(JL))**2 |
|---|
| 962 | ZSQUAR(JL,JK)=EXP(-ZVARIA(JL,JK)) |
|---|
| 963 | ELSE |
|---|
| 964 | ZVARIA(JL,JK)=0.0_JPRB |
|---|
| 965 | ZSQUAR(JL,JK)=1.0_JPRB |
|---|
| 966 | ENDIF |
|---|
| 967 | |
|---|
| 968 | !-- scaling a la Barker |
|---|
| 969 | IF (NINHOM ==2) THEN |
|---|
| 970 | ZTAU(JL,JSW,JK)=ZTAU(JL,JSW,JK)*ZSQUAR(JL,JK) |
|---|
| 971 | |
|---|
| 972 | !-- scaling a la Cairns et al. |
|---|
| 973 | ELSEIF (NINHOM == 3) THEN |
|---|
| 974 | ZVI=ZVARIA(JL,JK) |
|---|
| 975 | ZTAU(JL,JSW,JK) = ZTAU(JL,JSW,JK)/(1.0_JPRB+ZVI) |
|---|
| 976 | ZOMEGA(JL,JSW,JK)= ZOMEGA(JL,JSW,JK) & |
|---|
| 977 | & /(1.0_JPRB + ZVI*(1.0_JPRB-ZOMEGA(JL,JSW,JK) ) ) |
|---|
| 978 | ZCG(JL,JSW,JK) = ZCG(JL,JSW,JK) & |
|---|
| 979 | & *(1.0_JPRB+ZVI*(1.0_JPRB-ZOMEGA(JL,JSW,JK))) & |
|---|
| 980 | & /(1.0_JPRB+ZVI*(1.0_JPRB-ZOMEGA(JL,JSW,JK)*ZCG(JL,JSW,JK))) |
|---|
| 981 | ENDIF |
|---|
| 982 | ENDDO |
|---|
| 983 | ! JL=KIDIA |
|---|
| 984 | ! print 9261,JSW,JK,ZTAU(JL,JSW,JK),ZAVTO(JL),ZSQTO(JL),ZVARIA(JL,JK),ZSQUAR(JL,JK) |
|---|
| 985 | 9261 format(1x,'Varia1 ',2I3,7F10.4) |
|---|
| 986 | ENDDO |
|---|
| 987 | ENDDO |
|---|
| 988 | |
|---|
| 989 | |
|---|
| 990 | DO JRTM=1,16 |
|---|
| 991 | DO JK=1,KLEV |
|---|
| 992 | DO JL=KIDIA,KFDIA |
|---|
| 993 | ZSQUAR(JL,JK)=0.0_JPRB |
|---|
| 994 | ZVARIA(JL,JK)=1.0_JPRB |
|---|
| 995 | ENDDO |
|---|
| 996 | ENDDO |
|---|
| 997 | !-- range to be defined from Hogan & Illingworth |
|---|
| 998 | DO JK=1+JEXPLR,KLEV-JEXPLR |
|---|
| 999 | DO JL=KIDIA,KFDIA |
|---|
| 1000 | ! ZAVDP(JL)=0.0_JPRB |
|---|
| 1001 | ZAVTO(JL)=0.0_JPRB |
|---|
| 1002 | ZSQTO(JL)=0.0_JPRB |
|---|
| 1003 | ENDDO |
|---|
| 1004 | DO JKI=JK-JEXPLR,JK+JEXPLR |
|---|
| 1005 | IKI=KLEV+1-JKI |
|---|
| 1006 | DO JL=KIDIA,KFDIA |
|---|
| 1007 | ! ZAVDP(JL)=ZAVDP(JL)+PDP(JL,IKI)/RG |
|---|
| 1008 | ZAVTO(JL)=ZAVTO(JL)+ZTAUCLD(JL,JKI,JRTM) |
|---|
| 1009 | ENDDO |
|---|
| 1010 | ENDDO |
|---|
| 1011 | DO JL=KIDIA,KFDIA |
|---|
| 1012 | ! ZAVTO(JL)=ZAVTO(JL)/ZAVDP(JL) |
|---|
| 1013 | ZAVTO(JL)=ZAVTO(JL)/JXPLDN |
|---|
| 1014 | ENDDO |
|---|
| 1015 | DO JKI=JK-JEXPLR,JK+JEXPLR |
|---|
| 1016 | IKI=KLEV+1-JKI |
|---|
| 1017 | DO JL=KIDIA,KFDIA |
|---|
| 1018 | ! ZSQTO(JL)=ZSQTO(JL)+(ZTAUCLD(JL,JKI,JRTM)/PDP(JL,IKI)-ZAVTO(JL))**2 |
|---|
| 1019 | ZSQTO(JL)=ZSQTO(JL)+(ZTAUCLD(JL,JKI,JRTM)-ZAVTO(JL))**2 |
|---|
| 1020 | ENDDO |
|---|
| 1021 | ENDDO |
|---|
| 1022 | DO JL=KIDIA,KFDIA |
|---|
| 1023 | ZSQTO(JL)=SQRT(ZSQTO(JL)/(JXPLDN*(JXPLDN-1))) |
|---|
| 1024 | IF (ZAVTO(JL) > 0.0_JPRB) THEN |
|---|
| 1025 | ZVARIA(JL,JK)=(ZSQTO(JL)/ZAVTO(JL))**2 |
|---|
| 1026 | ZSQUAR(JL,JK)=EXP(-ZVARIA(JL,JK)) |
|---|
| 1027 | ELSE |
|---|
| 1028 | ZVARIA(JL,JK)=0.0_JPRB |
|---|
| 1029 | ZSQUAR(JL,JK)=1.0_JPRB |
|---|
| 1030 | ENDIF |
|---|
| 1031 | |
|---|
| 1032 | !-- scaling a la Barker |
|---|
| 1033 | IF (NINHOM ==2) THEN |
|---|
| 1034 | ZTAUCLD(JL,JK,JRTM)=ZTAUCLD(JL,JK,JRTM)*ZSQUAR(JL,JK) |
|---|
| 1035 | |
|---|
| 1036 | !-- scaling a la Cairns et al. |
|---|
| 1037 | ELSEIF (NINHOM == 3) THEN |
|---|
| 1038 | ZVI=ZVARIA(JL,JK) |
|---|
| 1039 | ZTAUCLD(JL,JK,JRTM)=ZTAUCLD(JL,JK,JRTM)/(1.0_JPRB+ZVI) |
|---|
| 1040 | ENDIF |
|---|
| 1041 | ENDDO |
|---|
| 1042 | ! JL=KIDIA |
|---|
| 1043 | ! print 9262,JRTM,JK,ZTAUCLD(JL,JK,JRTM),ZAVTO(JL),ZSQTO(JL),ZVARIA(JL,JK),ZSQUAR(JL,JK) |
|---|
| 1044 | 9262 format(1x,'Varia2 ',2I3,7F10.4) |
|---|
| 1045 | ENDDO |
|---|
| 1046 | ENDDO |
|---|
| 1047 | ENDIF |
|---|
| 1048 | |
|---|
| 1049 | |
|---|
| 1050 | |
|---|
| 1051 | ! ------------------------------------------------------------------ |
|---|
| 1052 | ! |
|---|
| 1053 | !* 2.7 DIFFUSIVITY FACTOR OR SATELLITE VIEWING ANGLE |
|---|
| 1054 | ! --------------------------------------------- |
|---|
| 1055 | |
|---|
| 1056 | DO JL = KIDIA,KFDIA |
|---|
| 1057 | ZVIEW(JL) = DIFF |
|---|
| 1058 | ENDDO |
|---|
| 1059 | |
|---|
| 1060 | ! ------------------------------------------------------------------ |
|---|
| 1061 | |
|---|
| 1062 | !* 3. CALL LONGWAVE RADIATION CODE |
|---|
| 1063 | ! ---------------------------- |
|---|
| 1064 | |
|---|
| 1065 | !* 3.1 FULL LONGWAVE RADIATION COMPUTATIONS |
|---|
| 1066 | ! ------------------------------------ |
|---|
| 1067 | |
|---|
| 1068 | !print *,'RADLSW: LPHYLIN, LRRTM',LPHYLIN, LRRTM |
|---|
| 1069 | IF (.NOT.LPHYLIN) THEN |
|---|
| 1070 | IF ( .NOT. LRRTM) THEN |
|---|
| 1071 | |
|---|
| 1072 | CALL LW & |
|---|
| 1073 | & ( KIDIA , KFDIA , KLON , KLEV , KMODE,& |
|---|
| 1074 | & PCCO2 , ZCLDLD, ZCLDLU,& |
|---|
| 1075 | & PDP , ZDT0 , ZEMIS , ZEMIW,& |
|---|
| 1076 | & ZPMB , POZON , ZTL,& |
|---|
| 1077 | & PAER , ZTAVE , ZVIEW , PQ,& |
|---|
| 1078 | & ZEMIT , PFLUX , PFLUC & |
|---|
| 1079 | & ) |
|---|
| 1080 | ! print *,'RADLSW: apres CALL LW' |
|---|
| 1081 | IF(LLDEBUG) THEN |
|---|
| 1082 | call writefield_phy('radlsw_flux1',PFLUX(:,1,:),klev+1) |
|---|
| 1083 | call writefield_phy('radlsw_flux2',PFLUX(:,2,:),klev+1) |
|---|
| 1084 | call writefield_phy('radlsw_fluc1',PFLUC(:,1,:),klev+1) |
|---|
| 1085 | call writefield_phy('radlsw_fluc2',PFLUC(:,2,:),klev+1) |
|---|
| 1086 | ENDIF |
|---|
| 1087 | |
|---|
| 1088 | ELSE |
|---|
| 1089 | |
|---|
| 1090 | !* 3.2 FULL LONGWAVE RADIATION COMPUTATIONS - RRTM |
|---|
| 1091 | ! ------------------------------------ ---- |
|---|
| 1092 | |
|---|
| 1093 | ! i) pass ZOZN (ozone mass mixing ratio) to RRTM; remove pressure |
|---|
| 1094 | ! weighting applied to POZON in driverMC (below) |
|---|
| 1095 | ! ii) pass ZEMIS and ZEMIW to RRTM; return ZEMIT from RRTM |
|---|
| 1096 | ! iii)pass ZTAUCLD, cloud optical depths (water+ice) to RRTM, |
|---|
| 1097 | ! computed from equations above |
|---|
| 1098 | ! iv) pass ECRT arrays to RRTM arrays in interface routine ECRTATM |
|---|
| 1099 | ! in module rrtm_ecrt.f |
|---|
| 1100 | |
|---|
| 1101 | DO JL = KIDIA,KFDIA |
|---|
| 1102 | DO JK = 1, KLEV |
|---|
| 1103 | ZOZN(JL,JK) = POZON(JL,JK)/PDP(JL,JK) |
|---|
| 1104 | ENDDO |
|---|
| 1105 | ENDDO |
|---|
| 1106 | |
|---|
| 1107 | ! print *,'RADLSW: avant CALL RRTM_RRTM_140GP,PAP=',PAP(1,:) |
|---|
| 1108 | CALL RRTM_RRTM_140GP & |
|---|
| 1109 | & ( KIDIA , KFDIA , KLON , KLEV,& |
|---|
| 1110 | & PAER , PAPH , PAP,& |
|---|
| 1111 | & PTS , PTH , PT,& |
|---|
| 1112 | & ZEMIS , ZEMIW,& |
|---|
| 1113 | & PQ , PCCO2 , ZOZN ,& |
|---|
| 1114 | & ZCLDSW , ZTAUCLD,& |
|---|
| 1115 | & PTAU_LW,& |
|---|
| 1116 | & ZEMIT , PFLUX , PFLUC , ZTCLEAR ) |
|---|
| 1117 | ! print *,'RADLSW: apres CALL RRTM_RRTM_140GP' |
|---|
| 1118 | |
|---|
| 1119 | ENDIF |
|---|
| 1120 | ELSE |
|---|
| 1121 | ZEMIT (:) = 0.0_JPRB |
|---|
| 1122 | PFLUX(:,:,:)= 0.0_JPRB |
|---|
| 1123 | PFLUC(:,:,:)= 0.0_JPRB |
|---|
| 1124 | ! print *,'RADLSW: ZEMIT,PFLUX et PFLUC = 0' |
|---|
| 1125 | ENDIF |
|---|
| 1126 | |
|---|
| 1127 | ! ------------------------------------------------------------------ |
|---|
| 1128 | |
|---|
| 1129 | !* 4. CALL SHORTWAVE RADIATION CODE |
|---|
| 1130 | ! ----------------------------- |
|---|
| 1131 | |
|---|
| 1132 | ZRMUZ=0.0_JPRB |
|---|
| 1133 | DO JL = KIDIA,KFDIA |
|---|
| 1134 | ZRMUZ = MAX (ZRMUZ, ZMU0(JL)) |
|---|
| 1135 | ENDDO |
|---|
| 1136 | |
|---|
| 1137 | IF (NSTEP == 0 .AND. LEDBUG .AND. ZMU0(KIDIA) > 0.0_JPRB) THEN |
|---|
| 1138 | WRITE(NULOUT,'(4E15.8)') PRII0,PCCO2,ZPSOL(KIDIA),ZMU0(KIDIA) |
|---|
| 1139 | WRITE(NULOUT,'("ZALBD ",6E15.8)') (ZALBD(KIDIA,JSW),JSW=1,NSW) |
|---|
| 1140 | WRITE(NULOUT,'("ZALBP ",6E15.8)') (ZALBP(KIDIA,JSW),JSW=1,NSW) |
|---|
| 1141 | WRITE(NULOUT,'("PQ ",10E12.5)') (PQ(KIDIA,JK),JK=1,KLEV) |
|---|
| 1142 | WRITE(NULOUT,'("PQS ",10E12.5)') (PQS(KIDIA,JK),JK=1,KLEV) |
|---|
| 1143 | WRITE(NULOUT,'("PDP ",10E12.5)') (PDP(KIDIA,JK),JK=1,KLEV) |
|---|
| 1144 | WRITE(NULOUT,'("ZPMB ",10E12.5)') (ZPMB(KIDIA,JK),JK=1,KLEV+1) |
|---|
| 1145 | WRITE(NULOUT,'("ZTAVE ",10E12.5)') (ZTAVE(KIDIA,JK),JK=1,KLEV) |
|---|
| 1146 | WRITE(NULOUT,'("ZCLDSW",10E12.5)') (ZCLDSW(KIDIA,JK),JK=1,KLEV) |
|---|
| 1147 | WRITE(NULOUT,'("ZTAU ",10E12.5)') ((ZTAU(KIDIA,JSW,JK),JK=1,KLEV),JSW=1,NSW) |
|---|
| 1148 | WRITE(NULOUT,'("ZCG ",10E12.5)') ((ZCG(KIDIA,JSW,JK),JK=1,KLEV),JSW=1,NSW) |
|---|
| 1149 | WRITE(NULOUT,'("ZOMEGA",10E12.5)') ((ZOMEGA(KIDIA,JSW,JK),JK=1,KLEV),JSW=1,NSW) |
|---|
| 1150 | WRITE(NULOUT,'("ZOZ ",10E12.5)') (ZOZ(KIDIA,JK),JK=1,KLEV) |
|---|
| 1151 | WRITE(NULOUT,'("PAER ",10E12.5)') ((PAER(KIDIA,JSW,JK),JK=1,KLEV),JSW=1,NSW) |
|---|
| 1152 | ENDIF |
|---|
| 1153 | |
|---|
| 1154 | IF (NSTEP == 0 .AND. LEDBUG .AND. ZMU0(KIDIA) > 0.0_JPRB) THEN |
|---|
| 1155 | WRITE(NULOUT,'(4E15.8)') PRII0,PCCO2,ZPSOL(KIDIA),ZMU0(KIDIA) |
|---|
| 1156 | WRITE(NULOUT,'("ZALBD ",6E15.8)') (ZALBD(KIDIA,JSW),JSW=1,NSW) |
|---|
| 1157 | WRITE(NULOUT,'("ZALBP ",6E15.8)') (ZALBP(KIDIA,JSW),JSW=1,NSW) |
|---|
| 1158 | WRITE(NULOUT,'("PQ ",10E12.5)') (PQ(KIDIA,JK),JK=1,KLEV) |
|---|
| 1159 | WRITE(NULOUT,'("PQS ",10E12.5)') (PQS(KIDIA,JK),JK=1,KLEV) |
|---|
| 1160 | WRITE(NULOUT,'("PDP ",10E12.5)') (PDP(KIDIA,JK),JK=1,KLEV) |
|---|
| 1161 | WRITE(NULOUT,'("ZPMB ",10E12.5)') (ZPMB(KIDIA,JK),JK=1,KLEV+1) |
|---|
| 1162 | WRITE(NULOUT,'("ZTAVE ",10E12.5)') (ZTAVE(KIDIA,JK),JK=1,KLEV) |
|---|
| 1163 | WRITE(NULOUT,'("ZCLDSW",10E12.5)') (ZCLDSW(KIDIA,JK),JK=1,KLEV) |
|---|
| 1164 | WRITE(NULOUT,'("ZTAU ",10E12.5)') ((ZTAU(KIDIA,JSW,JK),JK=1,KLEV),JSW=1,NSW) |
|---|
| 1165 | WRITE(NULOUT,'("ZCG ",10E12.5)') ((ZCG(KIDIA,JSW,JK),JK=1,KLEV),JSW=1,NSW) |
|---|
| 1166 | WRITE(NULOUT,'("ZOMEGA",10E12.5)') ((ZOMEGA(KIDIA,JSW,JK),JK=1,KLEV),JSW=1,NSW) |
|---|
| 1167 | WRITE(NULOUT,'("ZOZ ",10E12.5)') (ZOZ(KIDIA,JK),JK=1,KLEV) |
|---|
| 1168 | WRITE(NULOUT,'("PAER ",10E12.5)') ((PAER(KIDIA,JSW,JK),JK=1,KLEV),JSW=1,NSW) |
|---|
| 1169 | ENDIF |
|---|
| 1170 | CALL SW & |
|---|
| 1171 | & ( KIDIA , KFDIA , KLON , KLEV , KAER,& |
|---|
| 1172 | & PRII0 , PCCO2 , ZPSOL , ZALBD , ZALBP , PQ , PQS,& |
|---|
| 1173 | & ZMU0 , ZCG , ZCLDSW, PDP , ZOMEGA, ZOZ , ZPMB,& |
|---|
| 1174 | & ZTAU , ZTAVE , PAER,& |
|---|
| 1175 | & PFSDN , PFSUP , PFSCDN, PFSCUP,& |
|---|
| 1176 | & ZFSDNN, ZFSDNV, ZFSUPN, ZFSUPV,& |
|---|
| 1177 | & ZFCDNN, ZFCDNV, ZFCUPN, ZFCUPV,& |
|---|
| 1178 | & ZSUDU , ZUVDF , ZPARF ,ZPARCF, ZDIFFS, ZDIRFS, & |
|---|
| 1179 | & LRDUST,PPIZA_DST,PCGA_DST,PTAUREL_DST& |
|---|
| 1180 | & ) |
|---|
| 1181 | PFSDNV=ZFSDNV |
|---|
| 1182 | PFSDNN=ZFSDNN |
|---|
| 1183 | IF (SIZE(PSFSWDIR,2)>1) THEN |
|---|
| 1184 | PSFSWDIR= ZDIRFS |
|---|
| 1185 | PSFSWDIF= ZDIFFS |
|---|
| 1186 | ELSE |
|---|
| 1187 | PSFSWDIR (:,1) = ZFSDNV(:) + ZFSDNN(:) |
|---|
| 1188 | PSFSWDIF (:,:) = 0. |
|---|
| 1189 | ENDIF |
|---|
| 1190 | |
|---|
| 1191 | IF (NSTEP == 0 .AND. LEDBUG .AND. ZMU0(KIDIA) > 0.0_JPRB) THEN |
|---|
| 1192 | WRITE(NULOUT,'("ZFSDWN",10E12.5)') (ZFSDWN(KIDIA,JK),JK=1,KLEV) |
|---|
| 1193 | WRITE(NULOUT,'("ZFSUP ",10E12.5)') (ZFSUP (KIDIA,JK),JK=1,KLEV) |
|---|
| 1194 | WRITE(NULOUT,'("ZFCDWN",10E12.5)') (ZFCDWN(KIDIA,JK),JK=1,KLEV) |
|---|
| 1195 | WRITE(NULOUT,'("ZFCUP ",10E12.5)') (ZFCUP (KIDIA,JK),JK=1,KLEV) |
|---|
| 1196 | LEDBUG=.FALSE. |
|---|
| 1197 | ENDIF |
|---|
| 1198 | IF (NSTEP == 0 .AND. LEDBUG .AND. ZMU0(KIDIA) > 0.0_JPRB) THEN |
|---|
| 1199 | WRITE(NULOUT,'("ZFSDWN",10E12.5)') (ZFSDWN(KIDIA,JK),JK=1,KLEV) |
|---|
| 1200 | WRITE(NULOUT,'("ZFSUP ",10E12.5)') (ZFSUP (KIDIA,JK),JK=1,KLEV) |
|---|
| 1201 | WRITE(NULOUT,'("ZFCDWN",10E12.5)') (ZFCDWN(KIDIA,JK),JK=1,KLEV) |
|---|
| 1202 | WRITE(NULOUT,'("ZFCUP ",10E12.5)') (ZFCUP (KIDIA,JK),JK=1,KLEV) |
|---|
| 1203 | LEDBUG=.FALSE. |
|---|
| 1204 | ENDIF |
|---|
| 1205 | ! ------------------------------------------------------------------ |
|---|
| 1206 | |
|---|
| 1207 | !* 5. FILL UP THE MODEL NET LW AND SW RADIATIVE FLUXES |
|---|
| 1208 | ! ------------------------------------------------ |
|---|
| 1209 | |
|---|
| 1210 | DO JKL = 1 , KLEV+1 |
|---|
| 1211 | JK = KLEV+1 + 1 - JKL |
|---|
| 1212 | DO JL = KIDIA,KFDIA |
|---|
| 1213 | PFLS(JL,JKL) = ZFSDWN(JL,JK) - ZFSUP(JL,JK) |
|---|
| 1214 | PFLT(JL,JKL) = - PFLUX(JL,1,JK) - PFLUX(JL,2,JK) |
|---|
| 1215 | PFCS(JL,JKL) = ZFCDWN(JL,JK) - ZFCUP(JL,JK) |
|---|
| 1216 | PFCT(JL,JKL) = - PFLUC(JL,1,JK) - PFLUC(JL,2,JK) |
|---|
| 1217 | ENDDO |
|---|
| 1218 | ENDDO |
|---|
| 1219 | |
|---|
| 1220 | DO JL = KIDIA,KFDIA |
|---|
| 1221 | PFRSOD(JL)=ZFSDWN(JL,1) |
|---|
| 1222 | PEMIT (JL)=ZEMIT (JL) |
|---|
| 1223 | PSUDU (JL)=ZSUDU (JL) |
|---|
| 1224 | PUVDF (JL)=ZUVDF (JL) |
|---|
| 1225 | PPARF (JL)=ZPARF (JL) |
|---|
| 1226 | PPARCF(JL)=ZPARCF(JL) |
|---|
| 1227 | PTINCF(JL)=PRII0 * ZMU0(JL) |
|---|
| 1228 | ENDDO |
|---|
| 1229 | !print 9501,(PUVDF(JL),JL=KIDIA,KFDIA) |
|---|
| 1230 | 9501 format(1x,'RADLSW PUVDF: ',30f6.1) |
|---|
| 1231 | !print 9502,(PPARF(JL),JL=KIDIA,KFDIA) |
|---|
| 1232 | 9502 format(1x,'RADLSW PPARF: ',30f6.1) |
|---|
| 1233 | |
|---|
| 1234 | ! -------------------------------------------------------------- |
|---|
| 1235 | |
|---|
| 1236 | IF (LHOOK) CALL DR_HOOK('RADLSW',1,ZHOOK_HANDLE) |
|---|
| 1237 | END SUBROUTINE RADLSW |
|---|
| 1238 | |
|---|
| 1239 | |
|---|
| 1240 | |
|---|
| 1241 | |
|---|
| 1242 | |
|---|
| 1243 | |
|---|
| 1244 | |
|---|
| 1245 | |
|---|
| 1246 | |
|---|
| 1247 | |
|---|
| 1248 | |
|---|
| 1249 | |
|---|
| 1250 | |
|---|
| 1251 | |
|---|
| 1252 | |
|---|
| 1253 | |
|---|
| 1254 | |
|---|
| 1255 | |
|---|
| 1256 | |
|---|
| 1257 | |
|---|
| 1258 | |
|---|
| 1259 | |
|---|
| 1260 | |
|---|
| 1261 | |
|---|
| 1262 | |
|---|
| 1263 | |
|---|
| 1264 | |
|---|
| 1265 | |
|---|
| 1266 | |
|---|
| 1267 | |
|---|
| 1268 | |
|---|
| 1269 | |
|---|
| 1270 | |
|---|
| 1271 | |
|---|
