% 2Dgraph of the DEM 
clf
contour(lons_zoom,lats_zoom,dem_zoom,20); 
%pcolor(dem_zoom);
%shading interp
axis square

%load of the location file for the different stations

positions=load('tungu3.txt');
hold on
plot(positions(1,1),positions(1,2),'b^', 'MarkerFaceColor','b','MarkerSize',10);
text(positions(1,1)-0.001,positions(1,2)-0.0001,'TLAV') 

plot(positions(2,1),positions(2,2),'b^', 'MarkerFaceColor','b','MarkerSize',10);
text(positions(2,1)+0.0002,positions(2,2)-0.0001,'TREX') 

plot(positions(3,1),positions(3,2),'b^', 'MarkerFaceColor','b','MarkerSize',10);
text(positions(3,1)-0.001,positions(3,2)-0.0001,'TMAR') 

plot(positions(4,1),positions(4,2),'b^', 'MarkerFaceColor','b','MarkerSize',10);
text(positions(4,1)-0.001,positions(4,2)-0.0001,'TGRS') 

plot(positions(5,1),positions(5,2),'b^', 'MarkerFaceColor','b','MarkerSize',10);
text(positions(5,1)-0.001,positions(5,2)-0.0001,'TLEF') 

plot(positions(6,1),positions(6,2),'b^', 'MarkerFaceColor','b','MarkerSize',10);
text(positions(6,1)+0.0001,positions(6,2)-0.0001,'TBAC') 

plot(positions(7,1),positions(7,2),'b^', 'MarkerFaceColor','b','MarkerSize',10);
text(positions(7,1)-0.001,positions(7,2)-0.0001,'TINF') 

plot(positions(8,1),positions(8,2),'b^', 'MarkerFaceColor','b','MarkerSize',10);
text(positions(8,1)-0.001,positions(8,2)-0.0001,'TBUG') 

xlabel('longitude');
ylabel('latitude');
title('zoom in Seismic array');