NCCOOS Trac Projects: Top | Web | Platforms | Processing | Viz | Sprints | Sandbox | (Wind)

Changeset 232

Show
Ignore:
Timestamp:
07/17/09 14:27:12
Author:
gdusek
Message:

--

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • DPWavesProc/trunk/DPWavesProc/adcp_matlab/radialtouvw.m

    r229 r232  
    2626%set up sysinfo file 
    2727samplesInBurst=sysinfo(3,:); 
     28binsOut=sysinfo(8,:); 
    2829bin1height=sysinfo(9,:); 
    2930bin2height=sysinfo(10,:); 
    …… 
    3536roll=sysinfo(20,:); 
    3637 
     38 
    3739%set up pressure 
    3840press=pressure/1000; 
    …… 
    5355orbit=orbit/1000; 
    5456orbitnew=orbit; 
     57 
     58% calculate the total number of orbital bins output (usually 20) 
     59orbOut=binsOut*4; 
     60 
    5561%interpolate to take out any NaNs and QC for bad data in orbital data 
    56 std_orbit=ones(1,20); 
    57 for i=1:20 
     62std_orbit=ones(1,orbOut); 
     63for i=1:orbOut 
    5864     
    5965    %first take out any points outside of 4 std deviations 
    …… 
    6369end 
    6470%find the avg std deviation for each group 4 beams 
    65 for i=4:4:20 
     71for i=4:4:orbOut 
    6672    avgstd_orbit(i-3:i)=mean(std_orbit(i-3:i)); 
    6773end 
    6874 
    6975%now remove the points outside 4avg std dev and interp 
    70 for i=1:20 
     76for i=1:orbOut 
    7177    ibad_std=find(abs(orbit(:,i)) > 4*avgstd_orbit(i)); 
    7278    orbit(ibad_std,i)=NaN; 
    …… 
    131137%reshape the orbital matrix to 3 dimensions 
    132138 
    133 radial=reshape(orbitnew,samplesInBurst,4,5); 
     139radial=reshape(orbitnew,samplesInBurst,4,binsOut); 
    134140    
    135141%  Compute uvw velocities and change to earth coordinates 
    …… 
    140146error_vel = D1*(radial(:,1,:)+radial(:,2,:)-radial(:,3,:)-radial(:,4,:)); 
    141147    
    142 u=reshape(u,samplesInBurst,5); 
    143 v=reshape(v,samplesInBurst,5); 
    144 w=reshape(w,samplesInBurst,5); 
    145 error_vel=reshape(error_vel,samplesInBurst,5); 
     148u=reshape(u,samplesInBurst,binsOut); 
     149v=reshape(v,samplesInBurst,binsOut); 
     150w=reshape(w,samplesInBurst,binsOut); 
     151error_vel=reshape(error_vel,samplesInBurst,binsOut); 
    146152 
    147153[m,n] = size(u); 
    148 uno = ones(m,5); 
     154uno = ones(m,binsOut); 
    149155unew = u.*(uno*a) + v.*(uno*b) + w.*(uno*c); 
    150156vnew = u.*(uno*d) + v.*(uno*e) + w.*(uno*f); 

  • DPWavesProc/trunk/DPWavesProc/adcp_matlab/specmultiplot.m

    r231 r232  
    6060 
    6161    %set up data with radial velocities, freq and dir at default 
    62     [ID,SM,EP]=radialtouvw(pressure,range,orbit,sysinfo,1); 
     62    [ID,SM,EP]=radialtouvw(pressure,range,orbit,sysinfo,3); 
    6363 
    6464    % run diwasp to generate this spectrum with EMEP