Changeset 5589 for LMDZ6/branches/contrails/libf/phylmd/rrtm

Timestamp:

Mar 26, 2025, 6:05:40 PM (5 months ago)

Author:

aborella

Message:

Multiple changes:

• added new radiative diagnostics for contrails
• added ok_rad_contrail option to allow for a double call of RRTM (w/ and w/o contrails)
• transformed resuspension of snow into ice sedimentation (poprecip)
• some modifications in poprecip in line with the ones from EV
• cleaned sublimation of ice clouds in lmdz_lscp_condensation, option ok_ice_supersat
• aviation emissions can now be read with IOIPSL (in lon/lat mode)

File:

• LMDZ6/branches/contrails/libf/phylmd/rrtm/recmwf_aero.F90 (modified) (4 diffs)

LMDZ6/branches/contrails/libf/phylmd/rrtm/recmwf_aero.F90

@@ -40 +40 @@
      !..end
      & ok_ade, ok_aie, ok_volcan, flag_aerosol,flag_aerosol_strat,&
-     & flag_aer_feedback)
+     & flag_aer_feedback, &
+     !--AB contrails radiative effect
+     & ok_rad_contrail, PCLFR_NOCONT, PQIWP_NOCONT, PREF_ICE_NOCONT, &
+     & PTOPSWNOCONT, PSOLSWNOCONT, PTOPLWNOCONT, PSOLLWNOCONT)
   !--fin
 
@@ -266 +269 @@
   REAL(KIND=JPRB)   ,INTENT(OUT)   :: volmip_solsw(KPROMA) ! SW clear sky in the case of VOLMIP
   INTEGER, INTENT(IN)              :: flag_volc_surfstrat !--VOlMIP Modif
+  !--AB contrails radiative effect
+  LOGICAL           ,INTENT(IN)    :: ok_rad_contrail
+  REAL(KIND=JPRB)   ,INTENT(IN)    :: PCLFR_NOCONT(KPROMA,KLEV) 
+  REAL(KIND=JPRB)   ,INTENT(IN)    :: PQIWP_NOCONT(KPROMA,KLEV)
+  REAL(KIND=JPRB)   ,INTENT(IN)    :: PREF_ICE_NOCONT(KPROMA,KLEV)
+  REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTOPSWNOCONT(KPROMA), PSOLSWNOCONT(KPROMA) ! No contrails experiment forcing at TOA and surface (SW)
+  REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTOPLWNOCONT(KPROMA), PSOLLWNOCONT(KPROMA) ! No contrails experiment forcing at TOA and surface (LW)
 
   !     ==== COMPUTED IN RADITE ===
@@ -346 +356 @@
   REAL(KIND=JPRB) ::  LWUP0_AERO(KPROMA,KLEV+1,5)
   REAL(KIND=JPRB) ::  LWDN0_AERO(KPROMA,KLEV+1,5)
+  !--AB contrails radiative effect
+  REAL(KIND=JPRB) :: ZRCLC_NOCONT(KPROMA,KLEV), ZQIWP_NOCONT(KPROMA,KLEV)
+  REAL(KIND=JPRB) :: PREF_LIQ_NOCONT(KPROMA,KLEV)
+  REAL(KIND=JPRB) :: PPIZA_NOCONT(KPROMA,KLEV,NSW)
+  REAL(KIND=JPRB) :: PCGA_NOCONT(KPROMA,KLEV,NSW)
+  REAL(KIND=JPRB) :: PTAU_NOCONT(KPROMA,KLEV,NSW)
+  REAL(KIND=JPRB) :: PTAU_LW_NOCONT(KPROMA,KLEV,NLW)
 
 #include "radlsw.intfb.h"
@@ -673 +690 @@
   ENDIF ! .not. AEROSOLFEEDBACK_ACTIVE
 
