FINDLOG: search near file: ../data/ladcp/r0144_02.rbb READSERIAL_WH: parse file ../data/ladcp/r0144_02.log Can not find RDI serial number down looking serial number is : NaN up looking serial number is : NaN LOADRDI: Loading data loading down-data ../data/ladcp/r0144_02.rbb ISBB: BB-data raw data has 752+2 bytes per ensemble DETCTED BEAM coordinates try to make earth coordinates! begin to convert 11275 profiles from beam to earth read 11275 ensemble and 32 bins using increased bottom echo amplitudes to create bottom track use bin below maximum for bottom track created 1668 bottom track data from normal velocities removed 0 values because of low percent good found 908 finite bottom velocities removed 0 values because of high error velocity removed 0 bottom values because of high error velocity using whole profile since no start time was given removed 0 profiles due to larger tilt than 22 removed 8 profiles due to larger tilt difference than 4 multiply weight of bin 1 by 0.5 Outlier discarded 8771 bins down looking ==> load data took 210 seconds load NAV time series ../ladcp_output/A16N_nav.mat number of NAV scans: 437404 delta t : 10 seconds update start pos from:NaN°N NaN' NaN°E NaN' to: 2°S 59.5977' 25°W 0.5723' updated magnetic deviation and applied correction update end pos from:NaN°N NaN' NaN°E NaN' to: 2°S 59.9223' 25°W 0.1271' LOADCTDPROF: load CTD profile ../ladcp/ctd/prelim_ctd/rh203w1442.ctd CTD max depth : 5315 LOADCTD: load CTD time series ../ladcp/ctd/sio_ascts/144.02 number of CTD scans: 26850 delta t : 0.50002 seconds Adjusted CTD Time by 14.991 minutes max CTD pressure 5412 at 07-Aug-2003 13:06:31 max LADCP depth 5262 at 07-Aug-2003 13:04:52 WARNING WARNING bottom time difference between CTD/LADCP is -2 minute Changed Start Time : 07-Aug-2003 11:11:46 to 07-Aug-2003 11:31:07 Changed End Time : 07-Aug-2003 15:16:03 to 07-Aug-2003 15:08:25 CTD max depth : 5315 best lag : -9 pings ~ -12 seconds corr:0.99385 adjust ADCP time to CTD time and shift depth record adjust start/end position GETDPTHI: Depth from vertical velocity inverse method starting run 1 to LADCP depth take soundspeed from CTD profile Warning: Divide by zero. > In /usr/local/matlab61/toolbox/matlab/polyfun/interp1q.m at line 34 In /home/felixt/ladcp/m/getdpthi.m at line 147 In /home/felixt/ladcp/m/laproc.m at line 75 In /g1/felixt/A16N/scripts/runme.m at line 127 correct velocities for sound speed maximum depth from int W is :5314.8094 should be :5314.8094 Warning: Polynomial is badly conditioned. Remove repeated data points or try centering and scaling as described in HELP POLYFIT. > In /usr/local/matlab61/toolbox/matlab/polyfun/polyfit.m at line 74 In /home/felixt/ladcp/m/getdpthi.m at line 310 In /home/felixt/ladcp/m/laproc.m at line 75 In /g1/felixt/A16N/scripts/runme.m at line 127 bottom found at 5321 +/- 2 m removing 1465 values below bottom ==> get depth took 363 seconds PREPINV: prepare data for inversion, form Super-Ensembles average profiles over (p.avdz) 10 meter discarded 0 bottom tracks velocities because of dist > 300 reduce weight for larger tilts 0.5 at 10 degree will not correct for sound speed twice remove reference velocity and average ensembles Outlier discarded 1549 bins down looking reduced ensemble size is 887 LOADSADCP: load SADCP data file ../sadcp/00644v3.nc found 4 SADCP profiles ==> prep inversion took 115 seconds GETINV: compute best velocity profile Barotropic velocity error 0.0046019 [m/s] super ensemble velocity error 0.010926 [m/s] set velocity error to:0.02337 [m/s] vertical resolution (ps.dz) is 10 [m] preaveraged GPS ships vel 1 times smooth Ocean velocity profile found 8 ill constrained elements smooth CTD velocity profile found 3 ill constrained elements bottom inversion mean bottom track weight : 0.83205 weight for bottom track is (ps.botfac) 1 mean sadcp weight : 0.59415 weight for SADCP vel is (ps.sadcpfac) 1 add barotropic constrain weight: 5.0782 weight for drag is (ps.dragfac) 0.19804 ready for inversion length of d : 18730 (CTD vel) length of A1:887 (ocean vel) length of A2:530 Moore-Penrose inverse solve only down trace smooth Ocean velocity profile found 8 ill constrained elements smooth CTD velocity profile found 1 ill constrained elements Cholesky transform solve only up trace smooth Ocean velocity profile found 8 ill constrained elements smooth CTD velocity profile found 2 ill constrained elements Cholesky transform ==> inversion took 369 seconds GETSHEAR2: average shear profile over (ps.dz) 10 [m] maximum std (stdf) 2 of data minimum weight 0.1 of data will use 95 % of data increased error because of shear - inverse difference compute vertical diffusivity use dz:640 m long segments save results