GFS Operational Physics Documentation  Revision: 81451
radiation_aerosols.f File Reference

This file contains climatological atmospheric aerosol schemes for radiation computations. More...

Data Types

type  module_radiation_aerosols::gocart_index_type
 index for gocart aerosol species to be included in the calculations of aerosol optical properties (ext, ssa, asy) More...
 
type  module_radiation_aerosols::tracer_index_type
 index for gocart aerosols from prognostic tracer fields More...
 

Functions/Subroutines

subroutine module_radiation_aerosols::aer_property
 This subroutine maps the 5 degree global climatological aerosol data set onto model grids, and compute aerosol optical properties for SW and LW radiations. More...
 
subroutine radclimaer
 This subroutine computes aerosols optical properties in NSWLWBD bands. there are seven different vertical profile structures. in the troposphere, aerosol distribution at each grid point is composed from up to six components out of ten different substances.
 
subroutine module_radiation_aerosols::gocart_init
 The initialization program for gocart aerosols. More...
 
subroutine set_aerspc (raddt, fdaer)
 This subroutine determines merging coefficients ctaer; setup aerosol specification.
 
subroutine rd_gocart_luts
 This subroutine reads input gocart aerosol optical data from Mie code calculations.
 
subroutine optavg_grt
 This subroutine computes mean aerosols optical properties over each SW/LW radiation spectral band for each of the species components. This program follows GFDL's approach for thick cloud optical property in SW radiation scheme (2000).
 
subroutine rd_gocart_clim
 This subroutine: More...
 
subroutine module_radiation_aerosols::setgocartaer
 This subroutine computes SW + LW aerosol optical properties for gocart aerosol species (merged from fcst and clim fields). More...
 
subroutine map_aermr
 This subroutine maps input tracer fields (trcly) to local tracer array (aermr).
 
subroutine aeropt_grt
 This subroutine computes aerosols optical properties in NSWLWBD SW/LW bands. Aerosol distribution at each grid point is composed from up to NMXG aerosol species (from NUM_GRIDCOMP components).
 

Variables

integer, parameter, public module_radiation_aerosols::nf_aesw = 3
 num of output fields for SW rad
 
integer, parameter, public module_radiation_aerosols::nf_aelw = 3
 num of output fields for LW rad
 
integer, parameter, public module_radiation_aerosols::nlwstr = 1
 starting band number in ir region
 
integer, parameter, public module_radiation_aerosols::nspc = 5
 num of species for output aod (opnl)
 
integer, parameter, public module_radiation_aerosols::nspc1 = NSPC + 1
 total+species
 
real(kind=kind_phys), parameter module_radiation_aerosols::f_zero = 0.0
 
real(kind=kind_phys), parameter module_radiation_aerosols::f_one = 1.0
 
integer, save module_radiation_aerosols::nswbnd = NBDSW
 number of actual bands for sw aerosols; calculated according to laswflg setting
 
integer, save module_radiation_aerosols::nlwbnd = NBDLW
 number of actual bands for lw aerosols; calculated according to lalwflg and lalw1bd settings
 
integer, save module_radiation_aerosols::nswlwbd = NBDSW+NBDLW
 total number of bands for sw+lw aerosols
 
integer, parameter, public module_radiation_aerosols::nwvsol = 151
 num of wvnum regions where solar flux is constant
 
integer, parameter, public module_radiation_aerosols::nwvtot = 57600
 total num of wvnum included
 
integer, parameter, public module_radiation_aerosols::nwvtir = 4000
 total num of wvnum in ir range
 
integer, dimension(nwvsol), save module_radiation_aerosols::nwvns0
 number of wavenumbers in each region where the solar flux is constant
 
real(kind=kind_phys), dimension(nwvsol), save module_radiation_aerosols::s0intv
 solar flux \(w/m^2\) in each wvnumb region where it is constant
 
integer, parameter module_radiation_aerosols::minvyr = 1850
 lower limit (year) data available
 
integer, parameter module_radiation_aerosols::maxvyr = 1999
 upper limit (year) data available
 
integer, dimension(:,:,:), allocatable, save module_radiation_aerosols::ivolae
 monthly, 45-deg lat-zone aerosols data set in subroutine 'aer_init'
 
integer module_radiation_aerosols::kyrstr
 starting year of data in the input file
 
integer module_radiation_aerosols::kyrend
 ending year of data in the input file
 
integer module_radiation_aerosols::kyrsav
 the year of data in use in the input file
 
integer module_radiation_aerosols::kmonsav
 the month of data in use in the input file
 
integer, parameter module_radiation_aerosols::nxc = 5
 num of max componets in a profile
 
integer, parameter module_radiation_aerosols::nae = 7
 num of aerosols profile structures
 
