%read slow campbell
e=['D:\hp_data\P2_12206'];
clear st1;
clear st2;


st1=zeros(130,1);
st2=zeros(130,60);
sig=5.67e-8;

flist=fopen(e,'r');
stuff1=fgets(flist);
jam=0;
for i=1:60
st1=fscanf(flist,'%g%*c%g%*c%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%*c%g%*c%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%*c%g%*c%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%*c%g%*c%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%*c%g%*c%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%*c%g%*c%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%*c%g%*c%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%*c%g%*c%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%g%*c%g%*c%g%g',[130,1]);
                  %1     2   3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 1    2    3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 1     2   3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8  1    2   3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 1     2    3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 1    2   3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8  1    2   3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 1    2    3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 1    2    3 4   
	[nr1,nl1]=size(st1);
	if nr1==130
      st2(:,i)=st1; 
      jam=jam+1;
	else
   	break
	end;
end;
[nn lq]=size(st2);
dex=1:jam;
pir1=st2(58,:)+sig*(st2(42,:)+273.16).^4+4*sig*((st2(44,:)+273.16).^4-(st2(42,:)+273.16).^4);
pir2=st2(60,:)+sig*(st2(46,:)+273.16).^4+4*sig*((st2(48,:)+273.16).^4-(st2(46,:)+273.16).^4);
psp1=st2(62,:);
psp2=st2(64,:);
org=st2(66,:);
Tvais=st2(52,:);
Rhvais=st2(54,:);
Tsea=st2(50,:);


figure ;plot(dex,st2(42,1:jam),dex,st2(44,1:jam),dex,st2(46,1:jam),dex,st2(48,1:jam),dex,Tsea(1:jam),dex,Tvais(1:jam));
figure ;plot(dex,pir1(1:jam),dex,pir2(1:jam),dex,psp1(1:jam),dex,psp2(1:jam));
figure;plot (dex,org(1:jam)); 