+
+  !--Double call to radiative routine for contrails
+  !--The calculation are done again WITHOUT contrails
+  IF (ok_rad_contrail) THEN
+
+     !--The same base case is used
+     !--NB. here we could use pointers...
+     IF ( flag_aerosol .EQ. 0 ) THEN
+        PREF_LIQ_NOCONT(:,:) = PREF_LIQ_PI(:,:)
+        PPIZA_NOCONT(:,:,:) = PPIZA_ZERO(:,:,:)
+        PCGA_NOCONT(:,:,:) = PCGA_ZERO(:,:,:)
+        PTAU_NOCONT(:,:,:) = PTAU_ZERO(:,:,:)
+        PTAU_LW_NOCONT(:,:,:) = PTAU_LW_ZERO(:,:,:)
+     ELSEIF ( .not. ok_ade .AND. .not. ok_aie ) THEN
+        PREF_LIQ_NOCONT(:,:) = PREF_LIQ_PI(:,:)
+        PPIZA_NOCONT(:,:,:) = PPIZA_NAT(:,:,:)
+        PCGA_NOCONT(:,:,:) = PCGA_NAT(:,:,:)
+        PTAU_NOCONT(:,:,:) = PTAU_NAT(:,:,:)
+        PTAU_LW_NOCONT(:,:,:) = PTAU_LW_NAT(:,:,:)
+     ELSEIF ( .not. ok_ade .AND. ok_aie ) THEN
+        PREF_LIQ_NOCONT(:,:) = PREF_LIQ(:,:)
+        PPIZA_NOCONT(:,:,:) = PPIZA_NAT(:,:,:)
+        PCGA_NOCONT(:,:,:) = PCGA_NAT(:,:,:)
+        PTAU_NOCONT(:,:,:) = PTAU_NAT(:,:,:)
+        PTAU_LW_NOCONT(:,:,:) = PTAU_LW_NAT(:,:,:)
+     ELSEIF ( ok_ade .AND. .not. ok_aie ) THEN
+        PREF_LIQ_NOCONT(:,:) = PREF_LIQ_PI(:,:)
+        PPIZA_NOCONT(:,:,:) = PPIZA_TOT(:,:,:)
+        PCGA_NOCONT(:,:,:) = PCGA_TOT(:,:,:)
+        PTAU_NOCONT(:,:,:) = PTAU_TOT(:,:,:)
+        PTAU_LW_NOCONT(:,:,:) = PTAU_LW_TOT(:,:,:)
+     ELSEIF ( ok_ade .AND. ok_aie ) THEN
+        PREF_LIQ_NOCONT(:,:) = PREF_LIQ(:,:)
+        PPIZA_NOCONT(:,:,:) = PPIZA_TOT(:,:,:)
+        PCGA_NOCONT(:,:,:) = PCGA_TOT(:,:,:)
+        PTAU_NOCONT(:,:,:) = PTAU_TOT(:,:,:)
+        PTAU_LW_NOCONT(:,:,:) = PTAU_LW_TOT(:,:,:)
+     ENDIF
+
+     DO JK=1,KLEV
+         DO JL=IBEG,IEND
+             ZRCLC_NOCONT(JL,JK)=MAX( 0.0_JPRB ,MIN( 1.0_JPRB ,PCLFR_NOCONT(JL,JK)))
+             IF (ZRCLC_NOCONT(JL,JK) > REPCLC) THEN
+                 ZQIWP_NOCONT(JL,JK)=PQIWP_NOCONT(JL,JK)
+             ELSE
+                 ZQIWP_NOCONT(JL,JK)=REPH2O*ZRCLC_NOCONT(JL,JK)
+             ENDIF
+         ENDDO
+     ENDDO
+
+     CALL RADLSW (&
+          & IBEG  , IEND   , KPROMA  , KLEV  , KMODE , NAER,&
+          & ZRII0 ,&
+          & ZRAER , PALBD  , PALBP   , PAPRS , ZRPR  ,&
+          & ZCCNL , ZCCNO  ,&
+          & PCCO2 , ZRCLC_NOCONT  , PDP     , PEMIS , ZEMIW ,PSLM    , ZRMU0 , ZPQO3,&
+          & ZQ    , ZQIWP_NOCONT  , ZQLWP   , ZQS   , ZQRAIN,ZQRAINT ,&
+          & PTH   , ZRTI   , PTS     , ZNBAS , ZNTOP ,&
+          & PREF_LIQ_NOCONT, PREF_ICE_NOCONT,&
+          & ZEMIT , ZFCT   , ZFLT    , ZFCS    , ZFLS  ,&
+          & ZFRSOD, ZSUDU  , ZUVDF   , ZPARF   , ZPARCF, ZTINCF, PSFSWDIR,&
+          & PSFSWDIF,PFSDNN, PFSDNV  ,&  
+          & LRDUST,PPIZA_NOCONT,PCGA_NOCONT,PTAU_NOCONT,PTAU_LW_NOCONT,PFLUX,PFLUC,&
+          & PFSDN , PFSUP  , PFSCDN  , PFSCUP )
+
+     ! save budgets LW and SW at TOA and surface
+     PTOPSWNOCONT(:) = PFSDN(:,KLEV+1) - PFSUP(:,KLEV+1)
+     PSOLSWNOCONT(:) = PFSDN(:,1) - PFSUP(:,1)
+     PTOPLWNOCONT(:) = PFLUX(:,1,KLEV+1) + PFLUX(:,2,KLEV+1)
+     PSOLLWNOCONT(:) = - PFLUX(:,1,1) - PFLUX(:,2,1)
+
+  ENDIF
+
   !*         4.2     TRANSFORM FLUXES TO MODEL HISTORICAL VARIABLES