Changeset 4072 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

Feb 1, 2022, 6:34:34 PM (3 years ago)

Author:

evignon

Message:

Mise à jour de la routine de nuages lscp
les principaux changements consistent en:

• des corrections de bug (déclaration et division

par zero) pour que la routine compile avec les modifications d'OB

• des restructurations pour que les fonctions de lscp_tools_mod

travaillent sur des tableaux et non des scalaires

• une séparation des coefficients d'évaporation de la neige et de la pluie

File:

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

LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90

@@ -559 +559 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: cldemi, cldfra, cldtau, fiwc, fl, re, flwc
 !$OMP THREADPRIVATE(cldemi, cldfra, cldtau, fiwc, fl, re, flwc)
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: qlth, qith
-!$OMP THREADPRIVATE(qlth, qith)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: qlth, qith, qsith
+!$OMP THREADPRIVATE(qlth, qith, qsith)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq, ref_ice, theta, zphi
 !$OMP THREADPRIVATE(ref_liq, ref_ice, theta, zphi)
@@ -943 +943 @@
       ALLOCATE(rain_lsc(klon))
       ALLOCATE(rain_num(klon))
-      ALLOCATE(qlth(klon,klev), qith(klon,klev))
+      ALLOCATE(qlth(klon,klev), qith(klon,klev), qsith(klon,klev))
       !
 #ifdef ISO
@@ -1004 +1004 @@
       ALLOCATE(ref_liq_pi(klon, klev), ref_ice_pi(klon, klev))
       ALLOCATE(zphi(klon, klev), zx_rh(klon, klev), zx_rhl(klon,klev), zx_rhi(klon,klev))
+      zx_rhl(:,:)=0.; zx_rhi(:,:)=0. ! because not always defined
       ALLOCATE(pmfd(klon, klev), pmfu(klon, klev))
 
@@ -1324 +1325 @@
       DEALLOCATE(rain_lsc)
       DEALLOCATE(rain_num)
-      DEALLOCATE(qlth, qith)
+      DEALLOCATE(qlth, qith, qsith)
 !
       DEALLOCATE(sens_x, sens_w)