%
%  plot some time series data from the geo adjustment
%  calcualtion done by geoadj.
%
%  by Jim Price, 1999.
%  
% 
        
clear all
close all

% a matlab program to set default graphics things

set(0,'DefaultTextFontSize',16)
set(0,'DefaultAxesFontSize',14)
set(0,'DefaultAxesLineWidth',1.4)
set(0,'DefaultLineLineWidth',1.2)

%  to draw a legend w/o a box

%hl = legend('product','source')
%set(hl,'box','off')
%set(hl,'Ycolor', [1 1 1])
%set(hl,'Xcolor', [1 1 1])


load basin10.dat

[m n] = size(basin10);
time = zeros(m,1);
u = time; v = time; h = time; t = time;

time(:) = basin10(:,1);
dt = time(2) - time(1)

u(:) = basin10(:,2);
v(:) = basin10(:,3);
h(:) = basin10(:,4);
t(:) = basin10(:,5);

figure(1)

 
rke = zeros(m,1);
enet = rke; rpe = rke;
             
rke(:) = basin10(:,6);
rpe(:) = basin10(:,8);
enet(:) = basin10(:,9);

subplot(2,1,1)
plot(time,rke,'--', time,rpe,'-.', time, enet) 
xlabel('time, days')
ylabel('energy')
grid 
hl = legend('kinetic','potential','total')
set(hl,'box','off')
set(hl,'Ycolor', [1 1 1])
set(hl,'Xcolor', [1 1 1])

figure(2)
subplot(3,1,1)
plot(time,u, time,v, '-.')
xlabel('time, days')
ylabel('currents, m/sec')
grid

subplot(3,1,2)
plot(time,t)
ylabel('displacement, m')
grid


u(:) = basin10(:,10);
v(:) = basin10(:,11);
t(:) = basin10(:,12);

figure(3)

subplot(3,1,1)
plot(time,u, time,v, '-.')
xlabel('time, days')
ylabel('currents, m/sec')
grid

subplot(3,1,2)
plot(time,t)
ylabel('displacement, m')
grid
 

u(:) = basin10(:,13);
v(:) = basin10(:,14);
t(:) = basin10(:,15);

figure(4)

subplot(3,1,1)
plot(time,u, time,v, '-.')
xlabel('time, days')
ylabel('currents, m/sec')
grid

subplot(3,1,2)
plot(time,t)
ylabel('displacement, m')
grid
 
u(:) = basin10(:,16);
v(:) = basin10(:,17);
t(:) = basin10(:,18);

figure(5)

subplot(3,1,1)
plot(time,u, time,v, '-.')
xlabel('time, days')
ylabel('currents, m/sec')
grid

subplot(3,1,2)
plot(time,t)
ylabel('displacement, m')
grid