Changeset 3082 for LMDZ6/trunk/libf/phylmd/rrtm

Timestamp:

Nov 21, 2017, 5:18:23 PM (7 years ago)

Author:

oboucher

Message:

Adding new diagnostics for clean (no aerosol) clear sky SW radiative fluxes

File:

• LMDZ6/trunk/libf/phylmd/rrtm/recmwf_aero.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/rrtm/recmwf_aero.F90

@@ -24 +24 @@
 !--end
  & PFLUX,PFLUC,&
- & PFSDN ,PFSUP , PFSCDN , PFSCUP,&
+ & PFSDN ,PFSUP , PFSCDN , PFSCUP, PFSCCDN, PFSCCUP, &
 !--OB diagnostics
  & PTOPSWADAERO,PSOLSWADAERO,&
@@ -110 +110 @@
 ! PFSCDN(KPROMA,KLEV+1)         ; SW clear sky flux down
 ! PFSCUP(KPROMA,KLEV+1)         ; SW clear sky flux up
+! PFSCCDN(KPROMA,KLEV+1)         ; SW clear sky clean (no aerosol) flux down
+! PFSCCUP(KPROMA,KLEV+1)         ; SW clear sky clean (no aerosol) flux up
 
 
@@ -152 +154 @@
 USE YOERDI   , ONLY : RRAE   ,REPCLC    ,REPH2O
 USE YOMARPHY , ONLY : LRDUST 
-USE phys_output_mod, ONLY : swaero_diag
+USE phys_output_mod, ONLY : swaerofree_diag, swaero_diag
 
 !-----------------------------------------------------------------------
@@ -241 +243 @@
 REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFSCDN(KPROMA,KLEV+1)  ! SW clear sky flux down
 REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFSCUP(KPROMA,KLEV+1)  ! SW clear sky flux up
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFSCCDN(KPROMA,KLEV+1)  ! SW clear sky clean (no aerosol) flux down
+REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFSCCUP(KPROMA,KLEV+1)  ! SW clear sky clean (no aerosol) flux up
 
 !     ==== COMPUTED IN RADITE ===
@@ -607 +611 @@
 
 ! case with no aerosols at all is also computed IF ACTIVEFEEDBACK_ACTIVE is false 
-!IF (swaero_diag .OR. .not. AEROSOLFEEDBACK_ACTIVE .OR. flag_aerosol .EQ. 0 ) THEN    
-IF (.not. AEROSOLFEEDBACK_ACTIVE .OR. flag_aerosol .EQ. 0 ) THEN    
+IF (.not. AEROSOLFEEDBACK_ACTIVE .OR. flag_aerosol .EQ. 0 .OR. swaerofree_diag) THEN    
 
 ! ZERO aerosol effect
@@ -757 +760 @@
   ENDIF
 
+! copy clear-sky clean (no aerosol) case
+  PFSCCUP(:,:) =   ZFSUP0_AERO(:,:,5)
+  PFSCCDN(:,:) =   ZFSDN0_AERO(:,:,5)
+
 ! The following allows to compute the forcing diagostics without
 ! letting the aerosol forcing act on the meteorology