% Cast 6 is Tow 5

load Tow_05.txt
height=Tow_05(:,1);
A=Tow_05(:,6)/100000;
PH=Tow_05(:,7);

j=find(~isnan(A) & A> 4000/100000 & ~isnan(PH));
figure(2); clf
[fit,par]=try_fit(PH(j),A(j),1);
kp=find(isnan(A) & ~isnan(PH));
A(kp)=polyval(par,PH(kp));
k=find(~isnan(A));

figure(1)
subplot(3,2,5)
plot(A(k),height(k),'--',A,height,'o');
axis([0 3 -50 50])
xlabel('                                       GC Area x 10^{-5}')
text(1.5,-40,'Tow 5')

