module_radiation_aerosols Module Reference

This module contains climatological atmospheric aerosol schemes for radiation computations.

Variables
character(40), parameter	vtagaer ='NCEP-Radiation_aerosols v5.2 Jan 2013 '
integer, parameter, public	nf_aesw = 3
	num of output fields for SW rad
integer, parameter, public	nf_aelw = 3
	num of output fields for LW rad
integer, parameter, public	nlwstr = 1
	starting band number in ir region
integer, parameter, public	nspc = 5
	num of species for output aod (opnl)
integer, parameter, public	nspc1 = NSPC + 1
	total+species
real(kind=kind_phys), parameter	f_zero = 0.0
real(kind=kind_phys), parameter	f_one = 1.0
integer, save	nswbnd = NBDSW
	number of actual bands for sw aerosols; calculated according to laswflg setting
integer, save	nlwbnd = NBDLW
	number of actual bands for lw aerosols; calculated according to lalwflg and lalw1bd settings
integer, save	nswlwbd = NBDSW+NBDLW
	total number of bands for sw+lw aerosols
logical, save	lalwflg = .true.
logical, save	laswflg = .true.
logical, save	lavoflg = .true.
logical, save	lmap_new = .true.
integer, parameter, public	nwvsol = 151
	num of wvnum regions where solar flux is constant
integer, parameter, public	nwvtot = 57600
	total num of wvnum included
integer, parameter, public	nwvtir = 4000
	total num of wvnum in ir range
integer, dimension(nwvsol), save	nwvns0
	number of wavenumbers in each region where the solar flux is constant
real(kind=kind_phys), dimension(nwvsol), save	s0intv
	solar flux \(w/m^2\) in each wvnumb region where it is constant
real(kind=kind_phys), dimension(nbdsw), save	wvn_sw1
real(kind=kind_phys), dimension(nbdsw), save	wvn_sw2
real(kind=kind_phys), dimension(nbdlw), save	wvn_lw1
real(kind=kind_phys), dimension(nbdlw), save	wvn_lw2
integer, parameter	minvyr = 1850
	lower limit (year) data available
integer, parameter	maxvyr = 1999
	upper limit (year) data available
integer, dimension(:,:,:), allocatable, save	ivolae
	monthly, 45-deg lat-zone aerosols data set in subroutine 'aer_init'
integer	kyrstr
	starting year of data in the input file
integer	kyrend
	ending year of data in the input file
integer	kyrsav
	the year of data in use in the input file
integer	kmonsav
	the month of data in use in the input file
integer, parameter	nxc = 5
	num of max componets in a profile
integer, parameter	nae = 7
	num of aerosols profile structures
integer, parameter	ndm = 5
	num of atmos aerosols domains
integer, parameter	imxae = 72
	num of lon-points in glb aeros data set
integer, parameter	jmxae = 37
	num of lat-points in glb aeros data set
integer, parameter	naerbnd =61
	num of bands for clim aer data (opac)
integer, parameter	nrhlev =8
	num of rh levels for rh-dep components
integer, parameter	ncm1 = 6
	num of rh independent aeros species
integer, parameter	ncm2 = 4
	num of rh dependent aeros species
integer, parameter	ncm = NCM1+NCM2
real(kind=kind_phys), dimension(nrhlev), save	rhlev
	predefined relative humidity levels
real(kind=kind_phys), dimension(ndm, nae), save	haer
	scale height of aerosols (km)
real(kind=kind_phys), dimension(ndm, nae), save	prsref
	ref pressure lev (sfc to toa) in mb (100Pa)
real(kind=kind_phys), dimension(ndm, nae), save	sigref
	ref sigma lev (sfc to toa)
real(kind=kind_phys), dimension(:,:), allocatable, save	extrhi
real(kind=kind_phys), dimension(:,:), allocatable, save	scarhi
real(kind=kind_phys), dimension(:,:), allocatable, save	ssarhi
real(kind=kind_phys), dimension(:,:), allocatable, save	asyrhi
real(kind=kind_phys), dimension(:,:,:), allocatable, save	extrhd
real(kind=kind_phys), dimension(:,:,:), allocatable, save	scarhd
real(kind=kind_phys), dimension(:,:,:), allocatable, save	ssarhd
real(kind=kind_phys), dimension(:,:,:), allocatable, save	asyrhd
real(kind=kind_phys), dimension(:), allocatable, save	extstra
topospheric aerosol profile distribution
real(kind=kind_phys), dimension(nxc, imxae, jmxae), save	cmixg
	aeros component mixing ratio
real(kind=kind_phys), dimension(2,imxae, jmxae), save	denng
	aeros number density
integer, dimension(nxc, imxae, jmxae), save	idxcg
	aeros component index
