%a=fopen('c:\neaqs\neaqs_sp.dat','r');
%c=0;
%while ~feof(a)
%   b=fgetl(a);
%   c=c+1;
   
%   if c>4490,
%       disp([c length(b)])
%       pause;
%   end
% end
load c:\neaqs\neaqs_sp.dat
SA=neaqs_sp;
clear neaqs_sp;
date=SA(:,1);
fson=SA(:,2:55);
Su=SA(:,56:109);
Sv=SA(:,110:163);
Sw=SA(:,164:217);
St=SA(:,218:271);
Sq=SA(:,272:325);
Suj=SA(:,326:379);
Svj=SA(:,380:433);
Swj=SA(:,434:487);
Cuw=SA(:,488:541);
Cvw=SA(:,542:595);
Ctw=SA(:,596:649);
Cqw=SA(:,650:703);
Cuwj=SA(:,704:757);
Cvwj=SA(:,758:811);
for ii=1:10;figure(ii);loglog(fson(1,:),fson(1,:).*Su(ii,:));title([num2str(ii)]);end;


