MATLAB Examples

Stanton2013Channelized

This example describes how to create an image similar to Figure 1 of Stanton et al.'s 2013 Science paper. The paper is titled Channelized Ice Melting in the Ocean Boundary Layer Beneath Pine Island Glacier, Antarctica.

Contents

Requirements

This example requires the following:

Figure 1a

Start by plotting ice shelves and grounded ice:

figure('pos',[100 100 600 500])
bedmap2 patchshelves
bedmap2('patchgl','facecolor',.9*[1 1 1],'rotation',90)

Add a graticule, zoom in, and label:

antmap('graticule','lats',-77.5:2.5:-70,'lons',-120:2.5:-80,...
     'linestyle','--','color',[.8 .8 .8],'frame','on')
mapzoom('pine island glacier',1200,'inset','se','insetsize',.4)
scarlabel({'Amundsen Sea' 'Pine Island Bay'},...
    'fontweight','bold','fontsize',15)

Figure 1b

Start by initializing a figure roughly the same size as in Stanton et al., then :

figure('position',[100 100 800 700])

% Create background image:
modismoa('pine island glacier',130,'contrast','lc')

Matlab only understands one colormap per figure, and we've already plotted a grayscale colormap with modismoa, so how do we overlay colordata? Fortunately, John Iversen has solved this problem with freezeColors. All we do is type freezeColors and overlay a semitransparent measures colormap:

freezeColors
measures('speed','pine island glacier','mapwidth',130,...
    'colorbar','off','alpha',.2,...
    'colormap',rgbmap('dark blue','medium blue','canary','light red',255))

Plot ASAID grounding line and rotate the figure:

asaid('gl','color','k','linewidth',3);

Overlay a graticule and add a scalebar.

antmap('graticule','lats',-78:.2:-73,'lons',-105:-95,...
    'linestyle','--','color',[.2 .2 .2])

scalebar('length',20,'location','se')

Author Info

Chad Greene wrote this example file. Chad claims absolutely no credit for the data, science, or creative display shown above.