fclose all;
way_rawdata='C:\Data';

read_parameters_rhb_GOMECC;
disp('read_son_day_rhb_07')

sonic=[];
plotit=1;
prtit=1;

jhold=0;
jax=1;
jaxx=1;
backchk=0;


for ddd=210;
    jd=num2str(ddd);
    if  ddd<100;
        jd=['0' num2str(ddd)];	
    end;
    for hhh=0:5;
        hr=num2str(hhh);
        if hhh<10;
           hr=['0' num2str(hhh)];	
        end;


        flist=fopen([way_rawdata, '\GOMECC\day' ,jd, '\P0_',jd,hr]);
        disp(['P0_', jd,hr]);
   
        if flist>0;
           dump=fgetl(flist);
           sonic_data=fscanf(flist,'%g,%g,%g,%g,%g',[5,inf]);    
           sonic_data=sonic_data'; 
     
        end;
    
        ss = size(sonic_data);
        if ss(1)==0;
            break;
        end;
    

        if outputsonic==1;
           Tsonic=sonic_data(:,5)/100;      %Temperature in degC
        else
           speedson=sonic_data(:,5)/50;
           Tsonic=(speedson.*speedson+(u.*u + w.*w)/2. + v.*v)/403-273.15;
        end;

        time=sonic_data(:,1)/1000;          %time in "s"
        u=sonic_data(:,2)/100;
        v=sonic_data(:,3)/100;
        w=sonic_data(:,4)/100;


        
        if rotationsonic==1
           C=strcmp(sonicmodel,'R3');
           if C==1                                              %means sonicmodel=R3
              U=u*cos(30/180*pi)-v*sin(30/180*pi);              %to North
              V=u*sin(30/180*pi)+v*cos(30/180*pi);              %to West
           else
              U=-u*cos(30/180*pi)-v*sin(30/180*pi);             %to North
              V=-u*sin(30/180*pi)+v*cos(30/180*pi);             %to West
           end;
        else
            U=u;
            V=v;
        end;

        if length(w)> 72000 | length(w)< 72000 
           time=interpft(time,72000);
           w=interpft(w,72000);
           v=interpft(V,72000);
           u=interpft(U,72000);
           Tsonic=interpft(Tsonic,72000);
        end;
        
        time=(1:71880);
        u=u(1:71880);
        v=v(1:71880);
        w=w(1:71880); 
        Tsonic=Tsonic(1:71880);
        W=[W;w];
        
        sonic=[time', u, v, w, Tsonic];
        save ([way_rawdata, '\GOMECC\Analysis\sonic_' , jd, hr, '.txt'], 'sonic', '-ascii')
    end;
end;


S=sqrt(u.*u+v.*v+w.*w);
dir=atan2(-V,U)*180/pi+180;		%jeff coordinates, +boward, +portward, 180 means wind from
dir=mod(dir+720,360);
np=jhold;

%if prtit==1;
%   prt_jas_son_rhb_07_GOMECC;
%end;

if plotit==1;
    %set some plot limits
    Tll=-25;
    hll=70;
    uul=10;
    vul=10;
    Tul=-0;
    hul=100;
    sul=12;
    Rlul=30;
    Rlll=-100;
    Rtul=300;
    Rtll=125;
    Rsul=1000;
    Rsll=0;

    %plot T profiles
    j_st=jadson(1);
    j_en=jadson(np);
    figure;plot(jadson(1:np),Tsonic(1:np));xlabel(['Jday', year]);ylabel('Temperature (C)');title('Sonic Temperature');
    figure;plot(jadson(1:np),U(1:np),jadson(1:np),V(1:np),jadson(1:np),w(1:np),jadson(1:np),S(1:np));xlabel(['Jday', year]);ylabel('Sonic Winds (m/s)');title('Relative Wind Components:  U=blue,  V=green,  W=red,   Speed=cyan');

    for i=1:np % Jeff Otten's mod to display direction -180 to +180.
        if dir(i) > 180
            dir(i) = dir(i) - 360;
        end;
    end;
    figure;plot(jadson(1:np),dir(1:np),'b', ... 		%data points
        [ddd+4/24 ddd+4/24],[-180 180],'k', ...		%4 hour watch dividers
        [ddd+8/24 ddd+8/24],[-180 180],'k', ...
        [ddd+12/24 ddd+12/24],[-180 180],'k', ...
        [ddd+16/24 ddd+16/24],[-180 180],'k', ...
        [ddd+20/24 ddd+20/24],[-180 180],'k', ...
        [ddd ddd+1],[45 45],'g', ...			%45 degree good zone
        [ddd ddd+1],[-45 -45],'g', ...
        [ddd ddd+1],[60 60],'r', ...			%60 degree questionable zone
        [ddd ddd+1],[-60 -60],'r'); ...
        xlabel('Jday (2006)'); ...
        ylabel('Sonic Wind Direction (deg)'); ...
        title('Relative Wind Direction'); ...
        axis([j_st j_en -180 180]);
        year=2007;
    figure;plot(jadson(1:np),sqrt(U(1:np).^2+V(1:np).^2+w(1:np).^2));
    xlabel(['Jday', year]);ylabel('Wind Speed m/s');title('Sonic Relative Wind Speed');%axis([j_st j_en 0 360]);
end;

clear w vul v uul u tzilch tzero sul st3 speedson prtit prt_it plotit nr1 npz np nl1 jl jk jj jhold jd jazx jaxx jax jam jadson jade j_st j_en j hul hr hll flist
clear f e dmp2 dmp dir backchk ans V U Tul Tsonic Tll S Rtul Rtll Rsul Rsll Rlul Rlll



