(:,:,1)); % the first "slice" of Z data
%  hold on
%  he = earthimage('bottom');
%  shading interp
%  view(2)
%  axis image off
%  cmocean balance
%  caxis([-1 1]*2)
%  shadem(-12) % tiny bit of hillshade
%  hs.ZData = hs.ZData-min(hs.ZData(:))+1; % raises above earthimage 
%  set(gca,'position',[0 0 1 1])
%  ax = axis; % gets axis limits 
%  
%  % Downsampled grid: 
%  skip = 15; 
%  Lonf = Lon(floor(skip/2):skip:end,floor(skip/2):skip:end); 
%  Latf = Lat(floor(skip/2):skip:end,floor(skip/2):skip:end); 
%