FINDLOG: search near file: ../data/ladcp/r0064_01.rbb
 READSERIAL_WH: parse file ../data/ladcp/r0064_01.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/r0064_01.rbb
ISBB: BB-data 
 raw data has 432+2 bytes per ensemble
 DETCTED BEAM coordinates try to make earth coordinates!
 begin to convert 11821 profiles from beam to earth
 read 11821 ensemble and 16 bins 
 using increased bottom echo amplitudes to create bottom track
 use bin below maximum for bottom track 
 created 801 bottom track data from normal velocities
 removed 0 values because of low percent good
 found 625 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 0 profiles due to larger tilt difference than 4
 multiply weight of bin 1 by 0.5
 Outlier discarded 5789 bins down looking
==> load data took 177 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:33°N 31.3446'   21°W 23.6243'
 updated magnetic deviation and applied correction 
 update end pos    from:NaN°N     NaN'  NaN°E     NaN'
                     to:33°N 29.9906'   21°W 24.4411'
LOADCTDPROF: load CTD profile ../ladcp/ctd/prelim_ctd/rh103w0641.ctd
 CTD max depth : 5314
LOADCTD: load CTD time series ../ladcp/ctd/sio_ascts/064.01
 number of CTD scans: 9563  delta t : 1.5 seconds
 Adjusted CTD Time by 7.9888 minutes
 max CTD pressure 5420 at 06-Jul-2003 23:19:25
 max LADCP depth 5216 at 06-Jul-2003 23:21:49
WARNING WARNING 
 bottom time difference between CTD/LADCP is 2 minute
 Changed Start Time : 06-Jul-2003 21:31:19  to 06-Jul-2003 21:44:50
 Changed End   Time : 07-Jul-2003 01:47:25  to 07-Jul-2003 01:36:24
 CTD max depth : 5315
 best lag : -100 pings ~ -130 seconds  corr:0.94353
 lag to big or correlation to weak ignore CTD time series
 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 :5314.9827
 should be                   :5314.9827
 use LADCP depth from constrained integrated W 
 starting run 2 to LADCP depth
 bottom:5260.858 @ 4477
 found 30 ill constrained elements
 maximum depth from int W is :5318.4079
 should be                   :5314.9827
 use LADCP depth from constrained integrated W 
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 5339 +/- 4 m
 removing 1275 values below bottom
 remove data at begining and end of cast
 adjust start and end time 
 find new start end position
==> get depth took 177 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 661 bins down looking
 reduced ensemble size is 884
LOADSADCP: load SADCP data file ../sadcp/00644v3.nc
 found 4 SADCP profiles 
==> prep inversion took 97 seconds
GETINV:  compute best velocity profile
 Barotropic velocity error 0.0043605 [m/s]
 super ensemble velocity error 0.0095268 [m/s]
 set velocity error to:0.015922 [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.70121
 weight for bottom track is (ps.botfac) 1
 mean sadcp weight : 0.52809
 weight for SADCP vel is (ps.sadcpfac) 1
 add barotropic constrain weight: 3.6514
 weight for drag is (ps.dragfac) 0.19805
 ready for inversion  length of d : 10092
           (CTD vel)  length of A1:884
         (ocean vel)  length of A2:532
Moore-Penrose inverse
 solve only down trace
 smooth Ocean velocity profile
 found 9 ill constrained elements
 smooth CTD velocity profile
 found 1 ill constrained elements
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 395 seconds
GETSHEAR2: average shear profile over (ps.dz) 10 [m]
 maximum std (stdf) 2 of data 
 minimum weight  0.1 of data 
 will use 89 % of data 
 increased error because of shear - inverse difference
 compute vertical diffusivity
 use dz:640 m long segments
 save results