Changeset 4758 for LMDZ6/branches/LMDZ_ECRad/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Dec 1, 2023, 10:09:29 PM (6 months ago)

Author:

idelkadi

Message:

• Update of Ecrad in LMDZ version 1.6.1 ​https://github.com/lguez/ecrad.git
• Continuation of work on the implementation of the double call of Ecrad in LMDZ

File:

• LMDZ6/branches/LMDZ_ECRad/libf/phylmd/phys_output_write_mod.F90 (modified) (4 diffs)

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/phys_output_write_mod.F90

@@ -231 +231 @@
     USE phys_output_ctrlout_mod, ONLY:  &
          o_soll0_s2,o_soll_s2,o_sols0_s2,o_sols_s2, &
-         o_topl0_s2,o_topl_s2,o_tops0_s2,o_tops_s2   
+         o_topl0_s2,o_topl_s2,o_tops0_s2,o_tops_s2, &
+         o_SWupTOA_s2,o_SWupTOAclr_s2,o_cloud_cover_sw, &
+         o_cloud_cover_sw_s2,o_SWdnTOA_s2,o_SWdnTOAclr_s2, &
+         o_LWupSFCclr_s2, o_LWdnSFCclr_s2, o_SWupSFC_s2, &
+         o_SWupSFCclr_s2, o_SWdnSFC_s2, o_SWdnSFCclr_s2, &
+         o_LWupSFC_s2, o_LWdnSFC_s2
 #endif
 
@@ -290 +295 @@
     USE phys_state_var_mod, ONLY: &
         sollw0_s2,sollw_s2,solsw0_s2,solsw_s2, &
-        toplw0_s2,toplw_s2,topsw0_s2,topsw_s2
+        toplw0_s2,toplw_s2,topsw0_s2,topsw_s2, &
+        swup0_s2,swup_s2,swdn_s2,swdn0_s2,sollwdownclr_s2, &
+        sollwdown_s2,lwdn0_s2
+    USE phys_output_var_mod, ONLY: cloud_cover_sw, &
+        cloud_cover_sw_s2
 #endif
+
 
 
@@ -1068 +1078 @@
 !AI 08 2023 Ecrad 3Deffect
 #ifdef CPP_ECRAD
+     CALL histwrite_phy(o_cloud_cover_sw, cloud_cover_sw)
      if (ok_3Deffect) then
         IF (vars_defined) THEN
@@ -1091 +1102 @@
        CALL histwrite_phy(o_topl_s2, toplw_s2)
        CALL histwrite_phy(o_topl0_s2, toplw0_s2)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swup_s2(:,klevp1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWupTOA_s2, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swup0_s2(:,klevp1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWupTOAclr_s2, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swdn_s2(:,klevp1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWdnTOA_s2, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swdn0_s2(:,klevp1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWdnTOAclr_s2, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:)=sollwdown_s2(:)-sollw_s2(:)
+       ENDIF
+       CALL histwrite_phy(o_LWupSFC_s2, zx_tmp_fi2d)
+       CALL histwrite_phy(o_LWdnSFC_s2, sollwdown_s2)
+
+       IF (vars_defined) THEN
+          sollwdownclr_s2(1:klon) = -1.*lwdn0_s2(1:klon,1)
+          zx_tmp_fi2d(1:klon)=sollwdownclr_s2(1:klon)-sollw0_s2(1:klon)
+       ENDIF
+       CALL histwrite_phy(o_LWupSFCclr_s2, zx_tmp_fi2d)
+       CALL histwrite_phy(o_LWdnSFCclr_s2, sollwdownclr_s2)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swup_s2(:,1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWupSFC_s2, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+          zx_tmp_fi2d(:) = swup0_s2(:,1)*swradcorr(:)
+       ENDIF
+       CALL histwrite_phy(o_SWupSFCclr_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)
+
      endif
 #endif