% script to plot blackcanyon city landslide displacements
%"For Plotting, BC ls data, 2000 to 2002"						
%	Benchmark Point number	2002 survey number	E2002	de (2000-2002)	N2002	dn (2000-2002)	H2002	dH (2000-2002)
%	1	1000	777.176	0.01	1083.28	-0.004	973.328	-0.013
%	2	1001	749.496	-0.005	1146.919	-0.011	973.304	-0.023
%	5	1002	763.049	-0.016	1188.812	-0.01	981.293	0.009
%	8	1003	752.193	-0.019	1220.75	-0.022	993.769	0.007
%	25	1004	693.471	0.011	1258.961	-0.008	1001.134	0.005
%	12	1009	768.876	-0.07	2015.189	-0.031	1167.954	0.042
%	27	1008	658.399	0.027	1405.225	0.014	1024.745	-0.002
%	28	1007	713.982	-0.001	1414.616	-0.005	1026.848	-0.01
%	10	1006	754.946	-0.024	1406.879	0.002	1026.907	-0.002
%	19	1005	792.594	-0.012	1401.524	-0.023	1026.734	0.001
%5cm   450 		0	2000 0.05	0	0  For a scale!

load ohtwodisp.txt

bm = ohtwodisp(:,1);
e = ohtwodisp(:,3);
de = ohtwodisp(:,4);
n = ohtwodisp(:,5);
dn = ohtwodisp(:,6);
h = ohtwodisp(:,7);
dh = ohtwodisp(:,8);

figure(1)
clf
plot(e,n, 'k+')
axis equal

for i = 1:length(bm)
   s = sprintf('%.f', bm(i));
   text(e(i)+20,n(i)+20,s);
   drawnow; 
end
hold on
h = QUIVER(e,n,-de,-dn);
%legend(h, 'vector', 2)

axis([400 1000 1000 2100])
text(500,2020,'cm')
xlabel('Easting (m)')
ylabel('Northing (m)')
title('2000 to 2002 displacements:  Black Canyon City Landslide')

print -depsc ohtwodisp.eps