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

Timestamp:

Oct 16, 2016, 11:17:32 AM (8 years ago)

Author:

oboucher

Message:

adding new diagnostics for evap/rain/snow weighted by surface types
this parallels what was done for wbilo (net flux)
new diagnostics are wevap_??? wrain_??? wsnow_???
where ??? is ter, oce, lic or sic.

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -59 +59 @@
          o_sens_srf, o_lat_srf, o_flw_srf, &
          o_fsw_srf, o_wbils_srf, o_wbilo_srf, &
+         o_wevap_srf, o_wrain_srf, o_wsnow_srf, &
          o_tke_srf, o_tke_max_srf,o_dltpbltke_srf, o_wstar, &
          o_l_mixmin,o_l_mix, &
@@ -218 +219 @@
          sissnow, runoff, albsol3_lic, evap_pot, &
          t2m, fluxt, fluxlat, fsollw, fsolsw, &
-         wfbils, wfbilo, cdragm, cdragh, cldl, cldm, &
+         wfbils, wfbilo, wfevap, wfrain, wfsnow, & 
+         cdragm, cdragh, cldl, cldm, &
          cldh, cldt, JrNt, cldljn, cldmjn, cldhjn, &
          cldtjn, cldq, flwp, fiwp, ue, ve, uq, vq, &
@@ -679 +681 @@
 
        DO nsrf = 1, nbsrf
+
           IF (vars_defined)             zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
           CALL histwrite_phy(o_pourc_srf(nsrf), zx_tmp_fi2d)
@@ -713 +716 @@
           IF (vars_defined)         zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
           CALL histwrite_phy(o_wbilo_srf(nsrf), zx_tmp_fi2d)
+          IF (vars_defined)         zx_tmp_fi2d(1 : klon) = wfevap( 1 : klon, nsrf)
+          CALL histwrite_phy(o_wevap_srf(nsrf), zx_tmp_fi2d)
+          IF (vars_defined)         zx_tmp_fi2d(1 : klon) = wfrain( 1 : klon, nsrf)
+          CALL histwrite_phy(o_wrain_srf(nsrf), zx_tmp_fi2d)
+          IF (vars_defined)         zx_tmp_fi2d(1 : klon) = wfsnow( 1 : klon, nsrf)
+          CALL histwrite_phy(o_wsnow_srf(nsrf), zx_tmp_fi2d)
 
           IF (iflag_pbl > 1) THEN