Changeset 5425 for LMDZ6/trunk/libf

Timestamp:

Dec 18, 2024, 1:12:33 PM (13 hours ago)

Author:

jyg

Message:

Correction des arguments d'appel de phytrac et de phytracr_spl dans
physiq_mod.F90 : remplacement de wdtrainA par wdtrainAS (avec, auparavant,
calcul de wdtrainAS=wdtrainA+wdtrainS). Le bug date de la mise en oeuvre
(le 10mai 2019 : svn3499) de l'ejection des precipitations liquides dans
le schema convectif, wdtrainA ayant alors change de signification.

De svn3499 a svn5412, le lessivage convectif est completement faux. A

partir de la presente release il y a a nouveau conservation des traceurs.

Location:

LMDZ6/trunk/libf/phylmd

Files:

• phys_local_var_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -512 +512 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: pmflxr, pmflxs
 !$OMP THREADPRIVATE(pmflxr, pmflxs)
-      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: wdtrainA, wdtrainS, wdtrainM
-!$OMP THREADPRIVATE(wdtrainA, wdtrainS, wdtrainM)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: wdtrainA, wdtrainS, wdtrainM, wdtrainAS
+!$OMP THREADPRIVATE(wdtrainA, wdtrainS, wdtrainM, wdtrainAS)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:)  :: da, mp
 !$OMP THREADPRIVATE(da, mp)
@@ -1125 +1125 @@
 !  Deep convective variables used in phytrac
       ALLOCATE(pmflxr(klon, klev+1), pmflxs(klon, klev+1))
-      ALLOCATE(wdtrainA(klon,klev),wdtrainS(klon,klev),wdtrainM(klon,klev))
+      ALLOCATE(wdtrainA(klon,klev),wdtrainS(klon,klev),wdtrainM(klon,klev),wdtrainAS(klon,klev))
       ALLOCATE(dnwd(klon, klev), upwd(klon, klev))
       ALLOCATE(ep(klon,klev))                          ! epmax_cape
@@ -1540 +1540 @@
 
       DEALLOCATE(pmflxr, pmflxs)
-      DEALLOCATE(wdtrainA, wdtrainS, wdtrainM)
+      DEALLOCATE(wdtrainA, wdtrainS, wdtrainM, wdtrainAS)
       DEALLOCATE(upwd, dnwd)
       DEALLOCATE(ep)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -301 +301 @@
        !    Deep convective variables used in phytrac
        pmflxr, pmflxs,  &
-       wdtrainA, wdtrainS, wdtrainM,  &
+       wdtrainA, wdtrainS, wdtrainM, wdtrainAS,  &
        upwd, dnwd, &
        ep,  &
@@ -3084 +3084 @@
     wdtrainS(:,:) = 0.
     wdtrainM(:,:) = 0.
+    wdtrainAS(:,:) = 0.
     upwd(:,:) = 0.
     dnwd(:,:) = 0.
@@ -5318 +5319 @@
     ENDIF
 
+    ! Merge wdtrainA and wdtrainS in the total source of precipitation due to
+    ! adiabatic updraughts.
+    !
+    wdtrainAS(:,:) = wdtrainA(:,:) + wdtrainS(:,:)
+
 IF (CPPKEY_DUST) THEN
     !  Avec SPLA, iflag_phytrac est forcé =1
@@ -5329 +5335 @@
                       da,phi,phi2,d1a,dam,mp,ep,sigd,sij,clw,elij,     &  ! I
                       epmlmMm,eplaMm,upwd,dnwd,itop_con,ibas_con,      &  ! I
-                      ev,wdtrainA,  wdtrainM,wght_cvfd,              &  ! I
+                      ev,wdtrainAS,  wdtrainM,wght_cvfd,              &  ! I
                       fm_therm, entr_therm, rneb,                      &  ! I
                       beta_prec_fisrt,beta_prec, & !I
@@ -5354 +5360 @@
          da,       phi,       mp,       upwd, &
          phi2,     d1a,       dam,      sij, wght_cvfd, &        !<<RomP+RL
-         wdtrainA, wdtrainM,  sigd,     clw,elij, &   !<<RomP
+         wdtrainAS, wdtrainM,  sigd,     clw,elij, &   !<<RomP
          ev,       ep,        epmlmMm,  eplaMm, &     !<<RomP
          dnwd,     aerosol_couple,      flxmass_w, &