hirano_flux function

Internal function to compute surface co2 flux at a given timepoint via Hirano et al 2005