Changeset 4535 for LMDZ6/trunk/libf/phylmd/lscp_mod.F90

Timestamp:

May 13, 2023, 2:12:05 PM (17 months ago)

Author:

evignon

Message:

poursuite de la replay-isation de lscp en vue de la session
de reecriture de lscp_mod en juin

File:

• LMDZ6/trunk/libf/phylmd/lscp_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/lscp_mod.F90

@@ -95 +95 @@
 USE lscp_ini_mod, ONLY : seuil_neb, ninter, iflag_evap_prec, t_coup, DDT0, ztfondue, rain_int_min
 USE lscp_ini_mod, ONLY : iflag_mpc_bl, ok_radocond_snow, a_tr_sca, ok_debug_autoconversion
-
+USE lscp_ini_mod, ONLY : iflag_cloudth_vert, iflag_rain_incloud_vol, iflag_t_glace, t_glace_min 
+USE lscp_ini_mod, ONLY : coef_eva, coef_eva_i,cld_tau_lsc, cld_tau_con, cld_lc_lsc, cld_lc_con
+USE lscp_ini_mod, ONLY : iflag_bergeron, iflag_fisrtilp_qsat
+USE lscp_ini_mod, ONLY : RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT, RD, RG
 
 IMPLICIT NONE
@@ -102 +105 @@
 ! VARIABLES DECLARATION
 !===============================================================================
-
-include "YOMCST.h"
-include "YOETHF.h"
-include "fisrtilp.h"
-include "nuage.h" 
 
 ! INPUT VARIABLES:
@@ -726 +724 @@
                   ! saturation may occur at a humidity different from qsat (gamma qsat), so gamma correction for dqs
                   zdqs(:) = gammasat(:)*zdqs(:)+zqs(:)*dgammasatdt(:)
-                  CALL icefrac_lscp(klon, zt(:),pplay(:,k)/paprs(:,1),zfice(:),dzfice(:))
+                  CALL icefrac_lscp(klon, zt(:), iflag_ice_thermo, pplay(:,k)/paprs(:,1),zfice(:),dzfice(:))
 
                   DO i=1,klon !todoan : check if loop in i is needed
@@ -821 +819 @@
 
         ! Partition function in stratiform clouds (will be overwritten in boundary-layer MPCs)
-        CALL icefrac_lscp(klon,zt(:),pplay(:,k)/paprs(:,1),zfice(:), dzfice(:))
+        CALL icefrac_lscp(klon,zt(:),iflag_ice_thermo,pplay(:,k)/paprs(:,1),zfice(:), dzfice(:))
 
         DO i=1, klon