% script to plot blackcanyon city landslide displacements
%"For Plotting, BC ls data, 1998 to 2000"						
%Benchmark Point number	E2000ROT	dE (1998-2000ROT)	N2000ROT	dN (1998-2000ROT)	H2000fixed	dH (1998-2000fixed)
%1	780.438148	-0.002147972	1091.50718	0.019820466	973.362	0.011
%2	755.1294686	-0.015468643	1156.125406	0.021593938	973.328	-0.015
%5	770.2205784	0.036421594	1197.486637	0.023363225	981.349	-0.003
%8	760.5567622	0.040237779	1229.7945	0.021500192	993.823	-0.011
%15						
%25	703.327382	0.018617984	1270.176421	0.036579357	1001.186	-0.001
%13						
%12	806.7305485	-0.001548464	2023.055888	-0.007887869	1168.043	0.004
%18	837.5026078	0.073392199	1732.651135	-0.050134531	1109.789	0.001
%27	673.7375789	0.054421059	1417.665269	0.038730673	1024.79	0.016
%28	729.602769	0.051230961	1424.964101	0.037898963	1026.885	0.016
%10	770.2278695	0.03313052	1415.716421	0.01057949	1026.952	0.011
%19	807.6616625	0.031337539	1408.939169	0.052831373	1026.782	0.006
%5cm   450 		0	2000 0.05	0	0  For a scale!

load bcdisplacements.txt

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

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('1998 to 2000 displacements:  Black Canyon City Landslide')

print -depsc bcdisplacements.eps