| 1 | ! ecrad_driver_config.F90 - Configure driver for offline ecRad radiation scheme |
|---|
| 2 | ! |
|---|
| 3 | ! (C) Copyright 2015- ECMWF. |
|---|
| 4 | ! |
|---|
| 5 | ! This software is licensed under the terms of the Apache Licence Version 2.0 |
|---|
| 6 | ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. |
|---|
| 7 | ! |
|---|
| 8 | ! In applying this licence, ECMWF does not waive the privileges and immunities |
|---|
| 9 | ! granted to it by virtue of its status as an intergovernmental organisation |
|---|
| 10 | ! nor does it submit to any jurisdiction. |
|---|
| 11 | ! |
|---|
| 12 | ! Author: Robin Hogan |
|---|
| 13 | ! Email: r.j.hogan@ecmwf.int |
|---|
| 14 | |
|---|
| 15 | module ecrad_driver_config |
|---|
| 16 | |
|---|
| 17 | use parkind1, only : jprb |
|---|
| 18 | |
|---|
| 19 | implicit none |
|---|
| 20 | |
|---|
| 21 | public |
|---|
| 22 | |
|---|
| 23 | ! Max length of "experiment" global attribute |
|---|
| 24 | integer, parameter :: NMaxStringLength = 2000 |
|---|
| 25 | |
|---|
| 26 | ! Maximum number of spectral diagnostics |
|---|
| 27 | integer, parameter :: NMaxSpectralDiag = 256 |
|---|
| 28 | |
|---|
| 29 | type driver_config_type |
|---|
| 30 | |
|---|
| 31 | ! Parallel settings |
|---|
| 32 | logical :: do_parallel |
|---|
| 33 | integer :: nblocksize ! Number of columns processed at once |
|---|
| 34 | |
|---|
| 35 | ! Override values from the radiation_override namelist (mostly |
|---|
| 36 | ! related to clouds): these will override any values in the |
|---|
| 37 | ! NetCDF data file (or scale them) |
|---|
| 38 | real(jprb) :: fractional_std_override |
|---|
| 39 | real(jprb) :: overlap_decorr_length_override |
|---|
| 40 | real(jprb) :: high_inv_effective_size_override = -1.0_jprb ! m-1 |
|---|
| 41 | real(jprb) :: middle_inv_effective_size_override = -1.0_jprb ! m-1 |
|---|
| 42 | real(jprb) :: low_inv_effective_size_override = -1.0_jprb ! m-1 |
|---|
| 43 | real(jprb) :: effective_size_scaling |
|---|
| 44 | real(jprb) :: sw_albedo_override |
|---|
| 45 | real(jprb) :: lw_emissivity_override |
|---|
| 46 | real(jprb) :: q_liq_scaling, q_ice_scaling |
|---|
| 47 | real(jprb) :: cloud_fraction_scaling |
|---|
| 48 | real(jprb) :: overlap_decorr_length_scaling |
|---|
| 49 | real(jprb) :: skin_temperature_override ! K |
|---|
| 50 | real(jprb) :: solar_irradiance_override ! W m-2 |
|---|
| 51 | real(jprb) :: solar_cycle_multiplier_override |
|---|
| 52 | real(jprb) :: cos_sza_override |
|---|
| 53 | real(jprb) :: cloud_inhom_separation_factor = 1.0_jprb |
|---|
| 54 | real(jprb) :: cloud_separation_scale_surface = -1.0_jprb |
|---|
| 55 | real(jprb) :: cloud_separation_scale_toa = -1.0_jprb |
|---|
| 56 | real(jprb) :: cloud_separation_scale_power = 1.0_jprb |
|---|
| 57 | real(jprb) :: h2o_scaling = 1.0_jprb |
|---|
| 58 | real(jprb) :: co2_scaling = 1.0_jprb |
|---|
| 59 | real(jprb) :: o3_scaling = 1.0_jprb |
|---|
| 60 | real(jprb) :: co_scaling = 1.0_jprb |
|---|
| 61 | real(jprb) :: ch4_scaling = 1.0_jprb |
|---|
| 62 | real(jprb) :: n2o_scaling = 1.0_jprb |
|---|
| 63 | real(jprb) :: o2_scaling = 1.0_jprb |
|---|
| 64 | real(jprb) :: cfc11_scaling = 1.0_jprb |
|---|
| 65 | real(jprb) :: cfc12_scaling = 1.0_jprb |
|---|
| 66 | real(jprb) :: hcfc22_scaling = 1.0_jprb |
|---|
| 67 | real(jprb) :: ccl4_scaling = 1.0_jprb |
|---|
| 68 | real(jprb) :: no2_scaling = 1.0_jprb |
|---|
| 69 | |
|---|
| 70 | ! Optional monotonically increasing wavelength bounds (m) for |
|---|
| 71 | ! shortwave spectral flux diagnostics, to be written to |
|---|
| 72 | ! sw_diagnostic_file_name |
|---|
| 73 | real(jprb) :: sw_diag_wavelength_bound(NMaxSpectralDiag+1) = -1.0_jprb |
|---|
| 74 | |
|---|
| 75 | ! Name of file to write shortwave spectral diagnostics to, but |
|---|
| 76 | ! only if sw_diag_wavelength_bound is populated via the namelist |
|---|
| 77 | character(len=NMaxStringLength) :: sw_diag_file_name = 'sw_diagnostics.nc' |
|---|
| 78 | |
|---|
| 79 | ! Number of shortwave diagnostics, worked out from first |
|---|
| 80 | ! unassigned value in sw_diag_wavelength_bound after reading |
|---|
| 81 | ! namelist |
|---|
| 82 | integer :: n_sw_diag |
|---|
| 83 | |
|---|
| 84 | ! Volume mixing ratios (m3 m-3) in model layers (or equivalently |
|---|
| 85 | ! mole fractions (mol mol-1)) are typically stored in the input |
|---|
| 86 | ! file with a name like co2_vmr, but the suffix can be overridden |
|---|
| 87 | ! by the user |
|---|
| 88 | character(len=32) :: vmr_suffix_str = '_vmr' |
|---|
| 89 | |
|---|
| 90 | ! Process a limited number of columns (iendcol=0 indicates to |
|---|
| 91 | ! process from istartcol up to the end) |
|---|
| 92 | integer :: istartcol, iendcol |
|---|
| 93 | |
|---|
| 94 | ! Save inputs in "inputs.nc" |
|---|
| 95 | logical :: do_save_inputs |
|---|
| 96 | |
|---|
| 97 | ! Save aerosol optical properties to "aerosol_optics.nc" |
|---|
| 98 | logical :: do_save_aerosol_optics |
|---|
| 99 | |
|---|
| 100 | ! Save aerosol optical properties to "hydrometeor_optics*.nc" |
|---|
| 101 | logical :: do_save_cloud_optics |
|---|
| 102 | |
|---|
| 103 | ! Save only net and surface/TOA fluxes, rather than up and down |
|---|
| 104 | logical :: do_save_net_fluxes |
|---|
| 105 | |
|---|
| 106 | ! Do we ignore the inv_inhom_effective_size variable and instead |
|---|
| 107 | ! assume the scale of cloud inhomogeneities is the same as the |
|---|
| 108 | ! scale of the clouds themselves? |
|---|
| 109 | logical :: do_ignore_inhom_effective_size = .false. |
|---|
| 110 | |
|---|
| 111 | ! Number of repeats (for benchmarking) |
|---|
| 112 | integer :: nrepeat |
|---|
| 113 | |
|---|
| 114 | ! Do we correct unphysical inputs (e.g. negative gas concentrations)? |
|---|
| 115 | logical :: do_correct_unphysical_inputs = .false. |
|---|
| 116 | |
|---|
| 117 | ! Do we write NetCDF4/HDF5 file format, needed for very large |
|---|
| 118 | ! files? |
|---|
| 119 | logical :: do_write_hdf5 = .false. |
|---|
| 120 | |
|---|
| 121 | ! Do we write fluxes in double precision? |
|---|
| 122 | logical :: do_write_double_precision = .false. |
|---|
| 123 | |
|---|
| 124 | ! Name of the experiment, to save in output file |
|---|
| 125 | character(len=NMaxStringLength) :: experiment_name = '' |
|---|
| 126 | |
|---|
| 127 | ! Control verbosity in driver routine: 0=none (no output to |
|---|
| 128 | ! standard output; write to standard error only if an error |
|---|
| 129 | ! occurs), 1=warning, 2=info, 3=progress, 4=detailed, 5=debug |
|---|
| 130 | integer :: iverbose |
|---|
| 131 | |
|---|
| 132 | contains |
|---|
| 133 | procedure :: read => read_config_from_namelist |
|---|
| 134 | |
|---|
| 135 | end type driver_config_type |
|---|
| 136 | |
|---|
| 137 | contains |
|---|
| 138 | |
|---|
| 139 | !--------------------------------------------------------------------- |
|---|
| 140 | ! This subroutine reads configuration data from a namelist file, and |
|---|
| 141 | ! anything that is not in the namelists will be set to default |
|---|
| 142 | ! values. If optional output argument "is_success" is present, then on |
|---|
| 143 | ! error (e.g. missing file) it will be set to .false.; if this |
|---|
| 144 | ! argument is missing then on error the program will be aborted. |
|---|
| 145 | subroutine read_config_from_namelist(this, file_name, is_success) |
|---|
| 146 | |
|---|
| 147 | use radiation_io, only : nulerr, radiation_abort |
|---|
| 148 | |
|---|
| 149 | class(driver_config_type), intent(inout) :: this |
|---|
| 150 | character(*), intent(in) :: file_name |
|---|
| 151 | logical, intent(out), optional :: is_success |
|---|
| 152 | |
|---|
| 153 | integer :: iosopen ! Status after calling open |
|---|
| 154 | |
|---|
| 155 | ! Override and scaling values |
|---|
| 156 | real(jprb) :: fractional_std |
|---|
| 157 | real(jprb) :: overlap_decorr_length |
|---|
| 158 | real(jprb) :: inv_effective_size |
|---|
| 159 | real(jprb) :: high_inv_effective_size |
|---|
| 160 | real(jprb) :: middle_inv_effective_size |
|---|
| 161 | real(jprb) :: low_inv_effective_size |
|---|
| 162 | real(jprb) :: effective_size_scaling |
|---|
| 163 | real(jprb) :: sw_albedo |
|---|
| 164 | real(jprb) :: lw_emissivity |
|---|
| 165 | real(jprb) :: q_liquid_scaling, q_ice_scaling |
|---|
| 166 | real(jprb) :: cloud_fraction_scaling |
|---|
| 167 | real(jprb) :: overlap_decorr_length_scaling |
|---|
| 168 | real(jprb) :: skin_temperature |
|---|
| 169 | real(jprb) :: cos_solar_zenith_angle |
|---|
| 170 | real(jprb) :: solar_irradiance_override |
|---|
| 171 | real(jprb) :: solar_cycle_multiplier_override |
|---|
| 172 | real(jprb) :: cloud_inhom_separation_factor |
|---|
| 173 | real(jprb) :: cloud_separation_scale_surface |
|---|
| 174 | real(jprb) :: cloud_separation_scale_toa |
|---|
| 175 | real(jprb) :: cloud_separation_scale_power |
|---|
| 176 | real(jprb) :: h2o_scaling |
|---|
| 177 | real(jprb) :: co2_scaling |
|---|
| 178 | real(jprb) :: o3_scaling |
|---|
| 179 | real(jprb) :: co_scaling |
|---|
| 180 | real(jprb) :: ch4_scaling |
|---|
| 181 | real(jprb) :: n2o_scaling |
|---|
| 182 | real(jprb) :: o2_scaling |
|---|
| 183 | real(jprb) :: cfc11_scaling |
|---|
| 184 | real(jprb) :: cfc12_scaling |
|---|
| 185 | real(jprb) :: hcfc22_scaling |
|---|
| 186 | real(jprb) :: ccl4_scaling |
|---|
| 187 | real(jprb) :: no2_scaling |
|---|
| 188 | real(jprb) :: sw_diag_wavelength_bound(NMaxSpectralDiag+1) |
|---|
| 189 | character(len=NMaxStringLength) :: sw_diag_file_name |
|---|
| 190 | character(len=32) :: vmr_suffix_str |
|---|
| 191 | character(len=NMaxStringLength) :: experiment_name |
|---|
| 192 | |
|---|
| 193 | ! Parallel settings |
|---|
| 194 | logical :: do_parallel |
|---|
| 195 | integer :: nblocksize |
|---|
| 196 | |
|---|
| 197 | logical :: do_save_inputs, do_save_aerosol_optics, do_save_net_fluxes, & |
|---|
| 198 | & do_save_cloud_optics, do_ignore_inhom_effective_size, & |
|---|
| 199 | & do_correct_unphysical_inputs, do_write_hdf5, & |
|---|
| 200 | & do_write_double_precision |
|---|
| 201 | integer :: nrepeat |
|---|
| 202 | |
|---|
| 203 | ! Process a limited number of columns (iendcol=0 indicates to |
|---|
| 204 | ! process from istartcol up to the end) |
|---|
| 205 | integer :: istartcol, iendcol |
|---|
| 206 | |
|---|
| 207 | ! Verbosity |
|---|
| 208 | integer :: iverbose |
|---|
| 209 | |
|---|
| 210 | ! Are we going to override the effective size? |
|---|
| 211 | logical :: do_override_eff_size |
|---|
| 212 | |
|---|
| 213 | ! Loop index |
|---|
| 214 | integer :: jdiag |
|---|
| 215 | |
|---|
| 216 | namelist /radiation_driver/ fractional_std, & |
|---|
| 217 | & overlap_decorr_length, inv_effective_size, sw_albedo, & |
|---|
| 218 | & high_inv_effective_size, middle_inv_effective_size, & |
|---|
| 219 | & low_inv_effective_size, cloud_inhom_separation_factor, & |
|---|
| 220 | & effective_size_scaling, cos_solar_zenith_angle, & |
|---|
| 221 | & lw_emissivity, q_liquid_scaling, q_ice_scaling, & |
|---|
| 222 | & istartcol, iendcol, solar_irradiance_override, & |
|---|
| 223 | & solar_cycle_multiplier_override, & |
|---|
| 224 | & cloud_fraction_scaling, overlap_decorr_length_scaling, & |
|---|
| 225 | & skin_temperature, do_parallel, nblocksize, iverbose, & |
|---|
| 226 | & nrepeat, do_save_inputs, do_ignore_inhom_effective_size, & |
|---|
| 227 | & do_save_aerosol_optics, do_save_net_fluxes, do_save_cloud_optics, & |
|---|
| 228 | & cloud_separation_scale_toa, cloud_separation_scale_surface, & |
|---|
| 229 | & cloud_separation_scale_power, do_correct_unphysical_inputs, & |
|---|
| 230 | & do_write_hdf5, h2o_scaling, co2_scaling, o3_scaling, co_scaling, & |
|---|
| 231 | & ch4_scaling, o2_scaling, cfc11_scaling, cfc12_scaling, & |
|---|
| 232 | & hcfc22_scaling, no2_scaling, n2o_scaling, ccl4_scaling, & |
|---|
| 233 | & vmr_suffix_str, experiment_name, do_write_double_precision, & |
|---|
| 234 | & sw_diag_wavelength_bound, sw_diag_file_name |
|---|
| 235 | |
|---|
| 236 | ! Default values |
|---|
| 237 | do_parallel = .true. |
|---|
| 238 | do_save_inputs = .false. |
|---|
| 239 | do_save_aerosol_optics = .false. |
|---|
| 240 | do_save_cloud_optics = .false. |
|---|
| 241 | do_save_net_fluxes = .false. |
|---|
| 242 | do_ignore_inhom_effective_size = .false. |
|---|
| 243 | nblocksize = 8 |
|---|
| 244 | |
|---|
| 245 | ! Negative values indicate no override will take place |
|---|
| 246 | fractional_std = -1.0_jprb |
|---|
| 247 | overlap_decorr_length = -1.0_jprb |
|---|
| 248 | inv_effective_size = -1.0_jprb |
|---|
| 249 | high_inv_effective_size = -1.0_jprb |
|---|
| 250 | middle_inv_effective_size = -1.0_jprb |
|---|
| 251 | low_inv_effective_size = -1.0_jprb |
|---|
| 252 | effective_size_scaling = -1.0_jprb |
|---|
| 253 | sw_albedo = -1.0_jprb |
|---|
| 254 | lw_emissivity = -1.0_jprb |
|---|
| 255 | q_liquid_scaling = -1.0_jprb |
|---|
| 256 | q_ice_scaling = -1.0_jprb |
|---|
| 257 | cloud_fraction_scaling = -1.0_jprb |
|---|
| 258 | overlap_decorr_length_scaling = -1.0_jprb |
|---|
| 259 | skin_temperature = -1.0_jprb |
|---|
| 260 | cos_solar_zenith_angle = -1.0_jprb |
|---|
| 261 | solar_irradiance_override = -1.0_jprb |
|---|
| 262 | solar_cycle_multiplier_override = -2.0e6_jprb |
|---|
| 263 | cloud_inhom_separation_factor = 1.0_jprb |
|---|
| 264 | cloud_separation_scale_toa = -1.0_jprb |
|---|
| 265 | cloud_separation_scale_surface = -1.0_jprb |
|---|
| 266 | cloud_separation_scale_power = 1.0_jprb |
|---|
| 267 | h2o_scaling = 1.0_jprb |
|---|
| 268 | co2_scaling = 1.0_jprb |
|---|
| 269 | o3_scaling = 1.0_jprb |
|---|
| 270 | co_scaling = 1.0_jprb |
|---|
| 271 | ch4_scaling = 1.0_jprb |
|---|
| 272 | n2o_scaling = 1.0_jprb |
|---|
| 273 | o2_scaling = 1.0_jprb |
|---|
| 274 | cfc11_scaling = 1.0_jprb |
|---|
| 275 | cfc12_scaling = 1.0_jprb |
|---|
| 276 | hcfc22_scaling = 1.0_jprb |
|---|
| 277 | ccl4_scaling = 1.0_jprb |
|---|
| 278 | no2_scaling = 1.0_jprb |
|---|
| 279 | vmr_suffix_str = '_vmr'; |
|---|
| 280 | iverbose = 2 ! Default verbosity is "warning" |
|---|
| 281 | istartcol = 0 |
|---|
| 282 | iendcol = 0 |
|---|
| 283 | nrepeat = 1 |
|---|
| 284 | do_correct_unphysical_inputs = .false. |
|---|
| 285 | do_write_hdf5 = .false. |
|---|
| 286 | do_write_double_precision = .false. |
|---|
| 287 | experiment_name = '' |
|---|
| 288 | sw_diag_wavelength_bound = this%sw_diag_wavelength_bound |
|---|
| 289 | sw_diag_file_name = this%sw_diag_file_name |
|---|
| 290 | |
|---|
| 291 | ! Open the namelist file and read the radiation_driver namelist |
|---|
| 292 | open(unit=10, iostat=iosopen, file=trim(file_name)) |
|---|
| 293 | if (iosopen /= 0) then |
|---|
| 294 | ! An error occurred |
|---|
| 295 | if (present(is_success)) then |
|---|
| 296 | is_success = .false. |
|---|
| 297 | ! We now continue the subroutine so that the default values |
|---|
| 298 | ! are placed in the config structure |
|---|
| 299 | else |
|---|
| 300 | write(nulerr,'(a,a,a)') '*** Error: namelist file "', & |
|---|
| 301 | & trim(file_name), '" not found' |
|---|
| 302 | call radiation_abort('Driver configuration error') |
|---|
| 303 | end if |
|---|
| 304 | else |
|---|
| 305 | ! Read the radiation_driver namelist, noting that it is not an |
|---|
| 306 | ! error if this namelist is not present, provided all the required |
|---|
| 307 | ! variables are present in the NetCDF data file instead |
|---|
| 308 | read(unit=10, nml=radiation_driver) |
|---|
| 309 | close(unit=10) |
|---|
| 310 | end if |
|---|
| 311 | |
|---|
| 312 | ! Copy namelist data into configuration object |
|---|
| 313 | this%do_parallel = do_parallel |
|---|
| 314 | this%do_save_inputs = do_save_inputs |
|---|
| 315 | this%do_save_aerosol_optics = do_save_aerosol_optics |
|---|
| 316 | this%do_save_cloud_optics = do_save_cloud_optics |
|---|
| 317 | this%do_save_net_fluxes = do_save_net_fluxes |
|---|
| 318 | this%do_ignore_inhom_effective_size = do_ignore_inhom_effective_size |
|---|
| 319 | this%nblocksize = nblocksize |
|---|
| 320 | this%iverbose = iverbose |
|---|
| 321 | this%nrepeat = nrepeat |
|---|
| 322 | if (istartcol < 1) then |
|---|
| 323 | this%istartcol = 1 |
|---|
| 324 | else |
|---|
| 325 | this%istartcol = istartcol |
|---|
| 326 | end if |
|---|
| 327 | if (iendcol < 1) then |
|---|
| 328 | this%iendcol = 0 |
|---|
| 329 | else |
|---|
| 330 | this%iendcol = iendcol |
|---|
| 331 | end if |
|---|
| 332 | |
|---|
| 333 | ! Set override values |
|---|
| 334 | this%fractional_std_override = fractional_std |
|---|
| 335 | this%overlap_decorr_length_override = overlap_decorr_length |
|---|
| 336 | |
|---|
| 337 | do_override_eff_size = .false. |
|---|
| 338 | if (inv_effective_size >= 0.0_jprb) then |
|---|
| 339 | this%high_inv_effective_size_override = inv_effective_size |
|---|
| 340 | this%middle_inv_effective_size_override = inv_effective_size |
|---|
| 341 | this%low_inv_effective_size_override = inv_effective_size |
|---|
| 342 | end if |
|---|
| 343 | if (high_inv_effective_size >= 0.0_jprb) then |
|---|
| 344 | this%high_inv_effective_size_override = high_inv_effective_size |
|---|
| 345 | do_override_eff_size = .true. |
|---|
| 346 | end if |
|---|
| 347 | if (middle_inv_effective_size >= 0.0_jprb) then |
|---|
| 348 | this%middle_inv_effective_size_override = middle_inv_effective_size |
|---|
| 349 | do_override_eff_size = .true. |
|---|
| 350 | end if |
|---|
| 351 | if (low_inv_effective_size >= 0.0_jprb) then |
|---|
| 352 | this%low_inv_effective_size_override = low_inv_effective_size |
|---|
| 353 | do_override_eff_size = .true. |
|---|
| 354 | end if |
|---|
| 355 | |
|---|
| 356 | if (do_override_eff_size & |
|---|
| 357 | & .and. (this%high_inv_effective_size_override < 0.0_jprb & |
|---|
| 358 | .or. this%middle_inv_effective_size_override < 0.0_jprb & |
|---|
| 359 | .or. this%low_inv_effective_size_override < 0.0_jprb)) then |
|---|
| 360 | write(nulerr,'(a)') '*** Error: inverse effective cloud size not specified for high, middle and low clouds"' |
|---|
| 361 | call radiation_abort('Driver configuration error') |
|---|
| 362 | end if |
|---|
| 363 | |
|---|
| 364 | this%effective_size_scaling = effective_size_scaling |
|---|
| 365 | this%sw_albedo_override = sw_albedo |
|---|
| 366 | this%lw_emissivity_override = lw_emissivity |
|---|
| 367 | this%q_liq_scaling = q_liquid_scaling |
|---|
| 368 | this%q_ice_scaling = q_ice_scaling |
|---|
| 369 | this%cloud_fraction_scaling = cloud_fraction_scaling |
|---|
| 370 | this%overlap_decorr_length_scaling = overlap_decorr_length_scaling |
|---|
| 371 | this%skin_temperature_override = skin_temperature |
|---|
| 372 | this%cos_sza_override = cos_solar_zenith_angle |
|---|
| 373 | this%solar_irradiance_override = solar_irradiance_override |
|---|
| 374 | this%solar_cycle_multiplier_override = solar_cycle_multiplier_override |
|---|
| 375 | this%cloud_inhom_separation_factor = cloud_inhom_separation_factor |
|---|
| 376 | this%cloud_separation_scale_toa = cloud_separation_scale_toa |
|---|
| 377 | this%cloud_separation_scale_surface = cloud_separation_scale_surface |
|---|
| 378 | this%cloud_separation_scale_power = cloud_separation_scale_power |
|---|
| 379 | this%do_correct_unphysical_inputs = do_correct_unphysical_inputs |
|---|
| 380 | this%do_write_hdf5 = do_write_hdf5 |
|---|
| 381 | this%do_write_double_precision = do_write_double_precision |
|---|
| 382 | this%h2o_scaling = h2o_scaling |
|---|
| 383 | this%co2_scaling = co2_scaling |
|---|
| 384 | this%o3_scaling = o3_scaling |
|---|
| 385 | this%co_scaling = co_scaling |
|---|
| 386 | this%ch4_scaling = ch4_scaling |
|---|
| 387 | this%n2o_scaling = n2o_scaling |
|---|
| 388 | this%o2_scaling = o2_scaling |
|---|
| 389 | this%cfc11_scaling = cfc11_scaling |
|---|
| 390 | this%cfc12_scaling = cfc12_scaling |
|---|
| 391 | this%hcfc22_scaling = hcfc22_scaling |
|---|
| 392 | this%ccl4_scaling = ccl4_scaling |
|---|
| 393 | this%no2_scaling = no2_scaling |
|---|
| 394 | this%vmr_suffix_str = trim(vmr_suffix_str) |
|---|
| 395 | this%experiment_name= trim(experiment_name) |
|---|
| 396 | |
|---|
| 397 | this%sw_diag_file_name = trim(sw_diag_file_name) |
|---|
| 398 | this%sw_diag_wavelength_bound = sw_diag_wavelength_bound |
|---|
| 399 | ! Work out number of shortwave diagnostics from first negative |
|---|
| 400 | ! wavelength bound, noting that the number of diagnostics is one |
|---|
| 401 | ! fewer than the number of valid bounds |
|---|
| 402 | do jdiag = 0,NMaxSpectralDiag |
|---|
| 403 | if (this%sw_diag_wavelength_bound(jdiag+1) < 0.0_jprb) then |
|---|
| 404 | this%n_sw_diag = max(0,jdiag-1) |
|---|
| 405 | exit |
|---|
| 406 | end if |
|---|
| 407 | end do |
|---|
| 408 | |
|---|
| 409 | end subroutine read_config_from_namelist |
|---|
| 410 | |
|---|
| 411 | end module ecrad_driver_config |
|---|
