%** plot_vertprof.m
%
%   script to plot vertical profiles of a few variables
%   NOTE: assumes plot_slices has been run, if not load appropriate data
%   set

Rn=86;      % row number of grid point to plot; 65 would be a good south #
Cn=42;      % column number of grid point to plot, near the center of the ice sheet
Cne=28;     % column number of grid point to plot, near the center of the ice sheet
%n=5;        % time step for plotting


%** preliminaries
%   start by calculating vertical coordinates in (x,y,z) from (x,y,zeta)
%   next, calculate pressure-dependent melt temperature at these depths
Tnot=273.15;            % triple point for water, K
g=9.81;                 % acceleration due to gravity
rho=910;                % density of ice, kg/m^3
Phi=8.71e-4 /(rho*g);  %pressure dep. of melt temp, K Pa^-1 (Payne, 1995)


zplot=double(h(n,Rn,Cn))*(1-zeta);
Tpmp=Tnot-910*9.81*(double(h(n,Rn,Cn))-zplot)*Phi;	    % pressure melting point

zplote=double(h(n,Rn,Cne))*(1-zeta);
Tpmpe=Tnot-910*9.81*(double(h(n,Rn,Cne))-zplote)*Phi;	    % pressure melting point

vr=[0 max([max(zplot) max(zplote)])];

%******
% figures
%******

figure(5)
clf
subplot(1,2,1)
plot(T(n,:,Rn,Cne), zplote)
hold on
plot(Tpmpe-Tnot, zplote, 'm-.')
axis([min(T(n,:,Rn,Cne))-1 0 vr])
ylabel('depth (m)')
xlabel('temperature (degrees C)')
title('site near margin')

subplot(1,2,2)
plot(T(n,:,Rn,Cn), zplot)
hold on
plot(Tpmp-Tnot, zplot, 'm-.')
axis([min(T(n,:,Rn,Cn))-1 0 vr])
xlabel('temperature (degrees C)')
title('site near center')


figure(6)
clf
subplot(1,2,1)
speed=sqrt(ux(n,:,Rn,Cne).^2+uy(n,:,Rn,Cne).^2);
plot(speed, zplote)
axis([0 max(speed)+1 vr])
ylabel('depth (m)')
xlabel('speed (m a^-^1)')
title('site near margin')

subplot(1,2,2)
speed=sqrt(ux(n,:,Rn,Cn).^2+uy(n,:,Rn,Cn).^2);
plot(speed, zplot)
axis([0 max(speed)+1 vr])
xlabel('speed (m a^-^1)')
title('site near center')