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

Timestamp:

Dec 23, 2024, 6:12:29 PM (40 hours ago)

Author:

aborella

Message:

First implementation of the contrails parameterisation
Lacks the emission of H2O + the impact on radiative transfer

File:

• LMDZ6/branches/contrails/libf/phylmd/phys_local_var_mod.F90 (modified) (3 diffs)

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

@@ -665 +665 @@
 
 !-- LSCP - aviation and contrails variables
-      REAL, SAVE, ALLOCATABLE :: Tcontr(:,:), qcontr(:,:), qcontr2(:,:)
-      !$OMP THREADPRIVATE(Tcontr, qcontr, qcontr2)
-      REAL, SAVE, ALLOCATABLE :: fcontrN(:,:), fcontrP(:,:)
-      !$OMP THREADPRIVATE(fcontrN, fcontrP)
+      REAL, SAVE, ALLOCATABLE :: rcont_seri(:,:), d_rcont_dyn(:,:)
+      !$OMP THREADPRIVATE(rcont_seri, d_rcont_dyn)
+      REAL, SAVE, ALLOCATABLE :: flight_dist(:,:), flight_h2o(:,:)
+      !$OMP THREADPRIVATE(flight_dist, flight_h2o)
+      REAL, SAVE, ALLOCATABLE :: Tcritcont(:,:), qcritcont(:,:)
+      !$OMP THREADPRIVATE(Tcritcont, qcritcont)
+      REAL, SAVE, ALLOCATABLE :: potcontfraP(:,:), potcontfraNP(:,:)
+      !$OMP THREADPRIVATE(potcontfraP, potcontfraNP)
       REAL, SAVE, ALLOCATABLE :: dcf_avi(:,:), dqi_avi(:,:), dqvc_avi(:,:)
       !$OMP THREADPRIVATE(dcf_avi, dqi_avi, dqvc_avi)
-      REAL, SAVE, ALLOCATABLE :: flight_dist(:,:), flight_h2o(:,:)
-      !$OMP THREADPRIVATE(flight_dist, flight_h2o)
 
 !-- LSCP - mixed phase clouds variables
@@ -1213 +1215 @@
 
 !-- LSCP - aviation and contrails variables
-      ALLOCATE(Tcontr(klon,klev), qcontr(klon,klev), qcontr2(klon,klev))
-      ALLOCATE(fcontrN(klon,klev), fcontrP(klon,klev))
+      ALLOCATE(rcont_seri(klon,klev), d_rcont_dyn(klon,klev))
+      ALLOCATE(flight_dist(klon,klev), flight_h2o(klon,klev))
+      ALLOCATE(Tcritcont(klon,klev), qcritcont(klon,klev))
+      ALLOCATE(potcontfraP(klon,klev), potcontfraNP(klon,klev))
       ALLOCATE(dcf_avi(klon,klev), dqi_avi(klon,klev), dqvc_avi(klon,klev))
-      ALLOCATE(flight_dist(klon,klev), flight_h2o(klon,klev))
 
 !-- LSCP - POPRECIP variables
@@ -1618 +1621 @@
 
 !-- LSCP - aviation and contrails variables
-      DEALLOCATE(Tcontr, qcontr, qcontr2)
-      DEALLOCATE(fcontrN, fcontrP)
+      DEALLOCATE(rcont_seri, d_rcont_dyn, flight_dist, flight_h2o)
+      DEALLOCATE(Tcritcont, qcritcont, potcontfraP, potcontfraNP)
       DEALLOCATE(dcf_avi, dqi_avi, dqvc_avi)
-      DEALLOCATE(flight_dist, flight_h2o)
 
 !-- LSCP - POPRECIP variables