subroutine, public rrtmg_lw::rrtmg_lw_run	(	real (kind=kind_phys), dimension(npts,nlay), intent(in)	plyr,
		real (kind=kind_phys), dimension(npts,nlp1), intent(in)	plvl,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	tlyr,
		real (kind=kind_phys), dimension(npts,nlp1), intent(in)	tlvl,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	qlyr,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	olyr,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	gasvmr_co2,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	gasvmr_n2o,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	gasvmr_ch4,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	gasvmr_o2,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	gasvmr_co,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	gasvmr_cfc11,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	gasvmr_cfc12,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	gasvmr_cfc22,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	gasvmr_ccl4,
		integer, dimension(npts), intent(in)	icseed,
		real (kind=kind_phys), dimension(npts,nlay,nbands), intent(in)	aeraod,
		real (kind=kind_phys), dimension(npts,nlay,nbands), intent(in)	aerssa,
		real (kind=kind_phys), dimension(npts), intent(in)	sfemis,
		real (kind=kind_phys), dimension(npts), intent(in)	sfgtmp,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	dzlyr,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	delpin,
		real (kind=kind_phys), dimension(npts), intent(in)	de_lgth,
		integer, intent(in)	npts,
		integer, intent(in)	nlay,
		integer, intent(in)	nlp1,
		logical, intent(in)	lprnt,
		real (kind=kind_phys), dimension(npts,nlay), intent(in)	cld_cf,
		logical, intent(in)	lslwr,
		real (kind=kind_phys), dimension(npts,nlay), intent(inout)	hlwc,
		type (topflw_type), dimension(npts), intent(inout)	topflx,
		type (sfcflw_type), dimension(npts), intent(inout)	sfcflx,
		real (kind=kind_phys), dimension(npts,nlay), intent(inout)	cldtau,
		real (kind=kind_phys), dimension(npts,nlay), intent(inout), optional	HLW0,
		real (kind=kind_phys), dimension(npts,nlay,nbands), intent(inout), optional	HLWB,
		type (proflw_type), dimension(npts,nlp1), intent(inout), optional	FLXPRF,
		real (kind=kind_phys), dimension(npts,nlay), intent(in), optional	cld_lwp,
		real (kind=kind_phys), dimension(npts,nlay), intent(in), optional	cld_ref_liq,
		real (kind=kind_phys), dimension(npts,nlay), intent(in), optional	cld_iwp,
		real (kind=kind_phys), dimension(npts,nlay), intent(in), optional	cld_ref_ice,
		real (kind=kind_phys), dimension(npts,nlay), intent(in), optional	cld_rwp,
		real (kind=kind_phys), dimension(npts,nlay), intent(in), optional	cld_ref_rain,
		real (kind=kind_phys), dimension(npts,nlay), intent(in), optional	cld_swp,
		real (kind=kind_phys), dimension(npts,nlay), intent(in), optional	cld_ref_snow,
		real (kind=kind_phys), dimension(npts,nlay), intent(in), optional	cld_od,
		character(len=*), intent(out)	errmsg,
		integer, intent(out)	errflg
	)

1. Change random number seed value for each radiation invocation (isubclw =1 or 2).
2. Read surface emissivity.
3. Prepare atmospheric profile for use in rrtm.
4. Set absorber amount for h2o, co2, and o3.
5. Set up column amount for rare gases n2o,ch4,o2,co,ccl4,cf11,cf12, cf22, convert from volume mixing ratio to molec/cm2 based on coldry (scaled to 1.0e-20).
6. Set aerosol optical properties.
7. Read cloud optical properties.
8. Compute precipitable water vapor for diffusivity angle adjustments.
9. Compute column amount for broadening gases.
10. Compute diffusivity angle adjustments.
11. For cloudy atmosphere, call cldprop() to set cloud optical properties.
12. Calling setcoef() to compute various coefficients needed in radiative transfer calculations.
13. Call taumol() to calculte the gaseous optical depths and Plank fractions for each longwave spectral band.
14. Call the radiative transfer routine based on cloud scheme selection. Compute the upward/downward radiative fluxes, and heating rates for both clear or cloudy atmosphere.
    - call rtrn(): clouds are assumed as randomly overlaping in a vertical column
    - call rtrnmr(): clouds are assumed as in maximum-randomly overlaping in a vertical column;
    - call rtrnmc(): clouds are treated with the mcica stochastic approach.
15. Save outputs.

References a0, a1, a2, amdo3, amdw, cldprop(), physcons::con_amd, physcons::con_amw, physcons::con_avgd, physcons::con_g, eps, f_one, f_zero, physparam::ilwcliq, physparam::ilwrgas, physparam::iovrlw, ipsdlw0, physparam::isubclw, physparam::ivflip, lflxprf, lhlw0, lhlwb, module_radlw_parameters::maxgas, module_radlw_parameters::nbands, rtrn(), rtrnmc(), rtrnmr(), semiss0, setcoef(), and taumol().

Here is the call graph for this function: