Changeset 4853 for LMDZ6/trunk/libf/phylmd/ecrad/radiation/radiation_config.F90

Timestamp:

Mar 19, 2024, 3:34:21 PM (2 months ago)

Author:

idelkadi

Message:

Ecrad update in LMDZ / Implementation of Ecrad double call in LMDZ

• version 1.6.1 (commit 210d7911380f53a788c3cad73b3cf9b4e022ef87)
• interface routines between lmdz and ecrad grouped in a new "lmdz" directory
• double call method redesigned so as to go through the Ecrad initialization and configuration part only once for the entire simulation
• clean-up in the read.F routine: delete unitules arguments
• modification of default gas model in namelist (default: ECCKD)

File:

• LMDZ6/trunk/libf/phylmd/ecrad/radiation/radiation_config.F90 (modified) (18 diffs)

LMDZ6/trunk/libf/phylmd/ecrad/radiation/radiation_config.F90

@@ -34 +34 @@
 !
 
+#include "ecrad_config.h"
+
 module radiation_config
 
@@ -202 +204 @@
     ! more random numbers being generated?  This is the default on NEC
     ! SX.
-#ifdef __SX__
+#ifdef DWD_VECTOR_OPTIMIZATIONS
     logical :: use_vectorizable_generator = .true.
 #else
@@ -276 +278 @@
 
     ! Codes describing the gas model
-    integer :: i_gas_model = IGasModelIFSRRTMG
+    integer :: i_gas_model_sw = IGasModelIFSRRTMG
+    integer :: i_gas_model_lw = IGasModelIFSRRTMG
 
     ! Optics if i_gas_model==IGasModelMonochromatic.
@@ -702 +705 @@
     character(511) :: ssi_override_file_name
     character(63)  :: liquid_model_name, ice_model_name, gas_model_name
+    character(63)  :: sw_gas_model_name, lw_gas_model_name
     character(63)  :: sw_solver_name, lw_solver_name, overlap_scheme_name
     character(63)  :: sw_entrapment_name, sw_encroachment_name, cloud_pdf_shape_name
@@ -709 +713 @@
          &     .false.,.false.,.false.,.false.,.false.,.false.]
     integer :: i_aerosol_type_map(NMaxAerosolTypes) ! More than 256 is an error
-
+    
     logical :: do_nearest_spectral_sw_albedo
     logical :: do_nearest_spectral_lw_emiss
@@ -716 +720 @@
     integer :: i_sw_albedo_index(NMaxAlbedoIntervals)
     integer :: i_lw_emiss_index (NMaxAlbedoIntervals)
+    integer :: i_gas_model
 
     integer :: iunit ! Unit number of namelist file
@@ -726 +731 @@
          &  do_surface_sw_spectral_flux, do_lw_derivatives, do_toa_spectral_flux, &
          &  do_lw_aerosol_scattering, do_lw_cloud_scattering, &
-         &  n_regions, directory_name, gas_model_name, &
+         &  n_regions, directory_name, gas_model_name, sw_gas_model_name, lw_gas_model_name, &
          &  ice_optics_override_file_name, liq_optics_override_file_name, &
          &  aerosol_optics_override_file_name, cloud_pdf_override_file_name, &
@@ -815 +820 @@
     encroachment_scaling = -1.0_jprb
     gas_model_name = '' !DefaultGasModelName
+    sw_gas_model_name = '' !DefaultGasModelName
+    lw_gas_model_name = '' !DefaultGasModelName
     liquid_model_name = '' !DefaultLiquidModelName
     ice_model_name = '' !DefaultIceModelName
@@ -1014 +1021 @@
          &            'ice_model_name', this%i_ice_model)
 
