FINDLOG: search near file: ../k002d/K002D000.000
 READSERIAL_WH: parse file ../k002d/K002D000.LOG
 down looking serial number is : 3898
 FINDLOG: search near file: ../k002u/K002U000.000
 READSERIAL_WH: parse file ../k002u/K002U000.LOG
   up looking serial number is : 754
LOADRDI: Loading data
 loading down-data ../k002d/K002D000.000
ISBB: BB-data 
 raw data has 733+2 bytes per ensemble
 WARNING  found 11 ensembles with bad checksum 
 read 474 ensemble and 25 bins 
 using RDI bottom track
 removed 0 values because of low percent good
 loading up-data ../k002u/K002U000.000
 raw data has 653+2 bytes per ensemble
 WARNING  found 12 ensembles with bad checksum 
 read 479 ensemble and 21 bins 
 removed 0 values because of low percent good
 try to shift timeseries by lag: -1 correlation: 0.87503
 best lag not obvious!  use time to match up-down looking ADCP
 mean lag is 0.97679 ensembles
 number of joint ensembles is : 474
 found 202 finite bottom velocities
 found 8 NaN bottom error velocities and discarded them
 removed 364 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 24 profiles due to larger tilt difference than 4
 multiply weight of bin 1 by 0.5
 Outlier discarded 400 bins down looking
 Outlier discarded 40 bins up looking
==> load data took 5 seconds
LOADNAV: load NAV time series ../ctdnav/00201.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: 1751  delta t : 0.43198 seconds
 update start pos  from:NaN°N     NaN'  NaN°E     NaN'
                     to:43°N 10.7358'   50°W 40.8942'
 updated magnetic deviation and applied correction 
 update end pos    from:NaN°N     NaN'  NaN°E     NaN'
                     to:43°N 10.8762'   50°W 40.9908'
LOADCTD: load CTD time series ../ctdnav/00201.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: 1751  delta t : 0.43198 seconds
 max CTD pressure 73 at 26-Sep-2003 01:26:30
 max LADCP depth 55 at 26-Sep-2003 01:26:17
 Changed Start Time : 26-Sep-2003 01:13:13  to 26-Sep-2003 01:23:53
 Changed End   Time : 26-Sep-2003 01:40:33  to 26-Sep-2003 01:31:25
 CTD max depth : 72
 best lag : 0 pings ~ 0 seconds  corr:0.84345
 adjust ADCP time to CTD time and shift depth record 
 adjust start/end position
 WARNING found 99% of  hbot=0  WARNING
 call own bottom track routine to get bottom
 using increased bottom echo amplitudes to create bottom track
 use bin below maximum for bottom track 
 created 103 bottom track data from normal velocities
GETDPTHI: Depth from vertical velocity inverse method
 starting run 1 to LADCP depth
 take soundspeed from CTD profile 
 correct velocities for sound speed 
 shallow station no surface detection 
 maximum depth from int W is :71.6504
 should be                   :71.6504
  bottom found at 86 +/- 1 m
 removing 1006 values below bottom
==> get depth took 4 seconds
PREPINV: prepare data for inversion, form Super-Ensembles
 average profiles over (p.avens) 1 ensembles
 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 -28.8583 deg
 mean compass offset based on tilt is -32.9562 deg
 rot up2down use mean up/down compass
 will not correct for sound speed twice
 avdz=NAN  => No pre-averaging done !!!
 mean compass offset based on tilt is -32.9562 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
 avdz=NAN  => No pre-averaging done !!!
LOADSADCP: load SADCP data file ../sadcp/A20_sadcp.mat
 found 1 SADCP profiles 
==> prep inversion took 3 seconds
GETINV:  compute best velocity profile
 Barotropic velocity error 0.13274 [m/s]
 super ensemble velocity error 0.064821 [m/s]
 set velocity error to:0.1011 [m/s]
 vertical resolution (ps.dz) is 10 [m]
 preaveraged GPS ships vel 1 times 
 smooth Ocean velocity profile
 smooth CTD velocity profile
 found 1 ill constrained elements
 bottom inversion 
 mean bottom track weight : 0.66067
 weight for bottom track is (ps.botfac) 1
 mean sadcp weight : 1
 weight for SADCP vel is (ps.sadcpfac) 1
 add barotropic constrain weight: 0.76163
 weight for drag is (ps.dragfac) 0.007162
 ready for inversion  length of d : 1074
           (CTD vel)  length of A1:125
         (ocean vel)  length of A2:8
Moore-Penrose inverse
 solve only down trace
 smooth Ocean velocity profile
 smooth CTD velocity profile
Cholesky transform
 solve only up trace
 smooth Ocean velocity profile
 found 1 ill constrained elements
 smooth CTD velocity profile
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 84 % of data 
 compute vertical diffusivity
 use dz:640 m long segments
 save results