Changeset 2077 for LMDZ5/trunk/libf/phylmd/nuage.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/nuage.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/nuage.F90

@@ -5 +5 @@
     cldtaupi, re, fl)
   USE dimphy
-  USE microphys_mod ! cloud microphysics (JBM 3/14)
+  USE icefrac_lsc_mod ! cloud microphysics (JBM 3/14)
   IMPLICIT NONE
   ! ======================================================================
@@ -54 +54 @@
 
   INTEGER i, k
-  REAL zflwp, zradef, zfice, zmsac
+  REAL zflwp, zradef, zfice(klon), zmsac
 
   REAL radius, rad_chaud
@@ -98 +98 @@
 
   DO k = 1, klev
+     IF (iflag_t_glace.EQ.0) THEN
+       DO i = 1, klon
+        zfice(i) = 1.0 - (t(i,k)-t_glace_min_old)/(273.13-t_glace_min_old)
+        zfice(i) = min(max(zfice(i),0.0), 1.0)
+        zfice(i) = zfice(i)**exposant_glace_old
+       ENDDO
+     ELSE ! of IF (iflag_t_glace.EQ.0)
+! JBM: icefrac_lsc is now a function contained in microphys_mod
+!       zfice(i) = icefrac_lsc(t(i,k), t_glace_min, &
+!                           t_glace_max, exposant_glace)
+         CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1), &
+&               t_glace_min,t_glace_max,exposant_glace,zfice(:))
+     ENDIF
+
     DO i = 1, klon
       rad_chaud = rad_chau1
@@ -104 +118 @@
       pclc(i, k) = max(pclc(i,k), seuil_neb)
       zflwp = 1000.*pqlwp(i, k)/rg/pclc(i, k)*(paprs(i,k)-paprs(i,k+1))
-      IF (iflag_t_glace.EQ.0) THEN
-        zfice = 1.0 - (t(i,k)-t_glace_min_old)/(273.13-t_glace_min_old)
-        zfice = min(max(zfice,0.0), 1.0)
-        zfice = zfice**exposant_glace_old
-      ELSE ! of IF (iflag_t_glace.EQ.0)
-! JBM: icefrac_lsc is now a function contained in microphys_mod
-        zfice = icefrac_lsc(t(i,k), t_glace_min, &
-                            t_glace_max, exposant_glace)
-      ENDIF
 
       IF (ok_aie) THEN
@@ -150 +155 @@
           ! effective radius as re/fl
           !
-        fl(i, k) = pclc(i, k)*(1.-zfice)
+        fl(i, k) = pclc(i, k)*(1.-zfice(i))
         re(i, k) = rad_chaud*fl(i, k)
 
@@ -159 +164 @@
           ! cdnc.
         radius = max(1.1E6*((pqlwp(i,k)*pplay(i,k)/(rd*t(i,k)))/(4./3.*rpi* &
-          1000.*cdnc_pi(i,k)))**(1./3.), 3.)*(1.-zfice) + rad_froid*zfice
+          1000.*cdnc_pi(i,k)))**(1./3.), 3.)*(1.-zfice(i)) + rad_froid*zfice(i)
         cldtaupi(i, k) = 3.0/2.0*zflwp/radius
       END IF ! ok_aie
 
-      radius = rad_chaud*(1.-zfice) + rad_froid*zfice
-      coef = coef_chau*(1.-zfice) + coef_froi*zfice
+      radius = rad_chaud*(1.-zfice(i)) + rad_froid*zfice(i)
+      coef = coef_chau*(1.-zfice(i)) + coef_froi*zfice(i)
       pcltau(i, k) = 3.0/2.0*zflwp/radius
       pclemi(i, k) = 1.0 - exp(-coef*zflwp)