Changeset 4182 for LMDZ6/branches/LMDZ_ECRad

Timestamp:

Jun 20, 2022, 5:29:41 PM (3 years ago)

Author:

lguez

Message:

Use custom procedure to read aerosol optics file

Replace call to config%aerosol_optics%setup by call to an LMDZ
procedure, setup_aerosol_optics_lmdz, adapted to our aerosol optics
NetCDF file. The main differences between our file and the file
expected by config%aerosol_optics%setup are:

• Our file has groups.

• The SW bands in our file are ordered (as befits a NetCDF coordinate).

• Our file does not contain variables for monochromatic asymmetry factor, single-scattering albedo and lidar backscattering ratio. These variables are not used by ECRad but they are nevertheless read and required by config%aerosol_optics%setup.

Location:

LMDZ6/branches/LMDZ_ECRad/libf/phylmd

Files:

• ecrad/radiation_aerosol_optics.F90 (modified) (2 diffs)
• setup_aerosol_optics_lmdz_m.F90 (added)

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_aerosol_optics.F90

@@ -36 +36 @@
     use radiation_aerosol_optics_data, only : aerosol_optics_type
     use radiation_io,                  only : nulerr, radiation_abort
+    use setup_aerosol_optics_lmdz_m, only: setup_aerosol_optics_lmdz
 
     type(config_type), intent(inout) :: config
@@ -46 +47 @@
       ! Load data from file and prepare to map config%n_aerosol_types
       ! aerosol types
-      call config%aerosol_optics%setup(trim(config%aerosol_optics_file_name), &
-           &                        config%n_aerosol_types, iverbose=config%iverbosesetup)
+       call setup_aerosol_optics_lmdz(config%aerosol_optics, &
+            trim(config%aerosol_optics_file_name))
 
       ! Check agreement in number of bands