rom Climate Data Toolox

%%
% You may also add column-averaged water velocity vectors like this:

u = tmd_predict('Gr1kmTM_v1.nc',Lat,Lon,t,'u'); 
v = tmd_predict('Gr1kmTM_v1.nc',Lat,Lon,t,'v'); 

hold on
q = quiversc(lon,lat,u,v,'k'); % quiversc is in Climate Data Toolbox 
q.AutoScaleFactor = 1.5; % lengthens the arrows 

%% 5. Accessing tide model data 
% You can use |tmd_interp| to interpolate quantities from the model 
% data file, like ocean masks or constituent coefficients. Here we 
% get the water c