Changeset 5589 for LMDZ6/branches/contrails/libf/phylmd/phys_local_var_mod.F90

Timestamp:

Mar 26, 2025, 6:05:40 PM (3 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/phys_local_var_mod.F90 (modified) (6 diffs)

LMDZ6/branches/contrails/libf/phylmd/phys_local_var_mod.F90 ¶

@@ -683 +683 @@
       REAL, SAVE, ALLOCATABLE :: dcf_avi(:,:), dqi_avi(:,:), dqvc_avi(:,:)
       !$OMP THREADPRIVATE(dcf_avi, dqi_avi, dqvc_avi)
+      REAL, SAVE, ALLOCATABLE :: cldfra_nocont(:,:), cldtau_nocont(:,:), cldemi_nocont(:,:)
+      !$OMP THREADPRIVATE(cldfra_nocont, cldtau_nocont, cldemi_nocont)
+      REAL, SAVE, ALLOCATABLE :: cldh_nocont(:), contcov(:)
+      !$OMP THREADPRIVATE(cldh_nocont, contcov)
+      REAL, SAVE, ALLOCATABLE :: fiwp_nocont(:), fiwc_nocont(:,:), ref_ice_nocont(:,:)
+      !$OMP THREADPRIVATE(fiwp_nocont, fiwc_nocont, ref_ice_nocont)
+      REAL, SAVE, ALLOCATABLE :: topsw_nocont(:), solsw_nocont(:)
+      !$OMP THREADPRIVATE(topsw_nocont, solsw_nocont)
+      REAL, SAVE, ALLOCATABLE :: toplw_nocont(:), sollw_nocont(:)
+      !$OMP THREADPRIVATE(toplw_nocont, sollw_nocont)
 
 !-- LSCP - mixed phase clouds variables
@@ -717 +727 @@
       REAL, SAVE, ALLOCATABLE :: dqsfreez(:,:)
       !$OMP THREADPRIVATE(dqsfreez)
-      REAL, SAVE, ALLOCATABLE :: dcfres(:,:)
-      !$OMP THREADPRIVATE(dcfres)
-      REAL, SAVE, ALLOCATABLE :: dqsres(:,:)
-      !$OMP THREADPRIVATE(dqsres)
-      REAL, SAVE, ALLOCATABLE :: dqvcres(:,:)
-      !$OMP THREADPRIVATE(dqvcres)
+      REAL, SAVE, ALLOCATABLE :: dqised(:,:)
+      !$OMP THREADPRIVATE(dqised)
+      REAL, SAVE, ALLOCATABLE :: dcfsed(:,:)
+      !$OMP THREADPRIVATE(dcfsed)
+      REAL, SAVE, ALLOCATABLE :: dqvcsed(:,:)
+      !$OMP THREADPRIVATE(dqvcsed)
 
 ! variables for stratospheric aerosol
@@ -1238 +1248 @@
       ALLOCATE(contfra(klon,klev), dcontfra_cir(klon,klev))
       ALLOCATE(dcf_avi(klon,klev), dqi_avi(klon,klev), dqvc_avi(klon,klev))
+      ALLOCATE(cldfra_nocont(klon,klev), cldtau_nocont(klon,klev), cldemi_nocont(klon,klev))
+      ALLOCATE(cldh_nocont(klon), contcov(klon))
+      ALLOCATE(fiwp_nocont(klon), fiwc_nocont(klon,klev), ref_ice_nocont(klon,klev))
+      ALLOCATE(topsw_nocont(klon), solsw_nocont(klon))
+      ALLOCATE(toplw_nocont(klon), sollw_nocont(klon))
 
 !-- LSCP - POPRECIP variables
@@ -1244 +1259 @@
       ALLOCATE(dqrauto(klon,klev), dqrcol(klon,klev), dqrmelt(klon,klev), dqrfreez(klon,klev))
       ALLOCATE(dqsauto(klon,klev), dqsagg(klon,klev), dqsrim(klon,klev), dqsmelt(klon,klev), dqsfreez(klon,klev))
-      ALLOCATE(dcfres(klon,klev), dqsres(klon,klev), dqvcres(klon,klev))
+      ALLOCATE(dqised(klon,klev), dcfsed(klon,klev), dqvcsed(klon,klev))
 
 IF (CPPKEY_STRATAER) THEN
@@ -1646 +1661 @@
       DEALLOCATE(contfra, dcontfra_cir)
       DEALLOCATE(dcf_avi, dqi_avi, dqvc_avi)
+      DEALLOCATE(cldfra_nocont, cldtau_nocont, cldemi_nocont)
+      DEALLOCATE(cldh_nocont, contcov, fiwp_nocont, fiwc_nocont, ref_ice_nocont)
+      DEALLOCATE(topsw_nocont, solsw_nocont, toplw_nocont, sollw_nocont)
 
 !-- LSCP - POPRECIP variables
@@ -1652 +1670 @@
       DEALLOCATE(dqrauto, dqrcol, dqrmelt, dqrfreez)
       DEALLOCATE(dqsauto, dqsagg, dqsrim, dqsmelt, dqsfreez)
-      DEALLOCATE(dcfres, dqsres, dqvcres)
+      DEALLOCATE(dqised, dcfsed, dqvcsed)
 
 IF (CPPKEY_STRATAER) THEN