Changeset 2881 for LMDZ5/trunk/libf

Timestamp:

May 15, 2017, 7:39:40 PM (7 years ago)

Author:

jyg

Message:

Added test on fl_cor_ebil in ener_conserv.F90 for compatibility with older versions

File:

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

LMDZ5/trunk/libf/phylmd/ener_conserv.F90

@@ -24 +24 @@
 USE phys_local_var_mod, ONLY : d_t_eva,d_t_lsc,d_q_eva,d_q_lsc
 USE phys_output_var_mod, ONLY : bils_ec,bils_ech,bils_tke,bils_kinetic,bils_enthalp,bils_latent,bils_diss
+USE add_phys_tend_mod, ONLY : fl_cor_ebil 
 
 IMPLICIT none
@@ -59 +60 @@
    DO k = 1, klev
    DO i = 1, klon
-      ZRCPD = RCPD*(1.0+RVTMP2*(pqn(i,k)+pqln(i,k)+pqsn(i,k)))
-      d_t_ec(i,k)=0.5/ZRCPD &
- &      *(puo(i,k)**2+pvo(i,k)**2-pun(i,k)**2-pvn(i,k)**2)
-      ENDDO
-      ENDDO
+     IF (fl_cor_ebil .GT. 0) then
+       ZRCPD = RCPD*(1.0+RVTMP2*(pqn(i,k)+pqln(i,k)+pqsn(i,k)))
+     ELSE
+       ZRCPD = RCPD*(1.0+RVTMP2*pqn(i,k))
+     ENDIF
+     d_t_ec(i,k)=0.5/ZRCPD &
+ &     *(puo(i,k)**2+pvo(i,k)**2-pun(i,k)**2-pvn(i,k)**2)
+   ENDDO
+   ENDDO
 !-jld ec_conser