Module gplately.utils.crustal_production

Functions

def compute_crustal_production_rate(agegrid_fn: str, age_threshold: int = 3)

mask out all values greater than the age threshold in the age grid; find the total area which is covered by the remaining age values; divide the area by the age threshold to obtain crustal production rate in km2/year

Parameters

agegrid_fn : str

the file path to the age grid file(NetCDF)

age_threshold : int

use the age threshold to mask age grid and calculate the crustal production rate

Returns

a floating point number(crustal production rate in km2/year)