Changeset 5577 for LMDZ6/branches/contrails/libf/phylmd

Timestamp:

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

Author:

aborella

Message:

Bugfix for ice fallspeed velocity (clipping too restrictive)

Location:

LMDZ6/branches/contrails/libf/phylmd

Files:

• lmdz_lscp_ini.f90 (modified) (3 diffs)
• lmdz_lscp_tools.f90 (modified) (2 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

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

@@ -17 +17 @@
     
     use lmdz_lscp_ini, only: iflag_vice, ffallv_con, ffallv_lsc
-    use lmdz_lscp_ini, only: cice_velo, dice_velo 
+    use lmdz_lscp_ini, only: cice_velo, dice_velo
+    use lmdz_lscp_ini, only: ok_bug_ice_fallspeed
 
     IMPLICIT NONE
@@ -46 +47 @@
 
         tempc=temp(i)-273.15 ! celcius temp
-        iwcg=MAX(iwc(i)*1000.,1E-3) ! iwc in g/m3. We set a minimum value to prevent from division by 0
+        IF ( ok_bug_ice_fallspeed ) THEN
+            iwcg=MAX(iwc(i)*1000.,1E-3) ! iwc in g/m3. We set a minimum value to prevent from division by 0
+        ELSE
+            ! AB the threshold is way too high, we reduce it
+            iwcg=MAX(iwc(i)*1000.,1E-10) ! iwc in g/m3. We set a minimum value to prevent from division by 0
+        ENDIF
         phpa=pres(i)/100.    ! pressure in hPa