clear;

h=subplot(1,1,1);set(h,'fontsize',15);

load -ASCII coastline.mat;

load data/katrina_winds.mat;
wind=sqrt(u.^2+v.^2);
[r,t,y]=contourf(xlon,xlat,wind,[0:2:50]);
set(t,'edgecolor','none');
hold on;
[d1 d2]=size(u);
quiver(xlon(1:7:d1,1:7:d2),xlat(1:7:d1,1:7:d2),u(1:7:d1,1:7:d2)/40,v(1:7:d1,1:7:d2)/40,0,'color',[0.7 0.7 0.7],'linewidth',1);
load hurr_dat/katrina_track.dat; 
o_lat=katrina_track(:,2);
o_lon=katrina_track(:,3);
for i=1:27
plot(o_lon(i),o_lat(i),'r*','linewidth',1);
end
plot(coastline(:,1),coastline(:,2),'w','linewidth',1);
X2=[-95 -95];
Y2=[19 19];
U2=[1 0];
V2=[0 0];
h=quiver(X2,Y2,U2,V2,0,'r');set(h,'linewidth',1);
h=text(-93.5,19,'40 m s^{-1}');set(h,'color','r');
title('Katina 10-m winds');
xlabel('Longitude');
ylabel('Latitude');
hold off;

print -djpeg100 exam7
