figure

ht=textread('ht.dat','%n');
ht=reshape(ht,17,34);

contourf(ht,10); hold on


axis equal;axis([1 34 1 17]);
title('Initial Height')

u=textread('u.dat','%n');
u=reshape(u,17,34);
v=textread('v.dat','%n');
v=reshape(v,17,34);
quiver(u,v) ; hold off

title('A Plot of the initial Wind Field')


figure
surf(ht);
axis([1 34 1 17 5000 5500]);
view(-25,75)
title('A Surface Plot of the intial height Field')