bove earthimage
set(gca,'position',[0 0 1 1])
ax = axis; % gets axis limits</pre><pre>% Downsampled grid:
skip = 15;
Lonf = Lon(floor(skip/2):skip:end,floor(skip/2):skip:end);
Latf = Lat(floor(skip/2):skip:end,floor(skip/2):skip:end);</pre><pre>% Land mask for the downsampled grid:
land = ~tmd_interp(fn,'mask',Latf,Lonf);</pre><pre>% 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);</