Changeset 2507 for LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90

Timestamp:

May 4, 2016, 4:36:48 PM (8 years ago)

Author:

jbmadeleine

Message:

Added the iflag_t_glace = 2 option that turns off the convertion of cloud
water to ice close to the surface when T< t_glace_max (273.15K). This option
is now obsolete because the liquid precipitation is now converted to ice when the
temperature is negative using the iflag_bergeron option.

File:

• LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90

@@ -14 +14 @@
   ! Doutriaux-Boucher & Quaas 2004, section 2.2.)
   !
-  ! (JBM 3/14 8/14)
+  ! (JBM 3/14 8/14 5/16)
   
+  USE print_control_mod, ONLY: lunout, prt_level
   INCLUDE "nuage.h"
 
@@ -33 +34 @@
 
   DO ip=1,np
-     www=(max(sig(ip)-sig0,0.))/(1.-sig0) ! w=1 at the surface and 0 for sig < sig0
+     IF (iflag_t_glace.EQ.1) THEN
+       ! Transition to ice close to surface for T<Tmax
+       ! w=1 at the surface and 0 for sig < sig0
+       www=(max(sig(ip)-sig0,0.))/(1.-sig0) 
+     ELSEIF (iflag_t_glace.GE.2) THEN
+       ! No convertion to ice close to surface
+       www = 0.
+     ENDIF
      tmin_tmp=www*t_glace_max+(1.-www)*t_glace_min
      icefrac_tmp= 1.0 - (temp(ip)-tmin_tmp) / (t_glace_max-tmin_tmp)