mar_extrap_daily.py

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

Source code

General Methods

SMBcorr.extrapolate_mar_daily(DIRECTORY, EPSG, VERSION, tdec, X, Y, VARIABLE='SMB', XNAME=None, YNAME=None, TIMENAME='TIME', SEARCH='BallTree', NN=10, POWER=2.0, SIGMA=1.5, FILL_VALUE=None, EXTRAPOLATE=False)[source]

Spatially extrapolates daily MAR surface mass balance products

Parameters:
DIRECTORY: str

Working data directory

EPSG: str or int

input coordinate reference system

VERSION: str

MAR Version

  • v3.5.2

  • v3.9

  • v3.10

  • v3.11

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’

MAR product to interpolate

  • SMB: Surface Mass Balance

  • PRECIP: Precipitation

  • SNOWFALL: Snowfall

  • RAINFALL: Rainfall

  • RUNOFF: Melt Water Runoff

  • SNOWMELT: Snowmelt

  • REFREEZE: Melt Water Refreeze

  • SUBLIM: Sublimation

XNAME: str or NoneType, default None

Name of the x-coordinate variable

YNAME: str or NoneType, default None

Name of the y-coordinate variable

TIMENAME: str or NoneType, default ‘TIME’

Name of the time variable

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

Directories

  • <path_to_mar>/MARv3.11.2/Greenland/7.5km_ERA5/

  • <path_to_mar>/MARv3.11.2/Greenland/10km_ERA5/

  • <path_to_mar>/MARv3.11.2/Greenland/15km_ERA5/

  • <path_to_mar>/MARv3.11.2/Greenland/20km_ERA5/

  • <path_to_mar>/MARv3.11.2/Greenland/20km_NCEP1/

  • <path_to_mar>/MARv3.11/Greenland/ERA_1958-2019-15km/daily_15km/

  • <path_to_mar>/MARv3.11/Greenland/NCEP1_1948-2020_20km/daily_20km/

  • <path_to_mar>/MARv3.10/Greenland/NCEP1_1948-2019_20km/daily_20km/

  • <path_to_mar>/MARv3.9/Greenland/ERA_1958-2018_10km/daily_10km/