Changeset 4368 for LMDZ6/branches/Ocean_skin/libf/phylmd/ecrad/radiation_config.F90

Timestamp:

Dec 6, 2022, 12:01:16 AM (22 months ago)

Author:

lguez

Message:

Sync latest trunk changes to Ocean_skin

Location:

LMDZ6/branches/Ocean_skin

Files:

• . (modified) (1 prop)
• libf/phylmd/ecrad/radiation_config.F90 (modified) (7 diffs)

LMDZ6/branches/Ocean_skin/libf/phylmd/ecrad/radiation_config.F90

@@ -609 +609 @@
 
     integer :: iunit ! Unit number of namelist file
+
+    logical :: lldeb_conf = .false.
 
     namelist /radiation/ do_sw, do_lw, do_sw_direct, &
@@ -844 +846 @@
     this%i_lw_emiss_index              = i_lw_emiss_index
 
+! AI mars 2022
+if (lldeb_conf) then
+print*,'**************PARAMETRES DE CONFIGURATION OFFLINE*******************'
+print*,'config%iverbosesetup   = ', iverbosesetup
+print*,'config%do_lw   = ', do_lw
+print*,'config%do_sw   = ', do_sw
+print*,'config%do_clear   = ', do_clear
+print*,'config%do_sw_direct   = ', do_sw_direct
+print*,'config%do_3d_effects   = ', do_3d_effects
+print*,'config%do_3d_lw_multilayer_effects   = ', do_3d_lw_multilayer_effects
+print*,'config%do_lw_side_emissivity   = ', do_lw_side_emissivity
+print*,'config%use_expm_everywhere   = ', use_expm_everywhere
+print*,'config%use_aerosols   = ', use_aerosols
+print*,'config%do_lw_cloud_scattering   = ', do_lw_cloud_scattering
+print*,'config%do_lw_aerosol_scattering   = ', do_lw_aerosol_scattering
+print*,'config%nregions   = ', n_regions
+print*,'config%do_surface_sw_spectral_flux   = ', do_surface_sw_spectral_flux
+print*,'config%do_sw_delta_scaling_with_gases   = ', &
+do_sw_delta_scaling_with_gases
+print*,'config%do_fu_lw_ice_optics_bug   = ', do_fu_lw_ice_optics_bug
+print*,'config%do_canopy_fluxes_sw   = ', do_canopy_fluxes_sw
+print*,'config%do_canopy_fluxes_lw   = ', do_canopy_fluxes_lw
+print*,'config%use_canopy_full_spectrum_sw   = ', use_canopy_full_spectrum_sw
+print*,'config%use_canopy_full_spectrum_lw   = ', use_canopy_full_spectrum_lw
+print*,'config%do_canopy_gases_sw   = ', do_canopy_gases_sw
+print*,'config%do_canopy_gases_lw   = ', do_canopy_gases_lw
+print*,'config%mono_lw_wavelength   = ', mono_lw_wavelength
+print*,'config%mono_lw_total_od   = ', mono_lw_total_od
+print*,'config%mono_sw_total_od   = ', mono_sw_total_od
+print*,'config%mono_lw_single_scattering_albedo   = ', &
+mono_lw_single_scattering_albedo
+print*,'config%mono_sw_single_scattering_albedo   = ', &
+mono_sw_single_scattering_albedo
+print*,'config%mono_lw_asymmetry_factor   = ', mono_lw_asymmetry_factor
+print*,'config%mono_sw_asymmetry_factor   = ', mono_sw_asymmetry_factor
+print*,'config%use_beta_overlap   = ', use_beta_overlap
+print*,'config%cloud_inhom_decorr_scaling   = ', cloud_inhom_decorr_scaling
+print*,'config%clear_to_thick_fraction   = ', clear_to_thick_fraction
+print*,'config%overhead_sun_factor   = ', overhead_sun_factor
+print*,'config%max_gas_od_3d   = ', max_gas_od_3d
+print*,'config%max_cloud_od   = ', max_cloud_od
+print*,'config%max_3d_transfer_rate   = ', max_3d_transfer_rate
+print*,'config%min_cloud_effective_size   = ', &
+max(1.0e-6_jprb,min_cloud_effective_size)
+print*,'config%overhang_factor   = ', encroachment_scaling
+
+print*,'config%directory_name  = ',directory_name
+print*,'config%cloud_pdf_override_file_name  = ',cloud_pdf_override_file_name
+print*,'config%liq_optics_override_file_name  = ',liq_optics_override_file_name
+print*,'config%ice_optics_override_file_name  = ',ice_optics_override_file_name
+print*,'config%aerosol_optics_override_file_name  = ', &
+aerosol_optics_override_file_name
+print*,'config%cloud_fraction_threshold  = ',cloud_fraction_threshold
+print*,'config%cloud_mixing_ratio_threshold  = ',cloud_mixing_ratio_threshold
+print*,'config%n_aerosol_types  = ',n_aerosol_types
+print*,'config%do_save_radiative_properties  = ',do_save_radiative_properties
+print*,'config%do_lw_derivatives  = ',do_lw_derivatives
+print*,'config%do_save_spectral_flux  = ',do_save_spectral_flux
+print*,'config%do_save_gpoint_flux  = ',do_save_gpoint_flux
+print*,'config%do_nearest_spectral_sw_albedo  = ',do_nearest_spectral_sw_albedo
+print*,'config%do_nearest_spectral_lw_emiss   = ',do_nearest_spectral_lw_emiss
+print*,'config%sw_albedo_wavelength_bound     = ',sw_albedo_wavelength_bound
+print*,'config%lw_emiss_wavelength_bound      = ',lw_emiss_wavelength_bound
+print*,'config%i_sw_albedo_index              = ',i_sw_albedo_index
+print*,'config%i_lw_emiss_index               = ',i_lw_emiss_index
+print*,'************************************************************************'
+endif
     if (do_save_gpoint_flux) then
       ! Saving the fluxes every g-point overrides saving as averaged
