Changeset 2903 for LMDZ5/trunk/libf/phylmd

Timestamp:

Jun 6, 2017, 12:24:43 PM (8 years ago)

Author:

fhourdin

Message:

Correction de Jean-Yves sur la conservation oubliée au cours
d'une commission ultérieure de Frédéric.
Permettant de garantir une continuité sur l'ancienne physique.

File:

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

LMDZ5/trunk/libf/phylmd/ener_conserv.F90

@@ -27 +27 @@
 USE phys_state_var_mod, ONLY : du_gwd_front,du_gwd_rando
 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
@@ -62 +64 @@
    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