subplot(2,1,1); cla;

%  Set up map axis and land
ax = axesm('mapprojection', 'stereo', 'maplatlimit', [20 90], ...
      'maplonlimit', [0 360], 'origin', [90 0], 'frame', 'off', ...
      'parallellabel', 'off', 'plinelocation', 20);
lm = worldlo('POpatch');
for i = 1:length(lm);
  [lm(i).lat, lm(i).long] = maptrimp(lm(i).lat, lm(i).long, ...
     [20 90], [-1000 1000]);
  lm(i).otherproperty = {'facecolor', 0.7*[1 1 1], 'edgecolor', 0.5*[1 1 1]};
end

%  For polar stereographic plots that wrap around the globe, include the
%  following two lines to get rid of the empty space at lon=0.
lon2 = [lon(:); lon(1)+360];
data2 = [data data(:,1)];

%  Now, contour the results
[cn, hn] = contourm(lat, lon2, data2, [-50:5:-5], '--k', 'linewidth', 1);
[cz, hz] = contourm(lat, lon2, data2, [0 0], '-k', 'linewidth', 2);
[cp, hp] = contourm(lat, lon2, data2, [5:5:50], '-k', 'linewdith', 1);

gridm on;
tightmap