Wesley-1989

Wesley (1989) gas-phase dry deposition model. Computes dry deposition velocities v_* [m/s] and rates k_* [1/s] for 132 atmospheric trace gas species. Each species is characterized by its GasData (Dh2oPerDx, Hstar, Fo) constants from Wesely (1989) Table 2 and extended GEOS-Chem species data. The deposition velocity v_i = ifelse(lev==1,1,0) / (Ra + Rb_i + Rc_i) where Ra is the aerodynamic resistance (shared), Rb_i = 5*Sc_i^(2/3)/u_star is the quasi-laminar sublayer resistance, and Rc_i is the Wesley (1989) surface resistance computed inline with species constants substituted. Rain and dew are hardcoded false in the upstream DryDepositionGas builder. Migrated from AtmosphericDeposition.jl@e4468c79 src/dry_deposition.jl:323.

Wesely (1989) bulk surface (canopy) resistance to gaseous dry deposition. Computes the total surface resistance r_c [s/m] along Wesely's three parallel mass-transfer pathways — the upper-canopy leaf cuticles, the lower-canopy structures, and the ground — from the gas species' Henry's-law solubility and reactivity, solar irradiation, surface air temperature, terrain slope, season, land-use category, and surface wetness (rain/dew). Migrated from AtmosphericDeposition.jl@e4468c79 src/wesley1989.jl:360 (the `WesleySurfaceResistance` builder function); implements Wesely (1989) Eq. 2 with the stomatal, mesophyll, cuticle, lower-canopy and ground sub-resistances of Eqs. 3-14 and the Walmsley & Wesely (1996) corrections.