clear all

lon=load('valdez.bathymetry.long');
lat=load('valdez.bathymetry.latg');
HPV_slide=load('valdez.HPV_1964');
SBM_slide=load('valdez.SBM_1964');
bathymetry=load('valdez.1964.bathymetry');

Hotel_lon=-146.26701;
Hotel_lat=61.115823;

NavigationalLight_lon=-146.652749;
NavigationalLight_lat=61.081169;


figure(1)
subplot(2,1,1)
pcolor(lon, lat, HPV_slide); shading flat  %Thickness, meters
hold on
contour(lon, lat, bathymetry,[0 0],'k')    %Bathymetry on which the slide is moving, meters
vh=plot(Hotel_lon, Hotel_lat, '^r','MarkerFaceColor','r');
nl=plot(NavigationalLight_lon, NavigationalLight_lat, '^g','MarkerFaceColor','g');
legend([vh, nl], {'Valdez Hotel', 'Navigational Light'})
hold off
title('Thickness of HPV64 slide, m')

subplot(2,1,2)
pcolor(lon, lat, SBM_slide); shading flat  %Thickness, meters
hold on
contour(lon, lat, bathymetry,[0 0],'k')    %Bathymetry on which the slide is moving, meters
vh=plot(Hotel_lon, Hotel_lat, '^r','MarkerFaceColor','r');
nl=plot(NavigationalLight_lon, NavigationalLight_lat, '^g','MarkerFaceColor','g');
legend([vh, nl], {'Valdez Hotel', 'Navigational Light'})
hold off
title('Thickness of SBM64 slide, m')


InundationLine_1964=shaperead('InundationLine.shp');
InundationLineIslands1964=shaperead('InundationLineIslands.shp');
InundationLineFirstWave=shaperead('InundationLineFirstWave.shp');
McKinley=shaperead('McKinley_St.shp');

figure(2)
[c,h]=contour(lon, lat, bathymetry,[0 0],'k')    %Bathymetry on which the slide is moving, meter
hold on
ic=plot(InundationLine_1964.X, InundationLine_1964.Y,'-y','LineWidth',2) 
for i=1:7
    plot(InundationLineIslands1964(i).X, InundationLineIslands1964(i).Y,'-y','LineWidth',2) 
end
for i=1:3
    fw=plot(InundationLineFirstWave(i).X, InundationLineFirstWave(i).Y,'-r','LineWidth',2);
end
mk=plot(McKinley.X, McKinley.Y,'-k','LineWidth',2); 
vh=plot(Hotel_lon, Hotel_lat, '^r','MarkerFaceColor','r');
hold off
legend([h, ic, fw, mk, vh], {'Shoreline after the HPV64 slide', 'Observed extent of inundation', 'Debris line from the 1st wave', 'McKinley Street','Valdez Hotel'})
axis([-146.31 -146.23 61.105 61.125])
set(gca,'PlotBoxAspectRatio',[1 cosd(61.1) 1])
title('Post-eq bathymetry and extent of inundation')