Changeset 5675 for LMDZ6/trunk/libf

Timestamp:

May 23, 2025, 2:36:20 PM (4 weeks ago)

Author:

lguez

Message:

Allow saving input of ECRad

Setting do_save_inputs=t saves input of ECRad at the first time step
of physics, in a format expected by the ECRad offline program. This is
useful in particular to test ECRad offline with LMDZ aerosol types.

Location:

LMDZ6/trunk/libf/phylmd/ecrad/lmdz

Files:

• radiation_scheme_mod.f90 (modified) (4 diffs)
• setup_config_from_lmdz.f90 (modified) (5 diffs)

LMDZ6/trunk/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
@@ -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)
@@ -658 +662 @@
 
 IF (LHOOK) CALL DR_HOOK('RADIATION_SCHEME',1,ZHOOK_HANDLE)
+debut_ecrad=.false.
 
 END SUBROUTINE RADIATION_SCHEME

LMDZ6/trunk/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