Changeset 2077 for LMDZ5/trunk/libf/phylmd/newmicro.F90

Timestamp:

Jul 3, 2014, 3:38:05 PM (10 years ago)

Author:

fhourdin

Message:

Modification relative à la phase mixte des nuages :

1. on fait tendre t_glace_min vers t_glace_max (en principe 0°C) linéairerment entre p/ps = 0.8 et p/ps=1.
2. Passage de tableau à la routine icefrac_lsc en remplacement d'une fonction scalaire.
3. Changement de nom pour le module (le module icefrac_lsc_mod.F90 contenant maintenant la routine icefrac_lsc).

Mofication concerning the mixte phase of clouds

1. t_glace_min -> t_glace_max when p/ps : 0.8 -> 1
2. passing arrays instead of scalars to icefrac_lsc (which becomes a routine rather than a function).
3. Changing the name microphys_mod.F90 -> icerfrac_lsc_mod.F90

File:

• LMDZ5/trunk/libf/phylmd/newmicro.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/newmicro.F90

@@ -12 +12 @@
     reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra
   USE phys_state_var_mod, ONLY: rnebcon, clwcon
-  USE microphys_mod ! cloud microphysics (JBM 3/14)
+  USE icefrac_lsc_mod ! cloud microphysics (JBM 3/14)
   IMPLICIT NONE
   ! ======================================================================
@@ -202 +202 @@
   ELSE ! of IF (iflag_t_glace.EQ.0)
     DO k = 1, klev
+        CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1), &
+&               t_glace_min,t_glace_max,exposant_glace,zfice(:,k))
+ 
+
+        ! JBM: icefrac_lsc is now a function contained in microphys_mod
+!       zfice(i, k) = icefrac_lsc(t(i,k), t_glace_min, &
+!                                 t_glace_max, exposant_glace)
       DO i = 1, klon
         ! -layer calculation
@@ -207 +214 @@
         zrho(i, k) = pplay(i, k)/t(i, k)/rd ! kg/m3
         dh(i, k) = rhodz(i, k)/zrho(i, k) ! m
-        ! JBM: icefrac_lsc is now a function contained in microphys_mod
-        zfice(i, k) = icefrac_lsc(t(i,k), t_glace_min, &
-                                  t_glace_max, exposant_glace)
         ! -IM Total Liquid/Ice water content
         xflwc(i, k) = (1.-zfice(i,k))*pqlwp(i, k)