Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/ecrad/lmdz/radiation_setup.f90

Timestamp:

Jul 28, 2025, 7:23:15 PM (7 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/ecrad/lmdz/radiation_setup.f90 (modified) (5 diffs)

LMDZ6/branches/contrails/libf/phylmd/ecrad/lmdz/radiation_setup.f90

@@ -1 +1 @@
 MODULE RADIATION_SETUP
 
-! RADIATION_SETUP - Setting up modular radiation scheme
-!
-! AUTHOR
-! ------
-!   Robin Hogan, ECMWF
-!   Original: 2015-09-16
-!
-! MODIFICATIONS
-! -------------
-!   Abderrahmane Idelkadi LMD, juillet 2023
-!
-!-----------------------------------------------------------------------
+  ! RADIATION_SETUP - Setting up modular radiation scheme
+  !
+  ! AUTHOR
+  ! ------
+  !   Robin Hogan, ECMWF
+  !   Original: 2015-09-16
+  !
+  ! MODIFICATIONS
+  ! -------------
+  !   Abderrahmane Idelkadi LMD, juillet 2023
+  !
+  !-----------------------------------------------------------------------
 
   USE PARKIND1,         ONLY : JPRB
@@ -52 +52 @@
   ! to print lots of information during the setup stage (default is
   ! no).
-! AI At the end of the routine, the parameters are read in namelist
-!    
+  ! AI At the end of the routine, the parameters are read in namelist
+  !
   SUBROUTINE SETUP_RADIATION_SCHEME(LOUTPUT,file_name,rad_config,driver_config)
 
-!    USE radiation_config, ONLY : config_type, &
-!        &                       ISolverMcICA, ISolverSpartacus, &
-!        &                       ILiquidModelSlingo, ILiquidModelSOCRATES, &
-!        &                       IIceModelFu, IIceModelBaran, &
-!        &                       IOverlapExponentialRandom
+    !    USE radiation_config, ONLY : config_type, &
+    !        &                       ISolverMcICA, ISolverSpartacus, &
+    !        &                       ILiquidModelSlingo, ILiquidModelSOCRATES, &
+    !        &                       IIceModelFu, IIceModelBaran, &
+    !        &                       IOverlapExponentialRandom
     USE mod_phys_lmdz_para
-      
+
     USE YOMHOOK,  ONLY : LHOOK, DR_HOOK
     USE YOMLUN,   ONLY : NULOUT, NULERR
     USE YOESRTWN, ONLY : NMPSRTM
 
- !   USE radiation_interface,      ONLY : setup_radiation
- !   USE setup_config_from_lmdz,   ONLY : driver_config_type
+    !   USE radiation_interface,      ONLY : setup_radiation
+    !   USE setup_config_from_lmdz,   ONLY : driver_config_type
 
     ! Whether or not to provide information on the radiation scheme
@@ -90 +90 @@
 
     ! Configure verbosity of setup of radiation scheme
- 
+
     print*,'********** Dans radiation_setup *****************'
 
-! AI ATTENTION (parameters read in namelist file)
-!    file_name="namelist_ecrad"
-    call rad_config%read(file_name=file_name)
-    call driver_config%read(file_name)
-
-    ! Use configuration data to set-up radiation scheme, including
-    ! reading scattering datafiles
-    CALL setup_radiation(rad_config)
+    rad_config%i_aerosol_type_map(1:13) = (/ &
+         &  -1, &
+         &  -2, &
+         &  -3, &
+         &  -4, &
+         &  -5, &
+         &  -6, &
+         &  -7, &
+         &   1, &
+         &   2, &
+         &   3, &
+         &  -8, &
+         &  -9, &
+         &   4 /)
 
     ! Populate the mapping between the 14 RRTM shortwave bands and the
@@ -108 +114 @@
     ! results...
     ! Note that NMPSRTM(:)=(/  6, 6, 5, 5, 5, 5, 5, 4, 4, 3, 2, 2, 1, 6 /)
-! AI (6 albedo SW bands) 
+    ! AI (6 albedo SW bands)
     call rad_config%define_sw_albedo_intervals(6, &
-       &  [0.25e-6_jprb, 0.44e-6_jprb, 0.69e-6_jprb, &
-       &   1.19e-6_jprb, 2.38e-6_jprb], [1,2,3,4,5,6])
+         &  [0.25e-6_jprb, 0.44e-6_jprb, 0.69e-6_jprb, &
+         &   1.19e-6_jprb, 2.38e-6_jprb], [1,2,3,4,5,6])
     ! Likewise between the 16 RRTM longwave bands and the 2 emissivity
     ! inputs (info taken from rrtm_ecrt_140gp_mcica.F90) representing
     ! outside and inside the window region of the spectrum
-!     rad_config%i_emiss_from_band_lw = (/ 1,1,1,1,1,2,2,2,1,1,1,1,1,1,1,1 /)
-! AI ATTENTION ?????
-!!    call rad_config%define_lw_emiss_intervals(3, &
-!!         &  (/ 8.0e-6_jprb,13.0e-6_jprb /),  (/ 1,2,1 /))
+    !     rad_config%i_emiss_from_band_lw = (/ 1,1,1,1,1,2,2,2,1,1,1,1,1,1,1,1 /)
+    ! AI ATTENTION ?????
+    !!    call rad_config%define_lw_emiss_intervals(3, &
+    !!         &  (/ 8.0e-6_jprb,13.0e-6_jprb /),  (/ 1,2,1 /))
 
-!    ! Get spectral weightings for UV and PAR
+    ! AI ATTENTION (parameters read in namelist file)
+    !    file_name="namelist_ecrad"
+    call rad_config%read(file_name=file_name)
+    call driver_config%read(file_name)
+    call rad_config%print(iverbose = 2)
+
+    ! Use configuration data to set-up radiation scheme, including
+    ! reading scattering datafiles
+    CALL setup_radiation(rad_config)
+
+    !    ! Get spectral weightings for UV and PAR
     call rad_config%get_sw_weights(0.2e-6_jprb, 0.4415e-6_jprb, &
          &  NWEIGHT_UV, IBAND_UV, WEIGHT_UV, 'ultraviolet')
@@ -127 +143 @@
          &  'photosynthetically active radiation, PAR')
 
-     rad_config%i_aerosol_type_map(1:13) = (/ &
-           &  -1, & 
-           &  -2, & 
-           &  -3, & 
-           &  -4, & 
-           &  -5, & 
-           &  -6, &
-           &  -7, &
-           &   1, &
-           &   2, &
-           &   3, &
-           &  -8, &
-           &  -9, &
-           &   4 /)
-
-
     IF (LHOOK) CALL DR_HOOK('RADIATION_SETUP:SETUP_RADIATION_SCHEME',1,ZHOOK_HANDLE)