Changeset 4647 for LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_scheme.F90

Timestamp:

Aug 7, 2023, 1:42:51 PM (14 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/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_scheme.F90 (modified) (8 diffs)

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_scheme.F90

@@ -9 +9 @@
 !             3. Configuration a partir de namelist
 !             4. frac_std = 0.75      
+! Juillet 2023 :
+!             
 ! ============================================================================
 
@@ -14 +16 @@
 ! Inputs
      & (KIDIA, KFDIA, KLON, KLEV, KAEROSOL, NSW, &
-     &  IDAY, TIME, &
+     &  namelist_file, ok_3Deffect, IDAY, TIME, &
      &  PSOLAR_IRRADIANCE, &
      &  PMU0, PTEMPERATURE_SKIN, &
@@ -97 +99 @@
 USE mod_phys_lmdz_para
 USE setup_config_from_lmdz,   ONLY : driver_config_type
+!USE RADIATION_SETUP, ONLY : config_type, driver_config_type
 
 IMPLICIT NONE
@@ -220 +223 @@
 ! AI ATTENTION
 type(config_type),save         :: rad_config
+!!$OMP THREADPRIVATE(rad_config)
 type(driver_config_type),save  :: driver_config
+!!$OMP THREADPRIVATE(driver_config)
+!type(config_type)        :: rad_config
+!type(driver_config_type)  :: driver_config
 TYPE(single_level_type)   :: single_level
 TYPE(thermodynamics_type) :: thermodynamics
@@ -261 +268 @@
 
 ! Name of file names specified on command line
-character(len=512) :: file_name
+character(len=512) :: namelist_file
 
 logical :: loutput=.true.
@@ -269 +276 @@
 !$OMP THREADPRIVATE(debut_ecrad)
 integer, save :: itap_ecrad=1
+logical :: ok_3Deffect
 
 IF (LHOOK) CALL DR_HOOK('RADIATION_SCHEME',0,ZHOOK_HANDLE)
-print*,'Entree dans radiation_scheme'
 
 ! A.I juillet 2023 : 
 ! Initialisation dans radiation_setup au 1er passage dans Ecrad
 !$OMP MASTER
-if (debut_ecrad) then
-! AI appel radiation_setup
- call SETUP_RADIATION_SCHEME(loutput,rad_config,driver_config)
- debut_ecrad=.false. 
-endif 
+if (.not.ok_3Deffect) then
+  if (debut_ecrad) then
+   call SETUP_RADIATION_SCHEME(loutput,namelist_file,rad_config,driver_config)
+   debut_ecrad=.false. 
+  endif 
+else
+   call SETUP_RADIATION_SCHEME(loutput,namelist_file,rad_config,driver_config)
+endif
 !$OMP END MASTER
 !$OMP BARRIER 
@@ -287 +297 @@
 !AI juillet 2023 : verif des param de config :
 if (lprint_config) then
+ IF (is_master) THEN       
    print*,'Parametres de configuration de ecrad, etape ',itap_ecrad     
+   print*,'Entree dans radiation_scheme'
+   print*,'ok_3Deffect = ',ok_3Deffect
+   print*,'Fichier namelist = ',namelist_file
+
    print*,'do_sw, do_lw, do_sw_direct, do_3d_effects = ', &
            rad_config%do_sw, rad_config%do_lw, rad_config%do_sw_direct, rad_config%do_3d_effects
@@ -362 +377 @@
 
            itap_ecrad=itap_ecrad+1
+   ENDIF          
 endif