FINDLOG: search near file: ../k040/A0813001.612
 READSERIAL_WH: parse file ../k040/RECOVER.LOG
 down looking serial number is : 150
 FINDLOG: search near file: ../k040s/K040S000.000
 READSERIAL_WH: parse file ../k040s/K040S000.LOG
   up looking serial number is : 150
LOADRDI: Loading data
 loading down-data ../k040/A0813001.612
ISBB: BB-data 
 raw data has 432+2 bytes per ensemble
 read 11384 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 1072 bottom track data from normal velocities
 removed 0 values because of low percent good
 loading up-data ../k040s/K040S000.000
 raw data has 733+2 bytes per ensemble
 read 4285 ensemble and 25 bins 
 removed 0 values because of low percent good
 WARNING WARNING WARNING WARNING WARNING WARNING WARNING WARNING 
 warning ping rate not equal in down instrument
  min down ping rate :0.99999  max down ping rate :1.6
 warning mean ping rate not equal between both instruments 
  mean down ping rate :1.3  mean up ping rate :3.5
 find best time match of up-looking ADCP to down looking ADCP
 up instrument is different by -2674 ensembles
 WARNING WARNING WARNING WARNING WARNING WARNING WARNING WARNING 
 try to shift timeseries by lag: -4 correlation: 0.93719
 shift ADCP timeseries by -4 ensembles
 number of joint ensembles is : 11380
 found 572 finite bottom velocities
 removed 456 values because of high error velocity
 removed 0 bottom values because of high error velocity
** WARNING   WARNING  WARNING  WARNING  WARNING  WARNING **  
**  found  6  large horizontal velocities
** WARNING  check maximum velocity setting on CMD-file   ** 
 using whole profile since no start time was given
 removed 0 profiles due to larger tilt than 22
 removed 65 profiles due to larger tilt difference than 4
 multiply weight of bin 1 by 0.5
 Outlier discarded 5842 bins down looking
 Outlier discarded 4422 bins up looking
==> load data took 230 seconds
LOADNAV: load NAV time series ../ctdnav/04001.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: 28032  delta t : 0.86397 seconds
 update start pos  from:NaN°N     NaN'  NaN°E     NaN'
                     to:24°N 11.1780'   52°W 20.5086'
 updated magnetic deviation and applied correction 
 update end pos    from:NaN°N     NaN'  NaN°E     NaN'
                     to:24°N 12.3942'   52°W 19.8888'
LOADCTD: load CTD time series ../ctdnav/04001.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: 28032  delta t : 0.86397 seconds
 max CTD pressure 5513 at 08-Oct-2003 10:55:14
 max LADCP depth 5326 at 08-Oct-2003 10:54:28
WARNING WARNING 
 bottom time difference between CTD/LADCP is -1 minute
 Changed Start Time : 08-Oct-2003 09:12:22  to 08-Oct-2003 09:21:41
 Changed End   Time : 08-Oct-2003 13:18:54  to 08-Oct-2003 13:08:35
 CTD max depth : 5408
 best lag : 5 pings ~ 6 seconds  corr:0.98121
 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 :5408.2403
 should be                   :5408.2403
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 5421 +/- 0 m
 removing 771 values below bottom
==> get depth took 882 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
 mean compass offset based on velocity is 84.6155 deg
 mean compass offset based on tilt is 87.8984 deg
 rot up2down use mean up/down compass
 will not correct for sound speed twice
 remove reference velocity and average ensembles 
 Outlier discarded 761 bins down looking
 Outlier discarded 500 bins up looking
 reduced ensemble size is 898
 mean compass offset based on tilt is 87.8984 deg
 remove first guess ocean velocity from raw data
 adjusted for velocity offset in up and down looking ADCP
 remove first guess ocean velocity from raw data
 rot up2down use mean up/down compass
 rotated earlier, use difference 
 will not correct for sound speed twice
 remove reference velocity and average ensembles 
 Outlier discarded 759 bins down looking
 Outlier discarded 479 bins up looking
 reduced ensemble size is 898
LOADSADCP: load SADCP data file ../sadcp/A20_sadcp.mat
 found 15 SADCP profiles 
==> prep inversion took 398 seconds
GETINV:  compute best velocity profile
 Barotropic velocity error 0.0044072 [m/s]
 super ensemble velocity error 0.029424 [m/s]
 set velocity error to:0.028559 [m/s]
 vertical resolution (ps.dz) is 10 [m]
 preaveraged GPS ships vel 1 times 
 smooth Ocean velocity profile
 found 1 ill constrained elements
 smooth CTD velocity profile
 bottom inversion 
 mean bottom track weight : 0.77422
 weight for bottom track is (ps.botfac) 1
 mean sadcp weight : 0.77389
 weight for SADCP vel is (ps.sadcpfac) 1
 add barotropic constrain weight: 6.4801
 weight for drag is (ps.dragfac) 0.19822
 ready for inversion  length of d : 13926
           (CTD vel)  length of A1:898
         (ocean vel)  length of A2:540
Moore-Penrose inverse
 sorry dz not constant loop 11314 elements
 solve only down trace
 smooth Ocean velocity profile
 found 1 ill constrained elements
 smooth CTD velocity profile
Cholesky transform
 solve only up trace
 smooth Ocean velocity profile
 smooth CTD velocity profile
Cholesky transform
==> inversion took 532 seconds
GETSHEAR2: average shear profile over (ps.dz) 10 [m]
 maximum std (stdf) 2 of data 
 minimum weight  0.1 of data 
 will use 59 % of data 
 increased error because of shear - inverse difference
 compute vertical diffusivity
 use dz:640 m long segments
 save results