%** plot_change.m
%
%   script to plot change in a few variables, using slice locations 
%   from plot_slices.m
%   NOTE: assumes plot_slices has been run, if not, plot appropriate data
%   set


%** surface elevation

figure(f+15)
clf
plot(x1(C1n:C2n)/1e3, squeeze(b(1,R1n,C1n:C2n)), 'k-')
hold on
plot(x1(C1n:C2n)/1e3, squeeze(b(1,R1n,C1n:C2n))+squeeze(h(1,R1n,C1n:C2n)), 'm*-')
for m=3:N
    plot(x1(C1n:C2n)/1e3, squeeze(b(m,R1n,C1n:C2n))+squeeze(h(m,R1n,C1n:C2n)), 'c-')
end
plot(x1(C1n:C2n)/1e3, squeeze(b(N+1,R1n,C1n:C2n))+squeeze(h(N+1,R1n,C1n:C2n)), 'b*-')
title('surface elevation over course of model run (N)')
xlabel('distance (km)')
ylabel('elevation (m)')



figure(f+16)
clf
plot(x1(C1s:C2s)/1e3, squeeze(b(1,R1s,C1s:C2s)), 'k-')
hold on
plot(x1(C1s:C2s)/1e3, squeeze(b(1,R1s,C1s:C2s))+squeeze(h(1,R1s,C1s:C2s)), 'm*-')
for m=3:N
    plot(x1(C1s:C2s)/1e3, squeeze(b(m,R1s,C1s:C2s))+squeeze(h(m,R1s,C1s:C2s)), 'c-')
end
plot(x1(C1s:C2s)/1e3, squeeze(b(N+1,R1s,C1s:C2s))+squeeze(h(N+1,R1s,C1s:C2s)), 'b*-')
title('surface elevation over course of model run (S)')
xlabel('distance (km)')
ylabel('elevation (m)')




%** surface speed

figure(f+17)
clf
speed=sqrt(squeeze(ux(2,1,R1n,C1n:C2n)).^2+squeeze(uy(2,1,R1n,C1n:C2n)).^2);
plot(x1(C1n:C2n)/1e3, speed, 'm*-')
hold on
for m=3:N
    speed=sqrt(squeeze(ux(m,1,R1n,C1n:C2n)).^2+squeeze(uy(m,1,R1n,C1n:C2n)).^2);
    plot(x1(C1n:C2n)/1e3, speed, 'c-')
end
speed=sqrt(squeeze(ux(N+1,1,R1n,C1n:C2n)).^2+squeeze(uy(N+1,1,R1n,C1n:C2n)).^2);
plot(x1(C1n:C2n)/1e3, speed, 'b*-')
title('surface speed over course of model run (N)')
axis([x1(C1n)/1e3 x1(C2n)/1e3 0 3e3])
xlabel('distance (km)')
ylabel('speed (m a^-^1)')


figure(f+18)
clf
speed=sqrt(squeeze(ux(2,1,R1s,C1s:C2s)).^2+squeeze(uy(2,1,R1s,C1s:C2s)).^2);
plot(x1(C1s:C2s)/1e3, speed, 'm*-')
hold on
for m=3:N
    speed=sqrt(squeeze(ux(m,1,R1s,C1s:C2s)).^2+squeeze(uy(m,1,R1s,C1s:C2s)).^2);
    plot(x1(C1s:C2s)/1e3, speed, 'c-')
end
speed=sqrt(squeeze(ux(N+1,1,R1s,C1s:C2s)).^2+squeeze(uy(N+1,1,R1s,C1s:C2s)).^2);
plot(x1(C1s:C2s)/1e3, speed, 'b*-')
title('surface speed over course of model run (N)')
axis([x1(C1s)/1e3 x1(C2s)/1e3 0 3e3])
xlabel('distance (km)')
ylabel('speed (m a^-^1)')



%**  temperature at a point
Rn=83;      % row number of grid point to plot; 56 would be a good south #
Cn=42;      % column number of grid point to plot

figure(f+19)
clf
zplot=double(h(1,Rn,Cn))*(1-zeta);
plot(T(1,:,Rn,Cn), zplot, 'm-')
hold on
for m=3:N
    zplot=double(h(m,Rn,Cn))*(1-zeta);
    plot(T(m,:,Rn,Cn), zplot, 'c-')
end
zplot=double(h(N+1,Rn,Cn))*(1-zeta);
plot(T(N+1,:,Rn,Cn), zplot, 'b-')
ylabel('depth, m')
xlabel('degrees C')
title(['temperature over time at row ' int2str(Rn) ', column ' int2str(Cn)])




%** change map

diffh=squeeze(h(N+1,:,:)-h(1,:,:));

figure(f+21)
clf
pcolor(x1/1e3, y1/1e3, double(diffh)), shading flat, axis image
colorbar
title('total thickness change, end - start (m)')
xlabel('km')
ylabel('km')