% REFSCRIPT2 - view stats between Ed490 and OC-100 wavelngths from satlatic surface references

reffile = 'refstats2.txt';
fid = fopen(reffile,'wt');
fprintf(fid,'%%mean(Ed412) std(Ed412) mean(Ed443) std(Ed443) mean(Ed490) mean(Ed490) mean(Ed510) std(Ed510) mean(Ed555) std(Ed555) mean(Ed670) std(Ed670) mean(Ed683) mean(Ed670)\n');
f = getfname('H:\n98\gp3-98-ka\SatlanticProfiler\*ref');
for i = 1:15
   x = qreadref(['H:\n98\gp3-98-ka\SatlanticProfiler\' f(i,:)]);
   [r c] = size(x)
   m = mean(x);
   s = std(x)./r;
   for j = 1:c
      fprintf(fid,'%12.7f%12.7f', m(j), s(j));
   end
   fprintf(fid,'\n');
end
fclose(fid);
fprintf(1,'Data written to file: %s\n\n',reffile);

x = load(reffile);
Color = 'bcgymr';
%Color = [0 0 1; 0 1 1; 0 1 0; 1 1 0; 1 0 1; 1 0 0;];
Ed = ['Ed412';'Ed443';'Ed510';'Ed555';'Ed670';'Ed683'];
n = 0;
for i = [1 3 7 9 11 13]
   n = n + 1;
   h = errorbar(x(:,5),x(:,i),x(:,i+1),'o');
   set(h,'Color',Color(n))
   if i == 1
      xlabel('Ed490')
      ylabel('Ed \itx')
      hold('on')
      YLim = get(gca,'YLim');
   end
   [a r] = regress_linear(x(:,5),x(:,i));
   plot(x(:,5),(a(1).*x(:,5)+a(2)),'Color',Color(n))

h = text(5,180 - i*3,[Ed(n,:) ' = ' num2str(a(1)) '*Ed490 + ' num2str(a(2)) ' r^2 = ' num2str(r)]);
set(h,'Color',Color(n));
end