tang_2003_flux function

Internal function to compute surface co2 flux at a given timepoint via Tang et al 2003