% interp_profiles.m

zr=[-50:1:50]';
m=length(zr);
n=11;

load cez
zmat(:,9:n)=zmat(:,9:n)+10;

Cint=NaN*ones(m,n);
Eint=cint;

for i=1:n
    k=find(~isnan(Cmat(:,i)));
    Cint(:,i)=interp1(zmat(k,i),Cmat(k,i),zr);
    Eint(:,i)=interp1(zmat(k,i),Emat(k,i),zr);
end
    
[J,N]=sum_nan(Cint');

[Cbar,Ebar,nbar]=mean_r(Cint,Eint,J);

figure(1); clf
plot(Cbar,zr)
ke=5:10:95;
xerrbar(Cbar(ke),zr(ke),Ebar(ke));
ylabel('z  (m)')
xlabel('Mean C  (fM)')
text(10,-35,'Mean Profile')
