racmo_extrap_firn_height.py

  • Interpolates and extrapolates firn heights to times and coordinates using inverse distance weighting

Source code

General Methods

SMBcorr.extrapolate_racmo_firn(base_dir, EPSG, MODEL, tdec, X, Y, VARIABLE='zs', SEARCH='BallTree', NN=10, POWER=2.0, SIGMA=1.5, FILL_VALUE=None, REFERENCE=False)[source]

Spatially extrapolates RACMO firn heights

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

SEARCH: str, default ‘BallTree’

nearest-neighbor search algorithm

NN: int, default 10

number of nearest-neighbor points to use

POWER: int or float, default 2.0

Inverse distance weighting power

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