MATLAB Examples

albmap_interp documentation

albmap_interp loads data from Le Brocq et al.'s ALBMAP datasets.

See also albmap_data.



zi = albmap_interp(variable,lati,loni)
zi = albmap_interp(variable,xi,yi)
zi = albmap_interp(...,'method',InterpMethod)


zi = albmap_interp(variable,lati,loni) gives the value(s) of a variable at the specified geolocations lati,loni. The variable can be:

  • 'topg' (m) bedrock topography
  • 'lsrf' (m) lower ice surface elevation
  • 'usrf' (m) 'upper ice surface elevation
  • 'mask' grounded/floating mask
  • 'firn' (m) firn correction
  • 'temp' (deg C) surface air temperature
  • 'accr' (m) accumulation - Van de berg et al
  • 'acca' (m) accumulation - Arthern et al
  • 'ghffm' (mW / m^2) geothermal heat flux - Fox Maule et al
  • 'ghfsr' (mW / m^2) geothermal heat flux - Shapiro & Ritzwoller
  • 'mask_plus' mask plus ice free areas & ice streams
  • 'glmask' grounding line uncertainty mask
  • 'bmask' bed/bathymetry modification mask
  • 'umask' thickness/surface modification mask
  • 'topg2' bedrock topography (incl. curvebed)
  • 'lsrf2' lower ice surface elevation (incl. curvebed)

zi = albmap_interp(variable,xi,yi) as above, but for ps71 coordinates xi,yi (meters).

zi = albmap_interp(...,'method',InterpMethod) specifies an interpolation method. Default is 'linear'.


Here's a grid centered on Fimbul Ice Shelf:

% 650x400 km wide, at 1 km spacing:
[lat,lon] = psgrid('fimbul ice shelf',[650 400],1);

The ice thickness at Fimbul is the difference between the upper and lower surfaces:

usrf = albmap_interp('usrf',lat,lon);
lsrf = albmap_interp('lsrf',lat,lon);

thickness = usrf-lsrf;

% Make ice transparent where thickness is zero:
thickness(thickness==0) = nan;


axis tight off

mapzoomps('nw','frame','off') % places an inset map.


If you use ALBMAP data, please cite the Le Brocq paper listed below. And if this function is useful for you, please do me a kindness and cite my Antarctic Mapping Tools paper.

Le Brocq, Anne M; Payne, Antony J; Vieli, Andreas (2010): An improved Antarctic dataset for high resolution numerical ice sheet models (ALBMAP v1). Earth System Science Data, 2(2), 247-260,

Le Brocq, Anne M; Payne, Antony J; Vieli, Andreas (2010): Antarctic dataset in NetCDF format. PANGAEA,

Greene, C. A., Gwyther, D. E., & Blankenship, D. D. Antarctic Mapping Tools for Matlab. Computers & Geosciences. 104 (2017) pp.151-157.

Author Info

This function and supporting documentation were written by Chad A. Greene, October 2018.