integer, parameter module_radiation_aerosols::ndm = 5
 num of atmos aerosols domains
 
integer, parameter module_radiation_aerosols::imxae = 72
 num of lon-points in glb aeros data set
 
integer, parameter module_radiation_aerosols::jmxae = 37
 num of lat-points in glb aeros data set
 
integer, parameter module_radiation_aerosols::naerbnd =61
 num of bands for clim aer data (opac)
 
integer, parameter module_radiation_aerosols::nrhlev =8
 num of rh levels for rh-dep components
 
integer, parameter module_radiation_aerosols::ncm1 = 6
 num of rh independent aeros species
 
integer, parameter module_radiation_aerosols::ncm2 = 4
 num of rh dependent aeros species
 
integer, parameter module_radiation_aerosols::ncm = NCM1+NCM2
 
real(kind=kind_phys), dimension(nrhlev), save module_radiation_aerosols::rhlev
 predefined relative humidity levels
 
real(kind=kind_phys), dimension(ndm, nae), save module_radiation_aerosols::haer
 scale height of aerosols (km)
 
real(kind=kind_phys), dimension(ndm, nae), save module_radiation_aerosols::prsref
 ref pressure lev (sfc to toa) in mb (100Pa)
 
real(kind=kind_phys), dimension(ndm, nae), save module_radiation_aerosols::sigref
 ref sigma lev (sfc to toa)
 
topospheric aerosol profile distribution
real(kind=kind_phys), dimension(nxc, imxae, jmxae), save module_radiation_aerosols::cmixg
 aeros component mixing ratio
 
real(kind=kind_phys), dimension(2,imxae, jmxae), save module_radiation_aerosols::denng
 aeros number density
 
integer, dimension(nxc, imxae, jmxae), save module_radiation_aerosols::idxcg
 aeros component index
 
integer, dimension( imxae, jmxae), save module_radiation_aerosols::kprfg
 aeros profile index
 
module variables for gocart aerosol optical properties
integer, parameter module_radiation_aerosols::kaerbnd =61
 num of bands for aer data (gocart)
 
integer, parameter module_radiation_aerosols::krhlev =36
 num of rh levels for rh-dep components
 
integer, save module_radiation_aerosols::kcm1 = 0
 num of rh indep aerosols (set in subr set_aerspc)
 
integer, save module_radiation_aerosols::kcm2 = 0
 num of rh dep aerosols (set in subr set_aerspc)
 
integer, save module_radiation_aerosols::kcm
 =KCM1+KCM2 (set in subr set_aerspc)
 
real(kind=kind_phys), dimension(krhlev) module_radiation_aerosols::rhlev_grt
 
integer, dimension(:), allocatable module_radiation_aerosols::iendwv_grt
 spectral band structure: ending wavenumber ( \(cm^-1\)) for each band
 
relative humidity independent aerosol optical properties:

species : dust (8 bins)

real(kind=kind_phys), dimension(:,:), allocatable module_radiation_aerosols::rhidext0_grt
 extinction coefficient
 
real(kind=kind_phys), dimension(:,:), allocatable module_radiation_aerosols::rhidssa0_grt
 single scattering albedo
 
real(kind=kind_phys), dimension(:,:), allocatable module_radiation_aerosols::rhidasy0_grt
 asymmetry parameter
 
relative humidity dependent aerosol optical properties:

species : soot, suso, waso, ssam, sscm

real(kind=kind_phys), dimension(:,:,:), allocatable module_radiation_aerosols::rhdpext0_grt
 extinction coefficient
 
real(kind=kind_phys), dimension(:,:,:), allocatable module_radiation_aerosols::rhdpssa0_grt
 single scattering albedo
 
real(kind=kind_phys), dimension(:,:,:), allocatable module_radiation_aerosols::rhdpasy0_grt
 asymmetry parameter
 

module variables for gocart aerosol clim data set

extinction coefficient for SW+LW spectral band

single scattering albedo for SW+LW band asymmetry parameter for SW+LW band

integer, parameter module_radiation_aerosols::imxg = 144
 num of lon-points in geos dataset
 
integer, parameter module_radiation_aerosols::jmxg = 91
 num of lat-points in geos dataset
 
integer, parameter module_radiation_aerosols::kmxg = 30
 num of vertical layers in geos dataset
 
integer, save module_radiation_aerosols::nmxg
 to be determined by set_aerspc
 
real(kind=kind_phys), parameter module_radiation_aerosols::dltx = 360.0 / float(IMXG)
 
real(kind=kind_phys), parameter module_radiation_aerosols::dlty = 180.0 / float(JMXG-1)
 
real(kind=kind_phys), dimension(:,:,:), allocatable, save module_radiation_aerosols::psclmg
 pressure in cb
 
