Changeset 2851 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Apr 14, 2017, 11:47:04 AM (8 years ago)

Author:

jyg

Message:

Fixing bug of previous commit

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -4185 +4185 @@
     d_t_ec(:,:)=0.
     forall (k=1: nbp_lev) exner(:, k) = (pplay(:, k)/paprs(:,1))**RKAPPA
-    CALL ener_conserv(klon,klev,pdtphys,u,v,t,qx(:,:,ivap), &
-         u_seri,v_seri,t_seri,q_seri,pbl_tke(:,:,is_ave)-tke0(:,:), &
+    CALL ener_conserv(klon,klev,pdtphys,u,v,t,qx(:,:,ivap),qx(:,:,iliq),qx(:,:,isol), &
+         u_seri,v_seri,t_seri,q_seri,ql_seri,qs_seri,pbl_tke(:,:,is_ave)-tke0(:,:), &
          zmasse,exner,d_t_ec)
     t_seri(:,:)=t_seri(:,:)+d_t_ec(:,:)