Changeset 2507 for LMDZ5/trunk/libf/phylmd

Timestamp:

May 4, 2016, 4:36:48 PM (9 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.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• fisrtilp.F90 (modified) (4 diffs)
• icefrac_lsc_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/fisrtilp.F90

@@ -623 +623 @@
                     if (iflag_t_glace.eq.0) then
                     zdelta = MAX(0.,SIGN(1.,t_glace_min_old-Tbef(i)))
-                    else if (iflag_t_glace.eq.1) then
+                    else if (iflag_t_glace.ge.1) then
                     zdelta = MAX(0.,SIGN(1.,t_glace_min-Tbef(i)))
                     endif
@@ -673 +673 @@
                  else
                  ! Iteration pour convergence avec qsat(T)
-                 if (iflag_t_glace.eq.1) then
+                 if (iflag_t_glace.ge.1) then
                  CALL icefrac_lsc(klon,zt(:),pplay(:,k)/paprs(:,1),zfice(:))
                  endif
@@ -687 +687 @@
                  endif
                  
-                 if (iflag_t_glace.eq.1) then
+                 if (iflag_t_glace.ge.1) then
                  dzfice(i)= exposant_glace * zfice(i)**(exposant_glace-1) / (t_glace_min - t_glace_max)
                  endif
@@ -811 +811 @@
         endif
      ELSE
-         if (iflag_t_glace.eq.1) then
+         if (iflag_t_glace.ge.1) then
             CALL icefrac_lsc(klon,zt(:),pplay(:,k)/paprs(:,1),zfice(:))
          endif

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)