Changeset 3999 for LMDZ6/trunk/libf/phylmd/newmicro.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/newmicro.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/newmicro.F90

@@ -1 +1 @@
 ! $Id$
 
-SUBROUTINE newmicro(flag_aerosol, ok_cdnc, bl95_b0, bl95_b1, paprs, pplay, t, pqlwp, pclc, &
+SUBROUTINE newmicro(flag_aerosol, ok_cdnc, bl95_b0, bl95_b1, paprs, pplay, t, pqlwp, picefra, pclc, &
     pcltau, pclemi, pch, pcl, pcm, pct, pctlwp, xflwp, xfiwp, xflwc, xfiwc, &
     mass_solu_aero, mass_solu_aero_pi, pcldtaupi, re, fl, reliq, reice, &
@@ -9 +9 @@
   USE phys_local_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, &
       reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra,  & 
-      zfice, dNovrN
+      zfice, dNovrN, ptconv
   USE phys_state_var_mod, ONLY: rnebcon, clwcon
   USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14)
   USE ioipsl_getin_p_mod, ONLY : getin_p
   USE print_control_mod, ONLY: lunout
+  USE lscp_tools_mod, only: icefrac_lscp
+
 
 
@@ -31 +33 @@
   ! pqlwp---input-R-eau liquide nuageuse dans l'atmosphere dans la partie
   ! nuageuse (kg/kg)
+  ! picefra--input-R-fraction de glace dans les nuages
   ! pclc----input-R-couverture nuageuse pour le rayonnement (0 a 1)
   ! mass_solu_aero-----input-R-total mass concentration for all soluble
@@ -58 +61 @@
   include "radepsi.h"
   include "radopt.h"
+  include "clesphys.h"
 
   ! choix de l'hypothese de recouvrement nuageuse via radopt.h (IM, 19.07.2016)
@@ -81 +85 @@
   REAL t(klon, klev)
   REAL pclc(klon, klev)
-  REAL pqlwp(klon, klev)
+  REAL pqlwp(klon, klev), picefra(klon,klev)
   REAL pcltau(klon, klev)
   REAL pclemi(klon, klev)
@@ -148 +152 @@
   ! jq-end
   ! IM cf. CR:parametres supplementaires
+  REAL dzfice(klon,klev)
   REAL zclear(klon)
   REAL zcloud(klon)
@@ -229 +234 @@
   ELSE ! of IF (iflag_t_glace.EQ.0)
     DO k = 1, klev
-        CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:,k))
- 
-
-        ! JBM: icefrac_lsc is now contained icefrac_lsc_mod
+
+! JBM: icefrac_lsc is now contained icefrac_lsc_mod
 !       zfice(i, k) = icefrac_lsc(t(i,k), t_glace_min, &
 !                                 t_glace_max, exposant_glace)
-      DO i = 1, klon
+
+      IF (ok_new_lscp) THEN
+          CALL icefrac_lscp(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:,k),dzfice(:,k))
+      ELSE
+          CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:,k))
+      ENDIF
+
+      DO i = 1, klon
+        
+        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
+            zfice(i,k)=picefra(i,k)
+        ENDIF
+
         ! -layer calculation
         rhodz(i, k) = (paprs(i,k)-paprs(i,k+1))/rg ! kg/m2