clear all
clf

dfname = ['fort.88_bc'];
fid=fopen(dfname);
fc1 = fread(fid,[119 4320],'float32'); status = fclose(fid);

dfname = ['fort.88_bca'];
fid=fopen(dfname);
fc2 = fread(fid,[119 2880],'float32'); status = fclose(fid);

% im1 = fread(fid,'float32'); status = fclose(fid);
% plot(im1(1,:),im1(50,:));

for i=1:119
   diff(i,:) = fc1(i,1441:4320) - fc2(i,:);
  end

fc2(2,:) = fc2(2,:) + 72;

subplot(2,1,1)

hold on

plot(fc1(2,:),fc1(10,:),'m-')
plot(fc2(2,:),fc2(10,:),'m-')
plot(fc1(2,:),fc1(20,:),'g-')
plot(fc2(2,:),fc2(20,:),'g-')
plot(fc1(2,:),fc1(50,:),'k-')
plot(fc2(2,:),fc2(50,:),'k-')
plot(fc1(2,:),fc1(51,:),'b-')
plot(fc2(2,:),fc2(51,:),'b-')
plot(fc1(2,:),fc1(100,:),'r-')
plot(fc2(2,:),fc2(100,:),'r-')

subplot(2,1,2)

hold on

plot(fc2(2,:),diff(10,:),'m-')
plot(fc2(2,:),diff(20,:),'g-')
plot(fc2(2,:),diff(50,:),'k-')
plot(fc2(2,:),diff(51,:),'b-')
plot(fc2(2,:),diff(100,:),'r-')

% 2880
