FINDLOG: search near file: ../data/ladcp/r0056_01.rbb READSERIAL_WH: parse file ../data/ladcp/r0056_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/r0056_01.rbb ISBB: BB-data raw data has 432+2 bytes per ensemble read 16985 ensemble and 16 bins using increased bottom echo amplitudes to create bottom track use bin below maximum for bottom track created 1336 bottom track data from normal velocities removed 0 values because of low percent good found 992 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 8491 bins down looking ==> load data took 246 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:37°N 30.9456' 19°W 59.9969' updated magnetic deviation and applied correction update end pos from:NaN°N NaN' NaN°E NaN' to:37°N 29.9499' 19°W 59.9954' LOADCTDPROF: load CTD profile ../ladcp/ctd/prelim_ctd/rh103w0561.ctd CTD max depth : 4837 LOADCTD: load CTD time series ../ladcp/ctd/sio_ascts/056.01 number of CTD scans: 14178 delta t : 1.5 seconds Adjusted CTD Time by 4.9863 minutes max CTD pressure 4929 at 04-Jul-2003 06:26:20 max LADCP depth 4787 at 04-Jul-2003 08:27:20 WARNING WARNING bottom time difference between CTD/LADCP is 121 minute maybe the hour was set wrong: shift CTD time by :0.083333 days Changed Start Time : 04-Jul-2003 04:44:07 to 04-Jul-2003 06:56:17 Changed End Time : 04-Jul-2003 10:52:06 to 04-Jul-2003 12:41:21 CTD max depth : 4836 best lag : -48 pings ~ -62 seconds corr:0.13451 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 use LADCP depth from constrained integrated W no surface reflections up-cast maximum depth from int W is :4836.4664 should be :4836.4664 no bottom found given maximum profile depth : 4836 extracted bottom depth : 4887 bottom depth error : 42 starting run 2 to LADCP depth bottom:1197.7272 @ 4203 found 386 ill constrained elements use LADCP depth from constrained integrated W maximum depth from int W is :4856.091 should be :4836.4664 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 321 In /home/felixt/ladcp/m/laproc.m at line 75 In /g1/felixt/A16N/scripts/runme.m at line 127 bottom found at 4902 +/- 4 m removing 271 values below bottom last LADCP depth is -3616 remove data at begining and end of cast adjust start and end time find new start end position ==> get depth took 231 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 691 LOADSADCP: load SADCP data file ../sadcp/00644v3.nc found 3 SADCP profiles ==> prep inversion took 64 seconds GETINV: compute best velocity profile Barotropic velocity error 0.0051199 [m/s] super ensemble velocity error 0.011604 [m/s] set velocity error to:0.017891 [m/s] vertical resolution (ps.dz) is 10 [m] preaveraged GPS ships vel 1 times smooth Ocean velocity profile found 6 ill constrained elements smooth CTD velocity profile found 4 ill constrained elements bottom inversion mean bottom track weight : 0.78893 weight for bottom track is (ps.botfac) 1 mean sadcp weight : 0.47078 weight for SADCP vel is (ps.sadcpfac) 1 add barotropic constrain weight: 3.4943 weight for drag is (ps.dragfac) 0.19691 ready for inversion length of d : 7962 (CTD vel) length of A1:691 (ocean vel) length of A2:488 Moore-Penrose inverse solve only down trace smooth Ocean velocity profile found 111 ill constrained elements smooth CTD velocity profile Cholesky transform solve only up trace smooth Ocean velocity profile found 16 ill constrained elements smooth CTD velocity profile found 3 ill constrained elements Cholesky transform ==> inversion took 189 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