-    ! Determine gas optics model
+    ! Determine gas optics model(s) - firstly try the generic gas_model_name
+    i_gas_model = -1
     call get_enum_code(gas_model_name, GasModelName, &
-         &            'gas_model_name', this%i_gas_model)
-
+         &            'gas_model_name', i_gas_model)
+    if (i_gas_model > -1) then
+      this%i_gas_model_sw = i_gas_model
+      this%i_gas_model_lw = i_gas_model
+    end if
+    ! ...then the band-specific values
+    call get_enum_code(sw_gas_model_name, GasModelName, &
+         &            'sw_gas_model_name', this%i_gas_model_sw)
+    call get_enum_code(lw_gas_model_name, GasModelName, &
+         &            'lw_gas_model_name', this%i_gas_model_lw)
+   
     ! Determine solvers
     call get_enum_code(sw_solver_name, SolverName, &
@@ -1055 +1072 @@
     end if
 
-    if (this%i_gas_model == IGasModelIFSRRTMG &
-         & .and. (this%use_general_cloud_optics &
-         &        .or. this%use_general_aerosol_optics)) then
-      if (this%do_sw .and. this%do_cloud_aerosol_per_sw_g_point) then
+    if (this%use_general_cloud_optics .or. this%use_general_aerosol_optics) then
+      if (this%do_sw .and. this%do_cloud_aerosol_per_sw_g_point &
+           &  .and. this%i_gas_model_sw == IGasModelIFSRRTMG) then
         write(nulout,'(a)') 'Warning: RRTMG SW only supports cloud/aerosol/surface optical properties per band, not per g-point'
         this%do_cloud_aerosol_per_sw_g_point = .false.
       end if
-      if (this%do_lw .and. this%do_cloud_aerosol_per_lw_g_point) then
+      if (this%do_lw .and. this%do_cloud_aerosol_per_lw_g_point &
+           &  .and. this%i_gas_model_lw == IGasModelIFSRRTMG) then
         write(nulout,'(a)') 'Warning: RRTMG LW only supports cloud/aerosol/surface optical properties per band, not per g-point'
         this%do_cloud_aerosol_per_lw_g_point = .false.
@@ -1128 +1145 @@
 
     ! If ecCKD gas optics model is being used set relevant file names
-    if (this%i_gas_model == IGasModelECCKD) then
+    if (this%i_gas_model_sw == IGasModelECCKD .or. this%i_gas_model_lw == IGasModelECCKD) then
 
       ! This gas optics model usually used with general cloud and
@@ -1138 +1155 @@
         write(nulout,'(a)') 'Warning: ecCKD gas optics model usually used with general aerosol optics'
       end if
+
+    end if
+
+    if (this%i_gas_model_sw == IGasModelECCKD) then
 
       if (len_trim(this%gas_optics_sw_override_file_name) > 0) then
@@ -1153 +1174 @@
       end if
 
+    end if
+
+    if (this%i_gas_model_lw == IGasModelECCKD) then
+
       if (len_trim(this%gas_optics_lw_override_file_name) > 0) then
         if (this%gas_optics_lw_override_file_name(1:1) == '/') then
@@ -1170 +1195 @@
 
     if (this%use_spectral_solar_cycle) then
-      if (this%i_gas_model /= IGasModelECCKD) then
+      if (this%i_gas_model_sw /= IGasModelECCKD) then
         write(nulerr,'(a)') '*** Error: solar cycle only available with ecCKD gas optics model'
         call radiation_abort('Radiation configuration error')
@@ -1278 +1303 @@
     end if
 
-    ! In the monochromatic case we need to override the liquid, ice
-    ! and aerosol models to ensure compatibility
-    if (this%i_gas_model == IGasModelMonochromatic) then
+    if (this%i_gas_model_sw == IGasModelMonochromatic .or. this%i_gas_model_lw == IGasModelMonochromatic) then
+
+      if (this%i_gas_model_sw /= this%i_gas_model_lw) then
+        write(nulerr,'(a,i0)') '*** Error: Monochromatic gas optics model must be used in shortwave and longwave'
+        call radiation_abort('Radiation configuration error')
+      end if
+    
+      ! In the monochromatic case we need to override the liquid, ice
+      ! and aerosol models to ensure compatibility
       this%i_liq_model = ILiquidModelMonochromatic
       this%i_ice_model = IIceModelMonochromatic
       this%use_aerosols = .false.
+      
     end if
 
@@ -1392 +1424 @@
       call print_logical('  Saving spectral flux profiles', &
            &   'do_save_spectral_flux', this%do_save_spectral_flux)
-      call print_enum('  Gas model is', GasModelName, 'i_gas_model', &
-           &          this%i_gas_model)
+      call print_enum('  Shortwave gas model is', GasModelName, 'i_gas_model_sw', &
+           &          this%i_gas_model_sw)
+      call print_enum('  Longwave gas model is', GasModelName, 'i_gas_model_lw', &
+           &          this%i_gas_model_lw)
       call print_logical('  Aerosols are', 'use_aerosols', this%use_aerosols)
       if (this%use_aerosols) then
@@ -1481 +1515 @@
              &          'i_solver_sw', this%i_solver_sw)
         
-        if (this%i_gas_model == IGasModelMonochromatic) then
+        if (this%i_gas_model_sw == IGasModelMonochromatic) then
           call print_real('  Shortwave atmospheric optical depth', &
                &   'mono_sw_total_od', this%mono_sw_total_od)
@@ -1502 +1536 @@
              &          this%i_solver_lw)
 
-        if (this%i_gas_model == IGasModelMonochromatic) then
+        if (this%i_gas_model_lw == IGasModelMonochromatic) then
           if (this%mono_lw_wavelength > 0.0_jprb) then
             call print_real('  Longwave effective wavelength (m)', &