k = 3;

switch k
  case 1 % Use m_map to plot Tasmania
    % tic
    clear all
    close all
    m_proj('lambert','lon',[145 152],'lat',[-45 -40]);
    m_tbase('contourf');
    m_grid('linestyle','none');
    colormap(jet);
    % toc
  case 2 % Use m_map to plot across dateline
    % tic
    clear all
    close all
    m_proj('lambert','lon',[160 200],'lat',[-45 -30]);
    m_tbase('contourf');
    m_grid('linestyle','none');
    colormap(jet);
    % toc
  case 3 % Use m_map to plot pole to pole
    % tic
    clear all
    close all
    m_proj('lambert','lon',[120 180],'lat',[-45 80]);
    m_tbase('contourf');
    m_grid('linestyle','none');
    colormap(jet);
    % toc
end
return
% Check by doing the plot directly

lon1 = 1741;
lon2 = 1825;
lat1 = 541;
lat2 = 601;
lon= getnc('terrainbase', 'lon', lon1, lon2);
lat = getnc('terrainbase', 'lat', lat1, lat2);
h = getnc('terrainbase', 'height', [lat1 lon1], [lat2 lon2]);
[xx, yy] = meshgrid(lon, lat);
figure
[cs, han] = contourf(xx, yy, h);
colorbarf('vert', cs)

