 
lat = linspace(75.6,71.9,length(t));
lon = linspace(-74.8,-57.4,length(t));

z = tmd_predict('Gr1kmTM_v1.nc',lat,lon,t); 

figure
geoscatter(lat,lon,20,z,'filled')
geolimits([68 78],[-85 -45])
geobasemap colorterrain
caxis([-1 1]*1.1)
cb = colorbar; 
ylabel(cb,'tide height (m)')
cmocean balance % optional colormap 

%% Map-view snapshot
% Suppose you just want a snapshot of the geographic distribution of 
% tide heights across a wide area at a single instant. Let's pick 4:15 pm
% on January 12, 2020, and 