Index: LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90
===================================================================
--- LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90	(revision 5576)
+++ LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_ini.f90	(revision 5577)
@@ -93,4 +93,7 @@
   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,4 +432,5 @@
     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,4 +533,5 @@
     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
Index: LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_tools.f90
===================================================================
--- LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_tools.f90	(revision 5576)
+++ LMDZ6/branches/contrails/libf/phylmd/lmdz_lscp_tools.f90	(revision 5577)
@@ -17,5 +17,6 @@
     
     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,5 +47,10 @@
 
         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