FINDLOG: search near file: ../k016d/K016D000.000 READSERIAL_WH: parse file ../k016d/K016D000.LOG down looking serial number is : 149 FINDLOG: search near file: ../k016u/K016U000.000 READSERIAL_WH: parse file ../k016u/K016U000.LOG up looking serial number is : 754 LOADRDI: Loading data loading down-data ../k016d/K016D000.000 ISBB: BB-data raw data has 733+2 bytes per ensemble read 3752 ensemble and 25 bins using RDI bottom track removed 0 values because of low percent good loading up-data ../k016u/K016U000.000 raw data has 653+2 bytes per ensemble read 3757 ensemble and 21 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 :3.55 max down ping rate :3.66 find best time match of up-looking ADCP to down looking ADCP up instrument is different by 0 ensembles WARNING WARNING WARNING WARNING WARNING WARNING WARNING WARNING try to shift timeseries by lag: 0 correlation: 0.99275 shift ADCP timeseries by 0 ensembles number of joint ensembles is : 3752 found 3001 finite bottom velocities found 9 NaN bottom error velocities and discarded them removed 409 values because of high error velocity removed 6 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 39 profiles due to larger tilt difference than 4 multiply weight of bin 1 by 0.5 Outlier discarded 3213 bins down looking Outlier discarded 2259 bins up looking ==> load data took 53 seconds LOADNAV: load NAV time series ../ctdnav/01601.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: 24916 delta t : 0.86397 seconds update start pos from:NaN°N NaN' NaN°E NaN' to:39°N 59.9718' 51°W 48.2082' updated magnetic deviation and applied correction update end pos from:NaN°N NaN' NaN°E NaN' to:40°N 1.0530' 51°W 50.9532' LOADCTD: load CTD time series ../ctdnav/01601.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: 24916 delta t : 0.86397 seconds max CTD pressure 4837 at 28-Sep-2003 11:13:53 max LADCP depth 4726 at 28-Sep-2003 11:13:28 Changed Start Time : 28-Sep-2003 09:43:37 to 28-Sep-2003 09:51:54 Changed End Time : 28-Sep-2003 13:25:08 to 28-Sep-2003 13:14:41 CTD max depth : 4746 best lag : 1 pings ~ 4 seconds corr:0.97356 adjust ADCP time to CTD time and shift depth record adjust start/end position WARNING found 100% 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 164 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 maximum depth from int W is :4746.2701 should be :4746.2701 bottom found at 4761 +/- 1 m removing 317 values below bottom ==> get depth took 657 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 -175.7267 deg mean compass offset based on tilt is -177.0622 deg rot up2down use mean up/down compass will not correct for sound speed twice remove reference velocity and average ensembles Outlier discarded 669 bins down looking Outlier discarded 530 bins up looking reduced ensemble size is 652 mean compass offset based on tilt is -177.0622 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 661 bins down looking Outlier discarded 529 bins up looking reduced ensemble size is 652 LOADSADCP: load SADCP data file ../sadcp/A20_sadcp.mat found 14 SADCP profiles ==> prep inversion took 150 seconds GETINV: compute best velocity profile Barotropic velocity error 0.0049314 [m/s] super ensemble velocity error 0.013362 [m/s] set velocity error to:0.026017 [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.74659 weight for bottom track is (ps.botfac) 1 mean sadcp weight : 0.52269 weight for SADCP vel is (ps.sadcpfac) 1 add barotropic constrain weight: 5.2758 weight for drag is (ps.dragfac) 0.19656 ready for inversion length of d : 15310 (CTD vel) length of A1:652 (ocean vel) length of A2:475 Moore-Penrose inverse 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 found 1 ill constrained elements smooth CTD velocity profile Cholesky transform ==> inversion took 132 seconds GETSHEAR2: average shear profile over (ps.dz) 10 [m] maximum std (stdf) 2 of data minimum weight 0.1 of data will use 75 % of data compute vertical diffusivity use dz:640 m long segments save results