% Plotting scalar and vector field

% load sample data

load wfs_data.octave

SST(SST==0) = NaN;

% enable overplotting

hold on

% draw temperature

pcolor(lon,lat,SST); colorbar

% draw velocity

quiver(lon2,lat2,u,v,'k')

% set title

title('Surface temperature and velocity')

% draw coastline

plot_coastline('gshhs_l.octave',[-90.519  -80.397],[24.271   30.816])

hold off

% change size of figure (in pixels) to get a correct 
% aspect ratio for print

set(gcf,'PaperPosition',[0 0 600 350]);