Changeset 4853 for LMDZ6/trunk/libf/phylmd/ecrad/radiation/radiation_spectral_definition.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_spectral_definition.F90 (modified) (6 diffs)

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

@@ -15 +15 @@
 !
 
+#include "ecrad_config.h"
+
 module radiation_spectral_definition
 
@@ -23 +25 @@
   public
 
-  real(jprb), parameter :: SolarReferenceTemperature = 5777.0_jprb ! K
+  real(jprb), parameter :: SolarReferenceTemperature       = 5777.0_jprb ! K
   real(jprb), parameter :: TerrestrialReferenceTemperature = 273.15_jprb ! K
 
@@ -85 +87 @@
   subroutine read_spectral_definition(this, file)
 
-    use easy_netcdf, only : netcdf_file
+#ifdef EASY_NETCDF_READ_MPI
+    use easy_netcdf_read_mpi, only : netcdf_file
+#else
+    use easy_netcdf,          only : netcdf_file
+#endif
     use yomhook,     only : lhook, dr_hook, jphook
 
@@ -132 +138 @@
 
   !---------------------------------------------------------------------
-  ! Store a simple band description by copying over the lower and
-  ! upper wavenumbers of each band
-  subroutine allocate_bands_only(this, wavenumber1, wavenumber2)
+  ! Store a simple band description by copying over the reference
+  ! temperature and the lower and upper wavenumbers of each band
+  subroutine allocate_bands_only(this, reference_temperature, wavenumber1, wavenumber2)
 
     use yomhook,     only : lhook, dr_hook, jphook
 
     class(spectral_definition_type), intent(inout) :: this
-    real(jprb),        dimension(:), intent(in)    :: wavenumber1, wavenumber2
+    real(jprb),                      intent(in)    :: reference_temperature    ! K
+    real(jprb),        dimension(:), intent(in)    :: wavenumber1, wavenumber2 ! cm-1
 
     real(jphook) :: hook_handle
@@ -152 +159 @@
     this%wavenumber1_band = wavenumber1
     this%wavenumber2_band = wavenumber2
-
+    this%reference_temperature = reference_temperature
+    
     if (lhook) call dr_hook('radiation_spectral_definition:allocate_bands_only',1,hook_handle)
 
@@ -167 +175 @@
     this%ng    = 0
     this%nband = 0
+    this%reference_temperature = -1.0_jprb
 
     if (allocated(this%wavenumber1))      deallocate(this%wavenumber1)