clear;
load data/topo;
load -ascii mat_data/coastline.mat;
load data/ike_track.dat; 
h=surf(XLONG,XLAT,HGT,'edgecolor','none'); 
hold on;
plot(ike_track(:,3),ike_track(:,2),'r*');
plot(coastline(:,1),coastline(:,2),'w');
axis([-110,30,1,50,0,10000]); 
%daspect([0.1 0.1 1])
view(25,60);
shading interp;
mycolor=[0.0:0.05:1;0:0.05:1;1:-0.05:0]';
colormap(mycolor);
camlight right; 
camproj perspective;

lighting flat
%colormap('gray');
hold off;
