"""
    britter_puff_factory(scenario)

Generates a Britter-McQuaid dispersion model on the given scenario and returns a
function giving the concentration of the form
c(x, y, z, t)

"""
function britter_puff_factory(scenario)
    # TODO
    error("Not Implemented")
end
