%% Tungurahua GPS_2_DEM Script
% Written by Carolyn Parcheta and Brando
% June 23, 2009 NMTech Field Course

%% Improvements

% 2. Make and inset map of TINF

%% Note:
% Crater to TBUG = 6.01 km
% Crater to TLEF = 5.95 km
% Crater to TBAG = 5.90 km
% Crater to TMAR = 5.81 km
% Crater to TGRS = 5.92 km
% Crater to TLAV = 5.96 km
% Crater to TREX = 5.77 km
% Crater to TINF = 11.47 km

%%
addpath /Users/geop572/Desktop/GEOP572_MAP_STUFF;
figure(1), clf
hold on
[dem,lats,lons] = readSRTM('TUNGURAHUA',150,150); 

contourf(lons,lats,dem, 30) % without specifying contours MATLAB default is 10
hold on
axis equal
xlabel('LONG')
ylabel('LAT')
title('SEISMIC AND INFRASONIC ARRAYS')
annotation('textarrow', [.85 .85],[.7 .85], 'string', 'N')
colormap(summer)
% hold on

plot(-78.44430, -1.36512,'ro', 'LineWidth',8)
text(-78.5, -1.38,'INFRASOUND ARRAY (TINF)','FontSize',20,'FontWeight','bold', 'Color', 'r')
% StD = 3.22 
% All StD in m

plot(-78.48888, -1.43382,'ro', 'LineWidth',8)
text(-78.55, -1.42,'SEISMIC ARRAY', 'FontSize',20,'FontWeight','bold', 'Color', 'r')
% StD = 73.55 

figure(2), clf
colormap(summer)
[dem,lats,lons] = readSRTM('TUNGURAHUA_2',5,5); 
contourf(lons,lats,dem, 30) % without specifying contours MATLAB default is 10
axis equal
xlabel('LONG')
ylabel('LAT')
title('SEISMIC ARRAY')
hold on


plot(-78.48888, -1.43382,'rx','LineWidth',2)
text(-78.4894, -1.4335,'TBUG','FontSize',13)
% StD = 73.55 

plot(-78.48843, -1.43427,'rx','LineWidth',2)
text(-78.4885, -1.434,'TLEF','FontSize',13)
% StD = 8.82

plot(-78.48664, -1.43430,'rx','LineWidth',2)
text(-78.4872, -1.434,'TREX','FontSize',13)
% StD = 1.92 

plot(-78.48749, -1.43515,'rx','LineWidth',2)
text(-78.488, -1.4355,'TMAR','FontSize',13)
% StD = 5.15

plot(-78.48760, -1.43428,'rx','LineWidth',2)
text(-78.488, -1.4345,'TBAG','FontSize',13)
% StD = 9.40

plot(-78.48765, -1.43337,'rx','LineWidth',2)
text(-78.488, -1.433,'TGRS','FontSize',13)
% StD = 3.28

plot(-78.48742, -1.43264,'rx','LineWidth',2)
text(-78.4875, -1.43225,'TLAV','FontSize',13)
% StD = 2.90