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

Timestamp:

Nov 21, 2017, 5:18:23 PM (7 years ago)

Author:

oboucher

Message:

Adding new diagnostics for clean (no aerosol) clear sky SW radiative fluxes

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -43 +43 @@
          o_snow, o_msnow, o_fsnow, o_evap, o_ep,o_epmax_diag, & ! epmax_cape
          o_tops, o_tops0, o_topl, o_topl0, &
-         o_SWupTOA, o_SWupTOAclr, o_SWdnTOA, &
+         o_SWupTOA, o_SWupTOAclr, o_SWupTOAcleanclr, o_SWdnTOA, &
          o_SWdnTOAclr, o_nettop, o_SWup200, &
          o_SWup200clr, o_SWdn200, o_SWdn200clr, &
@@ -49 +49 @@
          o_LWdn200clr, o_sols, o_sols0, &
          o_soll, o_radsol, o_soll0, o_SWupSFC, &
-         o_SWupSFCclr, o_SWdnSFC, o_SWdnSFCclr, &
+         o_SWupSFCclr, o_SWupSFCcleanclr, o_SWdnSFC, o_SWdnSFCclr, o_SWdnSFCcleanclr, &
          o_LWupSFC, o_LWdnSFC, o_LWupSFCclr, &
          o_LWdnSFCclr, o_bils, o_bils_diss, &
@@ -154 +154 @@
          o_du_gwd_front, o_dv_gwd_front, &
          o_east_gwstress, o_west_gwstress, &
-         o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, &
+         o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, o_rsucsaf, o_rsdcsaf, &
          o_rlucs, o_rldcs, o_tnt, o_tntr, &
          o_tntscpbl, o_tnhus, o_tnhusscpbl, &
@@ -210 +210 @@
          nday_rain, rain_con, snow_con, &
          topsw, toplw, toplw0, swup, swdn, &
-         topsw0, swup0, swdn0, SWup200, SWup200clr, &
+         topsw0, swupc0, swdnc0, swup0, swdn0, SWup200, SWup200clr, &
          SWdn200, SWdn200clr, LWup200, LWup200clr, &
          LWdn200, LWdn200clr, solsw, solsw0, sollw, &
@@ -225 +225 @@
          ratqs, entr_therm, zqasc, detr_therm, f0, &
          lwup, lwdn, lwup0, coefm, &
-         swupp, lwupp, swup0p, lwup0p, swdnp, lwdnp, &
-         swdn0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
+         swupp, lwupp, swupc0p, swup0p, lwup0p, swdnp, lwdnp, &
+         swdnc0p, swdn0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
          vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &
          vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &
@@ -708 +708 @@
 
        IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swupc0(:,klevp1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWupTOAcleanclr, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
           zx_tmp_fi2d(:) = swdn(:,klevp1)*swradcorr(:)
        ENDIF
@@ -746 +751 @@
 
        IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swupc0(:,1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWupSFCcleanclr, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
           zx_tmp_fi2d(:) = swdn(:,1)*swradcorr(:)
        ENDIF
@@ -754 +764 @@
        ENDIF
        CALL histwrite_phy(o_SWdnSFCclr, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swdnc0(:,1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWdnSFCcleanclr, zx_tmp_fi2d)
 
        IF (vars_defined) THEN
@@ -1738 +1753 @@
        CALL histwrite_phy(o_rsucs, zx_tmp_fi3d1)
        DO k=1, klevp1
+         zx_tmp_fi3d1(:,k)=swupc0(:,k)*swradcorr(:)
+       ENDDO
+       CALL histwrite_phy(o_rsucsaf, zx_tmp_fi3d1)
+       DO k=1, klevp1
          zx_tmp_fi3d1(:,k)=swdn0(:,k)*swradcorr(:)
        ENDDO
        CALL histwrite_phy(o_rsdcs, zx_tmp_fi3d1)
+       DO k=1, klevp1
+         zx_tmp_fi3d1(:,k)=swdnc0(:,k)*swradcorr(:)
+       ENDDO
+       CALL histwrite_phy(o_rsdcsaf, zx_tmp_fi3d1)
 
        CALL histwrite_phy(o_rlu, lwup)