ly time array for 25 hours:
t = datetime('mar 29, 2017'):hours(1):datetime('mar 30, 2017');</pre><pre>% Get water column thickness for the grid:
wct = tmd_interp(fn,'wct',Lat,Lon);</pre><pre>% Predict tides:
Z = tmd_predict(fn,Lat,Lon,t,'h');</pre><h2 id="2">Animate</h2><pre>figure('position',[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 b