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.



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:

     'linestyle','--','color',[.8 .8 .8],'frame','on')
mapzoom('pine island glacier',1200,'inset','se','insetsize',.4)
scarlabel({'Amundsen Sea' 'Pine Island Bay'},...

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:

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

Plot ASAID grounding line and rotate the figure:


Overlay a graticule and add a scalebar.

    'linestyle','--','color',[.2 .2 .2])


Author Info

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