Changeset 3071 for LMDZ6/trunk/libf/phylmd/ener_conserv.F90

Timestamp:

Nov 14, 2017, 4:17:16 PM (7 years ago)

Author:

jyg

Message:

Correction in ener_conserv.F90: bils_tke is
set to zero if iflag_pbl <= 1.

Cosmetic correction in physiq_mod.F90: pbl_tke is
initialized to zero if iflag_pbl <= 1 and
klon_glo=1.

File:

• LMDZ6/trunk/libf/phylmd/ener_conserv.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/ener_conserv.F90

@@ -222 +222 @@
       DO k=1,klev
         bils_ec(:)=bils_ec(:)-d_t_ec(:,k)*masse(:,k)
-        bils_tke(:)=bils_tke(:)+0.5*(dtke(:,k)+dtke(:,k+1))*masse(:,k)
         bils_diss(:)=bils_diss(:)-d_t_diss(:,k)*masse(:,k)
         bils_kinetic(:)=bils_kinetic(:)+masse(:,k)* &
@@ -235 +234 @@
       ENDDO
       bils_ec(:)=rcpd*bils_ec(:)/pdtphys
-      bils_tke(:)=bils_tke(:)/pdtphys
       bils_diss(:)=rcpd*bils_diss(:)/pdtphys
       bils_kinetic(:)= 0.5*bils_kinetic(:)/pdtphys
       bils_enthalp(:)=rcpd*bils_enthalp(:)/pdtphys
       bils_latent(:)=rlvtt*bils_latent(:)/pdtphys
+!jyg<
+      IF (iflag_pbl > 1) THEN
+        DO k=1,klev
+          bils_tke(:)=bils_tke(:)+0.5*(dtke(:,k)+dtke(:,k+1))*masse(:,k)
+        ENDDO
+        bils_tke(:)=bils_tke(:)/pdtphys
+      ENDIF  ! (iflag_pbl > 1)
+!>jyg
 
 IF (iflag_ener_conserv>=1) THEN