mar_extrap_seasonal.py

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

  • Seasonal files are climatology files for each day of the year

Source code

General Methods

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

Spatially extrapolate daily climatologies of MAR 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

RANGE: list

Start and end year of seasonal

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

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/