Changeset 3106 for LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Dec 3, 2017, 10:03:26 PM (7 years ago)

Author:

oboucher

Message:

Adding clear-sky clean (no aerosol) LW radiative fluxes at TOA and SFC as diagnostic

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -51 +51 @@
          o_SWupSFCclr, o_SWupSFCcleanclr, o_SWdnSFC, o_SWdnSFCclr, o_SWdnSFCcleanclr, &
          o_LWupSFC, o_LWdnSFC, o_LWupSFCclr, &
-         o_LWdnSFCclr, o_bils, o_bils_diss, &
+         o_LWdnSFCclr, o_LWupTOAcleanclr, o_LWdnSFCcleanclr, o_bils, o_bils_diss, &
          o_bils_ec,o_bils_ech, o_bils_tke, o_bils_kinetic, &
          o_bils_latent, o_bils_enthalp, o_sens, &
@@ -216 +216 @@
          LWdn200, LWdn200clr, solsw, solsw0, sollw, &
          radsol, swradcorr, sollw0, sollwdown, sollw, gustiness, &
-         sollwdownclr, lwdn0, ftsol, ustar, u10m, &
+         sollwdownclr, lwdnc0, lwdn0, ftsol, ustar, u10m, &
          v10m, pbl_tke, wake_delta_pbl_TKE, &
          delta_tsurf, &
@@ -226 +226 @@
          rnebcon, wo, falb1, albsol2, coefh, clwcon0, &
          ratqs, entr_therm, zqasc, detr_therm, f0, &
-         lwup, lwdn, lwup0, coefm, &
-         swupp, lwupp, swupc0p, swup0p, lwup0p, swdnp, lwdnp, &
-         swdnc0p, swdn0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
+         lwup, lwdn, lwupc0, lwup0, coefm, &
+         swupp, lwupp, swupc0p, swup0p, lwupc0p, lwup0p, swdnp, lwdnp, &
+         swdnc0p, swdn0p, lwdnc0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
          vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &
          vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &
@@ -787 +787 @@
        CALL histwrite_phy(o_LWupSFCclr, zx_tmp_fi2d)
        CALL histwrite_phy(o_LWdnSFCclr, sollwdownclr)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = lwupc0(:,klevp1)
+       ENDIF
+       CALL histwrite_phy(o_LWupTOAcleanclr, zx_tmp_fi2d)
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = -1.*lwdnc0(:,1)
+       ENDIF
+       CALL histwrite_phy(o_LWdnSFCcleanclr, zx_tmp_fi2d)
+
        CALL histwrite_phy(o_bils, bils)
        CALL histwrite_phy(o_bils_diss, bils_diss)