racmo_extrap_daily.py

  • Interpolates and extrapolates daily RACMO products to times and coordinates using inverse distance weighting

Source code

General Methods

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

Spatially extrapolates daily RACMO products

Parameters:
base_dir: str

Working data directory

EPSG: str or int

input coordinate reference system

MODEL: str

Daily model outputs to interpolate

  • FGRN055: 5.5km Greenland 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 ‘smb’

RACMO product to interpolate

  • smb: Surface Mass Balance

  • hgtsrf: Change of Surface Height

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

EXTRAPOLATE: bool, default False

Create a regression model to extrapolate in time