@@ -850 +919 @@
       ! save anything
       this%do_save_spectral_flux = .true.
+      print*,'config%do_save_spectral_flux = .true.'
     end if
 
@@ -855 +925 @@
     call get_enum_code(liquid_model_name, LiquidModelName, &
          &            'liquid_model_name', this%i_liq_model)
+    print*,'config%i_liq_model =', this%i_liq_model
 
     ! Determine ice optics model
     call get_enum_code(ice_model_name, IceModelName, &
          &            'ice_model_name', this%i_ice_model)
-
+    print*,'config%i_ice_model =', this%i_ice_model
     ! Determine gas optics model
     call get_enum_code(gas_model_name, GasModelName, &
          &            'gas_model_name', this%i_gas_model)
+    print*,'config%%i_gas_model = ', this%i_gas_model
 
     ! Determine solvers
     call get_enum_code(sw_solver_name, SolverName, &
          &            'sw_solver_name', this%i_solver_sw)
+    print*,'config%i_solver_sw = ', this%i_solver_sw
     call get_enum_code(lw_solver_name, SolverName, &
          &            'lw_solver_name', this%i_solver_lw)
-
+    print*,'config%i_solver_lw = ', this%i_solver_lw
     if (len_trim(sw_encroachment_name) > 1) then
       call get_enum_code(sw_encroachment_name, EncroachmentName, &
@@ -877 +950 @@
       call get_enum_code(sw_entrapment_name, EntrapmentName, &
            &             'sw_entrapment_name', this%i_3d_sw_entrapment)
+      print*,'config%i_3d_sw_entrapment = ', this%i_3d_sw_entrapment
     end if
 
@@ -882 +956 @@
     call get_enum_code(overlap_scheme_name, OverlapName, &
          &             'overlap_scheme_name', this%i_overlap_scheme)
-    
+    print*,'config%i_overlap_scheme = ', this%i_overlap_scheme
     ! Determine cloud PDF shape 
     call get_enum_code(cloud_pdf_shape_name, PdfShapeName, &
          &             'cloud_pdf_shape_name', this%i_cloud_pdf_shape)
-
+    print*,'config%i_cloud_pdf_shape = ', this%i_cloud_pdf_shape
     this%i_aerosol_type_map = 0
     if (this%use_aerosols) then
       this%i_aerosol_type_map(1:n_aerosol_types) &
            &  = i_aerosol_type_map(1:n_aerosol_types)
+      print*,'config%i_aerosol_type_map = ', this%i_aerosol_type_map
     end if
 
@@ -900 +975 @@
       this%do_clouds = .false.
     end if
+    print*,'config%do_clouds = ', this%do_clouds
 
     ! Normal subroutine exit