function plot_profile(i)
% plot_profile(i); Ledwell; 9/6/05

s=int2str(i);
if length(s)==1
    s=['0',s];
end

eval(['load Tow_',s,'.txt'])
eval(['X=Tow_',s,';'])
height=X(:,1);
A=X(:,6);
PH=Cast6(:,7);

j=find(~isnan(A) & A> 4000 & ~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); clf
scl=1E-05;
plot(scl*A(k),height(k),'--',scl*A,height,'o');
axis([0 3 -50 50])
xlabel('GC Area x 10^{-5}')
text(1.5,-40,['Tow ',s])

