Changeset 5238 for LMDZ6/trunk

Timestamp:

Oct 14, 2024, 5:21:07 PM (2 months ago)

Author:

Laurent Fairhead

Message:

Temporary fix to correct divergence of results between revision 5210 and 5211 until we understand why
LMDZ workshop (EV, AI, LG, LF)

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_lscp.F90 (modified) (2 diffs)
• lmdz_lscp_condensation.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/lmdz_lscp.F90

@@ -106 +106 @@
 USE lmdz_lscp_tools, ONLY : icefrac_lscp, icefrac_lscp_turb
 USE lmdz_lscp_tools, ONLY : fallice_velocity, distance_to_cloud_top
-USE lmdz_lscp_condensation, ONLY : condensation_lognormal, condensation_ice_supersat
+USE lmdz_lscp_condensation, ONLY : condensation_lognormal_test, condensation_ice_supersat
 USE lmdz_lscp_poprecip, ONLY : poprecip_precld, poprecip_postcld
 
@@ -971 +971 @@
                   !--generalised lognormal condensation scheme (Bony and Emanuel 2001)
 
-                   CALL condensation_lognormal( &
-                       klon, Tbef, zq, zqs, gammasat, ratqs(:,k), &
-                       keepgoing, rneb(:,k), zqn, qvc)
+                   CALL condensation_lognormal_test( &
+                       klon, klev, k, keepgoing, ratqs, gammasat, zq, zqs, &
+                       zqn, qvc, rneb)
 
 

LMDZ6/trunk/libf/phylmd/lmdz_lscp_condensation.F90

@@ -8 +8 @@
 
 CONTAINS
+
+!*********************************************************************************
+SUBROUTINE condensation_lognormal_test(klon,klev,k,keepgoing,ratqs, gammasat,zq, zqs, zqn, qvc, rneb)
+
+
+IMPLICIT NONE
+
+INTEGER, INTENT(IN)   :: klon,klev,k
+LOGICAL,  INTENT(IN)   , DIMENSION(klon) :: keepgoing     ! .TRUE. if a new condensation loop should be computed
+REAL, INTENT(IN), DIMENSION(klon,klev)   :: ratqs
+REAL, INTENT(IN), DIMENSION(klon)        :: zq, zqs,gammasat
+REAL, INTENT(INOUT), DIMENSION(klon)     :: zqn, qvc
+REAL, INTENT(INOUT), DIMENSION(klon,klev) :: rneb
+
+REAL, DIMENSION(klon) :: zpdf_sig,zpdf_k,zpdf_delta             ! lognormal parameters
+REAL, DIMENSION(klon) :: Zpdf_a,zpdf_b,zpdf_e1,zpdf_e2          ! lognormal intermediate variables
+INTEGER               :: i
+
+
+DO i=1,klon !todoan : check if loop in i is needed
+
+    IF (keepgoing(i)) THEN
+
+    zpdf_sig(i)=ratqs(i,k)*zq(i)
+    zpdf_k(i)=-sqrt(log(1.+(zpdf_sig(i)/zq(i))**2))
+    zpdf_delta(i)=log(zq(i)/(gammasat(i)*zqs(i)))
+    zpdf_a(i)=zpdf_delta(i)/(zpdf_k(i)*sqrt(2.))
+    zpdf_b(i)=zpdf_k(i)/(2.*sqrt(2.))
+    zpdf_e1(i)=zpdf_a(i)-zpdf_b(i)
+    zpdf_e1(i)=sign(min(ABS(zpdf_e1(i)),5.),zpdf_e1(i))
+    zpdf_e1(i)=1.-erf(zpdf_e1(i))
+    zpdf_e2(i)=zpdf_a(i)+zpdf_b(i)
+    zpdf_e2(i)=sign(min(ABS(zpdf_e2(i)),5.),zpdf_e2(i))
+    zpdf_e2(i)=1.-erf(zpdf_e2(i))
+
+    IF (zpdf_e1(i).LT.1.e-10) THEN
+        rneb(i,k)=0.
+        zqn(i)=zqs(i)
+        !--AB grid-mean vapor in the cloud - we assume saturation adjustment
+        qvc(i) = 0.
+    ELSE
+        rneb(i,k)=0.5*zpdf_e1(i)
+        zqn(i)=zq(i)*zpdf_e2(i)/zpdf_e1(i)
+        !--AB grid-mean vapor in the cloud - we assume saturation adjustment
+        qvc(i) = rneb(i,k) * zqs(i)
+    ENDIF
+
+    ENDIF ! keepgoing
+ENDDO ! loop on klon
+
+
+END SUBROUTINE condensation_lognormal_test
+
 
 !**********************************************************************************