FINDLOG: search near file: ../k078/A1716001.612
 FINDLOG: give up
 FINDLOG: search near file: ../k078/A1716001.612
 READSERIAL_BB: parse file ../k078/K078
 down looking serial number is : 2612
   up looking serial number is : NaN
LOADRDI: Loading data
 loading down-data ../k078/A1716001.612
ISBB: BB-data 
 raw data has 432+2 bytes per ensemble
 read 2106 ensemble and 16 bins 
 using increased bottom echo amplitudes to create bottom track
 use bin below maximum for bottom track 
 created 388 bottom track data from normal velocities
 removed 0 values because of low percent good
 found 388 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 10 profiles due to larger tilt difference than 4
 multiply weight of bin 1 by 0.5
 Outlier discarded 903 bins down looking
==> load data took 14 seconds
LOADNAV: load NAV time series ../ctdnav/07801.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: 4041  delta t : 0.43198 seconds
 update start pos  from:NaN°N     NaN'  NaN°E     NaN'
                     to: 7°N 25.2342'   53°W 20.0520'
 updated magnetic deviation and applied correction 
 update end pos    from:NaN°N     NaN'  NaN°E     NaN'
                     to: 7°N 25.5960'   53°W 19.8768'
LOADCTD: load CTD time series ../ctdnav/07801.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: 4041  delta t : 0.43198 seconds
 max CTD pressure 471 at 17-Oct-2003 16:31:01
 max LADCP depth 448 at 17-Oct-2003 16:30:41
 Changed Start Time : 17-Oct-2003 16:10:29  to 17-Oct-2003 16:21:06
 Changed End   Time : 17-Oct-2003 16:56:05  to 17-Oct-2003 16:48:47
 CTD max depth : 467
 best lag : 41 pings ~ 53 seconds  corr:0.98114
 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 156
  In /home/felixt/ladcp/m/laproc.m at line 75
  In /home/felixt/A20/scripts/runme.m at line 85
 correct velocities for sound speed 
 maximum depth from int W is :467.3236
 should be                   :467.3236
  bottom found at 478 +/- 1 m
 removing 237 values below bottom
==> get depth took 23 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 48 bins down looking
 reduced ensemble size is 49
LOADSADCP: load SADCP data file ../sadcp/A20_sadcp.mat
 found 3 SADCP profiles 
==> prep inversion took 5 seconds
GETINV:  compute best velocity profile
 Barotropic velocity error 0.036123 [m/s]
 super ensemble velocity error 0.010771 [m/s]
 set velocity error to:0.042511 [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 4 ill constrained elements
 bottom inversion 
 mean bottom track weight : 0.7819
 weight for bottom track is (ps.botfac) 1
 mean sadcp weight : 1.1632
 weight for SADCP vel is (ps.sadcpfac) 1
 add barotropic constrain weight: 1.1768
 weight for drag is (ps.dragfac) 0.0459
 ready for inversion  length of d : 775
           (CTD vel)  length of A1:49
         (ocean vel)  length of A2:47
Moore-Penrose inverse
 solve only down trace
 smooth Ocean velocity profile
 found 10 ill constrained elements
 smooth CTD velocity profile
Cholesky transform
 solve only up trace
 smooth Ocean velocity profile
 found 9 ill constrained elements
 smooth CTD velocity profile
 found 2 ill constrained elements
Cholesky transform
==> inversion took 2 seconds
GETSHEAR2: average shear profile over (ps.dz) 10 [m]
 maximum std (stdf) 2 of data 
 minimum weight  0.1 of data 
 will use 90 % of data 
 increased error because of shear - inverse difference
 compute vertical diffusivity
 use dz:640 m long segments
 save results