source: LMDZ6/branches/LMDZ_ECRad/libf/phylmd/rrtm/recmwf_aero.F90 @ 5134

Last change on this file since 5134 was 3989, checked in by oboucher, 3 years ago

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