Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/ecrad/lmdz

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_scheme_mod.f90 (modified) (8 diffs)
• libf/phylmd/ecrad/lmdz/radiation_setup.f90 (modified) (5 diffs)
• libf/phylmd/ecrad/lmdz/setup_config_from_lmdz.f90 (modified) (5 diffs)

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

@@ -65 +65 @@
 
 USE mod_phys_lmdz_para
+use geometry_mod, only: longitude_deg, latitude_deg
 
 IMPLICIT NONE
@@ -160 +161 @@
 
 ! Direct component of surface flux into horizontal plane
-REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR(KLON)
-REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR_CLEAR(KLON)
+REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR(KLON,KLEV+1)
+REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR_CLEAR(KLON,KLEV+1)
 ! As PFLUX_DIR but into a plane perpendicular to the sun
 REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR_INTO_SUN(KLON)
@@ -260 +261 @@
   if (debut_ecrad) then
    call SETUP_RADIATION_SCHEME(loutput,namelist_file,rad_config,driver_config)
-   debut_ecrad=.false. 
   endif 
 !$OMP END MASTER
@@ -541 +541 @@
 !    endif
 !   endif
+if (debut_ecrad .and. driver_config%do_save_inputs) &
+     call save_inputs('inputs.nc', rad_config, single_level, thermodynamics, &
+     gas, cloud, aerosol, lat = latitude_deg, lon = longitude_deg)
+
 CALL radiation(KLON, KLEV, KIDIA, KFDIA, rad_config, &
      &  single_level, thermodynamics, gas, cloud, aerosol, flux)
@@ -585 +589 @@
 !PFLUX_SW_UP_CLEAR_SURF(KIDIA:KFDIA) = flux%sw_up_clear(KIDIA:KFDIA,KLEV+1)
 !PFLUX_LW_UP_CLEAR_SURF(KIDIA:KFDIA) = flux%lw_up_clear(KIDIA:KFDIA,KLEV+1)
-PFLUX_DIR(KIDIA:KFDIA) = flux%sw_dn_direct(KIDIA:KFDIA,KLEV+1)
-PFLUX_DIR_CLEAR(KIDIA:KFDIA) = flux%sw_dn_direct_clear(KIDIA:KFDIA,KLEV+1)
+! Direct component of flux into horizontal plane
+PFLUX_DIR(KIDIA:KFDIA,:) = flux%sw_dn_direct(KIDIA:KFDIA,:)
+PFLUX_DIR_CLEAR(KIDIA:KFDIA,:) = flux%sw_dn_direct_clear(KIDIA:KFDIA,:)
 PFLUX_DIR_INTO_SUN(KIDIA:KFDIA) = 0.0_JPRB
 WHERE (PMU0(KIDIA:KFDIA) > EPSILON(1.0_JPRB))
-  PFLUX_DIR_INTO_SUN(KIDIA:KFDIA) = PFLUX_DIR(KIDIA:KFDIA) / PMU0(KIDIA:KFDIA)
+! Direct Surface component of flux into a plane perpendicular to the sun        
+  PFLUX_DIR_INTO_SUN(KIDIA:KFDIA) = PFLUX_DIR(KIDIA:KFDIA,KLEV+1) / PMU0(KIDIA:KFDIA)
 END WHERE
 ! Top-of-atmosphere downwelling flux
@@ -656 +662 @@
 
 IF (LHOOK) CALL DR_HOOK('RADIATION_SCHEME',1,ZHOOK_HANDLE)
+debut_ecrad=.false.
 
 END SUBROUTINE RADIATION_SCHEME
@@ -849 +856 @@
 
 ! Direct component of surface flux into horizontal plane
-REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR(KLON)
-REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR_CLEAR(KLON)
+REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR(KLON,KLEV+1)
+REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR_CLEAR(KLON,KLEV+1)
 ! As PFLUX_DIR but into a plane perpendicular to the sun
 REAL(KIND=JPRB),  INTENT(OUT) :: PFLUX_DIR_INTO_SUN(KLON)
@@ -1283 +1290 @@
 !PFLUX_SW_UP_CLEAR_SURF(KIDIA:KFDIA) = flux%sw_up_clear(KIDIA:KFDIA,KLEV+1)
 !PFLUX_LW_UP_CLEAR_SURF(KIDIA:KFDIA) = flux%lw_up_clear(KIDIA:KFDIA,KLEV+1)
-PFLUX_DIR(KIDIA:KFDIA) = flux%sw_dn_direct(KIDIA:KFDIA,KLEV+1)
-PFLUX_DIR_CLEAR(KIDIA:KFDIA) = flux%sw_dn_direct_clear(KIDIA:KFDIA,KLEV+1)
+PFLUX_DIR(KIDIA:KFDIA,:) = flux%sw_dn_direct(KIDIA:KFDIA,:)
+PFLUX_DIR_CLEAR(KIDIA:KFDIA,:) = flux%sw_dn_direct_clear(KIDIA:KFDIA,:)
 PFLUX_DIR_INTO_SUN(KIDIA:KFDIA) = 0.0_JPRB
 WHERE (PMU0(KIDIA:KFDIA) > EPSILON(1.0_JPRB))
-  PFLUX_DIR_INTO_SUN(KIDIA:KFDIA) = PFLUX_DIR(KIDIA:KFDIA) / PMU0(KIDIA:KFDIA)
+  PFLUX_DIR_INTO_SUN(KIDIA:KFDIA) = PFLUX_DIR(KIDIA:KFDIA,KLEV+1) / PMU0(KIDIA:KFDIA)
 END WHERE
 ! Top-of-atmosphere downwelling flux

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)
 

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

@@ -27 +27 @@
      real(jprb) :: high_decorrelation_length = 2000.0_jprb
 
+     ! Save inputs in "inputs.nc"
+     logical :: do_save_inputs
  contains
  procedure :: read => read_config_from_namelist
@@ -57 +59 @@
     real(jprb) :: mid_decorrelation_length 
     real(jprb) :: high_decorrelation_length 
+    logical :: do_save_inputs
 
     namelist /radiation_driver/ ok_effective_size, ok_separation, &
@@ -63 +66 @@
          &  high_inv_effective_size, middle_inv_effective_size, low_inv_effective_size, &
          &  cloud_inhom_separation_factor, cloud_separation_scale_surface, &
-         &  cloud_separation_scale_toa, cloud_separation_scale_power 
+         &  cloud_separation_scale_toa, cloud_separation_scale_power, &
+         do_save_inputs
 
     ok_effective_size = .false.
@@ -80 +84 @@
     mid_decorrelation_length = 2000.0_jprb
     high_decorrelation_length = 2000.0_jprb
+    do_save_inputs = .false.
 
     ! Open the namelist file and read the radiation_driver namelist
@@ -118 +123 @@
     this%middle_inv_effective_size = middle_inv_effective_size
     this%low_inv_effective_size = low_inv_effective_size
+    this%do_save_inputs = do_save_inputs
 
   end subroutine read_config_from_namelist