clear all
clf

subplot(2,1,1)

load fort.87
load fc1.87_bc
fort(:,1) = fort(:,1) + 72;
% fc1(:,1) = fc1(:,1) + 72;
hold on
hp = plot(fc1(:,1),fc1(:,2),'k-');
% set(hp,'LineWidth',6.0)
hp = plot(fort(:,1),fort(:,2),'r-');
% set(hp,'LineWidth',1.0)

subplot(2,1,2)

diff1 = fc1(1441:4320,2) - fort(:,2);
hp = plot(fort(:,1),diff1,'k-');

