# 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.