Changeset 2908 for LMDZ5/trunk/libf

Timestamp:

Jun 12, 2017, 11:53:24 AM (8 years ago)

Author:

jyg

Message:

Corrections for energy conservation in the
convective scheme (activated by fl_cor_ebil>=2):
the change of Cp of air with moisture is now
better accounted for in the tendency due to
vertical advection.

File:

• LMDZ5/trunk/libf/phylmd/cv3_routines.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/cv3_routines.F90

@@ -3017 +3017 @@
 
     USE print_control_mod, ONLY: lunout, prt_level
+    USE add_phys_tend_mod, only : fl_cor_ebil
 
   IMPLICIT NONE
@@ -3253 +3254 @@
 
       IF ((0.01*grav*work(il)*am(il))>=delti) iflag(il) = 1 !consist vect
-      ft(il, 1) = ft(il, 1) + 0.01*grav*work(il)*am(il) * &
-                                   (t(il,2)-t(il,1)+(gz(il,2)-gz(il,1))/cpn(il,1))
+!jyg<
+        IF (fl_cor_ebil >= 2) THEN
+          ft(il, 1) = ft(il, 1) + 0.01*grav*work(il)*am(il) * &
+                    ((t(il,2)-t(il,1))*cpn(il,2)+gz(il,2)-gz(il,1))/cpn(il,1)
+        ELSE
+          ft(il, 1) = ft(il, 1) + 0.01*grav*work(il)*am(il) * &
+                    (t(il,2)-t(il,1)+(gz(il,2)-gz(il,1))/cpn(il,1))
+        ENDIF
+!>jyg
     END IF ! iflag
   END DO
@@ -3562 +3570 @@
 
 ! sature
-        ft(il, i) = ft(il, i) + 0.01*grav*dpinv * &
+!jyg<
+        IF (fl_cor_ebil >= 2) THEN
+          ft(il, i) = ft(il, i) + 0.01*grav*dpinv * &
+              ( amp1(il)*( (t(il,i+1)-t(il,i))*cpn(il,i+1) + gz(il,i+1)-gz(il,i))*cpinv - &
+                ad(il)*( (t(il,i)-t(il,i-1))*cpn(il,i-1) + gz(il,i)-gz(il,i-1))*cpinv)
+        ELSE
+          ft(il, i) = ft(il, i) + 0.01*grav*dpinv * &
                      (amp1(il)*(t(il,i+1)-t(il,i) + (gz(il,i+1)-gz(il,i))*cpinv) - &
                       ad(il)*(t(il,i)-t(il,i-1)+(gz(il,i)-gz(il,i-1))*cpinv))
+        ENDIF
+!>jyg