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

Timestamp:

Mar 30, 2009, 3:59:31 PM (16 years ago)

Author:

idelkadi

Message:

Correction bug (gravite=10) dans le schema de convection : cv3_routines
Utilisation du meme seuil dans calltherm et wake pour la vapeur d eau

File:

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

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

@@ -2968 +2968 @@
       do 503 il=1,ncum
       IF (iflag(il) .le. 1) THEN
+        if (cvflag_grav) then
+      ax=0.01*grav*ment(il,inb(il),inb(il))*(hp(il,inb(il))
+     : -h(il,inb(il))+t(il,inb(il))*(cpv-cpd)
+     : *(rr(il,inb(il))-qent(il,inb(il),inb(il))))
+     :  /(cpn(il,inb(il))*(ph(il,inb(il))-ph(il,inb(il)+1)))
+      ft(il,inb(il))=ft(il,inb(il))-ax
+      ft(il,inb(il)-1)=ft(il,inb(il)-1)+ax*cpn(il,inb(il))
+     :    *(ph(il,inb(il))-ph(il,inb(il)+1))/(cpn(il,inb(il)-1)
+     :    *(ph(il,inb(il)-1)-ph(il,inb(il))))
+
+      bx=0.01*grav*ment(il,inb(il),inb(il))*(qent(il,inb(il),inb(il))
+     :    -rr(il,inb(il)))/(ph(il,inb(il))-ph(il,inb(il)+1))
+      fr(il,inb(il))=fr(il,inb(il))-bx
+      fr(il,inb(il)-1)=fr(il,inb(il)-1)
+     :   +bx*(ph(il,inb(il))-ph(il,inb(il)+1))
+     :      /(ph(il,inb(il)-1)-ph(il,inb(il)))
+
+      cx=0.01*grav*ment(il,inb(il),inb(il))*(uent(il,inb(il),inb(il))
+     :       -u(il,inb(il)))/(ph(il,inb(il))-ph(il,inb(il)+1))
+      fu(il,inb(il))=fu(il,inb(il))-cx
+      fu(il,inb(il)-1)=fu(il,inb(il)-1)
+     :     +cx*(ph(il,inb(il))-ph(il,inb(il)+1))
+     :        /(ph(il,inb(il)-1)-ph(il,inb(il)))
+
+      dx=0.01*grav*ment(il,inb(il),inb(il))*(vent(il,inb(il),inb(il))
+     :      -v(il,inb(il)))/(ph(il,inb(il))-ph(il,inb(il)+1))
+      fv(il,inb(il))=fv(il,inb(il))-dx
+      fv(il,inb(il)-1)=fv(il,inb(il)-1)
+     :    +dx*(ph(il,inb(il))-ph(il,inb(il)+1))
+     :       /(ph(il,inb(il)-1)-ph(il,inb(il)))
+       else
        ax=0.1*ment(il,inb(il),inb(il))*(hp(il,inb(il))
      : -h(il,inb(il))+t(il,inb(il))*(cpv-cpd)
@@ -2997 +3028 @@
      :    +dx*(ph(il,inb(il))-ph(il,inb(il)+1))
      :       /(ph(il,inb(il)-1)-ph(il,inb(il)))
+       endif
       ENDIF    !iflag
 503   continue