FINDLOG: search near file: ../k045/A1003001.612 FINDLOG: give up FINDLOG: search near file: ../k045/A1003001.612 READSERIAL_BB: parse file ../k045/K045 down looking serial number is : 2612 up looking serial number is : NaN LOADRDI: Loading data loading down-data ../k045/A1003001.612 ISBB: BB-data raw data has 432+2 bytes per ensemble read 9595 ensemble and 16 bins using increased bottom echo amplitudes to create bottom track Warning: Divide by zero. > In /home/felixt/ladcp/m/localmax2.m at line 36 In /home/felixt/ladcp/m/loadrdi.m (updown) at line 614 In /home/felixt/ladcp/m/loadrdi.m at line 51 In /home/felixt/ladcp/m/laproc.m at line 32 In /home/felixt/A20/scripts/runme.m at line 85 use bin below maximum for bottom track created 386 bottom track data from normal velocities removed 0 values because of low percent good found 379 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 22 profiles due to larger tilt difference than 4 multiply weight of bin 1 by 0.5 Outlier discarded 4216 bins down looking ==> load data took 98 seconds LOADNAV: load NAV time series ../ctdnav/04501.ctd Warning: COLON arguments must be real scalars. > In /home/felixt/A20/scripts/loadnav.m at line 33 In /home/felixt/ladcp/m/laproc.m at line 43 In /home/felixt/A20/scripts/runme.m at line 85 number of NAV scans: 22663 delta t : 0.43198 seconds update start pos from:NaN°N NaN' NaN°E NaN' to:20°N 50.2404' 52°W 20.3946' updated magnetic deviation and applied correction update end pos from:NaN°N NaN' NaN°E NaN' to:20°N 51.0894' 52°W 20.7360' LOADCTD: load CTD time series ../ctdnav/04501.ctd Warning: COLON arguments must be real scalars. > In /home/felixt/A20/scripts/loadctd.m at line 46 In /home/felixt/ladcp/m/laproc.m at line 60 In /home/felixt/A20/scripts/runme.m at line 85 number of CTD scans: 22663 delta t : 0.43198 seconds max CTD pressure 4404 at 10-Oct-2003 01:08:07 max LADCP depth 4292 at 10-Oct-2003 01:07:23 WARNING WARNING bottom time difference between CTD/LADCP is -1 minute Changed Start Time : 09-Oct-2003 23:39:37 to 09-Oct-2003 23:52:21 Changed End Time : 10-Oct-2003 03:07:29 to 10-Oct-2003 02:55:46 CTD max depth : 4332 best lag : 5 pings ~ 6 seconds corr:0.98502 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 correct velocities for sound speed maximum depth from int W is :4332.0236 should be :4332.0236 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 330 In /home/felixt/ladcp/m/laproc.m at line 75 In /home/felixt/A20/scripts/runme.m at line 85 bottom found at 4343 +/- 1 m removing 73 values below bottom ==> get depth took 554 seconds PREPINV: prepare data for inversion, form Super-Ensembles average profiles over (p.avdz) 16 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 465 bins down looking reduced ensemble size is 478 LOADSADCP: load SADCP data file ../sadcp/A20_sadcp.mat found 13 SADCP profiles ==> prep inversion took 36 seconds GETINV: compute best velocity profile Barotropic velocity error 0.0054521 [m/s] super ensemble velocity error 0.011055 [m/s] set velocity error to:0.020924 [m/s] vertical resolution (ps.dz) is 10 [m] preaveraged GPS ships vel 1 times smooth Ocean velocity profile found 10 ill constrained elements smooth CTD velocity profile found 3 ill constrained elements bottom inversion mean bottom track weight : 0.817 weight for bottom track is (ps.botfac) 1 mean sadcp weight : 1.031 weight for SADCP vel is (ps.sadcpfac) 1 add barotropic constrain weight: 3.8378 weight for drag is (ps.dragfac) 0.19481 ready for inversion length of d : 5003 (CTD vel) length of A1:478 (ocean vel) length of A2:432 Moore-Penrose inverse solve only down trace smooth Ocean velocity profile found 19 ill constrained elements smooth CTD velocity profile Cholesky transform solve only up trace smooth Ocean velocity profile found 10 ill constrained elements smooth CTD velocity profile Cholesky transform ==> inversion took 65 seconds GETSHEAR2: average shear profile over (ps.dz) 10 [m] maximum std (stdf) 2 of data minimum weight 0.1 of data will use 83 % of data compute vertical diffusivity use dz:640 m long segments save results