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.

## 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.