Internal function to compute surface co2 flux at a given timepoint via De Jong and Schappert (1972)
Useful links