Changeset 5577 for LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90

Timestamp:

Mar 15, 2025, 11:33:38 PM (3 months ago)

Author:

aborella

Message:

Bugfix for ice fallspeed velocity (clipping too restrictive)

File:

• LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90 (modified) (3 diffs)

LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90

@@ -93 +93 @@
   LOGICAL, SAVE, PROTECTED :: ok_bug_phase_lscp=.true.      ! bug on phase partitioning after precipitation processes
   !$OMP THREADPRIVATE(ok_bug_phase_lscp)
+
+  LOGICAL, SAVE, PROTECTED :: ok_bug_ice_fallspeed=.true.   ! flag to activate the high clipping of iwc when calculating ice  fallspeed velocity
+  !$OMP THREADPRIVATE(ok_bug_ice_fallspeed)
 
   LOGICAL, SAVE, PROTECTED :: reevap_ice=.false.            ! no liquid precip for T< threshold
@@ -429 +432 @@
     CALL getin_p('temp_nowater',temp_nowater)
     CALL getin_p('ok_bug_phase_lscp',ok_bug_phase_lscp)
+    CALL getin_p('ok_bug_ice_fallspeed',ok_bug_ice_fallspeed)
     ! for poprecip
     CALL getin_p('ok_poprecip',ok_poprecip)
@@ -529 +533 @@
     WRITE(lunout,*) 'lscp_ini, temp_nowater', temp_nowater
     WRITE(lunout,*) 'lscp_ini, ok_bug_phase_lscp', ok_bug_phase_lscp
+    WRITE(lunout,*) 'lscp_ini, ok_bug_ice_fallspeed', ok_bug_ice_fallspeed
     ! for poprecip
     WRITE(lunout,*) 'lscp_ini, ok_poprecip', ok_poprecip