
%  
%  % Create 0.2 degree resolution grid: 
%  [Lat,Lon] = cdtgrid(0.2); 
%  
%  % Create an hourly time array for 25 hours: 
%  t = datetime('mar 29, 2017'):hours(1):datetime('mar 30, 2017'); 
%  
%  % Get water column thickness for the grid: 
%  wct = tmd_interp(fn,'wct',Lat,Lon); 
%  
%  % Predict tides: 
%  Z = tmd_predict(fn,Lat,Lon,t,'h'); 
%  U = tmd_predict(fn,Lat,Lon,t,'U'); 
%  V = tmd_predict(fn,Lat,Lon,t,'V'); 
% 
% 
%  figure('position',[10 10 560 280],'color','k')
%  hs = surf(Lon,Lat,-wct,Z