puff (version 0.1.0)

gpuff: Gaussian Puff Concentration Calculation

Description

Calculates the concentration of an emission event at a specified location and time due to a Gaussian puff.

This function uses wind speed and direction components, advection adjustments, and stability class calculations to accurately measure the dispersion of a puff in the atmosphere.

Usage

gpuff(Q, stab_class, x_p, y_p, x_r_vec, y_r_vec, z_r_vec, total_dist, H, U)

Value

Numeric. Pollutant concentration at the specified (x, y, z) locations and time `t`.

Arguments

Q

Numeric. Mass per puff. E.g., for 100 puffs/hour of a 100 kg/hr emission, put 1 kg of mass into each puff.

stab_class

Character vector. Stability class ("A" to "F").

x_p

Numeric. Puff position in the X direction.

y_p

Numeric. Puff position in the Y direction.

x_r_vec

Numeric vector. The x-coordinate (east-west) where the concentration is calculated.

y_r_vec

Numeric vector. The y-coordinate (north-south) where the concentration is calculated.

z_r_vec

Numeric vector. The z-coordinate (height) where the concentration is calculated.

total_dist

Numeric. The total distance the puff has traveled from the source in m.

H

Numeric. Source height.

U

Numeric. Wind speed in m/s.

Examples

out <- gpuff(Q = 1, stab_class = "D", x_p = 0, y_p = 0,
  x_r_vec = 100, y_r_vec = 0, z_r_vec = 2,
  total_dist = 100, H = 2, U = 5
)