load d:\temp\gasex\fson.prn
load d:\temp\gasex\fson_ls.prn
Bdate=SA(:,1);
BSu=SA(:,2:36);
BSv=SA(:,37:71);
BSw=SA(:,72:106);
BSt=SA(:,107:141);
BSq=SA(:,142:176);
BSoh=SA(:,177:211);
BSoc=SA(:,212:246);
BSlh=SA(:,247:281);
BSlc=SA(:,282:316);
BSsc=SA(:,317:339);
BCuw=SA(:,340:374);
BCvw=SA(:,375:409);
BCtw=SA(:,410:444);
BCqw=SA(:,445:479);
BCohw=SA(:,480:514);
BCocw=SA(:,515:549);
BClhw=SA(:,550:584);
BClcw=SA(:,585:619);
Adate=SB(:,1);
ASu=SB(:,2:36);
ASv=SB(:,37:71);
ASw=SB(:,72:106);
ASt=SB(:,107:141);
ASq=SB(:,142:176);
ASoh=SB(:,177:211);
ASoc=SB(:,212:246);
ASlh=SB(:,247:287);
ASlc=SB(:,288:328);
ASsc=SB(:,329:351);
ACuw=SB(:,352:386);
ACvw=SB(:,387:421);
ACtw=SB(:,422:456);
ACqw=SB(:,457:491);
ACohw=SB(:,492:526);
ACocw=SB(:,527:561);
AClhw=SB(:,562:602);
AClcw=SB(:,603:643);

%for ii=1:length(ASoc);figure(1);loglog(fson(1:35),fson(1:35).*ASoc(ii,1:35),fson(11:35),((fson(11:35).^(-2/3)))./(fson(19)^(-2/3)).*fson(19).*ASoc(ii,19));title([num2str(ii)]);pause;end;
for ii=1:length(BSoc);if ii>935,figure(1);loglog(fson(1:35),fson(1:35).*BSoc(ii,1:35),fson(2:35),fson(2:35).*BSlc(ii,2:35),fson(11:35),((fson(11:35).^(-2/3)))./(fson(19)^(-2/3)).*fson(19).*BSoc(ii,19));title([num2str(ii)]);pause;end;end;
%for ii=1:length(ASoh);figure(1);loglog(fson(1:35),fson(1:35).*ASoh(ii,1:35),fson(1:35),fson(1:35).*ASlh(ii,1:35),fson(1:35),fson(1:35).*ASq(ii,1:35),fson(11:35),((fson(11:35).^(-2/3)))./(fson(19)^(-2/3)).*fson(19).*ASoh(ii,19));title([num2str(ii)]);pause;end;
%for ii=1:length(BSoh);figure(1);loglog(fson(1:35),fson(1:35).*BSoh(ii,1:35),fson(1:35),fson(1:35).*BSlh(ii,1:35),fson(1:35),fson(1:35).*BSq(ii,1:35),fson(11:35),((fson(11:35).^(-2/3)))./(fson(19)^(-2/3)).*fson(19).*BSoh(ii,19));title([num2str(ii)]);pause;end;
%for ii=1:length(BSoh);if ii>2122,figure(1);loglog(fson(1:35),fson(1:35).*BSoh(ii,1:35),fson(2:35),fson(2:35).*BSlh(ii,2:35),fson(2:35),fson(2:35).*BSq(ii,2:35),fson(11:35),((fson(11:35).^(-2/3)))./(fson(19)^(-2/3)).*fson(19).*BSoh(ii,19));title([num2str(ii)]);pause;end;end;
%for ii=1:length(BSoh);if ii==1328      ,figure(1);loglog(fson(1:35),fson(1:35).*BSoh(ii,1:35),fson(2:35),fson(2:35).*BSlh(ii,2:35),fson(2:35),fson(2:35).*BSq(ii,2:35),fson(11:35),((fson(11:35).^(-2/3)))./(fson(19)^(-2/3)).*fson(19).*BSoh(ii,19));end;end;


