Changeset 1322 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cv3_routines.F

Timestamp:

Mar 12, 2010, 11:54:11 AM (14 years ago)

Author:

Laurent Fairhead

Message:

Improvements concerning wake parametrisation (from JYG, NR, IT, with more to come).
Alp_offset is read in form physiq.def file

---

Améliorations à la paramétrisation des poches froides (de JYG, NR, IT, d'autres
sont à venir)
Alp_offset est rajouté à la liste des paramètres lus dans physiq.def

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/cv3_routines.F (modified) (8 diffs)

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

@@ -74 +74 @@
 c      dtcrit = -5.0
 c      tau    = 3000.
-cc      tau = 1800.
-c     tau= 2800.
-      tau=8000.
+      tau = 1800.
+cc      tau=8000.
       beta   = 1.0 - delt/tau
       alpha1 = 1.5e-3
@@ -2094 +2093 @@
 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)
@@ -2100 +2099 @@
       if(c6.gt.0.0)then
        revap=0.5*(-b6+sqrt(b6*b6+4.*c6))
-       water(il,i)=revap*revap      !equation de conservation
+cjyg    Dans sa formulation originale, Emanuel calcule l'evaporation par:
+cc             evap(il,i)=sigt*afac*revap
+c     ce qui n'est pas correct. Dans cv_routines, la formulation a été modifiee.
+c     Ici,l'evaporation evap est simplement calculee par l'equation de
+c     conservation.
+       water(il,i)=revap*revap
       else
+cjyg----   Correction : si c6 <= 0, water(il,i)=0.
        water(il,i) = 0.
       endif
@@ -2345 +2350 @@
       real esum(nloc), fsum(nloc), gsum(nloc), hsum(nloc)
       real th_wake(nloc,nd)
-      real alpha_qpos(nloc)
+      real alpha_qpos(nloc),alpha_qpos1(nloc)
       real qcond(nloc,nd), nqcond(nloc,nd), wa(nloc,nd)  ! cld
       real siga(nloc,nd), sax(nloc,nd), mac(nloc,nd)      ! cld
@@ -3050 +3055 @@
        do il=1,ncum
         IF (iflag(il) .le. 1) THEN
+        IF (cvflag_grav) then
+        ex=0.01*grav*ment(il,inb(il),inb(il))
+     :      *(traent(il,inb(il),inb(il),j)-tra(il,inb(il),j))
+     :      /(ph(i l,inb(il))-ph(il,inb(il)+1))
+        ftra(il,inb(il),j)=ftra(il,inb(il),j)-ex
+        ftra(il,inb(il)-1,j)=ftra(il,inb(il)-1,j)
+     :       +ex*(ph(il,inb(il))-ph(il,inb(il)+1))
+     :          /(ph(il,inb(il)-1)-ph(il,inb(il)))
+        else
         ex=0.1*ment(il,inb(il),inb(il))
      :      *(traent(il,inb(il),inb(il),j)-tra(il,inb(il),j))
@@ -3057 +3071 @@
      :       +ex*(ph(il,inb(il))-ph(il,inb(il)+1))
      :          /(ph(il,inb(il)-1)-ph(il,inb(il)))
+        ENDIF   !cvflag grav
         ENDIF    !iflag
        enddo
@@ -3122 +3137 @@
 c        in order to ensure moisture positivity
       DO il = 1,ncum
+      alpha_qpos(il)=1.
        IF (iflag(il) .le. 1) THEN
-        alpha_qpos(il) = max(1. , -delt*fr(il,1)/
+        if (fr(il,1) .le. 0.) then
+            alpha_qpos(il) = max(alpha_qpos(il) ,
+     :     (-delt*fr(il,1))/
      :     (s_wake(il)*rr_wake(il,1)+(1.-s_wake(il))*rr(il,1)))
+        end if
        ENDIF
       ENDDO
@@ -3130 +3149 @@
        DO il = 1,ncum
         IF (iflag(il) .le. 1) THEN
-        alpha_qpos(il) = max(alpha_qpos(il) , -delt*fr(il,i)/
+          IF (fr(il,i) .le. 0.) THEN
+           alpha_qpos1(il)=max(1. , (-delt*fr(il,i))/
      :     (s_wake(il)*rr_wake(il,i)+(1.-s_wake(il))*rr(il,i)))
+             IF (alpha_qpos1(il) .ge. alpha_qpos(il))
+     :           alpha_qpos(il)=alpha_qpos1(il)
+          ENDIF
         ENDIF
        ENDDO