Changeset 4059 for LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

Timestamp:

Jan 21, 2022, 3:50:54 PM (2 years ago)

Author:

oboucher

Message:

Audran Borella's parametrisation for ice supersaturation
activated with flag_ice_sursat (FALSE by default)

File:

• LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -16 +16 @@
       REAL, SAVE, ALLOCATABLE :: u_seri(:,:), v_seri(:,:)
       !$OMP THREADPRIVATE(u_seri, v_seri)
+      REAL, SAVE, ALLOCATABLE :: rneb_seri(:,:)
+      !$OMP THREADPRIVATE(rneb_seri)
+      REAL, SAVE, ALLOCATABLE :: d_rneb_dyn(:,:)
+      !$OMP THREADPRIVATE(d_rneb_dyn)
       REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:),l_mix(:,:,:),tke_dissip(:,:,:),wprime(:,:,:)
       !$OMP THREADPRIVATE(l_mixmin, l_mix, tke_dissip,wprime)
@@ -473 +477 @@
 !$OMP THREADPRIVATE(zn2mout)
 
+      REAL, SAVE, ALLOCATABLE :: qclr(:,:)
+      !$OMP THREADPRIVATE(qclr)
+      REAL, SAVE, ALLOCATABLE :: qcld(:,:)
+      !$OMP THREADPRIVATE(qcld)
+      REAL, SAVE, ALLOCATABLE :: qss(:,:)
+      !$OMP THREADPRIVATE(qss)
+      REAL, SAVE, ALLOCATABLE :: qvc(:,:)
+      !$OMP THREADPRIVATE(qvc)
+      REAL, SAVE, ALLOCATABLE :: rnebclr(:,:)
+      !$OMP THREADPRIVATE(rnebclr)
+      REAL, SAVE, ALLOCATABLE :: rnebss(:,:)
+      !$OMP THREADPRIVATE(rnebss)
+      REAL, SAVE, ALLOCATABLE :: gamma_ss(:,:)
+      !$OMP THREADPRIVATE(gamma_ss)
+      REAL, SAVE, ALLOCATABLE :: N1_ss(:,:)
+      !$OMP THREADPRIVATE(N1_ss)
+      REAL, SAVE, ALLOCATABLE :: N2_ss(:,:)
+      !$OMP THREADPRIVATE(N2_ss)
+      REAL, SAVE, ALLOCATABLE :: drneb_sub(:,:)
+      !$OMP THREADPRIVATE(drneb_sub)
+      REAL, SAVE, ALLOCATABLE :: drneb_con(:,:)
+      !$OMP THREADPRIVATE(drneb_con)
+      REAL, SAVE, ALLOCATABLE :: drneb_tur(:,:)
+      !$OMP THREADPRIVATE(drneb_tur)
+      REAL, SAVE, ALLOCATABLE :: drneb_avi(:,:)
+      !$OMP THREADPRIVATE(drneb_avi)
+      REAL, SAVE, ALLOCATABLE :: zqsatl(:,:)
+      !$OMP THREADPRIVATE(zqsatl)
+      REAL, SAVE, ALLOCATABLE :: zqsats(:,:)
+      !$OMP THREADPRIVATE(zqsats)
+      REAL, SAVE, ALLOCATABLE :: Tcontr(:,:)
+      !$OMP THREADPRIVATE(Tcontr)
+      REAL, SAVE, ALLOCATABLE :: qcontr(:,:)
+      !$OMP THREADPRIVATE(qcontr)
+      REAL, SAVE, ALLOCATABLE :: qcontr2(:,:)
+      !$OMP THREADPRIVATE(qcontr2)
+      REAL, SAVE, ALLOCATABLE :: fcontrN(:,:)
+      !$OMP THREADPRIVATE(fcontrN)
+      REAL, SAVE, ALLOCATABLE :: fcontrP(:,:)
+      !$OMP THREADPRIVATE(fcontrP)
+
 #ifdef CPP_StratAer
 !
@@ -828 +873 @@
 
       ALLOCATE(zn2mout(klon,6))
+
+! Supersaturation
+      ALLOCATE(rneb_seri(klon,klev))
+      ALLOCATE(d_rneb_dyn(klon,klev))
+      ALLOCATE(qclr(klon,klev), qcld(klon,klev), qss(klon,klev), qvc(klon,klev))
+      ALLOCATE(rnebclr(klon,klev), rnebss(klon,klev), gamma_ss(klon,klev))
+      ALLOCATE(N1_ss(klon,klev), N2_ss(klon,klev))
+      ALLOCATE(drneb_sub(klon,klev), drneb_con(klon,klev), drneb_tur(klon,klev), drneb_avi(klon,klev))
+      ALLOCATE(zqsatl(klon,klev), zqsats(klon,klev))
+      ALLOCATE(Tcontr(klon,klev), qcontr(klon,klev), qcontr2(klon,klev), fcontrN(klon,klev), fcontrP(klon,klev))
 
 #ifdef CPP_StratAer
@@ -1117 +1172 @@
       DEALLOCATE(zn2mout)
 
+! Supersaturation
+      DEALLOCATE(rneb_seri)
+      DEALLOCATE(d_rneb_dyn)
+      DEALLOCATE(qclr, qcld, qss, qvc)
+      DEALLOCATE(rnebclr, rnebss, gamma_ss)
+      DEALLOCATE(N1_ss, N2_ss)
+      DEALLOCATE(drneb_sub, drneb_con, drneb_tur, drneb_avi)
+      DEALLOCATE(zqsatl, zqsats)
+      DEALLOCATE(Tcontr, qcontr, qcontr2, fcontrN, fcontrP)
+
 #ifdef CPP_StratAer
 ! variables for strat. aerosol CK