Changeset 4808

Timestamp:

Feb 8, 2024, 4:25:36 PM (3 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

Location:

LMDZ6/trunk

Files:

• DefLists/field_def_lmdz.xml (modified) (2 diffs)
• libf/phylmd/phys_output_write_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/DefLists/field_def_lmdz.xml

@@ -825 +825 @@
         <field id="rsdcs"    long_name="SW CS downward radiation"    unit="W m-2" />
         <field id="rlucs"    long_name="LW CS upward radiation"    unit="W m-2" />
+        <field id="rsucs_s2"    long_name="SW CS upward radiation"    unit="W m-2" />
+        <field id="rsdcs_s2"    long_name="SW CS downward radiation"    unit="W m-2" />
+        <field id="rlucs_s2"    long_name="LW CS upward radiation"    unit="W m-2" />   
         <field id="rldcs4co2"    long_name="Downwelling CS LW 4xCO2 atmosphere"    unit="W/m2" />
         <field id="rlucs4co2"    long_name="Upwelling CS LW 4xCO2 atmosphere"    unit="W/m2" />
@@ -842 +845 @@
         <field id="rlu"    long_name="LW upward radiation"    unit="W m-2" />
         <field id="rld"    long_name="LW downward radiation"    unit="W m-2" />
+        <field id="rldcs_s2" long_name="LW CS downward radiation"    unit="W m-2" />
+        <field id="rsu_s2"    long_name="SW upward radiation"    unit="W m-2" />
+        <field id="rsd_s2"    long_name="SW downward radiation"    unit="W m-2" />
+        <field id="rlu_s2"    long_name="LW upward radiation"    unit="W m-2" />
+        <field id="rld_s2"    long_name="LW downward radiation"    unit="W m-2" />      
         <field id="kz"    long_name="Kz melange"    unit="m2/s" enabled="FALSE" />
         <field id="evu"    long_name="Eddy viscosity coefficient for Momentum Variables"    unit="m2 s-1" />

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