CDF
lat lon date name zbot tim z I z_sadcp 4 BAR_ref_U YV BAR_ref_V ?z%i
BAR_ref_descr ![Navigation][Bottom-track][SADCP]
BAR_ref_error ?x<
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 ?\, GEN_LADCP_ensemble_time_std_sec ?ږw GEN_LADCP_station GEN_Magnetic_deviation_deg EU*R GEN_Matlab_version 9.3.0.713579 (R2017b) GEN_Ocean_depth_m b 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 @Z GEN_Profile_end_decimal_day @ZozC GEN_Profile_end_latitude P_+j GEN_Profile_end_longitude @> 3En GEN_Profile_max_depth_m [ GEN_Profile_start_decimal_day @ZN GEN_Profile_start_latitude P_M GEN_Profile_start_longitude @> 2| 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 @n
=p@n
=p@n
=p@n
=p LADCP_dn_btrk_u_bias y0Y LADCP_dn_btrk_u_std ?#e% LADCP_dn_btrk_v_bias b i. LADCP_dn_btrk_v_std ?,e8B 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 ?ٗs 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 @[AC LADCP_dn_xmit_pings LADCP_dn_xmit_vol @e{8rؒ LADCP_up_beam_range @Z
G{@Z
G{@Z
G{@\
G{ 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 ?J' 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 LADCP_up_xmit_pings LADCP_up_xmit_vol @cu/9 LOG_Inverse_log 2ALDEO LADCP software: Version IX_14beta2
################ [014] step 1: LOAD LADCP DATA #######################
LOADING CHECKPOINT checkpoints/014_0
RE-LOADING PER-CAST PARAMETERS
LOADRDI:
loading down-data ../Data/PD0/014DL.100
DETECTED BEAM coordinates: rotating to EARTH coordinates
31360 3-beam solutions calculated (21% of total)
read 7884 ensembles with 30 bins each
DETECTED BEAM bottom track coordinates!
0 3-beam solutions calculated (0% of total)
loading up-data ../Data/PD0/014UL.100
DETECTED BEAM coordinates: rotating to EARTH coordinates
7873 3-beam solutions calculated (11% of total)
read 7886 ensemble and 25 bins
WARNING WARNING WARNING WARNING WARNING WARNING WARNING WARNING
Warning: non-constant ping rate in downlooker data (staggered pinging?)
min down ping rate :1.35 max down ping rate :1.61
find best time match of up-looking ADCP to down looking ADCP
up instrument is different by 1 ensembles
WARNING WARNING WARNING WARNING WARNING WARNING WARNING WARNING
bestlag removed 288 spikes
try to shift timeseries by lag: 1 correlation: 0.99073
shift ADCP timeseries by 1 ensembles
number of joint ensembles is : 7883
use target strength of up looking to find surface
found 504 finite RDI bottom velocities
removed 16 values because of error velocity > 0.5 m/s
removed 192 values because of horizontal speed > 2.5 m/s
** found 46 (0.6% 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 51 profiles due to tilt derivative > 4 degrees
WARNING one beam might be weak !
weak down looking beam 1
Outlier discarded 8011 bins down looking
Outlier discarded 40 bottom track
Outlier discarded 5053 bins up looking
==> STEP 1 TOOK 7.7 seconds
SAVING CHECKPOINT checkpoints/014_1
################ [014] step 2: FIX LADCP-DATA PROBLEMS ###############
==> STEP 2 TOOK 0.0 seconds
################ [014] step 3: LOAD GPS DATA #########################
LOADNAV: load NAV time series ../Data/CTD/014.1Hz
number of NAV scans: 10278 delta t : 0.99999 seconds
corrected for magnetic declination of -42.7 deg
==> STEP 3 TOOK 0.1 seconds
################ [014] step 4: GET BOTTOM-TRACK DATA #################
GETBTRACK creates own bottom track in addition to RDI
in: p.btrk_mode 3 and p.btrk_used 1
using increased bottom echo amplitudes to create bottom track
localmax2: found 4241 valid values
found 88 bottom depth below btrk_range 50
use 0.5 bins below maximum target strength for own bottom track velocity
created 2058 bottom distances keeping original
removed 476 bottom track profiles W_btrk - W_ref difference > 0.05
boutlier removed 96 bottom track velocities
created 1244 bottom track velocities keeping original
out: p.btrk_mode 3 and p.btrk_used 1
==> STEP 4 TOOK 0.4 seconds
################ [014] step 5: LOAD CTD PROFILE ######################
==> STEP 5 TOOK 0.0 seconds
################ [014] step 6: LOAD CTD TIME SERIES ##################
LOADCTD: load CTD time series ../Data/CTD/014.1Hz
read 10278 CTD scans; median delta_t = 1.00 seconds
interpolated to 10278 CTD scans; delta_t = 1.00 seconds
90% LADCP depth 4219 at 18-Apr-2019 21:52:05
90% CTD pressure 4319 at 18-Apr-2019 21:38:11
Changed Start Time : 18-Apr-2019 20:22:02 to 18-Apr-2019 20:39:05
Changed End Time : 18-Apr-2019 23:33:47 to 18-Apr-2019 23:24:49
median CTD time difference 0.99999 s
removed 0 pressure spikes
use 1 time base for W_ctd
CTD max depth : 4699
bestlag removed 17 spikes
lag: 14 correlation: 0.96883
bestlag removed 19 spikes
lag: 14 correlation: 0.97822
bestlag removed 17 spikes
lag: 14 correlation: 0.97683
bestlag removed 23 spikes
lag: 14 correlation: 0.98047
bestlag removed 25 spikes
lag: 14 correlation: 0.98376
bestlag removed 26 spikes
lag: 14 correlation: 0.98282
bestlag removed 34 spikes
lag: 14 correlation: 0.9834
bestlag removed 16 spikes
lag: 14 correlation: 0.97699
bestlag removed 13 spikes
lag: 14 correlation: 0.98511
bestlag removed 18 spikes
lag: 14 correlation: 0.98514
bestlag removed 20 spikes
lag: 14 correlation: 0.97835
bestlag removed 20 spikes
lag: 14 correlation: 0.97455
bestlag removed 26 spikes
lag: 14 correlation: 0.97246
bestlag removed 19 spikes
lag: 14 correlation: 0.96967
bestlag removed 18 spikes
lag: 14 correlation: 0.97953
bestlag removed 16 spikes
lag: 14 correlation: 0.97988
bestlag removed 15 spikes
lag: 14 correlation: 0.97222
bestlag removed 17 spikes
lag: 14 correlation: 0.97395
bestlag removed 22 spikes
lag: 14 correlation: 0.98149
median lag 14
most popular lag 14
best correlated lag 14
BESTTLAG: lag is: 14 for which 106% of 18 lags agree
best lag W: 14 CTD scans ~ -14 seconds corr:0.98514
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:65°S 29.9652' 30°E 0.0464'
update end pos from:NaN°N NaN' NaN°E NaN'
to:65°S 29.9640' 30°E 0.0468'
==> STEP 6 TOOK 1.3 seconds
################ [014] 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: -16.636 std: 10.5041
maximum depth from int W is :4699
should be :4699
[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 4706 +/- 1 m
correct bin length for sound speed
removing 2155 values below bottom
==> STEP 7 TOOK 1.2 seconds
################ [014] step 8: APPLY PITCH/ROLL CORRECTIONS ##########
==> STEP 8 TOOK 0.0 seconds
################ [014] 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 129912 weights to NaN
side-lobe contamination : set 262 weights to NaN
==> STEP 9 TOOK 0.4 seconds
################ [014] step 10: FORM SUPER ENSEMBLES #################
PREPINV: prepare data for inversion, form Super-Ensembles
average profiles over (p.avdz) 7.7059 meter
reduce weight for larger tilts 0.5 at 10 degree
mean heading offset from compasses = -1.0114 deg
mean heading offset from pitch/roll = 3.0768 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 781 bins up looking
found 3 bottom track std==0 set to 0.1 m/s
found 35 finite bottom track ensembles
discarded 2 bottom tracks velocities because of wstd > 0.22487
removed 4 non finite super ensembles
set 154 weight values to nan because super ensemble std =0
set 1748 values to minimum super ensemble std 0.077539
reduced profile length = 923 super-ensemble bins
==> STEP 10 TOOK 2.9 seconds
################ [014] step 11: REMOVE SUPER-ENSEMBLE OUTLIERS #######
==> STEP 11 TOOK 1.7 seconds
################ [014] step 12: RE-FORM SUPER ENSEMBLES ##############
mean heading offset from pitch/roll = 3.0768 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 781 bins up looking
found 3 bottom track std==0 set to 0.1 m/s
found 35 finite bottom track ensembles
discarded 2 bottom tracks velocities because of wstd > 0.22487
removed 4 non finite super ensembles
set 154 weight values to nan because super ensemble std =0
set 1748 values to minimum super ensemble std 0.077539
reduced profile length = 923 super-ensemble bins
==> STEP 12 TOOK 3.4 seconds
################ [014] step 13: (RE-)LOAD SADCP DATA #################
LOADSADCP: load SADCP data file ../Data/SADCP/I6S_SADCP.mat
found 11 SADCP profiles
==> STEP 13 TOOK 0.2 seconds
################ [014] step 14: CALCULATE INVERSE SOLUTION ###########
GETINV: compute best velocity profile
Barotropic velocity error 0.0060338 [m/s]
super ensemble velocity error NaN [m/s]
set velocity error to:0.046463 [m/s]
vertical resolution (ps.dz) is 8 [m]
use super ensemble std based weights normalized by 0.046463 m/s
preaveraged GPS ships vel 1 times
remove 1219 constaints below minimum weight
smooth Ocean velocity profile
found 1 ill constrained elements will smooth
smooth CTD velocity profile
no smoothness constraint applied
bottom inversion
33 bottom track ctd-vel weights of about : 0.53125
weight for bottom track is (ps.botfac) 1
mean sadcp weight : 1.4161
weight for SADCP vel is (ps.sadcpfac) 1
lainbaro: 2% of profile have no useful data
normalized barotropic constrain weight: 6.783
mean individual ctd velocity weight : 405.7722
ready for inversion length of d: 8592
(CTD vel) length of A1: 922
(ocean vel) length of A2: 585
Moore-Penrose inverse
solve only down trace
smooth Ocean velocity profile
found 1 ill constrained elements will smooth
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
no smoothness constraint applied
Moore-Penrose inverse w/o errors
CHECKINV all values are given in [m/s]
Velocity profile error: 0.047 should be about noise: 0.033
Check bottom track rms: 0.060 should be smaller than 0.150 / 1.000
Check SADCP rms: 0.009 should be smaller than 0.035 / 1.000
GPS-LADCP ship spd diff: 0.001 should be smaller than 0.006 / 1.000
use reference layer W
CHECKBTRK: check bottom track against U_ctd solution
profiles within give acceptable range: 415
U bias :-0.006 [m/s] std: 0.114 [m/s]
V bias :-0.002 [m/s] std: 0.106 [m/s]
W bias : 0.000 [m/s] std: 0.039 [m/s]
W slope fact :0.0063 [1/m] lower W below bottom
W diff :0.1908 [m/s] ping to ping w difference
H std : 4.5 [m] large means bottom is rough/sloped
Tilt mean :1.3 +/- 0.5 [^o]
==> STEP 14 TOOK 6.7 seconds
################ [014] 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.063549 to 0.9957
Removed 9 % of data with lowest weights from shear calculation.
New weights range from 0.58894 to 0.9957
Inversion average error : 0.044549 m/s
> Increasing error estimate because of elevated shear - inverse difference
> by a factor of 2.032
> std of difference between regular and shear profile : 0.13579 m/s
Final average error : 0.090525 m/s
==> STEP 15 TOOK 1.4 seconds
################ [014] step 16: PLOT RESULTS & SHOW WARNINGS #########
do not know calibration of this instrument make a guess:
Battery Voltage is 56.7 V
LADCP WARNINGS
weak down looking beam 1
LADCP processing warnings:
> Increasing error estimate because of elevated shear - inverse difference
==> STEP 16 TOOK 0.6 seconds
################ [014] step 17: SAVE OUTPUT ##########################
save results
ambiguity @ avdz @Ѳ avens avpercent d barofac barvelerr ?x< battery @L[ beamangle blen_d blen_u blnk_d blnk_u botfac
bottomdist
btrk_below ? btrk_mode
btrk_range , 2
btrk_rough @ btrk_tilt_mean ?p>5
btrk_tilt_std ?ు*1 btrk_ts
btrk_u_bias y0Y
btrk_u_std ?#e% btrk_used btrk_v_bias b i.
btrk_v_std ?,e8B btrk_w_bias ?4կx
v
btrk_w_std ?^
btrk_wdiff ?l h btrk_weight_nblen btrk_wlim ? btrk_wstd ?ȁi checkpoints checkpoints/014 cm_save cruise_id I6S ctd ../Data/CTD/014.1Hz ctd_endtime AB|߬ ctd_fields_per_line ctd_header_lines 7 ctd_pressure_field ctd_salinity_field
ctd_starttime ABn$V ctd_temperature_field
ctd_time_base ctd_time_field ctddepth ctdmaxlag ctdmaxlagnp X
ctdmincorr ?} ctdprof ctdtime ctdtimoff ?ö cut
debug dist_d @ p
> dist_u @ W
=p dn_range @n
=p@n
=p@n
=p@n
=p down_sn down_up dragfac drot EU*R
dt_profile @l 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/014DL.100 ladcpr_CTD_depth_std 0XO@%%i ladcpup ../Data/PD0/014UL.100 lat P_U\R lon @> 2
magdec_source maxbinrange maxdepth @[!) name I6S station #14 (V5) nav ../Data/CTD/014.1Hz nav_end P@ =ȴ9W@> ? nav_error nav_fields_per_line nav_header_lines 7
nav_lat_field
nav_lon_field nav_start P@ =Xp@> ?êeP
nav_time_base nav_time_field navdata
navtime_av ?Vll nbin_d nbin_u nbins 7 nping_total nt
offsetup2down orig outlier outlier_n
oversample pg_save pglim pose P@ =ȴ9W@> ? poss P@ =Xp@> ?êeP res V5/014
rotup2down sadcp ../Data/SADCP/I6S_SADCP.mat
sadcp_dtok sadcpfac serial_cpu_d " , serial_cpu_u N 8 shear single_ping_accuracy smallfac smoofac software 'LDEO LADCP software: Version IX_14beta2 solve soundcorr
std_weight superens_std_min ?ٗs surfdist tilt_weight
tiltcor tiltmax time_end 2
time_start ' timoff tint @^ٜ@^(F( ts_att_dn ?+ ts_att_up ?Q ts_save
ts_signal_min ucorr ?߄Mn up_dn_comp_off .ܪ up_dn_looker
up_dn_pit_off ԁ up_dn_pit_rol_comp_off @O;dc
up_dn_rol_off ?]/ up_range @Z
G{@Z
G{@Z
G{@\
G{ up_sn uship >sEv vcorr ?r velerr ?5e vlim @ vship ?-K!'R warn )LADCP WARNINGS
weak down looking beam 1 warnings i LADCP processing warnings:
> Increasing error estimate because of elevated shear - inverse difference
warnp i LADCP processing warnings:
> Increasing error estimate because of elevated shear - inverse difference wbslope ?y?AHO?3>* weighbin1 weightmin ? weightpower whoami A.M. Thurnherr wizr
wlim ?ə xdisp ?Щw
xmc @[AC@f xmv @e{8rؒ@cu/9 xmv_min ydisp @yH zbottom @b49 zbottomerror ?$\a0 zpar @#6@[!)@#x!U} 2 lat long_name Latitude units Degree North m lon long_name Longitude units Degree East n date long_name Date units Y M D H M S n name long_name Cast ID n zbot long_name Bottom Referenced Profile Depth units m \ n0 tim long_name Station Time Series units Julian Days n z long_name Depth units Meters $ \ z_sadcp long_name SADCP Profile Depth units m ctd_s long_name CTD profile salinity units psu $ P ctd_t long_name CTD profile temperature units Degree C $ t ensemble_vel_err long_name ADCP ensemble velocity error units m/s $ nvel long_name !LADCP number of ensembles per bin p long_name Pressure units dBar $ P range long_name ADCP total range of data units m $ t range_do long_name ADCP down looking range of data units m $ ǘ range_up long_name ADCP up looking range of data units m $ м shiplat long_name Latitude units Degree North h shiplon long_name Longitude units Degree East h H tim_hour long_name Station Time Series units Hour of Day h ts long_name !ADCP echo amplitude profile bin 1 units dB $ ts_out long_name )ADCP echo amplitude profile last down bin units dB $ <