Changeset 1286 for LMDZ4/trunk/libf/phylmd/newmicro.F

Timestamp:

Dec 17, 2009, 2:47:10 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Les parametres definissant la transition eau glace dans les nuages sont
maintenant lus dans physiq.def. Les valeurs par defaut donnent les memes
resultats que precedemment. JLD

File:

• LMDZ4/trunk/libf/phylmd/newmicro.F (modified) (3 diffs)

LMDZ4/trunk/libf/phylmd/newmicro.F

@@ -73 +73 @@
       REAL coef, coef_froi, coef_chau
       PARAMETER (coef_chau=0.13, coef_froi=0.09)
-      REAL seuil_neb, t_glace
-      PARAMETER (seuil_neb=0.001, t_glace=273.0-15.0)
+      REAL seuil_neb
+      PARAMETER (seuil_neb=0.001)
       INTEGER nexpo ! exponentiel pour glace/eau
       PARAMETER (nexpo=6)
@@ -143 +143 @@
          DO k = 1, klev
             DO i = 1, klon
-               zfice2(i,k) = 1.0 - (t(i,k)-t_glace) / (273.13-t_glace)
+c               zfice2(i,k) = 1.0 - (t(i,k)-t_glace) / (273.13-t_glace)
+               zfice2(i,k) = 1.0 - (t(i,k)-t_glace_min) / 
+     &                             (t_glace_max-t_glace_min)
                zfice2(i,k) = MIN(MAX(zfice2(i,k),0.0),1.0)
 c     IM Total Liquid/Ice water content                                    
@@ -339 +341 @@
                   
                   zfice1 = MIN(
-     &                 MAX( 1.0 - (t(i,k)-t_glace) / (273.13-t_glace)
-     &                 ,0.0),1.0)**nexpo
+     &                 MAX( 1.0 - (t(i,k)-t_glace_min) / 
+     &                    (t_glace_max-t_glace_min),0.0),1.0)**nexpo
                   
                   radius = rad_chaud * (1.-zfice1) + rad_froid * zfice1