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

Timestamp:

Oct 27, 2017, 6:14:12 PM (7 years ago)

Author:

musat

Message:

Set up lwoff configuration : add ok_lwoff=y to activate it.
Add solbnd CMIP6 output.
Enable 3D outputs of SW/LW total sky and clear sky.
IM

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -119 +119 @@
          o_swtoacf_zero, o_swsrfcf_zero, &
          o_topswai, o_solswai, o_toplwai, o_sollwai, o_scdnc, &
-         o_cldncl, o_reffclws, o_reffclwc, o_stratomask,&
+         o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,&
          o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
          o_lcc3dstra, o_reffclwtop, o_ec550aer, &
@@ -389 +389 @@
     REAL, DIMENSION (klon,klev) :: zx_tmp_fi3d, zpt_conv
     REAL, DIMENSION (klon,klev+1) :: zx_tmp_fi3d1
+    REAL, DIMENSION (klon,NSW) :: zx_tmp_fi3dsp
     CHARACTER (LEN=4)              :: bb2
     INTEGER, DIMENSION(nbp_lon*nbp_lat)  :: ndex2d
@@ -1383 +1384 @@
        CALL histwrite_phy(o_zfull,zphi/RG)
 
+!solbnd begin
+#ifdef CPP_RRTM
+      IF (iflag_rrtm.EQ.1) THEN
+       IF (vars_defined) THEN
+        DO ISW=1, NSW
+          zx_tmp_fi3dsp(:,ISW) = swdn(:,klevp1)*swradcorr(:)*RSUN(ISW)
+        ENDDO
+        CALL histwrite_phy(o_solbnd, zx_tmp_fi3dsp)
+       ENDIF
+      ENDIF
+#endif
+
+!solbnd end
+
        IF (flag_aerosol_strat.EQ.2) THEN
          CALL histwrite_phy(o_stratomask, stratomask)