Changeset 3999 for LMDZ6/trunk/libf/phylmd/nuage.F90

Timestamp:

Nov 5, 2021, 12:40:08 PM (3 years ago)

Author:

evignon

Message:

commission de la nouvelle routine de condensation
grande echelle simplifiee (lscp, version epuree de fistilp)
et du schema de nuages de phase mixte (en developpement)

La routine lscp n'est active que sous flag
ok_new_lscp=y

File:

• LMDZ6/trunk/libf/phylmd/nuage.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/nuage.F90

@@ -1 +1 @@
 ! $Id$
 
-SUBROUTINE nuage(paprs, pplay, t, pqlwp, pclc, pcltau, pclemi, pch, pcl, pcm, &
+SUBROUTINE nuage(paprs, pplay, t, pqlwp,picefra, pclc, pcltau, pclemi, pch, pcl, pcm, &
     pct, pctlwp, ok_aie, mass_solu_aero, mass_solu_aero_pi, bl95_b0, bl95_b1, &
     cldtaupi, re, fl)
   USE dimphy
+  USE lscp_tools_mod, only: icefrac_lscp
   USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14)
+  USE phys_local_var_mod, ONLY: ptconv
   IMPLICIT NONE
   ! ======================================================================
@@ -14 +16 @@
   ! t-------input-R-temperature
   ! pqlwp---input-R-eau liquide nuageuse dans l'atmosphere (kg/kg)
+  ! picefra--inout-R-fraction de glace dans les nuages (-)
   ! pclc----input-R-couverture nuageuse pour le rayonnement (0 a 1)
   ! ok_aie--input-L-apply aerosol indirect effect or not
@@ -36 +39 @@
   include "YOMCST.h"
   include "nuage.h" ! JBM 3/14
+  include "clesphys.h"
 
   REAL paprs(klon, klev+1), pplay(klon, klev)
@@ -41 +45 @@
 
   REAL pclc(klon, klev)
-  REAL pqlwp(klon, klev)
+  REAL pqlwp(klon, klev), picefra(klon,klev)
   REAL pcltau(klon, klev), pclemi(klon, klev)
 
@@ -89 +93 @@
 
   REAL cldtaupi(klon, klev) ! pre-industrial cloud opt thickness for diag
+  REAl dzfice(klon)
   ! jq-end
 
@@ -106 +111 @@
 !       zfice(i) = icefrac_lsc(t(i,k), t_glace_min, &
 !                           t_glace_max, exposant_glace)
-         CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:))
+        IF (ok_new_lscp) THEN
+            CALL icefrac_lscp(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:),dzfice(:))
+        ELSE
+            CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:))
+
+        ENDIF
+
+        IF ((.NOT. ptconv(i,k)) .AND. ok_new_lscp .AND. ok_icefra_lscp) THEN
+        ! EV: take the ice fraction directly from the lscp code
+        ! consistent only for non convective grid points
+        ! critical for mixed phase clouds
+            DO i=1,klon
+            zfice(i)=picefra(i,k)
+            ENDDO
+        ENDIF
+
+
      ENDIF