[2005] | 1 | ! |
---|
| 2 | ! $Id: recmwf_aero.F90 5294 2024-10-29 18:35:00Z aborella $ |
---|
| 3 | ! |
---|
[2004] | 4 | !OPTIONS XOPT(NOEVAL) |
---|
| 5 | SUBROUTINE RECMWF_AERO (KST, KEND, KPROMA, KTDIA , KLEV,& |
---|
[4875] | 6 | & KMODE,& |
---|
| 7 | & PALBD , PALBP , PAPRS , PAPRSF , PCCO2 , PCLFR,& |
---|
| 8 | & PQO3 , PAER , PDP , PEMIS , PMU0,& |
---|
| 9 | & PQ , PQS , PQIWP , PQLWP , PSLM , PT , PTS,& |
---|
| 10 | & PREF_LIQ, PREF_ICE,& |
---|
| 11 | !--OB |
---|
| 12 | & PREF_LIQ_PI, PREF_ICE_PI,& |
---|
| 13 | !--fin |
---|
| 14 | & PEMTD , PEMTU , PTRSO,& |
---|
| 15 | & PTH , PCTRSO, PCEMTR, PTRSOD,& |
---|
| 16 | & PLWFC, PLWFT, PSWFC, PSWFT, PSFSWDIR, PSFSWDIF,& |
---|
| 17 | & PFSDNN, PFSDNV,& |
---|
| 18 | & PPIZA_TOT,PCGA_TOT,PTAU_TOT, & |
---|
| 19 | !--OB |
---|
| 20 | & PPIZA_NAT,PCGA_NAT,PTAU_NAT, & |
---|
| 21 | !--fin OB |
---|
| 22 | !--C.Kleinschmitt |
---|
| 23 | & PTAU_LW_TOT, PTAU_LW_NAT, & |
---|
| 24 | !--end |
---|
| 25 | & PFLUX,PFLUC,& |
---|
| 26 | & PFSDN ,PFSUP , PFSCDN , PFSCUP, PFSCCDN, PFSCCUP, PFLCCDN, PFLCCUP,& |
---|
| 27 | !--OB diagnostics |
---|
| 28 | & PTOPSWADAERO,PSOLSWADAERO,& |
---|
| 29 | & PTOPSWAD0AERO,PSOLSWAD0AERO,& |
---|
| 30 | & PTOPSWAIAERO,PSOLSWAIAERO,& |
---|
| 31 | & PTOPSWCFAERO,PSOLSWCFAERO,& |
---|
| 32 | & PSWADAERO,& !--NL |
---|
| 33 | !--LW diagnostics CK |
---|
| 34 | & PTOPLWADAERO,PSOLLWADAERO,& |
---|
| 35 | & PTOPLWAD0AERO,PSOLLWAD0AERO,& |
---|
| 36 | & PTOPLWAIAERO,PSOLLWAIAERO,& |
---|
| 37 | & PLWADAERO,& !--NL |
---|
| 38 | !--ajout volmip |
---|
| 39 | & volmip_solsw, flag_volc_surfstrat,& |
---|
| 40 | !..end |
---|
| 41 | & ok_ade, ok_aie, ok_volcan, flag_aerosol,flag_aerosol_strat,& |
---|
| 42 | & flag_aer_feedback) |
---|
| 43 | !--fin |
---|
[2004] | 44 | |
---|
[4875] | 45 | !**** *RECMWF* - METEO-FRANCE RADIATION INTERFACE TO ECMWF RADIATION SCHEME |
---|
[2004] | 46 | |
---|
[4875] | 47 | ! PURPOSE. |
---|
| 48 | ! -------- |
---|
| 49 | ! SIMPLE INTERFACE TO RADLSW (NO INTERPOLATION) |
---|
[2004] | 50 | |
---|
[4875] | 51 | !** INTERFACE. |
---|
| 52 | ! ---------- |
---|
[2004] | 53 | |
---|
[4875] | 54 | ! EXPLICIT ARGUMENTS : |
---|
| 55 | ! -------------------- |
---|
| 56 | ! KST : START INDEX OF DATA IN KPROMA-LONG VECTOR |
---|
| 57 | ! KEND : END INDEX OF DATA IN KPROMA-LONG VECTOR |
---|
| 58 | ! KPROMA : VECTOR LENGTH |
---|
| 59 | ! KTDIA : INDEX OF TOP LEVEL FROM WHICH COMPUTATIONS ARE ACTIVE |
---|
| 60 | ! KLEV : NUMBER OF LEVELS |
---|
| 61 | ! PAER : (KPROMA,KLEV ,6) ; OPTICAL THICKNESS OF THE AEROSOLS |
---|
| 62 | ! PALBD : (KPROMA,NSW) ; DIFFUSE ALBEDO IN THE 2 SW INTERVALS |
---|
| 63 | ! PALBP : (KPROMA,NSW) ; PARALLEL ALBEDO IN THE 2 SW INTERVALS |
---|
| 64 | ! PAPRS : (KPROMA,KLEV+1) ; HALF LEVEL PRESSURE |
---|
| 65 | ! PAPRSF : (KPROMA,KLEV ) ; FULL LEVEL PRESSURE |
---|
| 66 | ! PCCO2 : ; CONCENTRATION IN CO2 (PA/PA) |
---|
| 67 | ! PCLFR : (KPROMA,KLEV ) ; CLOUD FRACTIONAL COVER |
---|
| 68 | ! PQO3 : (KPROMA,KLEV ) ; OZONE MIXING RATIO (MASS) |
---|
| 69 | ! PDP : (KPROMA,KLEV) ; LAYER PRESSURE THICKNESS |
---|
| 70 | ! PEMIS : (KPROMA) ; SURFACE EMISSIVITY |
---|
| 71 | ! PMU0 : (KPROMA) ; SOLAR ANGLE |
---|
| 72 | ! PQ : (KPROMA,KLEV ) ; SPECIFIC HUMIDITY PA/PA |
---|
| 73 | ! PQS : (KPROMA,KLEV ) ; SATURATION SPECIFIC HUMIDITY PA/PA |
---|
| 74 | ! PQIWP : (KPROMA,KLEV ) ; ICE WATER KG/KG |
---|
| 75 | ! PQLWP : (KPROMA,KLEV ) ; LIQUID WATER KG/KG |
---|
| 76 | ! PSLM : (KPROMA) ; LAND-SEA MASK |
---|
| 77 | ! PT : (KPROMA,KLEV) ; FULL LEVEL TEMPERATURE |
---|
| 78 | ! PTS : (KPROMA) ; SURFACE TEMPERATURE |
---|
| 79 | ! PPIZA_TOT : (KPROMA,KLEV,NSW); Single scattering albedo of total aerosol |
---|
| 80 | ! PCGA_TOT : (KPROMA,KLEV,NSW); Assymetry factor for total aerosol |
---|
| 81 | ! PTAU_TOT: (KPROMA,KLEV,NSW) ; Optical depth of total aerosol |
---|
| 82 | ! PREF_LIQ (KPROMA,KLEV) ; Liquid droplet radius (um) - present-day |
---|
| 83 | ! PREF_ICE (KPROMA,KLEV) ; Ice crystal radius (um) - present-day |
---|
| 84 | !--OB |
---|
| 85 | ! PREF_LIQ_PI (KPROMA,KLEV) ; Liquid droplet radius (um) - pre-industrial |
---|
| 86 | ! PREF_ICE_PI (KPROMA,KLEV) ; Ice crystal radius (um) - pre-industrial |
---|
| 87 | ! ok_ade---input-L- apply the Aerosol Direct Effect or not? |
---|
| 88 | ! ok_aie---input-L- apply the Aerosol Indirect Effect or not? |
---|
| 89 | ! ok_volcan-input-L- activate volcanic diags (SW heat & LW cool rate, SW & LW flux) |
---|
| 90 | ! flag_aerosol-input-I- aerosol flag from 0 to 7 |
---|
| 91 | ! flag_aerosol_strat-input-I- use stratospheric aerosols flag (T/F) |
---|
| 92 | ! flag_aer_feedback-input-I- use aerosols radiative effect flag (T/F) |
---|
| 93 | ! PPIZA_NAT : (KPROMA,KLEV,NSW); Single scattering albedo of natural aerosol |
---|
| 94 | ! PCGA_NAT : (KPROMA,KLEV,NSW); Assymetry factor for natural aerosol |
---|
| 95 | ! PTAU_NAT: (KPROMA,KLEV,NSW) ; Optical depth of natural aerosol |
---|
| 96 | ! PTAU_LW_TOT (KPROMA,KLEV,NLW); LW Optical depth of total aerosols |
---|
| 97 | ! PTAU_LW_NAT (KPROMA,KLEV,NLW); LW Optical depth of natural aerosols |
---|
| 98 | !--fin OB |
---|
[2004] | 99 | |
---|
[4875] | 100 | ! ==== OUTPUTS === |
---|
| 101 | ! PEMTD (KPROMA,KLEV+1) ; TOTAL DOWNWARD LONGWAVE EMISSIVITY |
---|
| 102 | ! PEMTU (KPROMA,KLEV+1) ; TOTAL UPWARD LONGWAVE EMISSIVITY |
---|
| 103 | ! PTRSO (KPROMA,KLEV+1) ; TOTAL SHORTWAVE TRANSMISSIVITY |
---|
| 104 | ! PTH (KPROMA,KLEV+1) ; HALF LEVEL TEMPERATURE |
---|
| 105 | ! PCTRSO(KPROMA,2) ; CLEAR-SKY SHORTWAVE TRANSMISSIVITY |
---|
| 106 | ! PCEMTR(KPROMA,2) ; CLEAR-SKY NET LONGWAVE EMISSIVITY |
---|
| 107 | ! PTRSOD(KPROMA) ; TOTAL-SKY SURFACE SW TRANSMISSITY |
---|
| 108 | ! PLWFC (KPROMA,2) ; CLEAR-SKY LONGWAVE FLUXES |
---|
| 109 | ! PLWFT (KPROMA,KLEV+1) ; TOTAL-SKY LONGWAVE FLUXES |
---|
| 110 | ! PSWFC (KPROMA,2) ; CLEAR-SKY SHORTWAVE FLUXES |
---|
| 111 | ! PSWFT (KPROMA,KLEV+1) ; TOTAL-SKY SHORTWAVE FLUXES |
---|
| 112 | ! Ajout flux LW et SW montants et descendants, et ciel clair (MPL 19.12.08) |
---|
| 113 | ! PFLUX (KPROMA,2,KLEV+1) ; LW total sky flux (1=up, 2=down) |
---|
| 114 | ! PFLUC (KPROMA,2,KLEV+1) ; LW clear sky flux (1=up, 2=down) |
---|
| 115 | ! PFSDN(KPROMA,KLEV+1) ; SW total sky flux down |
---|
| 116 | ! PFSUP(KPROMA,KLEV+1) ; SW total sky flux up |
---|
| 117 | ! PFSCDN(KPROMA,KLEV+1) ; SW clear sky flux down |
---|
| 118 | ! PFSCUP(KPROMA,KLEV+1) ; SW clear sky flux up |
---|
| 119 | ! PFSCCDN(KPROMA,KLEV+1) ; SW clear sky clean (no aerosol) flux down |
---|
| 120 | ! PFSCCUP(KPROMA,KLEV+1) ; SW clear sky clean (no aerosol) flux up |
---|
| 121 | ! PFLCCDN(KPROMA,KLEV+1) ; LW clear sky clean (no aerosol) flux down |
---|
| 122 | ! PFLCCUP(KPROMA,KLEV+1) ; LW clear sky clean (no aerosol) flux up |
---|
[2004] | 123 | |
---|
| 124 | |
---|
[4875] | 125 | ! IMPLICIT ARGUMENTS : NONE |
---|
| 126 | ! -------------------- |
---|
[2004] | 127 | |
---|
[4875] | 128 | ! METHOD. |
---|
| 129 | ! ------- |
---|
| 130 | ! SEE DOCUMENTATION |
---|
[2004] | 131 | |
---|
[4875] | 132 | ! EXTERNALS. |
---|
| 133 | ! ---------- |
---|
[2004] | 134 | |
---|
[4875] | 135 | ! REFERENCE. |
---|
| 136 | ! ---------- |
---|
| 137 | ! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS |
---|
[2004] | 138 | |
---|
[4875] | 139 | ! AUTHORS. |
---|
| 140 | ! -------- |
---|
| 141 | ! ORIGINAL BY B. RITTER *ECMWF* 83-10-13 |
---|
| 142 | ! REWRITING FOR IFS BY J.-J. MORCRETTE 94-11-15 |
---|
| 143 | ! 96-11: Ph. Dandin. Meteo-France |
---|
| 144 | ! REWRITING FOR DM BY J.PH. PIEDELIEVRE 1998-07 |
---|
| 145 | ! Duplication of RFMR to use present (cy25) ECMWF radiation scheme : Y. Bouteloup 09-2003 |
---|
| 146 | ! Use of 6 aerosols & introduce NSW : F. Bouyssel 09-2004 |
---|
| 147 | ! 04-11-18 : 4 New arguments for AROME : Y. Seity |
---|
| 148 | ! 2005-10-10 Y. Seity : 3 optional arguments for dust optical properties |
---|
| 149 | ! JJMorcrette 20060721 PP of clear-sky PAR and TOA incident solar radiation (ECMWF) |
---|
| 150 | ! Olivier Boucher: added LMD radiation diagnostics 2014-03 |
---|
[2004] | 151 | |
---|
[4875] | 152 | !----------------------------------------------------------------------- |
---|
[2004] | 153 | |
---|
[4875] | 154 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
---|
| 155 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
---|
| 156 | USE YOEAERD , ONLY : RCAEROS |
---|
| 157 | USE YOMCST , ONLY : RMD ,RMO3 |
---|
| 158 | USE YOMPHY3 , ONLY : RII0 |
---|
| 159 | USE YOERAD , ONLY : NLW, NAER, RCCNLND ,RCCNSEA |
---|
| 160 | USE YOERAD , ONLY : NAER, RCCNLND ,RCCNSEA |
---|
| 161 | USE YOERDU , ONLY : REPSCQ |
---|
| 162 | USE YOMGEM , ONLY : NGPTOT |
---|
| 163 | USE YOERDI , ONLY : RRAE ,REPCLC ,REPH2O |
---|
| 164 | USE YOMARPHY , ONLY : LRDUST |
---|
| 165 | USE phys_output_mod, ONLY : swaerofree_diag, swaero_diag |
---|
[5294] | 166 | ! Temporary fix waiting for cleaner interface (or not) |
---|
| 167 | USE clesphys_mod_h, ONLY: NSW |
---|
| 168 | |
---|
[4875] | 169 | !----------------------------------------------------------------------- |
---|
[2004] | 170 | |
---|
[4875] | 171 | !* 0.1 ARGUMENTS. |
---|
| 172 | ! ---------- |
---|
[2004] | 173 | |
---|
[4875] | 174 | IMPLICIT NONE |
---|
[5294] | 175 | !! INCLUDE "clesphys.h" |
---|
[2004] | 176 | |
---|
[4875] | 177 | INTEGER(KIND=JPIM),INTENT(IN) :: KPROMA |
---|
| 178 | INTEGER(KIND=JPIM),INTENT(IN) :: KLEV |
---|
| 179 | INTEGER(KIND=JPIM),INTENT(IN) :: KST |
---|
| 180 | INTEGER(KIND=JPIM),INTENT(IN) :: KEND |
---|
| 181 | INTEGER(KIND=JPIM) :: KTDIA ! Argument NOT used |
---|
| 182 | INTEGER(KIND=JPIM),INTENT(IN) :: KMODE |
---|
| 183 | REAL(KIND=JPRB) ,INTENT(IN) :: PALBD(KPROMA,NSW) |
---|
| 184 | REAL(KIND=JPRB) ,INTENT(IN) :: PALBP(KPROMA,NSW) |
---|
| 185 | REAL(KIND=JPRB) ,INTENT(IN) :: PAPRS(KPROMA,KLEV+1) |
---|
| 186 | REAL(KIND=JPRB) ,INTENT(IN) :: PAPRSF(KPROMA,KLEV) |
---|
| 187 | REAL(KIND=JPRB) ,INTENT(IN) :: PCCO2 |
---|
| 188 | REAL(KIND=JPRB) ,INTENT(IN) :: PCLFR(KPROMA,KLEV) |
---|
| 189 | REAL(KIND=JPRB) ,INTENT(IN) :: PQO3(KPROMA,KLEV) |
---|
| 190 | REAL(KIND=JPRB) ,INTENT(IN) :: PAER(KPROMA,KLEV,6) |
---|
| 191 | REAL(KIND=JPRB) ,INTENT(IN) :: PDP(KPROMA,KLEV) |
---|
| 192 | REAL(KIND=JPRB) ,INTENT(IN) :: PEMIS(KPROMA) |
---|
| 193 | REAL(KIND=JPRB) ,INTENT(IN) :: PMU0(KPROMA) |
---|
| 194 | REAL(KIND=JPRB) ,INTENT(IN) :: PQ(KPROMA,KLEV) |
---|
| 195 | REAL(KIND=JPRB) ,INTENT(IN) :: PQS(KPROMA,KLEV) |
---|
| 196 | REAL(KIND=JPRB) ,INTENT(IN) :: PQIWP(KPROMA,KLEV) |
---|
| 197 | REAL(KIND=JPRB) ,INTENT(IN) :: PQLWP(KPROMA,KLEV) |
---|
| 198 | REAL(KIND=JPRB) ,INTENT(IN) :: PSLM(KPROMA) |
---|
| 199 | REAL(KIND=JPRB) ,INTENT(IN) :: PT(KPROMA,KLEV) |
---|
| 200 | REAL(KIND=JPRB) ,INTENT(IN) :: PTS(KPROMA) |
---|
| 201 | REAL(KIND=JPRB) ,INTENT(IN) :: PPIZA_TOT(KPROMA,KLEV,NSW) |
---|
| 202 | REAL(KIND=JPRB) ,INTENT(IN) :: PCGA_TOT(KPROMA,KLEV,NSW) |
---|
| 203 | REAL(KIND=JPRB) ,INTENT(IN) :: PTAU_TOT(KPROMA,KLEV,NSW) |
---|
| 204 | !--OB |
---|
| 205 | REAL(KIND=JPRB) ,INTENT(IN) :: PPIZA_NAT(KPROMA,KLEV,NSW) |
---|
| 206 | REAL(KIND=JPRB) ,INTENT(IN) :: PCGA_NAT(KPROMA,KLEV,NSW) |
---|
| 207 | REAL(KIND=JPRB) ,INTENT(IN) :: PTAU_NAT(KPROMA,KLEV,NSW) |
---|
| 208 | REAL(KIND=JPRB) :: PPIZA_ZERO(KPROMA,KLEV,NSW) |
---|
| 209 | REAL(KIND=JPRB) :: PCGA_ZERO(KPROMA,KLEV,NSW) |
---|
| 210 | REAL(KIND=JPRB) :: PTAU_ZERO(KPROMA,KLEV,NSW) |
---|
| 211 | !--fin |
---|
| 212 | !--C.Kleinschmitt |
---|
| 213 | REAL(KIND=JPRB) :: PTAU_LW_ZERO(KPROMA,KLEV,NLW) |
---|
| 214 | REAL(KIND=JPRB) ,INTENT(IN) :: PTAU_LW_TOT(KPROMA,KLEV,NLW) |
---|
| 215 | REAL(KIND=JPRB) ,INTENT(IN) :: PTAU_LW_NAT(KPROMA,KLEV,NLW) |
---|
| 216 | !--end |
---|
| 217 | REAL(KIND=JPRB) ,INTENT(IN) :: PREF_LIQ(KPROMA,KLEV) |
---|
| 218 | REAL(KIND=JPRB) ,INTENT(IN) :: PREF_ICE(KPROMA,KLEV) |
---|
| 219 | !--OB |
---|
| 220 | REAL(KIND=JPRB) ,INTENT(IN) :: PREF_LIQ_PI(KPROMA,KLEV) |
---|
| 221 | REAL(KIND=JPRB) ,INTENT(IN) :: PREF_ICE_PI(KPROMA,KLEV) |
---|
| 222 | LOGICAL, INTENT(in) :: ok_ade, ok_aie ! switches whether to use aerosol direct (indirect) effects or not |
---|
| 223 | LOGICAL, INTENT(in) :: ok_volcan ! produce volcanic diags (SW/LW heat flux and rate) |
---|
| 224 | INTEGER, INTENT(in) :: flag_aerosol ! takes value 0 (no aerosol) or 1 to 6 (aerosols) |
---|
| 225 | LOGICAL, INTENT(in) :: flag_aerosol_strat ! use stratospheric aerosols |
---|
| 226 | LOGICAL, INTENT(in) :: flag_aer_feedback ! use aerosols radiative feedback |
---|
| 227 | REAL(KIND=JPRB) ,INTENT(out) :: PTOPSWADAERO(KPROMA), PSOLSWADAERO(KPROMA) ! Aerosol direct forcing at TOA and surface |
---|
| 228 | REAL(KIND=JPRB) ,INTENT(OUT) :: PTOPSWAD0AERO(KPROMA), PSOLSWAD0AERO(KPROMA) ! Aerosol direct forcing at TOA and surface |
---|
| 229 | REAL(KIND=JPRB) ,INTENT(OUT) :: PTOPSWAIAERO(KPROMA), PSOLSWAIAERO(KPROMA) ! ditto, indirect |
---|
| 230 | REAL(KIND=JPRB) ,INTENT(OUT) :: PTOPSWCFAERO(KPROMA,3), PSOLSWCFAERO(KPROMA,3) !--do we keep this ? |
---|
| 231 | !--fin |
---|
| 232 | !--NL |
---|
| 233 | REAL(KIND=JPRB) ,INTENT(OUT) :: PSWADAERO(KPROMA, KLEV+1) ! SW Aerosol direct forcing |
---|
| 234 | REAL(KIND=JPRB) ,INTENT(OUT) :: PLWADAERO(KPROMA, KLEV+1) ! LW Aerosol direct forcing |
---|
| 235 | !--CK |
---|
| 236 | REAL(KIND=JPRB) ,INTENT(out) :: PTOPLWADAERO(KPROMA), PSOLLWADAERO(KPROMA) ! LW Aerosol direct forcing at TOA + surface |
---|
| 237 | REAL(KIND=JPRB) ,INTENT(OUT) :: PTOPLWAD0AERO(KPROMA), PSOLLWAD0AERO(KPROMA) ! LW Aerosol direct forcing at TOA + surface |
---|
| 238 | REAL(KIND=JPRB) ,INTENT(OUT) :: PTOPLWAIAERO(KPROMA), PSOLLWAIAERO(KPROMA) ! LW Aer. indirect forcing at TOA + surface |
---|
| 239 | !--end |
---|
| 240 | REAL(KIND=JPRB) ,INTENT(OUT) :: PEMTD(KPROMA,KLEV+1) |
---|
| 241 | REAL(KIND=JPRB) ,INTENT(OUT) :: PEMTU(KPROMA,KLEV+1) |
---|
| 242 | REAL(KIND=JPRB) ,INTENT(OUT) :: PTRSO(KPROMA,KLEV+1) |
---|
| 243 | REAL(KIND=JPRB) ,INTENT(INOUT) :: PTH(KPROMA,KLEV+1) |
---|
| 244 | REAL(KIND=JPRB) ,INTENT(OUT) :: PCTRSO(KPROMA,2) |
---|
| 245 | REAL(KIND=JPRB) ,INTENT(OUT) :: PCEMTR(KPROMA,2) |
---|
| 246 | REAL(KIND=JPRB) ,INTENT(OUT) :: PTRSOD(KPROMA) |
---|
| 247 | REAL(KIND=JPRB) ,INTENT(OUT) :: PLWFC(KPROMA,2) |
---|
| 248 | REAL(KIND=JPRB) ,INTENT(OUT) :: PLWFT(KPROMA,KLEV+1) |
---|
| 249 | REAL(KIND=JPRB) ,INTENT(OUT) :: PSWFC(KPROMA,2) |
---|
| 250 | REAL(KIND=JPRB) ,INTENT(OUT) :: PSWFT(KPROMA,KLEV+1) |
---|
| 251 | REAL(KIND=JPRB) ,INTENT(OUT) :: PSFSWDIR(KPROMA,NSW) |
---|
| 252 | REAL(KIND=JPRB) ,INTENT(OUT) :: PSFSWDIF(KPROMA,NSW) |
---|
| 253 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSDNN(KPROMA) |
---|
| 254 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSDNV(KPROMA) |
---|
| 255 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFLUX(KPROMA,2,KLEV+1) ! LW total sky flux (1=up, 2=down) |
---|
| 256 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFLUC(KPROMA,2,KLEV+1) ! LW clear sky flux (1=up, 2=down) |
---|
| 257 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSDN(KPROMA,KLEV+1) ! SW total sky flux down |
---|
| 258 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSUP(KPROMA,KLEV+1) ! SW total sky flux up |
---|
| 259 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSCDN(KPROMA,KLEV+1) ! SW clear sky flux down |
---|
| 260 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSCUP(KPROMA,KLEV+1) ! SW clear sky flux up |
---|
| 261 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSCCDN(KPROMA,KLEV+1) ! SW clear sky clean (no aerosol) flux down |
---|
| 262 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFSCCUP(KPROMA,KLEV+1) ! SW clear sky clean (no aerosol) flux up |
---|
| 263 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFLCCDN(KPROMA,KLEV+1) ! LW clear sky clean (no aerosol) flux down |
---|
| 264 | REAL(KIND=JPRB) ,INTENT(OUT) :: PFLCCUP(KPROMA,KLEV+1) ! LW clear sky clean (no aerosol) flux up |
---|
| 265 | !--ajout VOLMIP |
---|
| 266 | REAL(KIND=JPRB) ,INTENT(OUT) :: volmip_solsw(KPROMA) ! SW clear sky in the case of VOLMIP |
---|
| 267 | INTEGER, INTENT(IN) :: flag_volc_surfstrat !--VOlMIP Modif |
---|
[2004] | 268 | |
---|
[4875] | 269 | ! ==== COMPUTED IN RADITE === |
---|
| 270 | ! ------------------------------------------------------------------ |
---|
| 271 | !* 0.2 LOCAL ARRAYS. |
---|
| 272 | ! ------------- |
---|
| 273 | REAL(KIND=JPRB) :: ZRAER (KPROMA,6,KLEV) |
---|
| 274 | REAL(KIND=JPRB) :: ZRCLC (KPROMA,KLEV) |
---|
| 275 | REAL(KIND=JPRB) :: ZRMU0 (KPROMA) |
---|
| 276 | REAL(KIND=JPRB) :: ZRPR (KPROMA,KLEV) |
---|
| 277 | REAL(KIND=JPRB) :: ZRTI (KPROMA,KLEV) |
---|
| 278 | REAL(KIND=JPRB) :: ZQLWP (KPROMA,KLEV ) , ZQIWP (KPROMA,KLEV ) |
---|
[2004] | 279 | |
---|
[4875] | 280 | REAL(KIND=JPRB) :: ZPQO3 (KPROMA,KLEV) |
---|
| 281 | REAL(KIND=JPRB) :: ZQOZ (NGPTOT,KLEV) |
---|
| 282 | REAL(KIND=JPRB) :: ZQS (KPROMA,KLEV) |
---|
| 283 | REAL(KIND=JPRB) :: ZQ (KPROMA,KLEV) |
---|
| 284 | REAL(KIND=JPRB) :: ZEMTD (KPROMA,KLEV+1) |
---|
| 285 | REAL(KIND=JPRB) :: ZEMTU (KPROMA,KLEV+1) |
---|
| 286 | REAL(KIND=JPRB) :: ZTRSOC (KPROMA,2) |
---|
| 287 | REAL(KIND=JPRB) :: ZEMTC (KPROMA,2) |
---|
[2004] | 288 | |
---|
[4875] | 289 | REAL(KIND=JPRB) :: ZNBAS (KPROMA) |
---|
| 290 | REAL(KIND=JPRB) :: ZNTOP (KPROMA) |
---|
| 291 | REAL(KIND=JPRB) :: ZQRAIN (KPROMA,KLEV) |
---|
| 292 | REAL(KIND=JPRB) :: ZQRAINT(KPROMA,KLEV) |
---|
| 293 | REAL(KIND=JPRB) :: ZCCNL (KPROMA) |
---|
| 294 | REAL(KIND=JPRB) :: ZCCNO (KPROMA) |
---|
[2004] | 295 | |
---|
[4875] | 296 | ! output of radlsw |
---|
[2004] | 297 | |
---|
[4875] | 298 | REAL(KIND=JPRB) :: ZEMIT (KPROMA) |
---|
| 299 | REAL(KIND=JPRB) :: ZFCT (KPROMA,KLEV+1) |
---|
| 300 | REAL(KIND=JPRB) :: ZFLT (KPROMA,KLEV+1) |
---|
| 301 | REAL(KIND=JPRB) :: ZFCS (KPROMA,KLEV+1) |
---|
| 302 | REAL(KIND=JPRB) :: ZFLS (KPROMA,KLEV+1) |
---|
| 303 | REAL(KIND=JPRB) :: ZFRSOD (KPROMA),ZSUDU(KPROMA) |
---|
| 304 | REAL(KIND=JPRB) :: ZPARF (KPROMA),ZUVDF(KPROMA),ZPARCF(KPROMA),ZTINCF(KPROMA) |
---|
[2004] | 305 | |
---|
[4875] | 306 | INTEGER(KIND=JPIM) :: IBEG, IEND, JK, JL |
---|
[2004] | 307 | |
---|
[4875] | 308 | REAL(KIND=JPRB) :: ZCRAE, ZRII0, ZEMIW(KPROMA) |
---|
| 309 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
---|
[2004] | 310 | |
---|
[4875] | 311 | !---aerosol radiative diagnostics |
---|
| 312 | ! Key to define the aerosol effect acting on climate |
---|
| 313 | ! OB: AEROSOLFEEDBACK_ACTIVE is now a LOGICAL |
---|
| 314 | ! TRUE: fluxes use natural and/or anthropogenic aerosols according to ok_ade and ok_aie, DEFAULT |
---|
| 315 | ! FALSE: fluxes use no aerosols (case 1) |
---|
| 316 | ! to be used only for maintaining bit reproducibility with aerosol diagnostics activated |
---|
| 317 | LOGICAL :: AEROSOLFEEDBACK_ACTIVE ! now externalized from .def files |
---|
[2004] | 318 | |
---|
[4875] | 319 | !OB - Fluxes including aerosol effects |
---|
| 320 | ! | direct effect |
---|
| 321 | !ind effect | no aerosol NATural TOTal |
---|
| 322 | !standard | 5 |
---|
| 323 | !natural (PI) | 1 3 |
---|
| 324 | !total (PD) | 2 4 |
---|
| 325 | ! so we need which case when ? |
---|
| 326 | ! if flag_aerosol is on |
---|
| 327 | ! ok_ade and ok_aie = 4-2, 4-3 and 4 to proceed |
---|
| 328 | ! ok_ade and not ok_aie = 3-1 and 3 to proceed |
---|
| 329 | ! not ok_ade and ok_aie = 2-1 and 2 to proceed |
---|
| 330 | ! not ok_ade and not ok_aie = 1 to proceed |
---|
| 331 | ! therefore the cases have the following corresponding switches |
---|
| 332 | ! 1 = not ok_ade and not ok_aie OR not ok_ade and ok_aie and swaero_diag OR ok_ade and not ok_aie and swaero_diag |
---|
| 333 | ! 2 = not ok_ade and ok_aie OR ok_aie and ok_ade and swaero_diag |
---|
| 334 | ! 3 = ok_ade and not ok_aie OR ok_aie and ok_ade and swaero_diag |
---|
| 335 | ! 4 = ok_ade and ok_aie |
---|
| 336 | ! 5 = no aerosol feedback wanted or no aerosol at all |
---|
| 337 | ! if they are called in this order then the correct call is used to proceed |
---|
[2004] | 338 | |
---|
[4875] | 339 | REAL(KIND=JPRB) :: ZFSUP_AERO(KPROMA,KLEV+1,5) |
---|
| 340 | REAL(KIND=JPRB) :: ZFSDN_AERO(KPROMA,KLEV+1,5) |
---|
| 341 | REAL(KIND=JPRB) :: ZFSUP0_AERO(KPROMA,KLEV+1,5) |
---|
| 342 | REAL(KIND=JPRB) :: ZFSDN0_AERO(KPROMA,KLEV+1,5) |
---|
| 343 | !--LW (CK): |
---|
| 344 | REAL(KIND=JPRB) :: LWUP_AERO(KPROMA,KLEV+1,5) |
---|
| 345 | REAL(KIND=JPRB) :: LWDN_AERO(KPROMA,KLEV+1,5) |
---|
| 346 | REAL(KIND=JPRB) :: LWUP0_AERO(KPROMA,KLEV+1,5) |
---|
| 347 | REAL(KIND=JPRB) :: LWDN0_AERO(KPROMA,KLEV+1,5) |
---|
[2004] | 348 | |
---|
| 349 | #include "radlsw.intfb.h" |
---|
| 350 | |
---|
[4875] | 351 | IF (LHOOK) CALL DR_HOOK('RECMWF_AERO',0,ZHOOK_HANDLE) |
---|
| 352 | IBEG=KST |
---|
| 353 | IEND=KEND |
---|
[2004] | 354 | |
---|
[4875] | 355 | AEROSOLFEEDBACK_ACTIVE = flag_aer_feedback !NL: externalize aer feedback |
---|
[3412] | 356 | |
---|
| 357 | |
---|
[4875] | 358 | !* 1. PREPARATORY WORK |
---|
| 359 | ! ---------------- |
---|
| 360 | !--OB |
---|
| 361 | ! 1.0 INITIALIZATIONS |
---|
| 362 | ! -------------- |
---|
[2004] | 363 | |
---|
[4875] | 364 | ZFSUP_AERO (:,:,:)=0. |
---|
| 365 | ZFSDN_AERO (:,:,:)=0. |
---|
| 366 | ZFSUP0_AERO(:,:,:)=0. |
---|
| 367 | ZFSDN0_AERO(:,:,:)=0. |
---|
[2004] | 368 | |
---|
[4875] | 369 | LWUP_AERO (:,:,:)=0. |
---|
| 370 | LWDN_AERO (:,:,:)=0. |
---|
| 371 | LWUP0_AERO(:,:,:)=0. |
---|
| 372 | LWDN0_AERO(:,:,:)=0. |
---|
[2146] | 373 | |
---|
[4875] | 374 | PTAU_ZERO(:,:,:) =1.e-15 |
---|
| 375 | PPIZA_ZERO(:,:,:)=1.0 |
---|
| 376 | PCGA_ZERO(:,:,:) =0.0 |
---|
[2004] | 377 | |
---|
[4875] | 378 | PTAU_LW_ZERO(:,:,:) =1.e-15 |
---|
[2004] | 379 | |
---|
[2146] | 380 | |
---|
[4875] | 381 | !* 1.1 LOCAL CONSTANTS |
---|
| 382 | ! --------------- |
---|
[2004] | 383 | |
---|
[4875] | 384 | ZRII0=RII0 |
---|
| 385 | ZCRAE=RRAE*(RRAE+2.0_JPRB) |
---|
[2004] | 386 | |
---|
[4875] | 387 | !* 2.1 FULL-LEVEL QUANTITIES |
---|
[2004] | 388 | |
---|
[4875] | 389 | ZRPR =PAPRSF |
---|
[2004] | 390 | |
---|
[4875] | 391 | DO JK=1,KLEV |
---|
| 392 | DO JL=IBEG,IEND |
---|
| 393 | ! ZPQO3(JL,JK)=PQO3(JL,JK)*PDP(JL,JK)*RMD/RMO3 |
---|
| 394 | ZPQO3(JL,JK)=PQO3(JL,JK)*PDP(JL,JK) |
---|
| 395 | ZRCLC(JL,JK)=MAX( 0.0_JPRB ,MIN( 1.0_JPRB ,PCLFR(JL,JK))) |
---|
| 396 | IF (ZRCLC(JL,JK) > REPCLC) THEN |
---|
| 397 | ZQLWP(JL,JK)=PQLWP(JL,JK) |
---|
| 398 | ZQIWP(JL,JK)=PQIWP(JL,JK) |
---|
| 399 | ELSE |
---|
| 400 | ZQLWP(JL,JK)=REPH2O*ZRCLC(JL,JK) |
---|
| 401 | ZQIWP(JL,JK)=REPH2O*ZRCLC(JL,JK) |
---|
| 402 | ENDIF |
---|
| 403 | ZQRAIN(JL,JK)=0. |
---|
| 404 | ZQRAINT(JL,JK)=0. |
---|
| 405 | ZRTI(JL,JK) =PT(JL,JK) |
---|
| 406 | ZQS (JL,JK)=MAX(2.0_JPRB*REPH2O,PQS(JL,JK)) |
---|
| 407 | ZQ (JL,JK)=MAX(REPH2O,MIN(PQ(JL,JK),ZQS(JL,JK)*(1.0_JPRB-REPH2O))) |
---|
| 408 | ZEMIW(JL)=PEMIS(JL) |
---|
| 409 | ENDDO |
---|
[2004] | 410 | ENDDO |
---|
| 411 | |
---|
[4875] | 412 | IF (NAER == 0) THEN |
---|
| 413 | ZRAER=RCAEROS |
---|
| 414 | ELSE |
---|
| 415 | DO JK=1,KLEV |
---|
| 416 | DO JL=IBEG,IEND |
---|
| 417 | ZRAER(JL,1,JK)=PAER(JL,JK,1) |
---|
| 418 | ZRAER(JL,2,JK)=PAER(JL,JK,2) |
---|
| 419 | ZRAER(JL,3,JK)=PAER(JL,JK,3) |
---|
| 420 | ZRAER(JL,4,JK)=PAER(JL,JK,4) |
---|
| 421 | ZRAER(JL,5,JK)=RCAEROS |
---|
| 422 | ZRAER(JL,6,JK)=PAER(JL,JK,6) |
---|
| 423 | ENDDO |
---|
| 424 | ENDDO |
---|
| 425 | ENDIF |
---|
| 426 | |
---|
| 427 | !* 2.2 HALF-LEVEL QUANTITIES |
---|
| 428 | |
---|
| 429 | DO JK=2,KLEV |
---|
| 430 | DO JL=IBEG,IEND |
---|
| 431 | PTH(JL,JK)=& |
---|
| 432 | & (PT(JL,JK-1)*PAPRSF(JL,JK-1)*(PAPRSF(JL,JK)-PAPRS(JL,JK))& |
---|
| 433 | & +PT(JL,JK)*PAPRSF(JL,JK)*(PAPRS(JL,JK)-PAPRSF(JL,JK-1)))& |
---|
| 434 | & *(1.0_JPRB/(PAPRS(JL,JK)*(PAPRSF(JL,JK)-PAPRSF(JL,JK-1)))) |
---|
| 435 | ENDDO |
---|
[2004] | 436 | ENDDO |
---|
| 437 | |
---|
[4875] | 438 | !* 2.3 QUANTITIES AT BOUNDARIES |
---|
[2004] | 439 | |
---|
| 440 | DO JL=IBEG,IEND |
---|
[4875] | 441 | PTH(JL,KLEV+1)=PTS(JL) |
---|
| 442 | PTH(JL,1)=PT(JL,1)-PAPRSF(JL,1)*(PT(JL,1)-PTH(JL,2))& |
---|
| 443 | & /(PAPRSF(JL,1)-PAPRS(JL,2)) |
---|
| 444 | ZNBAS(JL)=1. |
---|
| 445 | ZNTOP(JL)=1. |
---|
| 446 | ZCCNL(JL)=RCCNLND |
---|
| 447 | ZCCNO(JL)=RCCNSEA |
---|
[2004] | 448 | ENDDO |
---|
| 449 | |
---|
[4875] | 450 | !* 3.1 SOLAR ZENITH ANGLE IS EARTH'S CURVATURE |
---|
| 451 | ! CORRECTED |
---|
[2004] | 452 | |
---|
[4875] | 453 | ! CCMVAL: on impose ZRMU0=PMU0 MPL 25032010 |
---|
| 454 | ! 2eme essai en 3D MPL 20052010 |
---|
| 455 | !DO JL=IBEG,IEND |
---|
| 456 | ! ZRMU0(JL)=PMU0(JL) |
---|
| 457 | !ENDDO |
---|
| 458 | !!!!! A REVOIR MPL 20091201: enleve cette correction pour comparer a AR4 |
---|
| 459 | DO JL=IBEG,IEND |
---|
| 460 | IF (PMU0(JL) > 1.E-10_JPRB) THEN |
---|
| 461 | ZRMU0(JL)=RRAE/(SQRT(PMU0(JL)**2+ZCRAE)-PMU0(JL)) |
---|
| 462 | ELSE |
---|
| 463 | ZRMU0(JL)= RRAE/SQRT(ZCRAE) |
---|
| 464 | ENDIF |
---|
| 465 | ENDDO |
---|
[2004] | 466 | |
---|
[4875] | 467 | !* 4.1 CALL TO ACTUAL RADIATION SCHEME |
---|
| 468 | ! |
---|
| 469 | !----now we make multiple calls to the radiation according to which |
---|
| 470 | !----aerosol flags are on |
---|
[2004] | 471 | |
---|
[4875] | 472 | IF (flag_aerosol .GT. 0 .OR. flag_aerosol_strat) THEN |
---|
[2004] | 473 | |
---|
[4875] | 474 | !--Case 1 |
---|
| 475 | IF ( ( .not. ok_ade .AND. .not. ok_aie ) .OR. & |
---|
| 476 | & ( .not. ok_ade .AND. ok_aie .AND. swaero_diag ) .OR. & |
---|
| 477 | & ( ok_ade .AND. .not. ok_aie .AND. swaero_diag ) ) THEN |
---|
[2004] | 478 | |
---|
[4875] | 479 | ! natural aerosols for direct and indirect effect |
---|
| 480 | ! PI cloud optical properties |
---|
| 481 | ! use PREF_LIQ_PI and PREF_ICE_PI |
---|
| 482 | ! use NAT aerosol optical properties |
---|
| 483 | ! store fluxes in index 1 |
---|
[2004] | 484 | |
---|
[4875] | 485 | CALL RADLSW (& |
---|
| 486 | & IBEG , IEND , KPROMA , KLEV , KMODE , NAER,& |
---|
| 487 | & ZRII0 ,& |
---|
| 488 | & ZRAER , PALBD , PALBP , PAPRS , ZRPR ,& |
---|
| 489 | & ZCCNL , ZCCNO ,& |
---|
| 490 | & PCCO2 , ZRCLC , PDP , PEMIS , ZEMIW ,PSLM , ZRMU0 , ZPQO3,& |
---|
| 491 | & ZQ , ZQIWP , ZQLWP , ZQS , ZQRAIN,ZQRAINT ,& |
---|
| 492 | & PTH , ZRTI , PTS , ZNBAS , ZNTOP ,& |
---|
| 493 | & PREF_LIQ_PI, PREF_ICE_PI,& |
---|
| 494 | & ZEMIT , ZFCT , ZFLT , ZFCS , ZFLS ,& |
---|
| 495 | & ZFRSOD, ZSUDU , ZUVDF , ZPARF , ZPARCF, ZTINCF, PSFSWDIR,& |
---|
| 496 | & PSFSWDIF,PFSDNN, PFSDNV ,& |
---|
| 497 | & LRDUST,PPIZA_NAT,PCGA_NAT,PTAU_NAT,PTAU_LW_NAT,PFLUX,PFLUC,& |
---|
| 498 | & PFSDN , PFSUP , PFSCDN , PFSCUP ) |
---|
[2004] | 499 | |
---|
[4875] | 500 | !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4 |
---|
| 501 | ZFSUP0_AERO(:,:,1) = PFSCUP(:,:) |
---|
| 502 | ZFSDN0_AERO(:,:,1) = PFSCDN(:,:) |
---|
[2004] | 503 | |
---|
[4875] | 504 | ZFSUP_AERO(:,:,1) = PFSUP(:,:) |
---|
| 505 | ZFSDN_AERO(:,:,1) = PFSDN(:,:) |
---|
[2004] | 506 | |
---|
[4875] | 507 | LWUP0_AERO(:,:,1) = PFLUC(:,1,:) |
---|
| 508 | LWDN0_AERO(:,:,1) = PFLUC(:,2,:) |
---|
[2004] | 509 | |
---|
[4875] | 510 | LWUP_AERO(:,:,1) = PFLUX(:,1,:) |
---|
| 511 | LWDN_AERO(:,:,1) = PFLUX(:,2,:) |
---|
[2004] | 512 | |
---|
[4875] | 513 | ENDIF |
---|
[2146] | 514 | |
---|
[4875] | 515 | !--Case 2 |
---|
| 516 | IF ( ( .not. ok_ade .AND. ok_aie ) .OR. & |
---|
| 517 | & ( ok_ade .AND. ok_aie .AND. swaero_diag ) ) THEN |
---|
[2146] | 518 | |
---|
[4875] | 519 | ! natural aerosols for direct indirect effect |
---|
| 520 | ! use NAT aerosol optical properties |
---|
| 521 | ! PD cloud optical properties |
---|
| 522 | ! use PREF_LIQ and PREF_ICE |
---|
| 523 | ! store fluxes in index 2 |
---|
[2004] | 524 | |
---|
[4875] | 525 | CALL RADLSW (& |
---|
| 526 | & IBEG , IEND , KPROMA , KLEV , KMODE , NAER,& |
---|
| 527 | & ZRII0 ,& |
---|
| 528 | & ZRAER , PALBD , PALBP , PAPRS , ZRPR ,& |
---|
| 529 | & ZCCNL , ZCCNO ,& |
---|
| 530 | & PCCO2 , ZRCLC , PDP , PEMIS , ZEMIW ,PSLM , ZRMU0 , ZPQO3,& |
---|
| 531 | & ZQ , ZQIWP , ZQLWP , ZQS , ZQRAIN,ZQRAINT ,& |
---|
| 532 | & PTH , ZRTI , PTS , ZNBAS , ZNTOP ,& |
---|
| 533 | & PREF_LIQ, PREF_ICE,& |
---|
| 534 | & ZEMIT , ZFCT , ZFLT , ZFCS , ZFLS ,& |
---|
| 535 | & ZFRSOD, ZSUDU , ZUVDF , ZPARF , ZPARCF, ZTINCF, PSFSWDIR,& |
---|
| 536 | & PSFSWDIF,PFSDNN, PFSDNV ,& |
---|
| 537 | & LRDUST,PPIZA_NAT,PCGA_NAT,PTAU_NAT,PTAU_LW_NAT,PFLUX,PFLUC,& |
---|
| 538 | & PFSDN , PFSUP , PFSCDN , PFSCUP ) |
---|
[2004] | 539 | |
---|
[4875] | 540 | !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4 |
---|
| 541 | ZFSUP0_AERO(:,:,2) = PFSCUP(:,:) |
---|
| 542 | ZFSDN0_AERO(:,:,2) = PFSCDN(:,:) |
---|
[2004] | 543 | |
---|
[4875] | 544 | ZFSUP_AERO(:,:,2) = PFSUP(:,:) |
---|
| 545 | ZFSDN_AERO(:,:,2) = PFSDN(:,:) |
---|
[2004] | 546 | |
---|
[4875] | 547 | LWUP0_AERO(:,:,2) = PFLUC(:,1,:) |
---|
| 548 | LWDN0_AERO(:,:,2) = PFLUC(:,2,:) |
---|
[2004] | 549 | |
---|
[4875] | 550 | LWUP_AERO(:,:,2) = PFLUX(:,1,:) |
---|
| 551 | LWDN_AERO(:,:,2) = PFLUX(:,2,:) |
---|
[2004] | 552 | |
---|
[4875] | 553 | ENDIF ! ok_aie |
---|
[2146] | 554 | |
---|
[4875] | 555 | !--Case 3 |
---|
| 556 | IF ( ( ok_ade .AND. .not. ok_aie ) .OR. & |
---|
| 557 | & ( ok_ade .AND. ok_aie .AND. swaero_diag ) ) THEN |
---|
[2146] | 558 | |
---|
[4875] | 559 | ! direct effect of total aerosol activated |
---|
| 560 | ! TOT aerosols for direct effect |
---|
| 561 | ! PI cloud optical properties |
---|
| 562 | ! use PREF_LIQ_PI and PREF_ICE_PI |
---|
| 563 | ! STORE fluxes in index 3 |
---|
[2004] | 564 | |
---|
[4875] | 565 | CALL RADLSW (& |
---|
| 566 | & IBEG , IEND , KPROMA , KLEV , KMODE , NAER,& |
---|
| 567 | & ZRII0 ,& |
---|
| 568 | & ZRAER , PALBD , PALBP , PAPRS , ZRPR ,& |
---|
| 569 | & ZCCNL , ZCCNO ,& |
---|
| 570 | & PCCO2 , ZRCLC , PDP , PEMIS , ZEMIW ,PSLM , ZRMU0 , ZPQO3,& |
---|
| 571 | & ZQ , ZQIWP , ZQLWP , ZQS , ZQRAIN,ZQRAINT ,& |
---|
| 572 | & PTH , ZRTI , PTS , ZNBAS , ZNTOP ,& |
---|
| 573 | & PREF_LIQ_PI, PREF_ICE_PI,& |
---|
| 574 | & ZEMIT , ZFCT , ZFLT , ZFCS , ZFLS ,& |
---|
| 575 | & ZFRSOD, ZSUDU , ZUVDF , ZPARF , ZPARCF, ZTINCF, PSFSWDIR,& |
---|
| 576 | & PSFSWDIF,PFSDNN, PFSDNV ,& |
---|
| 577 | & LRDUST,PPIZA_TOT,PCGA_TOT,PTAU_TOT,PTAU_LW_TOT,PFLUX,PFLUC,& |
---|
| 578 | & PFSDN , PFSUP , PFSCDN , PFSCUP ) |
---|
[2004] | 579 | |
---|
[4875] | 580 | !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4 |
---|
| 581 | ZFSUP0_AERO(:,:,3) = PFSCUP(:,:) |
---|
| 582 | ZFSDN0_AERO(:,:,3) = PFSCDN(:,:) |
---|
[2004] | 583 | |
---|
[4875] | 584 | ZFSUP_AERO(:,:,3) = PFSUP(:,:) |
---|
| 585 | ZFSDN_AERO(:,:,3) = PFSDN(:,:) |
---|
[2004] | 586 | |
---|
[4875] | 587 | LWUP0_AERO(:,:,3) = PFLUC(:,1,:) |
---|
| 588 | LWDN0_AERO(:,:,3) = PFLUC(:,2,:) |
---|
[2004] | 589 | |
---|
[4875] | 590 | LWUP_AERO(:,:,3) = PFLUX(:,1,:) |
---|
| 591 | LWDN_AERO(:,:,3) = PFLUX(:,2,:) |
---|
[2146] | 592 | |
---|
[4875] | 593 | ENDIF !-end ok_ade |
---|
[2146] | 594 | |
---|
[4875] | 595 | !--Case 4 |
---|
| 596 | IF (ok_ade .and. ok_aie) THEN |
---|
[2004] | 597 | |
---|
[4875] | 598 | ! total aerosols for direct indirect effect |
---|
| 599 | ! use TOT aerosol optical properties |
---|
| 600 | ! PD cloud optical properties |
---|
| 601 | ! use PREF_LIQ and PREF_ICE |
---|
| 602 | ! store fluxes in index 4 |
---|
[2004] | 603 | |
---|
[4875] | 604 | CALL RADLSW (& |
---|
| 605 | & IBEG , IEND , KPROMA , KLEV , KMODE , NAER,& |
---|
| 606 | & ZRII0 ,& |
---|
| 607 | & ZRAER , PALBD , PALBP , PAPRS , ZRPR ,& |
---|
| 608 | & ZCCNL , ZCCNO ,& |
---|
| 609 | & PCCO2 , ZRCLC , PDP , PEMIS , ZEMIW ,PSLM , ZRMU0 , ZPQO3,& |
---|
| 610 | & ZQ , ZQIWP , ZQLWP , ZQS , ZQRAIN,ZQRAINT ,& |
---|
| 611 | & PTH , ZRTI , PTS , ZNBAS , ZNTOP ,& |
---|
| 612 | & PREF_LIQ, PREF_ICE,& |
---|
| 613 | & ZEMIT , ZFCT , ZFLT , ZFCS , ZFLS ,& |
---|
| 614 | & ZFRSOD, ZSUDU , ZUVDF , ZPARF , ZPARCF, ZTINCF, PSFSWDIR,& |
---|
| 615 | & PSFSWDIF,PFSDNN, PFSDNV ,& |
---|
| 616 | & LRDUST,PPIZA_TOT,PCGA_TOT,PTAU_TOT,PTAU_LW_TOT,PFLUX,PFLUC,& |
---|
| 617 | & PFSDN , PFSUP , PFSCDN , PFSCUP ) |
---|
[2004] | 618 | |
---|
[4875] | 619 | !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4 |
---|
| 620 | ZFSUP0_AERO(:,:,4) = PFSCUP(:,:) |
---|
| 621 | ZFSDN0_AERO(:,:,4) = PFSCDN(:,:) |
---|
[2004] | 622 | |
---|
[4875] | 623 | ZFSUP_AERO(:,:,4) = PFSUP(:,:) |
---|
| 624 | ZFSDN_AERO(:,:,4) = PFSDN(:,:) |
---|
[2004] | 625 | |
---|
[4875] | 626 | LWUP0_AERO(:,:,4) = PFLUC(:,1,:) |
---|
| 627 | LWDN0_AERO(:,:,4) = PFLUC(:,2,:) |
---|
[2004] | 628 | |
---|
[4875] | 629 | LWUP_AERO(:,:,4) = PFLUX(:,1,:) |
---|
| 630 | LWDN_AERO(:,:,4) = PFLUX(:,2,:) |
---|
[2146] | 631 | |
---|
[4875] | 632 | ENDIF ! ok_ade .and. ok_aie |
---|
[2146] | 633 | |
---|
[4875] | 634 | ENDIF !--if flag_aerosol GT 0 OR flag_aerosol_strat |
---|
[2004] | 635 | |
---|
[4875] | 636 | ! case with no aerosols at all is also computed IF ACTIVEFEEDBACK_ACTIVE is false |
---|
| 637 | IF (.not. AEROSOLFEEDBACK_ACTIVE .OR. flag_aerosol .EQ. 0 .OR. swaerofree_diag) THEN |
---|
[2004] | 638 | |
---|
[4875] | 639 | ! ZERO aerosol effect |
---|
| 640 | ! ZERO aerosol optical depth |
---|
| 641 | ! STANDARD cloud optical properties |
---|
| 642 | ! STORE fluxes in index 5 |
---|
[2004] | 643 | |
---|
[4875] | 644 | CALL RADLSW (& |
---|
| 645 | & IBEG , IEND , KPROMA , KLEV , KMODE , NAER,& |
---|
| 646 | & ZRII0 ,& |
---|
| 647 | & ZRAER , PALBD , PALBP , PAPRS , ZRPR ,& |
---|
| 648 | & ZCCNL , ZCCNO ,& |
---|
| 649 | & PCCO2 , ZRCLC , PDP , PEMIS , ZEMIW ,PSLM , ZRMU0 , ZPQO3,& |
---|
| 650 | & ZQ , ZQIWP , ZQLWP , ZQS , ZQRAIN,ZQRAINT ,& |
---|
| 651 | & PTH , ZRTI , PTS , ZNBAS , ZNTOP ,& |
---|
| 652 | !--this needs to be changed to fixed cloud optical properties |
---|
| 653 | & PREF_LIQ_PI, PREF_ICE_PI,& |
---|
| 654 | & ZEMIT , ZFCT , ZFLT , ZFCS , ZFLS ,& |
---|
| 655 | & ZFRSOD, ZSUDU , ZUVDF , ZPARF , ZPARCF, ZTINCF, PSFSWDIR,& |
---|
| 656 | & PSFSWDIF,PFSDNN, PFSDNV ,& |
---|
| 657 | & LRDUST,PPIZA_ZERO,PCGA_ZERO,PTAU_ZERO, PTAU_LW_ZERO,PFLUX,PFLUC,& |
---|
| 658 | & PFSDN , PFSUP , PFSCDN , PFSCUP ) |
---|
[2004] | 659 | |
---|
[4875] | 660 | !* SAVE VARIABLES IN INTERIM VARIABLES A LA SW_AEROAR4 |
---|
| 661 | ZFSUP0_AERO(:,:,5) = PFSCUP(:,:) |
---|
| 662 | ZFSDN0_AERO(:,:,5) = PFSCDN(:,:) |
---|
[2004] | 663 | |
---|
[4875] | 664 | ZFSUP_AERO(:,:,5) = PFSUP(:,:) |
---|
| 665 | ZFSDN_AERO(:,:,5) = PFSDN(:,:) |
---|
[2004] | 666 | |
---|
[4875] | 667 | LWUP0_AERO(:,:,5) = PFLUC(:,1,:) |
---|
| 668 | LWDN0_AERO(:,:,5) = PFLUC(:,2,:) |
---|
[2004] | 669 | |
---|
[4875] | 670 | LWUP_AERO(:,:,5) = PFLUX(:,1,:) |
---|
| 671 | LWDN_AERO(:,:,5) = PFLUX(:,2,:) |
---|
[2146] | 672 | |
---|
[4875] | 673 | ENDIF ! .not. AEROSOLFEEDBACK_ACTIVE |
---|
[2146] | 674 | |
---|
[4875] | 675 | !* 4.2 TRANSFORM FLUXES TO MODEL HISTORICAL VARIABLES |
---|
[2004] | 676 | |
---|
[4875] | 677 | DO JK=1,KLEV+1 |
---|
| 678 | DO JL=IBEG,IEND |
---|
| 679 | PSWFT(JL,JK)=ZFLS(JL,JK)/(ZRII0*ZRMU0(JL)) |
---|
| 680 | PLWFT(JL,JK)=ZFLT(JL,JK) |
---|
| 681 | ENDDO |
---|
| 682 | ENDDO |
---|
[2004] | 683 | |
---|
[4875] | 684 | ZEMTD=PLWFT |
---|
| 685 | ZEMTU=PLWFT |
---|
| 686 | |
---|
[2004] | 687 | DO JL=IBEG,IEND |
---|
[4875] | 688 | ZTRSOC(JL, 1)=ZFCS(JL, 1)/(ZRII0*ZRMU0(JL)) |
---|
| 689 | ZTRSOC(JL, 2)=ZFCS(JL,KLEV+1)/(ZRII0*ZRMU0(JL)) |
---|
| 690 | ZEMTC (JL, 1)=ZFCT(JL, 1) |
---|
| 691 | ZEMTC (JL, 2)=ZFCT(JL,KLEV+1) |
---|
[2004] | 692 | ENDDO |
---|
| 693 | |
---|
[4875] | 694 | ! ------------ -- ------- -- ---- ----- |
---|
| 695 | !* 5.1 STORAGE OF TRANSMISSIVITY AND EMISSIVITIES |
---|
| 696 | !* IN KPROMA-LONG ARRAYS |
---|
[2004] | 697 | |
---|
[4875] | 698 | DO JK=1,KLEV+1 |
---|
| 699 | DO JL=IBEG,IEND |
---|
| 700 | PEMTD(JL,JK)=ZEMTD(JL,JK) |
---|
| 701 | PEMTU(JL,JK)=ZEMTU(JL,JK) |
---|
| 702 | PTRSO(JL,JK)=MAX(0.0_JPRB,MIN(1.0_JPRB,PSWFT(JL,JK))) |
---|
| 703 | ENDDO |
---|
[2004] | 704 | ENDDO |
---|
[4875] | 705 | DO JK=1,2 |
---|
| 706 | DO JL=IBEG,IEND |
---|
| 707 | PCEMTR(JL,JK)=ZEMTC (JL,JK) |
---|
| 708 | PCTRSO(JL,JK)=MAX( 0.0_JPRB,MIN(1.0_JPRB,ZTRSOC(JL,JK))) |
---|
| 709 | ENDDO |
---|
| 710 | ENDDO |
---|
[2004] | 711 | DO JL=IBEG,IEND |
---|
[4875] | 712 | PTRSOD(JL)=MAX(0.0_JPRB,MIN(1.0_JPRB,ZFRSOD(JL)/(ZRII0*ZRMU0(JL)))) |
---|
[2004] | 713 | ENDDO |
---|
| 714 | |
---|
[4875] | 715 | !* 7.3 RECONSTRUCT FLUXES FOR DIAGNOSTICS |
---|
[2004] | 716 | |
---|
| 717 | DO JL=IBEG,IEND |
---|
[4875] | 718 | IF (PMU0(JL) < 1.E-10_JPRB) ZRMU0(JL)=0.0_JPRB |
---|
[2004] | 719 | ENDDO |
---|
[4875] | 720 | DO JK=1,KLEV+1 |
---|
| 721 | DO JL=IBEG,IEND |
---|
| 722 | PLWFT(JL,JK)=PEMTD(JL,JK) |
---|
| 723 | PSWFT(JL,JK)=ZRMU0(JL)*ZRII0*PTRSO(JL,JK) |
---|
| 724 | ENDDO |
---|
[2004] | 725 | ENDDO |
---|
[4875] | 726 | DO JK=1,2 |
---|
| 727 | DO JL=IBEG,IEND |
---|
| 728 | PSWFC(JL,JK)=ZRMU0(JL)*ZRII0*PCTRSO(JL,JK) |
---|
| 729 | PLWFC(JL,JK)=PCEMTR(JL,JK) |
---|
| 730 | ENDDO |
---|
| 731 | ENDDO |
---|
[2004] | 732 | |
---|
[4875] | 733 | !* 8.0 DIAGNOSTICS |
---|
| 734 | !---Now we copy back the correct fields to proceed to the next timestep |
---|
[2004] | 735 | |
---|
[4875] | 736 | IF ( AEROSOLFEEDBACK_ACTIVE .AND. (flag_aerosol .GT. 0 .OR. flag_aerosol_strat) ) THEN |
---|
[2004] | 737 | |
---|
[4875] | 738 | IF ( ok_ade .and. ok_aie ) THEN |
---|
| 739 | PFSUP(:,:) = ZFSUP_AERO(:,:,4) |
---|
| 740 | PFSDN(:,:) = ZFSDN_AERO(:,:,4) |
---|
| 741 | PFSCUP(:,:) = ZFSUP0_AERO(:,:,4) |
---|
| 742 | PFSCDN(:,:) = ZFSDN0_AERO(:,:,4) |
---|
[2146] | 743 | |
---|
[4875] | 744 | PFLUX(:,1,:) = LWUP_AERO(:,:,4) |
---|
| 745 | PFLUX(:,2,:) = LWDN_AERO(:,:,4) |
---|
| 746 | PFLUC(:,1,:) = LWUP0_AERO(:,:,4) |
---|
| 747 | PFLUC(:,2,:) = LWDN0_AERO(:,:,4) |
---|
| 748 | ENDIF |
---|
[2004] | 749 | |
---|
[4875] | 750 | IF ( ok_ade .and. (.not. ok_aie) ) THEN |
---|
| 751 | PFSUP(:,:) = ZFSUP_AERO(:,:,3) |
---|
| 752 | PFSDN(:,:) = ZFSDN_AERO(:,:,3) |
---|
| 753 | PFSCUP(:,:) = ZFSUP0_AERO(:,:,3) |
---|
| 754 | PFSCDN(:,:) = ZFSDN0_AERO(:,:,3) |
---|
[2146] | 755 | |
---|
[4875] | 756 | PFLUX(:,1,:) = LWUP_AERO(:,:,3) |
---|
| 757 | PFLUX(:,2,:) = LWDN_AERO(:,:,3) |
---|
| 758 | PFLUC(:,1,:) = LWUP0_AERO(:,:,3) |
---|
| 759 | PFLUC(:,2,:) = LWDN0_AERO(:,:,3) |
---|
| 760 | ENDIF |
---|
[2004] | 761 | |
---|
[4875] | 762 | IF ( (.not. ok_ade) .and. ok_aie ) THEN |
---|
| 763 | PFSUP(:,:) = ZFSUP_AERO(:,:,2) |
---|
| 764 | PFSDN(:,:) = ZFSDN_AERO(:,:,2) |
---|
| 765 | PFSCUP(:,:) = ZFSUP0_AERO(:,:,2) |
---|
| 766 | PFSCDN(:,:) = ZFSDN0_AERO(:,:,2) |
---|
[2146] | 767 | |
---|
[4875] | 768 | PFLUX(:,1,:) = LWUP_AERO(:,:,2) |
---|
| 769 | PFLUX(:,2,:) = LWDN_AERO(:,:,2) |
---|
| 770 | PFLUC(:,1,:) = LWUP0_AERO(:,:,2) |
---|
| 771 | PFLUC(:,2,:) = LWDN0_AERO(:,:,2) |
---|
| 772 | ENDiF |
---|
[2004] | 773 | |
---|
[4875] | 774 | IF ((.not. ok_ade) .and. (.not. ok_aie)) THEN |
---|
| 775 | PFSUP(:,:) = ZFSUP_AERO(:,:,1) |
---|
| 776 | PFSDN(:,:) = ZFSDN_AERO(:,:,1) |
---|
| 777 | PFSCUP(:,:) = ZFSUP0_AERO(:,:,1) |
---|
| 778 | PFSCDN(:,:) = ZFSDN0_AERO(:,:,1) |
---|
[2146] | 779 | |
---|
[4875] | 780 | PFLUX(:,1,:) = LWUP_AERO(:,:,1) |
---|
| 781 | PFLUX(:,2,:) = LWDN_AERO(:,:,1) |
---|
| 782 | PFLUC(:,1,:) = LWUP0_AERO(:,:,1) |
---|
| 783 | PFLUC(:,2,:) = LWDN0_AERO(:,:,1) |
---|
| 784 | ENDIF |
---|
| 785 | |
---|
| 786 | ! The following allows to compute the forcing diagostics without |
---|
| 787 | ! letting the aerosol forcing act on the meteorology |
---|
| 788 | ! SEE logic above |
---|
| 789 | |
---|
| 790 | ELSE !--not AEROSOLFEEDBACK_ACTIVE |
---|
| 791 | |
---|
| 792 | PFSUP(:,:) = ZFSUP_AERO(:,:,5) |
---|
| 793 | PFSDN(:,:) = ZFSDN_AERO(:,:,5) |
---|
| 794 | PFSCUP(:,:) = ZFSUP0_AERO(:,:,5) |
---|
| 795 | PFSCDN(:,:) = ZFSDN0_AERO(:,:,5) |
---|
| 796 | |
---|
| 797 | PFLUX(:,1,:) = LWUP_AERO(:,:,5) |
---|
| 798 | PFLUX(:,2,:) = LWDN_AERO(:,:,5) |
---|
| 799 | PFLUC(:,1,:) = LWUP0_AERO(:,:,5) |
---|
| 800 | PFLUC(:,2,:) = LWDN0_AERO(:,:,5) |
---|
| 801 | |
---|
[2004] | 802 | ENDIF |
---|
| 803 | |
---|
[4875] | 804 | !--VolMIP Strat/Surf |
---|
| 805 | !--only ok_ade + ok_aie case treated |
---|
| 806 | IF (ok_ade.AND.ok_aie.AND.ok_volcan) THEN |
---|
| 807 | !--in this case the fluxes used for the heating rates come from case 4 but SW surface radiation is kept from case 2 |
---|
| 808 | IF (flag_volc_surfstrat.EQ.2) THEN ! STRAT HEATING |
---|
| 809 | volmip_solsw(:)= ZFSDN_AERO(:,1,2)-ZFSUP_AERO(:,1,2) |
---|
| 810 | ELSEIF (flag_volc_surfstrat.EQ.1) THEN ! SURF COOLING |
---|
| 811 | !--in this case the fluxes used for the heating rates come from case 2 but SW surface radiation is kept from case 4 |
---|
| 812 | PFSUP(:,:) = ZFSUP_AERO(:,:,2) |
---|
| 813 | PFSDN(:,:) = ZFSDN_AERO(:,:,2) |
---|
| 814 | PFSCUP(:,:) = ZFSUP0_AERO(:,:,2) |
---|
| 815 | PFSCDN(:,:) = ZFSDN0_AERO(:,:,2) |
---|
| 816 | PFLUX(:,1,:) = LWUP_AERO(:,:,2) |
---|
| 817 | PFLUX(:,2,:) = LWDN_AERO(:,:,2) |
---|
| 818 | PFLUC(:,1,:) = LWDN0_AERO(:,:,2) |
---|
| 819 | PFLUC(:,2,:) = LWDN0_AERO(:,:,2) |
---|
| 820 | volmip_solsw(:)= ZFSDN_AERO(:,1,4)-ZFSUP_AERO(:,1,4) |
---|
| 821 | ENDIF |
---|
| 822 | ENDIF |
---|
| 823 | !--End VolMIP Strat/Surf |
---|
[2004] | 824 | |
---|
[4875] | 825 | IF (swaerofree_diag) THEN |
---|
| 826 | ! copy shortwave clear-sky clean (no aerosol) case |
---|
| 827 | PFSCCUP(:,:) = ZFSUP0_AERO(:,:,5) |
---|
| 828 | PFSCCDN(:,:) = ZFSDN0_AERO(:,:,5) |
---|
| 829 | ! copy longwave clear-sky clean (no aerosol) case |
---|
| 830 | PFLCCUP(:,:) = LWUP0_AERO(:,:,5) |
---|
| 831 | PFLCCDN(:,:) = LWDN0_AERO(:,:,5) |
---|
| 832 | ENDIF |
---|
[2004] | 833 | |
---|
[4875] | 834 | !OB- HERE CHECK WITH MP IF BOTTOM AND TOP INDICES ARE OK !!!!!!!!!!!!!!!!!! |
---|
| 835 | ! net anthropogenic forcing direct and 1st indirect effect diagnostics |
---|
| 836 | ! requires a natural aerosol field read and used |
---|
| 837 | ! Difference of net fluxes from double call to radiation |
---|
| 838 | ! Will need to be extended to LW radiation -> done by CK (2014-05-23) |
---|
[2004] | 839 | |
---|
[4875] | 840 | IF (flag_aerosol .GT. 0 .OR. flag_aerosol_strat) THEN |
---|
[2146] | 841 | |
---|
[4875] | 842 | IF (ok_ade.AND.ok_aie) THEN |
---|
[2004] | 843 | |
---|
[4875] | 844 | ! direct anthropogenic forcing |
---|
| 845 | PSOLSWADAERO(:) = (ZFSDN_AERO(:,1,4) -ZFSUP_AERO(:,1,4)) -(ZFSDN_AERO(:,1,2) -ZFSUP_AERO(:,1,2)) |
---|
| 846 | PTOPSWADAERO(:) = (ZFSDN_AERO(:,KLEV+1,4) -ZFSUP_AERO(:,KLEV+1,4)) -(ZFSDN_AERO(:,KLEV+1,2) -ZFSUP_AERO(:,KLEV+1,2)) |
---|
| 847 | PSOLSWAD0AERO(:) = (ZFSDN0_AERO(:,1,4) -ZFSUP0_AERO(:,1,4)) -(ZFSDN0_AERO(:,1,2) -ZFSUP0_AERO(:,1,2)) |
---|
| 848 | PTOPSWAD0AERO(:) = (ZFSDN0_AERO(:,KLEV+1,4)-ZFSUP0_AERO(:,KLEV+1,4))-(ZFSDN0_AERO(:,KLEV+1,2)-ZFSUP0_AERO(:,KLEV+1,2)) |
---|
| 849 | IF(ok_volcan) THEN |
---|
| 850 | PSWADAERO(:,:) = (ZFSDN_AERO(:,:,4) -ZFSUP_AERO(:,:,4)) -(ZFSDN_AERO(:,:,2) -ZFSUP_AERO(:,:,2)) !--NL |
---|
| 851 | ENDIF |
---|
[3989] | 852 | |
---|
[4875] | 853 | ! indirect anthropogenic forcing |
---|
| 854 | PSOLSWAIAERO(:) = (ZFSDN_AERO(:,1,4) -ZFSUP_AERO(:,1,4)) -(ZFSDN_AERO(:,1,3) -ZFSUP_AERO(:,1,3)) |
---|
| 855 | PTOPSWAIAERO(:) = (ZFSDN_AERO(:,KLEV+1,4)-ZFSUP_AERO(:,KLEV+1,4))-(ZFSDN_AERO(:,KLEV+1,3)-ZFSUP_AERO(:,KLEV+1,3)) |
---|
[3083] | 856 | |
---|
[4875] | 857 | ! Cloud radiative forcing with natural aerosol for direct effect |
---|
| 858 | PSOLSWCFAERO(:,1) = (ZFSDN_AERO(:,1,2) -ZFSUP_AERO(:,1,2)) -(ZFSDN0_AERO(:,1,2) -ZFSUP0_AERO(:,1,2)) |
---|
| 859 | PTOPSWCFAERO(:,1) = (ZFSDN_AERO(:,KLEV+1,2)-ZFSUP_AERO(:,KLEV+1,2))-(ZFSDN0_AERO(:,KLEV+1,2)-ZFSUP0_AERO(:,KLEV+1,2)) |
---|
| 860 | ! Cloud radiative forcing with anthropogenic aerosol for direct effect |
---|
| 861 | PSOLSWCFAERO(:,2) = (ZFSDN_AERO(:,1,4) -ZFSUP_AERO(:,1,4)) -(ZFSDN0_AERO(:,1,4) -ZFSUP0_AERO(:,1,4)) |
---|
| 862 | PTOPSWCFAERO(:,2) = (ZFSDN_AERO(:,KLEV+1,4)-ZFSUP_AERO(:,KLEV+1,4))-(ZFSDN0_AERO(:,KLEV+1,4)-ZFSUP0_AERO(:,KLEV+1,4)) |
---|
| 863 | ! Cloud radiative forcing with no direct effect at all |
---|
| 864 | PSOLSWCFAERO(:,3) = 0.0 |
---|
| 865 | PTOPSWCFAERO(:,3) = 0.0 |
---|
[2004] | 866 | |
---|
[4875] | 867 | ! LW direct anthropogenic forcing |
---|
| 868 | PSOLLWADAERO(:) = (-LWDN_AERO(:,1,4) -LWUP_AERO(:,1,4)) -(-LWDN_AERO(:,1,2) -LWUP_AERO(:,1,2)) |
---|
| 869 | PTOPLWADAERO(:) = (-LWDN_AERO(:,KLEV+1,4) -LWUP_AERO(:,KLEV+1,4)) -(-LWDN_AERO(:,KLEV+1,2) -LWUP_AERO(:,KLEV+1,2)) |
---|
| 870 | PSOLLWAD0AERO(:) = (-LWDN0_AERO(:,1,4) -LWUP0_AERO(:,1,4)) -(-LWDN0_AERO(:,1,2) -LWUP0_AERO(:,1,2)) |
---|
| 871 | PTOPLWAD0AERO(:) = (-LWDN0_AERO(:,KLEV+1,4)-LWUP0_AERO(:,KLEV+1,4))-(-LWDN0_AERO(:,KLEV+1,2)-LWUP0_AERO(:,KLEV+1,2)) |
---|
| 872 | IF(ok_volcan) THEN |
---|
| 873 | PLWADAERO(:,:) = (-LWDN_AERO(:,:,4) -LWUP_AERO(:,:,4)) -(-LWDN_AERO(:,:,2) -LWUP_AERO(:,:,2)) !--NL |
---|
| 874 | ENDIF |
---|
[2004] | 875 | |
---|
[4875] | 876 | ! LW indirect anthropogenic forcing |
---|
| 877 | PSOLLWAIAERO(:) = (-LWDN_AERO(:,1,4) -LWUP_AERO(:,1,4)) -(-LWDN_AERO(:,1,3) -LWUP_AERO(:,1,3)) |
---|
| 878 | PTOPLWAIAERO(:) = (-LWDN_AERO(:,KLEV+1,4)-LWUP_AERO(:,KLEV+1,4))-(-LWDN_AERO(:,KLEV+1,3)-LWUP_AERO(:,KLEV+1,3)) |
---|
[2004] | 879 | |
---|
[3479] | 880 | ENDIF |
---|
[2004] | 881 | |
---|
[4875] | 882 | IF (ok_ade.AND..NOT.ok_aie) THEN |
---|
[2004] | 883 | |
---|
[4875] | 884 | ! direct anthropogenic forcing |
---|
| 885 | PSOLSWADAERO(:) = (ZFSDN_AERO(:,1,3) -ZFSUP_AERO(:,1,3)) -(ZFSDN_AERO(:,1,1) -ZFSUP_AERO(:,1,1)) |
---|
| 886 | PTOPSWADAERO(:) = (ZFSDN_AERO(:,KLEV+1,3) -ZFSUP_AERO(:,KLEV+1,3)) -(ZFSDN_AERO(:,KLEV+1,1) -ZFSUP_AERO(:,KLEV+1,1)) |
---|
| 887 | PSOLSWAD0AERO(:) = (ZFSDN0_AERO(:,1,3) -ZFSUP0_AERO(:,1,3)) -(ZFSDN0_AERO(:,1,1) -ZFSUP0_AERO(:,1,1)) |
---|
| 888 | PTOPSWAD0AERO(:) = (ZFSDN0_AERO(:,KLEV+1,3)-ZFSUP0_AERO(:,KLEV+1,3))-(ZFSDN0_AERO(:,KLEV+1,1)-ZFSUP0_AERO(:,KLEV+1,1)) |
---|
| 889 | IF(ok_volcan) THEN |
---|
| 890 | PSWADAERO(:,:) = (ZFSDN_AERO(:,:,3) -ZFSUP_AERO(:,:,3)) -(ZFSDN_AERO(:,:,1) -ZFSUP_AERO(:,:,1)) !--NL |
---|
| 891 | ENDIF |
---|
[2004] | 892 | |
---|
[4875] | 893 | ! indirect anthropogenic forcing |
---|
| 894 | PSOLSWAIAERO(:) = 0.0 |
---|
| 895 | PTOPSWAIAERO(:) = 0.0 |
---|
[2146] | 896 | |
---|
[4875] | 897 | ! Cloud radiative forcing with natural aerosol for direct effect |
---|
| 898 | PSOLSWCFAERO(:,1) = (ZFSDN_AERO(:,1,1) -ZFSUP_AERO(:,1,1)) -(ZFSDN0_AERO(:,1,1) -ZFSUP0_AERO(:,1,1)) |
---|
| 899 | PTOPSWCFAERO(:,1) = (ZFSDN_AERO(:,KLEV+1,1)-ZFSUP_AERO(:,KLEV+1,1))-(ZFSDN0_AERO(:,KLEV+1,1)-ZFSUP0_AERO(:,KLEV+1,1)) |
---|
| 900 | ! Cloud radiative forcing with anthropogenic aerosol for direct effect |
---|
| 901 | PSOLSWCFAERO(:,2) = (ZFSDN_AERO(:,1,3) -ZFSUP_AERO(:,1,3)) -(ZFSDN0_AERO(:,1,3) -ZFSUP0_AERO(:,1,3)) |
---|
| 902 | PTOPSWCFAERO(:,2) = (ZFSDN_AERO(:,KLEV+1,3)-ZFSUP_AERO(:,KLEV+1,3))-(ZFSDN0_AERO(:,KLEV+1,3)-ZFSUP0_AERO(:,KLEV+1,3)) |
---|
| 903 | ! Cloud radiative forcing with no direct effect at all |
---|
| 904 | PSOLSWCFAERO(:,3) = 0.0 |
---|
| 905 | PTOPSWCFAERO(:,3) = 0.0 |
---|
[2146] | 906 | |
---|
[4875] | 907 | ! LW direct anthropogenic forcing |
---|
| 908 | PSOLLWADAERO(:) = (-LWDN_AERO(:,1,3) -LWUP_AERO(:,1,3)) -(-LWDN_AERO(:,1,1) -LWUP_AERO(:,1,1)) |
---|
| 909 | PTOPLWADAERO(:) = (-LWDN_AERO(:,KLEV+1,3) -LWUP_AERO(:,KLEV+1,3)) -(-LWDN_AERO(:,KLEV+1,1) -LWUP_AERO(:,KLEV+1,1)) |
---|
| 910 | PSOLLWAD0AERO(:) = (-LWDN0_AERO(:,1,3) -LWUP0_AERO(:,1,3)) -(-LWDN0_AERO(:,1,1) -LWUP0_AERO(:,1,1)) |
---|
| 911 | PTOPLWAD0AERO(:) = (-LWDN0_AERO(:,KLEV+1,3)-LWUP0_AERO(:,KLEV+1,3))-(-LWDN0_AERO(:,KLEV+1,1)-LWUP0_AERO(:,KLEV+1,1)) |
---|
| 912 | IF(ok_volcan) THEN |
---|
| 913 | PLWADAERO(:,:) = (-LWDN_AERO(:,:,3) -LWUP_AERO(:,:,3)) -(-LWDN_AERO(:,:,1) -LWUP_AERO(:,:,1)) !--NL |
---|
| 914 | ENDIF |
---|
[2004] | 915 | |
---|
[4875] | 916 | ! LW indirect anthropogenic forcing |
---|
| 917 | PSOLLWAIAERO(:) = 0.0 |
---|
| 918 | PTOPLWAIAERO(:) = 0.0 |
---|
[2004] | 919 | |
---|
[3479] | 920 | ENDIF |
---|
[2004] | 921 | |
---|
[4875] | 922 | IF (.NOT.ok_ade.AND.ok_aie) THEN |
---|
[2004] | 923 | |
---|
[4875] | 924 | ! direct anthropogenic forcing |
---|
| 925 | PSOLSWADAERO(:) = 0.0 |
---|
| 926 | PTOPSWADAERO(:) = 0.0 |
---|
| 927 | PSOLSWAD0AERO(:) = 0.0 |
---|
| 928 | PTOPSWAD0AERO(:) = 0.0 |
---|
| 929 | IF(ok_volcan) THEN |
---|
| 930 | PSWADAERO(:,:) = 0.0 !--NL |
---|
| 931 | ENDIF |
---|
[2004] | 932 | |
---|
[4875] | 933 | ! indirect anthropogenic forcing |
---|
| 934 | PSOLSWAIAERO(:) = (ZFSDN_AERO(:,1,2) -ZFSUP_AERO(:,1,2)) -(ZFSDN_AERO(:,1,1) -ZFSUP_AERO(:,1,1)) |
---|
| 935 | PTOPSWAIAERO(:) = (ZFSDN_AERO(:,KLEV+1,2)-ZFSUP_AERO(:,KLEV+1,2))-(ZFSDN_AERO(:,KLEV+1,1)-ZFSUP_AERO(:,KLEV+1,1)) |
---|
[2146] | 936 | |
---|
[4875] | 937 | ! Cloud radiative forcing with natural aerosol for direct effect |
---|
| 938 | PSOLSWCFAERO(:,1) = (ZFSDN_AERO(:,1,2) -ZFSUP_AERO(:,1,2)) -(ZFSDN0_AERO(:,1,2) -ZFSUP0_AERO(:,1,2)) |
---|
| 939 | PTOPSWCFAERO(:,1) = (ZFSDN_AERO(:,KLEV+1,2)-ZFSUP_AERO(:,KLEV+1,2))-(ZFSDN0_AERO(:,KLEV+1,2)-ZFSUP0_AERO(:,KLEV+1,2)) |
---|
| 940 | ! Cloud radiative forcing with anthropogenic aerosol for direct effect |
---|
| 941 | PSOLSWCFAERO(:,2) = 0.0 |
---|
| 942 | PTOPSWCFAERO(:,2) = 0.0 |
---|
| 943 | ! Cloud radiative forcing with no direct effect at all |
---|
| 944 | PSOLSWCFAERO(:,3) = 0.0 |
---|
| 945 | PTOPSWCFAERO(:,3) = 0.0 |
---|
[2004] | 946 | |
---|
[4875] | 947 | ! LW direct anthropogenic forcing |
---|
| 948 | PSOLLWADAERO(:) = 0.0 |
---|
| 949 | PTOPLWADAERO(:) = 0.0 |
---|
| 950 | PSOLLWAD0AERO(:) = 0.0 |
---|
| 951 | PTOPLWAD0AERO(:) = 0.0 |
---|
| 952 | IF(ok_volcan) THEN |
---|
| 953 | PLWADAERO(:,:) = 0.0 !--NL |
---|
| 954 | ENDIF |
---|
[2004] | 955 | |
---|
[4875] | 956 | ! LW indirect anthropogenic forcing |
---|
| 957 | PSOLLWAIAERO(:) = (-LWDN_AERO(:,1,2) -LWUP_AERO(:,1,2)) -(-LWDN_AERO(:,1,1) -LWUP_AERO(:,1,1)) |
---|
| 958 | PTOPLWAIAERO(:) = (-LWDN_AERO(:,KLEV+1,2)-LWUP_AERO(:,KLEV+1,2))-(-LWDN_AERO(:,KLEV+1,1)-LWUP_AERO(:,KLEV+1,1)) |
---|
| 959 | |
---|
[3479] | 960 | ENDIF |
---|
[2004] | 961 | |
---|
[4875] | 962 | IF (.NOT.ok_ade.AND..NOT.ok_aie) THEN |
---|
[2004] | 963 | |
---|
[4875] | 964 | ! direct anthropogenic forcing |
---|
| 965 | PSOLSWADAERO(:) = 0.0 |
---|
| 966 | PTOPSWADAERO(:) = 0.0 |
---|
| 967 | PSOLSWAD0AERO(:) = 0.0 |
---|
| 968 | PTOPSWAD0AERO(:) = 0.0 |
---|
| 969 | IF(ok_volcan) THEN |
---|
| 970 | PSWADAERO(:,:) = 0.0 !--NL |
---|
| 971 | ENDIF |
---|
[2146] | 972 | |
---|
[4875] | 973 | ! indirect anthropogenic forcing |
---|
| 974 | PSOLSWAIAERO(:) = 0.0 |
---|
| 975 | PTOPSWAIAERO(:) = 0.0 |
---|
[2004] | 976 | |
---|
[4875] | 977 | ! Cloud radiative forcing with natural aerosol for direct effect |
---|
| 978 | PSOLSWCFAERO(:,1) = (ZFSDN_AERO(:,1,1) -ZFSUP_AERO(:,1,1)) -(ZFSDN0_AERO(:,1,1) -ZFSUP0_AERO(:,1,1)) |
---|
| 979 | PTOPSWCFAERO(:,1) = (ZFSDN_AERO(:,KLEV+1,1)-ZFSUP_AERO(:,KLEV+1,1))-(ZFSDN0_AERO(:,KLEV+1,1)-ZFSUP0_AERO(:,KLEV+1,1)) |
---|
| 980 | ! Cloud radiative forcing with anthropogenic aerosol for direct effect |
---|
| 981 | PSOLSWCFAERO(:,2) = 0.0 |
---|
| 982 | PTOPSWCFAERO(:,2) = 0.0 |
---|
| 983 | ! Cloud radiative forcing with no direct effect at all |
---|
| 984 | PSOLSWCFAERO(:,3) = 0.0 |
---|
| 985 | PTOPSWCFAERO(:,3) = 0.0 |
---|
[2004] | 986 | |
---|
[4875] | 987 | ! LW direct anthropogenic forcing |
---|
| 988 | PSOLLWADAERO(:) = 0.0 |
---|
| 989 | PTOPLWADAERO(:) = 0.0 |
---|
| 990 | PSOLLWAD0AERO(:) = 0.0 |
---|
| 991 | PTOPLWAD0AERO(:) = 0.0 |
---|
| 992 | IF(ok_volcan) THEN |
---|
| 993 | PLWADAERO(:,:) = 0.0 !--NL |
---|
| 994 | ENDIF |
---|
[2004] | 995 | |
---|
[4875] | 996 | ! LW indirect anthropogenic forcing |
---|
| 997 | PSOLLWAIAERO(:) = 0.0 |
---|
| 998 | PTOPLWAIAERO(:) = 0.0 |
---|
[2004] | 999 | |
---|
[3479] | 1000 | ENDIF |
---|
[2146] | 1001 | |
---|
[4875] | 1002 | ENDIF |
---|
[2004] | 1003 | |
---|
[4875] | 1004 | !IF (swaero_diag .OR. .NOT. AEROSOLFEEDBACK_ACTIVE) THEN |
---|
| 1005 | IF (.NOT. AEROSOLFEEDBACK_ACTIVE) THEN |
---|
| 1006 | ! Cloudforcing without aerosol at all |
---|
[2004] | 1007 | PSOLSWCFAERO(:,3) = (ZFSDN_AERO(:,1,5) -ZFSUP_AERO(:,1,5)) -(ZFSDN0_AERO(:,1,5) -ZFSUP0_AERO(:,1,5)) |
---|
| 1008 | PTOPSWCFAERO(:,3) = (ZFSDN_AERO(:,KLEV+1,5)-ZFSUP_AERO(:,KLEV+1,5))-(ZFSDN0_AERO(:,KLEV+1,5)-ZFSUP0_AERO(:,KLEV+1,5)) |
---|
[2146] | 1009 | |
---|
[4875] | 1010 | ENDIF |
---|
[2004] | 1011 | |
---|
[4875] | 1012 | IF (LHOOK) CALL DR_HOOK('RECMWF_AERO',1,ZHOOK_HANDLE) |
---|
[2004] | 1013 | END SUBROUTINE RECMWF_AERO |
---|