racmo_interp_firn_height.py

  • Interpolates and extrapolates firn heights to times and coordinates using surface splines

Source code

General Methods

SMBcorr.interpolate_racmo_firn(base_dir, EPSG, MODEL, tdec, X, Y, VARIABLE='zs', SIGMA=1.5, FILL_VALUE=None, REFERENCE=False)[source]

Reads and interpolates downscaled RACMO firn height products

Parameters:
base_dir: str

Working data directory

EPSG: str or int

input coordinate reference system

MODEL: str

RACMO firn model

  • FGRN055: 5.5km Greenland RACMO2.3p2

  • FGRN11: 11km Greenland RACMO2.3p2

  • XANT27: 27km Antarctic RACMO2.3p2

  • ASE055: 5.5km Amundsen Sea Embayment RACMO2.3p2

  • XPEN055: 5.5km Antarctic Peninsula RACMO2.3p2

tdec: float

time coordinates to interpolate in year-decimal

X: float

x-coordinates to interpolate

Y: float

y-coordinates to interpolate

VARIABLE: str, default ‘zs’

RACMO product to interpolate

  • zs: Firn height

  • FirnAir: Firn air content

SIGMA: float, default 1.5

Standard deviation for Gaussian kernel

FILL_VALUE: float or NoneType, default None

Output fill_value for invalid points

Default will use fill values from data file

REFERENCE: bool, default False

Calculate firn variables in reference to first field