Changeset 2799 for LMDZ5/trunk/libf/phylmd/reevap.F90

Timestamp:

Feb 24, 2017, 7:50:33 PM (7 years ago)

Author:

jyg

Message:

New diagnostics to verify energy conservation.
Some corrections to improve energy conservation.
Some bug correction.

New output variables and new flags are

introduced:

(1)New variables: d_xx_col is the budget over

each column of variable "xx".

(2) fl_ebil: 0 -> no diag; 1 -> diags activated
(3) fl_cor_ebil: 0 -> no correction ensuring

energy conservation; 1 -> first set of
corrections.

+ ok_bad_ecmwf_thermo: if true, a bug setting

RVTMP2 to 0 is fixed; for backward compatibility
the default is "false".

File:

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

LMDZ5/trunk/libf/phylmd/reevap.F90

@@ -2 +2 @@
    &         d_t_eva,d_q_eva,d_ql_eva,d_qs_eva)
 
-
+    ! flag to include modifications to ensure energy conservation (if flag >0)
+    USE add_phys_tend_mod, only : fl_cor_ebil 
+    
     IMPLICIT none
     !>======================================================================
@@ -25 +27 @@
     DO k = 1, klev  ! re-evaporation de l'eau liquide nuageuse
        DO i = 1, klon
+        if (fl_cor_ebil .GT. 0) then
+          zlvdcp=RLVTT/RCPD/(1.0+RVTMP2*(q_seri(i,k)+ql_seri(i,k)+qs_seri(i,k)))
+          zlsdcp=RLSTT/RCPD/(1.0+RVTMP2*(q_seri(i,k)+ql_seri(i,k)+qs_seri(i,k)))
+        else
           zlvdcp=RLVTT/RCPD/(1.0+RVTMP2*q_seri(i,k))
           !jyg<
@@ -30 +36 @@
           !                  A verifier !!!
           zlsdcp=RLSTT/RCPD/(1.0+RVTMP2*q_seri(i,k))
+        end if
           IF (iflag_ice_thermo .EQ. 0) THEN
              zlsdcp=zlvdcp