Changeset 4808 for LMDZ6/trunk/libf/phylmd

Timestamp:

Feb 8, 2024, 4:25:36 PM (10 months ago)

Author:

idelkadi

Message:

Following implementation of the double call of Ecrad:

• addition of output variables in the code
• addition of output variables in XML file for XIOS

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -238 +238 @@
          o_LWupSFCclr_s2, o_LWdnSFCclr_s2, o_SWupSFC_s2, &
          o_SWupSFCclr_s2, o_SWdnSFC_s2, o_SWdnSFCclr_s2, &
-         o_LWupSFC_s2, o_LWdnSFC_s2   
+         o_LWupSFC_s2, o_LWdnSFC_s2, , o_rlu_s2, o_rld_s2, &
+         o_rlucs_s2, o_rldcs_s2, o_rsu_s2, o_rsd_s2, &
+         o_rsucs_s2, o_rsdcs_s2   
 #endif
 
@@ -301 +303 @@
         toplw0_s2,toplw_s2,topsw0_s2,topsw_s2, &
         swup0_s2,swup_s2,swdn_s2,swdn0_s2,sollwdownclr_s2, &
-        sollwdown_s2,lwdn0_s2
+        sollwdown_s2,lwdn0_s2,lwup_s2,lwdn_s2,lwup0_s2,lwdn0_s2
     USE phys_output_var_mod, ONLY: cloud_cover_sw, &
         cloud_cover_sw_s2
@@ -1158 +1160 @@
        ENDIF
        CALL histwrite_phy(o_SWdnSFCclr_s2, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swdn_s2(:,1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWdnSFC_s2, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swdn0_s2(:,1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWdnSFCclr_s2, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+         DO k=1, klevp1
+           zx_tmp_fi3d1(:,k)=swup_s2(:,k)*swradcorr(:)
+         ENDDO
+       ENDIF
+       CALL histwrite_phy(o_rsu_s2, zx_tmp_fi3d1)
+
+       IF (vars_defined) THEN
+         DO k=1, klevp1
+           zx_tmp_fi3d1(:,k)=swdn_s2(:,k)*swradcorr(:)
+         ENDDO
+       ENDIF
+       CALL histwrite_phy(o_rsd_s2, zx_tmp_fi3d1)
+
+       IF (vars_defined) THEN
+         DO k=1, klevp1
+           zx_tmp_fi3d1(:,k)=swup0_s2(:,k)*swradcorr(:)
+         ENDDO
+       ENDIF
+       CALL histwrite_phy(o_rsucs_s2, zx_tmp_fi3d1)
+
+       IF (vars_defined) THEN
+         DO k=1, klevp1
+           zx_tmp_fi3d1(:,k)=swdn0_s2(:,k)*swradcorr(:)
+         ENDDO
+       ENDIF
+       CALL histwrite_phy(o_rsdcs_s2, zx_tmp_fi3d1)
+
+       CALL histwrite_phy(o_rlu_s2, lwup_s2)
+       CALL histwrite_phy(o_rld_s2, lwdn_s2)
+       CALL histwrite_phy(o_rlucs_s2, lwup0_s2)
+       CALL histwrite_phy(o_rldcs_s2, lwdn0_s2)
     ENDIF !ok_3Deffect
 #endif