Changeset 5727 for LMDZ6/trunk/libf

Timestamp:

Jun 27, 2025, 4:33:03 PM (5 weeks ago)

Author:

evignon

Message:

legere modification de lscp pour enlever un warning obsolete sur la parameterisation de phase

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_lscp_ini.f90 (modified) (1 diff)
• lmdz_lscp_main.f90 (modified) (3 diffs)
• lmdz_lscp_tools.f90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.f90

@@ -615 +615 @@
     ENDIF
 
-    IF ( (iflag_icefrac .GE. 1) .AND. (.NOT. ok_poprecip .AND. (iflag_evap_prec .LT. 4)) ) THEN
-      abort_message = 'in lscp, icefracturb works with poprecip or with precip evap option >=4'
-      CALL abort_physic (modname,abort_message,1)
-    ENDIF
-
-
 
     !AA Temporary initialisation

LMDZ6/trunk/libf/phylmd/lmdz_lscp_main.f90

@@ -824 +824 @@
                      ENDDO
                      CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbef, zp, paprs(:,k), paprs(:,k+1), wls(:,k), zqi_ini,    &
-                     ziflcld, znebprecipcld, qincloud, zcf, tke(:,k), tke_dissip(:,k), sursat_e, invtau_e, zqliq, zqvapcl, zqice, & 
+                     ziflcld, qincloud, zcf, tke(:,k), tke_dissip(:,k), sursat_e, invtau_e, zqliq, zqvapcl, zqice, & 
                      zficeenv, dzficeenv, cldfraliq(:,k),sigma2_icefracturb(:,k),mean_icefracturb(:,k))                     
                      DO i=1,klon
@@ -853 +853 @@
                      ENDDO
                      CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbefth, zp, paprs(:,k), paprs(:,k+1), wth(:,k), zqi_ini,  &
-                     zeroklon, znebprecipcld, qincloud, zcf, zeroklon, zeroklon, sursat_e, invtau_e, zqliqth, zqvapclth, zqiceth, &
+                     zeroklon, qincloud, zcf, zeroklon, zeroklon, sursat_e, invtau_e, zqliqth, zqvapclth, zqiceth, &
                      zficeth, dzficeth,cldfraliqth(:,k), sigma2_icefracturbth(:,k), mean_icefracturbth(:,k))
                      !Environment
@@ -871 +871 @@
                      ENDDO
                      CALL icefrac_lscp_turb(klon, dtime, pticefracturb, Tbef, zp, paprs(:,k), paprs(:,k+1), wls(:,k), zqi_ini,    &
-                     ziflcld, znebprecipcld, qincloud, zcf, tke(:,k), tke_dissip(:,k), sursat_e, invtau_e, zqliq, zqvapcl, zqice, &
+                     ziflcld, qincloud, zcf, tke(:,k), tke_dissip(:,k), sursat_e, invtau_e, zqliq, zqvapcl, zqice, &
                      zfice, dzfice, cldfraliq(:,k),sigma2_icefracturb(:,k), mean_icefracturb(:,k)) 
   

LMDZ6/trunk/libf/phylmd/lmdz_lscp_tools.f90

@@ -241 +241 @@
 
 
-SUBROUTINE ICEFRAC_LSCP_TURB(klon, dtime, pticefracturb, temp, pplay, paprsdn, paprsup, wvel, qice_ini, snowcld, snowfracld, qtot_incl, cldfra, tke,   &
+SUBROUTINE ICEFRAC_LSCP_TURB(klon, dtime, pticefracturb, temp, pplay, paprsdn, paprsup, wvel, qice_ini, snowcld, qtot_incl, cldfra, tke,   &
                              tke_dissip, sursat_e, invtau_e, qliq, qvap_cld, qice, icefrac, dicefracdT, cldfraliq, sigma2_icefracturb, mean_icefracturb)
 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
@@ -275 +275 @@
    REAL,      INTENT(IN),       DIMENSION(klon)    :: qice_ini            !--initial specific ice content gridbox-mean     [kg/kg]
    REAL,      INTENT(IN),       DIMENSION(klon)    :: snowcld             !--in-cloud snowfall flux                        [kg/m2/s]
-   REAL,      INTENT(IN),       DIMENSION(klon)    :: snowfracld          !--cloudy precip fraction                        [-]
    REAL,      INTENT(IN),       DIMENSION(klon)    :: sursat_e            !--environment supersaturation                   [-]
    REAL,      INTENT(IN),       DIMENSION(klon)    :: invtau_e            !--inverse time-scale of mixing with environment [s-1]
@@ -376 +375 @@
            ! we consider here the mean snowflake concentration in the mesh (not the in-cloud concentration)
            ! when poprecip is active, it will be worth testing considering the incloud fraction, dividing 
-           ! by snowfracld     
+           ! by znebprecipcld     
            ! qiceini_incl  = qice_ini(i) / cldfra1D + &
            !              gamma_snwretro * snowcld(i) * RG * dtime / ( paprsdn(i) - paprsup(i) )