integer, dimension( imxae, jmxae), save	kprfg
	aeros profile index
module variables for gocart aerosol optical properties
integer, parameter	kaerbndd =61
	num of bands for aer data (gocart)
integer, parameter	kaerbndi =56
	num of bands for aer data (gocart)
integer, parameter	krhlev =36
	num of rh levels for rh-dep components
integer, parameter	kcm1 = 5
	num of gocart rh indep aerosols
integer, parameter	kcm2 = 10
	num of gocart rh dep aerosols
integer, parameter	kcm = KCM1 + KCM2
	num of gocart aerosols
real(kind=kind_phys), dimension(:)	rhlev_grt
	num of bands for aer data (gocart)
real(kind=kind_phys), dimension(krhlev)	data
	num of bands for aer data (gocart)
relative humidity independent aerosol optical properties:
species: du001, du002, du003, du004, du005
real(kind=kind_phys), dimension(:,:), allocatable, save	extrhi_grt
real(kind=kind_phys), dimension(:,:), allocatable, save	extrhi_grt_550
real(kind=kind_phys), dimension(:,:), allocatable, save	scarhi_grt
real(kind=kind_phys), dimension(:,:), allocatable, save	ssarhi_grt
real(kind=kind_phys), dimension(:,:), allocatable, save	asyrhi_grt

relative humidity dependent aerosol optical properties:
species : ss001, ss002, ss003, ss004, ss005, so4, bcphobic, bcphilic, ocphobic, ocphilic
real(kind=kind_phys), dimension(:,:,:), allocatable, save	extrhd_grt
	extinction coefficient
real(kind=kind_phys), dimension(:,:,:), allocatable, save	extrhd_grt_550
	extinction coefficient
real(kind=kind_phys), dimension(:,:,:), allocatable, save	scarhd_grt
	extinction coefficient
real(kind=kind_phys), dimension(:,:,:), allocatable, save	ssarhd_grt
	extinction coefficient
real(kind=kind_phys), dimension(:,:,:), allocatable, save	asyrhd_grt
	extinction coefficient
integer, parameter	num_gc = 5
	gocart species
character *2, dimension(num_gc)	gridcomp
	extinction coefficient
integer, dimension(num_gc)	num_radius
	extinction coefficient
integer, dimension(num_gc)	radius_lower
	extinction coefficient
integer, dimension(num_gc)	trc_to_aod
	extinction coefficient
integer, dimension(ncm)	idxspc
	index conversion array:data idxspc / 1, 2, 1, 1, 1, 1, 3, 5, 5, 4 /
real(kind=kind_phys), parameter	wvn550 = 1.0e4/0.55
	the wavenumber ( \(cm^-1\)) of wavelength 550nm for diagnostic aod output
integer, save	nv_aod = 1
	the sw spectral band covering wvn550 (comp in aer_init)
integer	i550
	extinction coefficient
integer	id550
	extinction coefficient
subroutine, public	aer_init (nlay, me, iaermdl, iaerflg, lalw1bd, aeros_file, con_pi, con_t0c, con_c, con_boltz, con_plnk, errflg, errmsg)
subroutine	clim_aerinit (solfwv, eirfwv, me, aeros_file, errflg, errmsg)
	This subroutine is the opac-climatology aerosol initialization program to set up necessary parameters and working arrays.
subroutine, public	aer_update (iyear, imon, me, iaermdl, aeros_file, errflg, errmsg)
	This subroutine checks and updates time varying climatology aerosol data sets.
subroutine, public	setaer (prsi, prsl, prslk, tvly, rhlay, slmsk, tracer, aerfld, xlon, xlat, imax, nlay, nlp1, lsswr, lslwr, iaermdl, iaerflg, top_at_1, con_pi, con_rd, con_g, aerosw, aerolw, aerodp, ext550, errflg, errmsg)
	This subroutine computes aerosols optical properties.
subroutine	aer_property (prsi, prsl, prslk, tvly, rhlay, dz, hz, tracer, alon, alat, slmsk, laersw, laerlw, imax, nlay, nlp1, top_at_1, aerosw, aerolw, aerodp, errflg, errmsg)
	This subroutine maps the 5 degree global climatological aerosol data set onto model grids, and compute aerosol optical properties for SW and LW radiations.
subroutine	gocart_aerinit (solfwv, eirfwv, me, errflg, errmsg)
	This subroutine is the gocart aerosol initialization program to set up necessary parameters and working arrays.
subroutine	aer_property_gocart (prsi, prsl, prslk, tvly, rhlay, dz, hz, tracer, aerfld, alon, alat, slmsk, laersw, laerlw, con_rd, imax, nlay, nlp1, aerosw, aerolw, aerodp, ext550, errflg, errmsg)
	This subroutine compute aerosol optical properties for SW and LW radiations.