import math
# Given values
Q = 315 # emission rate in g/s
H = 295 * 0.3048 # stack height converted from feet to meters
u = 1.3 # wind speed at 10m in m/s
x = 2200 # distance from stack in meters
σ_y = 200 # dispersion coefficient in y-direction (adjust based on stability condition)
σ_z = 100 # dispersion coefficient in z-direction (adjust based on stability condition)
V_t = 0.1 # settling velocity in m/s (example value, adjust based on particle data)
# Compute ground-level concentration using Eqn 5-41
C_x_0_0 = (Q / (2 * math.pi * u * σ_y * σ_z)) * math.exp(- (H**2) / (2 * σ_z**2))
# Display result
print(f"Ground-level concentration of PM at x = {x} m: {C_x_0_0:.5f} g/m³")