Changeset 1276 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/cv3_routines.F

Timestamp:

Dec 7, 2009, 3:30:39 PM (15 years ago)

Author:

musat

Message:

Corrections conservation de l'eau dans la nouvelle physique
JYG/IM

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/cv3_routines.F (modified) (4 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/cv3_routines.F

@@ -72 +72 @@
 c      tau    = 3000.
 cc      tau = 1800.
+c     tau= 2800.
       tau=8000.
       beta   = 1.0 - delt/tau
@@ -2086 +2087 @@
 cc---end jyg---
 c
-c--------retour à la formulation originale d'Emanuel.
+c--------retour à la formulation originale d''Emanuel.
       b6=bfac*50.*sigd(il)*(ph(il,i)-ph(il,i+1))*sigt*afac
       c6=water(il,i+1)+bfac*wdtrain(il)
@@ -2092 +2093 @@
       if(c6.gt.0.0)then
        revap=0.5*(-b6+sqrt(b6*b6+4.*c6))
-       evap(il,i)=sigt*afac*revap
-       water(il,i)=revap*revap
+       water(il,i)=revap*revap      !equation de conservation
       else
+       water(il,i) = 0.
+      endif
 cJYG/IM : ci-dessous formulation originale de KE
 c      evap(il,i)=-evap(il,i+1)
@@ -2104 +2106 @@
 c        lorsqu'il y a evaporation totale de la pluie
 c
-       evap(il,i)= +(wdtrain(il)+sigd(il)*wt(il,i)*water(il,i+1))
+c       evap(il,i)= +(wdtrain(il)+sigd(il)*wt(il,i)*water(il,i+1)) !itlmd(jyg)
+c     :                 /(sigd(il)*(ph(il,i)-ph(il,i+1))*100.)
+c      end if  !itlmd(jyg)
+cjyg---   Dans tous les cas, evaporation = [tt ce qui entre dans la couche i]
+c                                    moins [tt ce qui sort de la couche i]
+       evap(il,i)=
+     :       (wdtrain(il)+sigd(il)*wt(il,i)*(water(il,i+1)-water(il,i)))
      :                 /(sigd(il)*(ph(il,i)-ph(il,i+1))*100.)
-      end if
+c
 ccc
 c    ***  calculate precipitating downdraft mass flux under     ***