CDF
lat lon date name zbot tim - z z_sadcp = BAR_ref_U ?6\B4 BAR_ref_V ?9/?˱
BAR_ref_descr ![Navigation][Bottom-track][SADCP]
BAR_ref_error ?qU5w
BAR_tide_U
BAR_tide_V BAR_tide_model [NA] GEN_Depth_source [measured P (CTD)] GEN_Editing_parm_descr GEN_Inverse_weight_bottom GEN_Inverse_weight_navigation GEN_Inverse_weight_smooth GEN_LADCP_cast GEN_LADCP_ensemble_time_mean_sec ?NF GEN_LADCP_ensemble_time_std_sec ?/ GEN_LADCP_station GEN_Magnetic_deviation_deg Db4g GEN_Matlab_version 9.3.0.713579 (R2017b) GEN_Ocean_depth_m 9 GEN_Percent_3beam GEN_Proc_general_comments GEN_Proc_methodology [inverse] GEN_Processing_date 11-Sep-2019 GEN_Processing_personnel A.M. Thurnherr GEN_Profile_bottom_decimal_day @\/ GEN_Profile_end_decimal_day @\#o\ GEN_Profile_end_latitude M
D GEN_Profile_end_longitude @>f1 GEN_Profile_max_depth_m / GEN_Profile_start_decimal_day @\R GEN_Profile_start_latitude L= GEN_Profile_start_longitude @> J GEN_Software_orig 'LDEO LADCP software: Version IX_14beta2 GEN_Sound_sp_calc [T-P] GEN_Velocity_Units m/s GEN_conf_general_comments INPUT_Bottom_track_profile_avail INPUT_CTD_time_series_avail INPUT_Nav_time_series_avail INPUT_Pegasus_profile_avail INPUT_SADCP_profile_avail LADCP_dn_beam_range @i
=p@i
=p@i
=p@i
=p LADCP_dn_btrk_u_bias ?o V^ LADCP_dn_btrk_u_std ?o LADCP_dn_btrk_v_bias ?V ̔ LADCP_dn_btrk_v_std ?o80M LADCP_dn_conf_bin_len_m LADCP_dn_conf_blank_intvl_m LADCP_dn_conf_bottom_trkr LADCP_dn_conf_coord_system [unknown] LADCP_dn_conf_general_comments LADCP_dn_conf_number_bins LADCP_dn_conf_number_pings LADCP_dn_conf_ping_stagr [NA] !LADCP_dn_conf_ping_trns_intvl_sec LADCP_dn_conf_pulse_len_m @ zG LADCP_dn_conf_single_ping_acc ?Y.Zϛ LADCP_dn_conf_vel_ambiguity @ LADCP_dn_hard_FV @I4zG LADCP_dn_hard_SN LADCP_dn_hard_TNO [convex4] LADCP_dn_hard_beam_ang_deg LADCP_dn_hard_comp_type [RDI] LADCP_dn_hard_freq_kHz LADCP_dn_hard_general_comments LADCP_dn_hard_type [RDI-150BB] LADCP_dn_xmit_cur @ZW䓊 LADCP_dn_xmit_pings + LADCP_dn_xmit_vol @ewLwbg LADCP_up_beam_range @ap
=@ap
=@ap
=@ap
= LADCP_up_compass [down-compass] LADCP_up_conf_bin_len_m LADCP_up_conf_blank_intvl_m LADCP_up_conf_coord_system [unknown] LADCP_up_conf_general_comments LADCP_up_conf_number_bins LADCP_up_conf_number_pings LADCP_up_conf_ping_stagr [NA] !LADCP_up_conf_ping_trns_intvl_sec LADCP_up_conf_pulse_len_m @ 8Q LADCP_up_conf_single_ping_acc ?8( LADCP_up_conf_vel_ambiguity @ LADCP_up_hard_FV @I333333 LADCP_up_hard_SN LADCP_up_hard_TNO [convex4] LADCP_up_hard_beam_ang_deg LADCP_up_hard_comp_type [RDI] LADCP_up_hard_freq_kHz , LADCP_up_hard_general_comments LADCP_up_hard_type [RDI-300WH] LADCP_up_xmit_cur @f_QyY LADCP_up_xmit_pings + LADCP_up_xmit_vol @b̷߬ LOG_Inverse_log 2LDEO LADCP software: Version IX_14beta2
################ [029] step 1: LOAD LADCP DATA #######################
LOADING CHECKPOINT checkpoints/029_0
RE-LOADING PER-CAST PARAMETERS
LOADRDI:
loading down-data ../Data/PD0/029DL.100
DETECTED BEAM coordinates: rotating to EARTH coordinates
21385 3-beam solutions calculated (9% of total)
read 11200 ensembles with 25 bins each
loading up-data ../Data/PD0/029UL.100
DETECTED BEAM coordinates: rotating to EARTH coordinates
10665 3-beam solutions calculated (8% of total)
read 11202 ensemble and 30 bins
WARNING WARNING WARNING WARNING WARNING WARNING WARNING WARNING
Warning: non-constant ping rate in downlooker data (staggered pinging?)
min down ping rate :1.36 max down ping rate :1.61
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
bestlag removed 458 spikes
try to shift timeseries by lag: 0 correlation: 0.98058
shift ADCP timeseries by 0 ensembles
number of joint ensembles is : 11200
removed 37 values because of error velocity > 0.5 m/s
removed 431 values because of horizontal speed > 2.5 m/s
** found 223 (2.0% of total) velocity measurements > 2.5 m/s
** WARNING check ambiguity velocity setting in CMD-file **
using whole profile since no start time was given
removed 4 profiles due to tilt > 22 degrees
removed 687 profiles due to tilt derivative > 4 degrees
removed 2 suspect non pinging? (low velocity gradient) ensembles from down-looker
WARNING one beam might be weak !
weak down looking beam 1
Outlier discarded 9889 bins down looking
Outlier discarded 0 bottom track
Outlier discarded 9095 bins up looking
==> STEP 1 TOOK 10.3 seconds
SAVING CHECKPOINT checkpoints/029_1
################ [029] step 2: FIX LADCP-DATA PROBLEMS ###############
==> STEP 2 TOOK 0.0 seconds
################ [029] step 3: LOAD GPS DATA #########################
LOADNAV: load NAV time series ../Data/CTD/029.1Hz
number of NAV scans: 14293 delta t : 0.99999 seconds
corrected for magnetic declination of -41.8 deg
==> STEP 3 TOOK 0.1 seconds
################ [029] step 4: GET BOTTOM-TRACK DATA #################
GETBTRACK creates own bottom track in addition to RDI
in: p.btrk_mode 3 and p.btrk_used 0
using increased bottom echo amplitudes to create bottom track
localmax2: found 7398 valid values
use 0.5 bins below maximum target strength for own bottom track velocity
created 4342 bottom distances
removed 852 bottom track profiles W_btrk - W_ref difference > 0.05
boutlier removed 220 bottom track velocities
created 3339 bottom track data from normal velocities
out: p.btrk_mode 3 and p.btrk_used 2
==> STEP 4 TOOK 0.7 seconds
################ [029] step 5: LOAD CTD PROFILE ######################
==> STEP 5 TOOK 0.0 seconds
################ [029] step 6: LOAD CTD TIME SERIES ##################
LOADCTD: load CTD time series ../Data/CTD/029.1Hz
read 14293 CTD scans; median delta_t = 1.00 seconds
interpolated to 14293 CTD scans; delta_t = 1.00 seconds
90% LADCP depth 4873 at 23-Apr-2019 11:08:12
90% CTD pressure 4990 at 23-Apr-2019 10:40:02
Changed Start Time : 23-Apr-2019 09:04:46 to 23-Apr-2019 09:37:04
Changed End Time : 23-Apr-2019 13:36:29 to 23-Apr-2019 13:26:54
median CTD time difference 0.99999 s
removed 0 pressure spikes
use 1 time base for W_ctd
CTD max depth : 5423
bestlag removed 19 spikes
lag: 44 correlation: 0.96739
bestlag removed 15 spikes
lag: 44 correlation: 0.97073
bestlag removed 18 spikes
lag: 44 correlation: 0.96104
lag: 44 correlation: 0.95893
bestlag removed 18 spikes
lag: 44 correlation: 0.96955
bestlag removed 11 spikes
lag: 44 correlation: 0.96439
bestlag removed 19 spikes
lag: 44 correlation: 0.97234
bestlag removed 23 spikes
lag: 44 correlation: 0.97272
bestlag removed 18 spikes
lag: 44 correlation: 0.97301
bestlag removed 40 spikes
lag: 44 correlation: 0.98803
bestlag removed 20 spikes
lag: 44 correlation: 0.98439
bestlag removed 22 spikes
lag: 44 correlation: 0.97834
bestlag removed 16 spikes
lag: 44 correlation: 0.98556
bestlag removed 17 spikes
lag: 44 correlation: 0.98228
bestlag removed 24 spikes
lag: 44 correlation: 0.98924
bestlag removed 28 spikes
lag: 44 correlation: 0.98671
bestlag removed 22 spikes
lag: 44 correlation: 0.98135
bestlag removed 16 spikes
lag: 44 correlation: 0.9798
bestlag removed 21 spikes
lag: 44 correlation: 0.98009
bestlag removed 25 spikes
lag: 44 correlation: 0.97962
bestlag removed 16 spikes
lag: 44 correlation: 0.97595
bestlag removed 16 spikes
lag: 44 correlation: 0.9708
bestlag removed 18 spikes
lag: 44 correlation: 0.9747
bestlag removed 19 spikes
lag: 44 correlation: 0.97127
bestlag removed 15 spikes
lag: 44 correlation: 0.96664
median lag 44
most popular lag 44
best correlated lag 44
BESTTLAG: lag is: 44 for which 104% of 24 lags agree
best lag W: 44 CTD scans ~ -44 seconds corr:0.98924
adjust ADCP time to CTD time and shift depth record
download SW routines to get more accurate soundspeed and N^2
made CTD profile from time series data.
update start pos from:NaN°N NaN' NaN°E NaN'
to:57°S 59.9724' 30°E 0.0684'
update end pos from:NaN°N NaN' NaN°E NaN'
to:58°S 0.4884' 30°E 0.4673'
==> STEP 6 TOOK 1.6 seconds
################ [029] step 7: FIND SURFACE & SEA BED ################
GETDPTHI: Depth from vertical velocity inverse method
starting run 1 to get LADCP depth
make soundspeed based on CTD pressure and temp
correct velocities for sound speed
use CTD time series depth, will not do depth inversion
LADCP minus CTD depth mean: -25.9165 std: 20.2066
maximum depth from int W is :5423
should be :5423
[Warning: Polynomial is badly conditioned. Add points with distinct X values, reduce
the degree of the polynomial, or try centering and scaling as described in HELP
POLYFIT.]
[> In polyfit (line 79)
In getdpthi (line 349)
In process_cast (line 298)]
bottom found at 5433 +/- 0 m
correct bin length for sound speed
removing 2927 values below bottom
==> STEP 7 TOOK 1.2 seconds
################ [029] step 8: APPLY PITCH/ROLL CORRECTIONS ##########
==> STEP 8 TOOK 0.0 seconds
################ [029] step 9: EDIT DATA #############################
bin masking : masking uplooker bin 1 because of zero blanking distance
bin masking : masking downlooker bin 1 because of zero blanking distance
bin masking : set 198047 weights to NaN
side-lobe contamination : set 476 weights to NaN
==> STEP 9 TOOK 0.4 seconds
################ [029] step 10: FORM SUPER ENSEMBLES #################
PREPINV: prepare data for inversion, form Super-Ensembles
average profiles over (p.avdz) 7.7497 meter
reduce weight for larger tilts 0.5 at 10 degree
mean heading offset from compasses = -0.71933 deg
mean heading offset from pitch/roll = -3.8124 deg
rot up2down use mean up/down compass
will not correct for sound speed twice
remove reference velocity and average ensembles
Outlier discarded 0 bins down looking
Outlier discarded 1046 bins up looking
found 3 bottom track std==0 set to 0.1 m/s
found 29 finite bottom track ensembles
discarded 1 bottom tracks velocities because of wstd > 0.20888
removed 4 non finite super ensembles
set 96 weight values to nan because super ensemble std =0
set 1806 values to minimum super ensemble std 0.071673
reduced profile length = 1070 super-ensemble bins
==> STEP 10 TOOK 3.5 seconds
################ [029] step 11: REMOVE SUPER-ENSEMBLE OUTLIERS #######
==> STEP 11 TOOK 2.3 seconds
################ [029] step 12: RE-FORM SUPER ENSEMBLES ##############
mean heading offset from pitch/roll = -3.8124 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 0 bins down looking
Outlier discarded 1046 bins up looking
found 3 bottom track std==0 set to 0.1 m/s
found 29 finite bottom track ensembles
discarded 1 bottom tracks velocities because of wstd > 0.20888
removed 4 non finite super ensembles
set 96 weight values to nan because super ensemble std =0
set 1806 values to minimum super ensemble std 0.071673
reduced profile length = 1070 super-ensemble bins
==> STEP 12 TOOK 4.2 seconds
################ [029] step 13: (RE-)LOAD SADCP DATA #################
LOADSADCP: load SADCP data file ../Data/SADCP/I6S_SADCP.mat
found 15 SADCP profiles
==> STEP 13 TOOK 0.2 seconds
################ [029] step 14: CALCULATE INVERSE SOLUTION ###########
GETINV: compute best velocity profile
Barotropic velocity error 0.004351 [m/s]
super ensemble velocity error NaN [m/s]
set velocity error to:0.041071 [m/s]
vertical resolution (ps.dz) is 8 [m]
use super ensemble std based weights normalized by 0.041071 m/s
preaveraged GPS ships vel 1 times
remove 1173 constaints below minimum weight
smooth Ocean velocity profile
found 1 ill constrained elements will smooth
smooth CTD velocity profile
found 2 ill constrained elements will smooth
bottom inversion
28 bottom track ctd-vel weights of about : 0.51119
weight for bottom track is (ps.botfac) 1
mean sadcp weight : 1.7328
weight for SADCP vel is (ps.sadcpfac) 1
lainbaro: 4% of profile have no useful data
normalized barotropic constrain weight: 6.7699
mean individual ctd velocity weight : 483.5219
ready for inversion length of d: 14389
(CTD vel) length of A1: 1069
(ocean vel) length of A2: 675
Moore-Penrose inverse
solve only down trace
smooth Ocean velocity profile
no smoothness constraint applied
smooth CTD velocity profile
no smoothness constraint applied
Moore-Penrose inverse w/o errors
solve only up trace
smooth Ocean velocity profile
found 2 ill constrained elements will smooth
smooth CTD velocity profile
found 2 ill constrained elements will smooth
Moore-Penrose inverse w/o errors
CHECKINV all values are given in [m/s]
Velocity profile error: 0.041 should be about noise: 0.029
Check bottom track rms: 0.065 should be smaller than 0.172 / 1.000
Check SADCP rms: 0.012 should be smaller than 0.035 / 1.000
GPS-LADCP ship spd diff: 0.001 should be smaller than 0.004 / 1.000
use reference layer W
CHECKBTRK: check bottom track against U_ctd solution
profiles within give acceptable range: 448
U bias : 0.025 [m/s] std: 0.108 [m/s]
V bias : 0.022 [m/s] std: 0.120 [m/s]
W bias : 0.014 [m/s] std: 0.031 [m/s]
W slope fact :0.0081 [1/m] lower W below bottom
W diff :0.2371 [m/s] ping to ping w difference
H std : 0.8 [m] large means bottom is rough/sloped
Tilt mean :1.6 +/- 0.5 [^o]
==> STEP 14 TOOK 7.8 seconds
################ [029] step 15: CALCULATE SHEAR SOLUTION #############
CALC_SHEAR3: calculate a baroclinic velocity profile based on shears only
Averaging shear profile over 16 m intervals
Maximum allowed std within calculation intervals : 2
Data deviating more from the median will be discarded.
Correlation-derived weights range from 0.013106 to 0.98867
Removed 10 % of data with lowest weights from shear calculation.
New weights range from 0.57146 to 0.98867
Inversion average error : 0.038348 m/s
Final average error : 0.038348 m/s
==> STEP 15 TOOK 2.0 seconds
################ [029] step 16: PLOT RESULTS & SHOW WARNINGS #########
do not know calibration of this instrument make a guess:
Battery Voltage is 56.7 V
LADCP WARNINGS
** found 223 (2.0% of total) velocity measurements > 2.5 m/s
weak down looking beam 1
LADCP processing warnings:
==> STEP 16 TOOK 0.6 seconds
################ [029] step 17: SAVE OUTPUT ##########################
save results
ambiguity @ avdz @X avens avpercent d barofac barvelerr ?qU5w battery @LUUK beamangle blen_d blen_u blnk_d blnk_u botfac
bottomdist
btrk_below ? btrk_mode
btrk_range , 2
btrk_rough ?(k| btrk_tilt_mean ?X\P*
btrk_tilt_std ?6" btrk_ts
btrk_u_bias ?o V^
btrk_u_std ?o btrk_used btrk_v_bias ?V ̔
btrk_v_std ?o80M btrk_w_bias ?/g
btrk_w_std ?
btrk_wdiff ?Wp btrk_weight_nblen btrk_wlim ? btrk_wstd ?ʼ/ checkpoints checkpoints/029 cm_save cruise_id I6S ctd ../Data/CTD/029.1Hz ctd_endtime ABǹ ctd_fields_per_line ctd_header_lines 7 ctd_pressure_field ctd_salinity_field
ctd_starttime ABK ctd_temperature_field
ctd_time_base ctd_time_field ctddepth ctdmaxlag ctdmaxlagnp X
ctdmincorr ?|8s ctdprof ctdtime ctdtimoff ? cut
debug dist_d @ p
> dist_u @ W
=p dn_range @i
=p@i
=p@i
=p@i
=p down_sn down_up dragfac drot Db4g
dt_profile @߀ dz dzbelow edit_mask_dn_bins
elim ?
fig3_colormap fig3_err_y_axis fix_compass getdepth guessbottom hbot_0 ignore_beam instid Anz@ Ap interp_ctd_times interp_missing_GPS interp_nav_times
ladcp_cast
ladcp_station
ladcpdepth ladcpdo ../Data/PD0/029DL.100 ladcpr_CTD_depth_std 9
XO@@44&C