real(kind=kind_phys), dimension(:,:,:,:), allocatable, save module_radiation_aerosols::dmclmg
 aerosol dry mass in g/m3 or aerosol mixing ration in mol/mol or Kg/Kg
 
real(kind=kind_phys), dimension(:), allocatable, save module_radiation_aerosols::geos_rlon
 geos-gocart longitude arrays
 
real(kind=kind_phys), dimension(:), allocatable, save module_radiation_aerosols::geos_rlat
 geos-gocart latitude arrays
 
character *4, save module_radiation_aerosols::gocart_climo = 'xxxx'
 control flag for gocart climo data set: xxxx as default; ver3 for geos3; ver4 for geos4; 0000 for unknown data
 
real(kind=kind_io4), dimension(:), allocatable module_radiation_aerosols::molwgt
 molecular wght of gocart aerosol species
 
logical, save module_radiation_aerosols::lgrtint = .true.
 logical parameter for gocart initialization control
 
logical, save module_radiation_aerosols::lckprnt = .false.
 logical parameter for gocart debug print control
 
real(kind=kind_phys), save module_radiation_aerosols::ctaer = f_zero
 merging coefficients for fcst/clim; determined from fdaer
 
logical, save module_radiation_aerosols::get_fcst = .true.
 option to get fcst gocart aerosol field
 
logical, save module_radiation_aerosols::get_clim = .true.
 option to get clim gocart aerosol field
 
integer, save module_radiation_aerosols::isoot
 index for rh dependent aerosol optical properties (2nd dimension for extrhd_grt, ssarhd_grt, and asyrhd_grt)
 
integer, save module_radiation_aerosols::iwaso
 
integer, save module_radiation_aerosols::isuso
 
integer, save module_radiation_aerosols::issam
 
integer, save module_radiation_aerosols::isscm
 
type(gocart_index_type), save module_radiation_aerosols::dm_indx
 
type(tracer_index_type), save module_radiation_aerosols::dmfcs_indx
 
integer, save module_radiation_aerosols::num_gridcomp = 0
 number of aerosol grid components
 
character, dimension(:), allocatable, save module_radiation_aerosols::gridcomp
 aerosol grid components
 
integer, parameter module_radiation_aerosols::max_num_gridcomp = 5
 default full-package setting
 
character *2, dimension(max_num_gridcomp) module_radiation_aerosols::max_gridcomp
 data max_gridcomp /'DU', 'BC', 'OC', 'SU', 'SS'/
 
integer, dimension(ncm) module_radiation_aerosols::idxspc
 index conversion array:data idxspc / 1, 2, 1, 1, 1, 1, 3, 5, 5, 4 /
 
real(kind=kind_phys), parameter module_radiation_aerosols::wvn550 = 1.0e4/0.55
 the wavenumber ( \(cm^-1\)) of wavelength 550nm for diagnostic aod output
 
integer, save module_radiation_aerosols::nv_aod = 1
 the sw spectral band covering wvn550 (comp in aer_init)
 
subroutine, public module_radiation_aerosols::aer_init
 The initialization program is to set up necessary parameters and working arrays. More...
 
subroutine wrt_aerlog
 This subroutine writes aerosol parameter configuration to run log file.
 
subroutine set_spectrum
 This subroutine defines the one wavenumber solar fluxes based on toa solar spectral distribution, and define the one wavenumber IR fluxes based on black-body emission distribution at a predefined temperature. More...
 
subroutine set_volcaer
 The initialization program for stratospheric volcanic aerosols.
 
subroutine module_radiation_aerosols::clim_aerinit
 This subroutine is the opac-climatology aerosol initialization program to set up necessary parameters and working arrays. More...
 
subroutine set_aercoef
 The initialization program for climatological aerosols. The program reads and maps the pre-tabulated aerosol optical spectral data onto corresponding SW radiation spectral bands. More...
 
subroutine optavg
 This subroutine computes mean aerosols optical properties over each SW radiation spectral band for each of the species components. This program follows GFDL's approach for thick cloud optical property in SW radiation scheme (2000).
 
subroutine, public module_radiation_aerosols::aer_update
 This subroutine checks and updates time varying climatology aerosol data sets. More...
 
subroutine trop_update
 This subroutine updates the monthly global distribution of aerosol profiles in five degree horizontal resolution.
 
subroutine volc_update
 This subroutine searches historical volcanic data sets to find and read in monthly 45-degree lat-zone band of optical depth.
 
subroutine, public module_radiation_aerosols::setaer
 This subroutine computes aerosols optical properties. More...
 

Detailed Description

This file contains climatological atmospheric aerosol schemes for radiation computations.

Definition in file radiation_aerosols.f.

Go to the source code of this file.