firn_air.py

code to deal with modeling firn-air diffusion

class firn_air.FirnAir(air_config, Gs, z, modeltime, Tz, rho, dz, gaschoice, bdot)

__init__(air_config, Gs, z, modeltime, Tz, rho, dz, gaschoice, bdot)

Initialize Firn Air class

diffusivity()

D_0 is CO2 in free air. gam_x is the diffusivity relative to CO2 D_x=D_0*gam_x is the free-air (or any) diffusity for that species diffu_full has units m^2/s

firn_air_diffusion(AirParams, iii)

Solve the diffusion equation. Calls solver.py

porosity()

Calculate the total, open, and closed porosity. co is close-off, cl is closed

firn_air.gasses(gaschoice, T, p_a, M_air)

Function to set up specifics for different gasses.