tion',[10 10 560 280],'color','k')
%  hs = surf(Lon,Lat,-wct,Z(:,:,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])
%  
%  txt = text(0.5,0.5,{'Tide Model Driver';'for MATLAB'},'horiz','center',...
%     'vert','middle','units','normalized','fontweight','bol