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

Timestamp:

Sep 7, 2023, 1:07:27 PM (9 months ago)

Author:

idelkadi

Message:

Implementation in the LMDZ code of the double call of the ECRAD radiative transfer code to estimate the 3D radiative effect of clouds.

• This double call of Ecrad is controlled by the ok_3Deffect logic key.
• If this key is enabled, 2 files of parameter configuration "namelists" for ECRAD are required at runtime: namelist_ecrad and namelist_ecrad_s2.
• If this key is deactivated, the configuration and initialization part (reading namelist and netcdf files) is performed only once during simulation (1st call to ECRAD). Otherwise, configuration and initialization are performed each time Ecrad is called.

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -228 +228 @@
          o_zxfluxt,o_zxfluxq
 
+#ifdef CPP_ECRAD
+    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   
+#endif
+
 #ifdef CPP_StratAer
     USE phys_output_ctrlout_mod, ONLY:  & 
@@ -279 +285 @@
 #endif
          dt_ns, delta_sst, dter, dser
+         
+! AI 08 2023 pour ECRAD 3Deffect
+#ifdef CPP_ECRAD
+    USE phys_state_var_mod, ONLY: &
+        sollw0_s2,sollw_s2,solsw0_s2,solsw_s2, &
+        toplw0_s2,toplw_s2,topsw0_s2,topsw_s2
+#endif
+
 
     USE phys_local_var_mod, ONLY: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, &
@@ -1051 +1065 @@
        ENDIF
 
+!AI 08 2023 Ecrad 3Deffect
+#ifdef CPP_ECRAD
+     if (ok_3Deffect) then
+        IF (vars_defined) THEN
+          zx_tmp_fi2d = solsw_s2*swradcorr
+       ENDIF
+       CALL histwrite_phy(o_sols_s2, zx_tmp_fi2d)
+       IF (vars_defined) THEN
+          zx_tmp_fi2d = solsw0_s2*swradcorr
+       ENDIF
+       CALL histwrite_phy(o_sols0_s2, zx_tmp_fi2d)
+       CALL histwrite_phy(o_soll_s2, sollw_s2)
+       CALL histwrite_phy(o_soll0_s2, sollw0_s2)
+       IF (vars_defined) THEN
+         zx_tmp_fi2d = topsw_s2*swradcorr
+       ENDIF
+       CALL histwrite_phy(o_tops_s2, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+         zx_tmp_fi2d = topsw0_s2*swradcorr
+       ENDIF
+       CALL histwrite_phy(o_tops0_s2, zx_tmp_fi2d)
+
+       CALL histwrite_phy(o_topl_s2, toplw_s2)
+       CALL histwrite_phy(o_topl0_s2, toplw0_s2)
+     endif
+#endif       
+
        CALL histwrite_phy(o_bils, bils)
        CALL histwrite_phy(o_bils_diss, bils_diss)