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

Timestamp:

Apr 1, 2026, 12:22:43 PM (5 days ago)

Author:

idelkadi

Message:

Correction to the output:
Direct radiation fluxes were output outside the condition for the double call of the Ecrad radiation code

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -1380 +1380 @@
 #ifdef CPP_ECRAD
      CALL histwrite_phy(o_cloud_cover_sw, cloud_cover_sw)
+     CALL histwrite_phy(o_ZFLUX_DIR, ZFLUX_DIR)
+     CALL histwrite_phy(o_ZFLUX_DIR_CLEAR, ZFLUX_DIR_CLEAR)
+     CALL histwrite_phy(o_ZFLUX_DIR_SUN, ZFLUX_DIR_SUN)
+
      IF (ok_3Deffect) THEN
        CALL histwrite_phy(o_cloud_cover_sw_s2, cloud_cover_sw_s2)
@@ -1490 +1494 @@
        CALL histwrite_phy(o_rlucs_s2, lwup0_s2)
        CALL histwrite_phy(o_rldcs_s2, lwdn0_s2)
-       CALL histwrite_phy(o_ZFLUX_DIR, ZFLUX_DIR)
        CALL histwrite_phy(o_ZFLUX_DIR_s2, ZFLUX_DIR_s2)
-       CALL histwrite_phy(o_ZFLUX_DIR_CLEAR, ZFLUX_DIR_CLEAR)
        CALL histwrite_phy(o_ZFLUX_DIR_CLEAR_s2, ZFLUX_DIR_CLEAR_s2)
-       CALL histwrite_phy(o_ZFLUX_DIR_SUN, ZFLUX_DIR_SUN)
        CALL histwrite_phy(o_ZFLUX_DIR_SUN_s2, ZFLUX_DIR_SUN_s2)
     ENDIF !ok_3Deffect