 
%  % Land mask for the downsampled grid: 
%  land = ~tmd_interp(fn,'mask',Latf,Lonf); 
% 
%  % Anti-aliased (by lowpass filtering) U and V: 
%  tmpu = interp2(Lon,Lat,filt2(U(:,:,1),1,skip*2,'lp'),Lonf,Latf); 
%  tmpv = interp2(Lon,Lat,filt2(V(:,:,1),1,skip*2,'lp'),Lonf,Latf); 
% 
%  % Ensure land pixels are NaN: 
%  tmpu(land) = nan; 
% 
%  % Plot arrows (with a Z component to get above surface topography): 
%  q = quiver3(Lonf,Latf,9*ones(size(Lonf)),tmpu,tmpv,zeros(size(Lonf)),'k'); 
%  q.LineWidth = 0