function h = MyMap()
% MyMap is a template for your own function
%
%

disp('Getting Matlab''s intrinsic topo Data')

load topo

map     = topo';

map(find(topo > 1.E35)) = NaN;

datmin = min(map(:));
datmax = max(map(:));

%----------------------------------------------------------------------
% Possible projections
%
% Equidistant Azimuthal         eqdazim
% Gnomonic                      gnomonic
% Lambert Azimuthal Equal-Area  eqaazim
% Othographic                   ortho
% Stereographic                 stereo
% Globe                         globe           (spherical)
%----------------------------------------------------------------------

clf; orient portrait;

load coast;                                     % GET MATLAB COAST DATA
[linelat,linelon] = maptriml(lat,long,[0 90],[-180 180]); % JUST NH LAND

elats   = [ -90 : 1 :  90 ];
elons   = [ 0.0 : 1 : 360 ]';
elatmat = ones(length(elons),1)*elats;  % Matrix of latitudes  of grid edges
elonmat = elons*ones(1,length(elats));  % Matrix of longitudes of grid edges

%  zmax  = 3000;
%  map   = thresh(map, zmax, 1);           % Replace anything above
%  map   = thresh(map,-zmax,-1);           % Replace anything below

h1 = axesm('MapProjection','stereo',      ...
         'Origin',[90 0 0],               ...
         'Grid','on',                     ...
         'MapLatLimit',[0 90],            ...
         'MapLonLimit',[-179.9 180],      ...
         'MLineLimit',[0 75],             ...
         'MLineException',[-90 0 90 180], ...
         'MLabelParallel',[0],            ...
         'Meridianlabel','on',            ...
         'Frame','Off');

axis off;                                    % Turns off "box"

plotm(linelat,linelon,'k');                  % OVERLAY LANDMASSES
s = surfacem(elatmat,elonmat,map);           % DRAW MAP
set(s,'facecolor','flat');                   % HELPS NaNs appear clear
title({'Elevation';'meters'},'FontSize',14);   % TITLE
colorbar

shading flat;
%  print -dpsc2 -painters mapdemo.ps;           % BEST WAY TO PRINT
