subroutine module_radiation_aerosols::setgocartaer ( real (kind=kind_phys), dimension(:), intent(in)  alon, real (kind=kind_phys), dimension(:), intent(in)  alat, real (kind=kind_phys), dimension(:,:), intent(in)  prslk, real (kind=kind_phys), dimension(:,:), intent(in)  rhlay, real (kind=kind_phys), dimension(:,:), intent(in)  dz, real (kind=kind_phys), dimension(:,:), intent(in)  hz, integer, intent(in)  NSWLWBD, real (kind=kind_phys), dimension(:,:), intent(in)  prsl, real (kind=kind_phys), dimension(:,:), intent(in)  tvly, real (kind=kind_phys), dimension(:,:,:), intent(in)  trcly, integer, intent(in)  IMAX, integer, intent(in)  NLAY, integer, intent(in)  NLP1, integer, intent(in)  ivflip, logical, intent(in)  lsswr, logical, intent(in)  lslwr, real (kind=kind_phys), dimension(:,:,:,:), intent(out)  aerosw, real (kind=kind_phys), dimension(:,:,:,:), intent(out)  aerolw  )

private

Parameters

alon	IMAX, longitude of given points in degree
alat	IMAX, latitude of given points in degree
prslk	(IMAX,NLAY), pressure in cb
rhlay	(IMAX,NLAY), layer mean relative humidity
dz	(IMAX,NLAY), layer thickness in m
hz	(IMAX,NLP1), level high in m
NSWLWBD	total number of sw+ir bands for aeros opt prop
prsl	(IMAX,NLAY), layer mean pressure in mb
tvly	(IMAX,NLAY), layer mean virtual temperature in K
trcly	(IMAX,NLAY,NTRAC), layer mean specific tracer in g/g
IMAX	horizontal dimension of arrays
NLAY,NLP1	vertical dimensions of arrays
ivflip	control flag for direction of vertical index =0: index from toa to surface =1: index from surface to toa
lsswr,lslwr	logical flag for sw/lw radiation calls
aerosw	(IMAX,NLAY,NBDSW,NF_AESW), aeros opt properties for SW (:,:,:,1): optical depth (:,:,:,2): single scattering albedo (:,:,:,3): asymmetry parameter
aerolw	(IMAX,NLAY,NBDLW,NF_AELW), aeros opt properties for LW (:,:,:,1): optical depth (:,:,:,2): single scattering albedo (:,:,:,3): asymmetry parameter

setgocartaer General Algorithm

1. Call map_aermr() to map input tracer array (trcly) to local tracer array (aermr).
2. Map geos-gocart climo (dmclmg) to gfs grids (dmclm).
3. Call aeropt_grt() to alculate sw/lw aerosol optical properties for the corresponding frequency bands.

References aeropt_grt(), physcons::con_amd, physcons::con_rd, ctaer, dmclmg, f_one, f_zero, funcphys::fpkap(), geos_rlat, geos_rlon, get_clim, get_fcst, gocart_climo, imxg, jmxg, kmxg, lalwflg, laswflg, map_aermr(), molwgt, module_radlw_parameters::nbdlw, module_radsw_parameters::nbdsw, nlwbnd, nmxg, and psclmg.

Referenced by setaer().

Here is the call graph for this function: