atmos: github.com/ctessum/atmos/seinfeld

package seinfeld

import "github.com/ctessum/atmos/seinfeld"

Constants ¶

const (
Rkcal = 1.9872041e-3 // [kcal K-1 mol-1] Universal gas constant
R     = 0.08205      // [atm L mol-1 K-1] Universal gas constant

)

func DryDepGas¶Uses

func DryDepGas(z, zo, ustar, L, T, rhoA, G, Θ float64,
gd *wesely1989.GasData,
iSeason wesely1989.SeasonCategory,
iLandUse wesely1989.LandUseCategory,
rain, dew, isSO2, isO3 bool) float64

Function DryDepGas calculates dry deposition velocity [m/s] for a gas species, where z is the height of the surface layer [m], zo is roughness length [m], ustar is friction velocity [m/s], L is Monin-Obukhov length [m], T is surface air temperature [K], rhoA is air density [kg/m3] gd is data about the gas species for surface resistance calculations, G is solar irradiation [W m-2], Θ is the slope of the local terrain [radians], iSeason and iLandUse are indexes for the season and land use, dew and rain indicate whether there is dew or rain on the ground, and isSO2 and isO3 indicate whether the gas species of interest is either SO2 or O3, respectively. Based on Seinfeld and Pandis (2006) equation 19.2.

func DryDepParticle¶Uses

func DryDepParticle(z, zo, ustar, L, Dp, T, P, ρParticle, ρAir float64, iSeason SeasonalCategory,
iLandUse LandUseCategory) float64

Function DryDepParticle calculates particle dry deposition velocity [m/s] where z is the height of the surface layer [m], zo is roughness length [m], ustar is friction velocity [m/s], L is Monin-Obukhov length [m], Dp is particle diameter [m], T is surface air temperature [K], P is pressure [Pa], ρParticle is particle density [kg/m3], ρAir is air density [kg/m3], and iSeason and iLandUse are indexes for the season and land use. Based on Seinfeld and Pandis (2006) equation 19.7.

func GasLiquidDistributionFactor¶Uses

func GasLiquidDistributionFactor(H, T, wL float64) float64

Seinfeld and Pandis equation 7.7. Calculates the ratio of gaseous-phase mass concentration to aqueous-phase mass concentration per unit volume of air. Inputs are Henry's law coefficient (H [M atm-1]), temperature (T [K]), cloud/fog liquid water mixing ratio (wL [vol water/vol air]).

func SulfurH2O2aqueousOxidationRate¶Uses

func SulfurH2O2aqueousOxidationRate(H2O2, pH, T, P, wL float64) float64

Seinfeld and Pandis equation 7.84. Calculates the aqueous oxdation rate of SO2 [1/s] by H2O2 when given the H2O2 mixing ratio (H2O2 [ppb]), the water droplet pH (pH [-log(M)]), temperature (T [K]), pressure (P [atm]), and cloud/fog liquid water mixing ratio (wL [vol water/vol air]).

func TemperatureAdjustRate(k298, EperR, T float64) (k float64)

Seinfeld and Pandis equation 7.5. Calculates a temperature adjusted reaction rate given reaction rate at 298K (k298), the heat of dissolution divided by the gas constant (EperR [K]), and the temperature (T [K]). Also works for adjusting Henry's law coefficients

type LandUseCategory¶Uses

type LandUseCategory int
const (
Evergreen LandUseCategory = iota //	0. Evergreen-needleleaf trees
Deciduous                        //	1. Deciduous broadleaf trees
Grass                            //	2. Grass
Desert                           //	3. Desert
Shrubs                           //	4. Shrubs and interrupted woodlands
)

type SeasonalCategory¶Uses

type SeasonalCategory int
const (
Midsummer    SeasonalCategory = iota //	0. Midsummer with lush vegetation
Autumn                               //	1. Autumn with cropland not harvested
LateAutumn                           //	2. Late autumn after frost, no snow
Winter                               //	3. Winter, snow on ground
Transitional                         //	4. transitional
)

Package seinfeld imports 3 packages (graph) and is imported by 1 packages. Updated 2016-07-19. Refresh now. Tools for package owners.