clear;

load data/topo;
load -ascii mat_data/coastline.mat;
load data/ike_track.dat;

surf(XLONG,XLAT,HGT/1000,'edgecolor','none'); 
hold on;
plot(ike_track(:,3),ike_track(:,2),'r*');
plot(coastline(:,1),coastline(:,2),'w');
axis([-110,30,1,50,0,10]);
view(25,60);
shading interp;
xlabel('longitude');
ylabel('latitude');
zlabel('Height (km)');
hold off;

