Changeset 2284 for LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

May 19, 2015, 6:33:30 PM (9 years ago)

Author:

jyg

Message:

Bug fixing concerning aerosol scavenging (bugs in
cvltr_scav and in lsc_scav) and aerosol
concentration output (the values at beginning,
instead of end, of time step were output).

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -215 +215 @@
          lcc, lcc3d, lcc3dcon, lcc3dstra, reffclwtop, &
          ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
-         ql_seri, zphi, u_seri, v_seri, omega, cldfra, &
+!jyg<
+!!         ql_seri, zphi, u_seri, v_seri, omega, cldfra, &
+         ql_seri, tr_seri, &
+         zphi, u_seri, v_seri, omega, cldfra, &
+!>jyg
          rneb, rnebjn, zx_rh, d_t_dyn, d_q_dyn, &
          d_u_dyn, d_v_dyn, d_t_con, d_t_ajsb, d_t_ajs, &
@@ -1452 +1456 @@
               IF (type_trac == 'lmdz' .OR. type_trac == 'repr') THEN
 
-             CALL histwrite_phy(o_trac(iq-nqo), qx(:,:,iq))
+!jyg<
+!!             CALL histwrite_phy(o_trac(iq-nqo), qx(:,:,iq))
+             CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
+!>jyg
              CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
              CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
@@ -1469 +1476 @@
              IF(vars_defined) THEN
                 DO k=1,klev
-                   zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*qx(:,k,iq)
+!jyg<
+!!                   zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*qx(:,k,iq)
+                   zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*tr_seri(:,k,iq-nqo)
+!>jyg
                 ENDDO
              ENDIF