disp('red_ts4')

happ=zeros(12822,3);
hapx=zeros(12000,2);
stux=zeros(8,1);
stuy=zeros(29,1);

fclose('all');
flist=fopen('c:\cwf\cope\jess_b.txt','r');
stuff=fgetl(flist);
%******************  setup read data loop  **********
nx=12822;				%# of lines of data 
jbg=1;
for ibg = 1:nx 			%major read loop
	stux=fscanf(flist,'%g%g%g%g%g%g%g%g',[8,1]);
	%*********************************************
		%*******   decode date  ************
   if stux(1)>10;
      happ(jbg,1)=stux(7);
      happ(jbg,2)=stux(1);
      happ(jbg,3)=stux(2);
      jbg=jbg+1;
   end;
end; %  data line loop
jx=jbg-1
	
fclose('all');
flist=fopen('c:\cwf\cope\jess_s.txt','r');
stuff=fgetl(flist);
nz=12000;				%# of lines of data  
jbg=1;
for ibg = 1:nz 			%major read loop
	stuz=fscanf(flist,'%g%g%g%',[3,1]);
	%*********************************************
		%*******   decode date  ************
	if stuz(3)>0;
      hapx(jbg,1)=stuz(1);
      hapx(jbg,2)=stuz(3);
      jbg=jbg+1;
	end;
end; %  data line loop
jz=jbg-1
fclose('all');
	figure;	
	plot(hapx(1:jz,1), hapx(1:jz,2));
   axis([260 272 16 19]);
   xlabel('Julian Day');
   ylabel('T(C)');
   hold on;
   plot(happ(1:jx,1), happ(1:jx,2), 'r');

  % figure;
  % plot(dtc(1,1:ny), qtot(1:ny), 'o');
   %axis([260 272 0 200]);
   %xlabel('Julian Day');
   %ylabel('Q_{tot} (W/m^2)');
   %figure;	
	%plot(dtc(1,1:ny), hl(1,1:ny), dtc(1,1:ny), hl(2,1:ny));
   %axis([260 272 -10 100]);
   %xlabel('Julian Day');
   %ylabel('Hl (W/m^2)');
   %figure;
   %plot(hl(1,1:ny), hl(2,1:ny), 'o');
   %axis([-10 100 -10 100]);
   %xlabel('Hl_{old} (W/m^2)');
   %ylabel('Hl_{new} (W/m^2)');
	 %figure;	
	%plot(dtc(1,1:ny), tu(1,1:ny), dtc(1,1:ny), tu(2,1:ny));
   %axis([260 272 -0 .6]);
   %xlabel('Julian Day');
   %ylabel('Tau (N/m^2)');
   %figure;
   %plot(tu(1,1:ny), tu(2,1:ny), 'o');
   %axis([-0 .6 -0 .6]);
   %xlabel('Tau_{cov} (W/m^2)');
   %ylabel('Tau_{new} (W/m^2)');

   
fclose('all');
%flist=fopen('c:\matlab5\cwf\progs\tact.txt','w');
%flist=1;
%for ibg = 1:ny 			%major write loop
%fprintf(flist,'%9.5f %6.2f\n',dtc);
%fclose('all');
	
	
	
	
	
