diff -r cc6c4309828a -r 2ccb81b7cea5 defaults.pl --- a/defaults.pl Wed May 25 12:14:29 2016 -0400 +++ b/defaults.pl Fri Aug 05 11:02:51 2016 -0400 @@ -1,9 +1,9 @@ #====================================================================== # D E F A U L T S . P L # doc: Tue Oct 11 17:11:21 2011 -# dlm: Tue Mar 29 07:23:24 2016 +# dlm: Mon Jun 6 22:07:19 2016 # (c) 2011 A.M. Thurnherr -# uE-Info: 74 39 NIL 0 0 72 0 2 4 NIL ofnI +# uE-Info: 179 18 NIL 0 0 72 0 2 4 NIL ofnI #====================================================================== # HISTORY: @@ -72,6 +72,10 @@ # Mar 19, 2016: - improved docu # Mar 29, 2016: - moved out dir creation to [LADCP_w_ocean] # - added opt_r support +# May 26, 2016: - added RDI_Coords::binMapping +# May 28, 2016: - added delta-residual filter (-r) +# Jun 1, 2016: - added $plotting_level +# Jun 2, 2016: - added $tilt_correction_* #====================================================================== # Output Log Files @@ -87,17 +91,35 @@ &antsCardOpt(\$opt_v,$ENV{VERB}); $opt_v = 1 unless numberp($opt_v); +#====================================================================== +# Output Plots +# - there are 3 plotting levels +# 0 : suppress all plots +# 1* : produce default plots; *DEFAULT +# 2 : produce default and diagnostic plots +# >2: produce all plots, including useless ones +#====================================================================== + +$plotting_level = 1; + #====================================================================== -# Data Input +# Input Data #====================================================================== # Set $opt_4 to 1 (or use the -4 option) to suppress 3-beam LADCP -# solutions +# solutions. This only has an effect for beam-coordinate data. #$opt_4 = 1; +# Set $RDI_Coords::binMapping to 'none' or use -d to disable +# linear bin interpolation (which is better than bin mapping). +# This only has an effect for beam-coordinate data. + +$RDI_Coords::binMapping = $opt_d ? 'none' : 'linterp'; + + # The following variables allow bias-correcting the attiude # sensors. # NB: heading is not used for vertical-velocity processing! @@ -125,10 +147,15 @@ # Data Editing #====================================================================== -# The following sets the max allowable rms residual w per ensemble; -# data from ensembles with larger rms residuals are discarded. +# The following sets the max allowable rms residual w per ensemble, as +# well as the max allowable difference between the two beam-pair +# residuals. Measurements that fail either of these tests are are +# discarded. The limiting values were chosed by inspection of +# log files and diagnostic plots of a few example profiles. In case of +# the delta-residual limit, histograms of this parameter show a very +# steep cutoff at 0.05 cm/s for both IWISE and 2016_I08S data. -&antsFloatOpt(\$opt_r,0.04); +&antsFloatOpt(\$opt_r,'0.06,0.06'); # By default, ensembles with uncertain time-lagging are discarded. @@ -147,25 +174,18 @@ &antsFloatOpt(\$opt_c,70); -# The following sets the default limit for instrument attitude -# (pitch/roll). -# -# The default value was established with IWISE profiles 004, 005, 045 -# and 049, which all show considerabe tilt-related discrepancies between -# the corresponding 2-beam solutions. The first and second pair of -# profiles were collected with 8 and 6m bins, respectively, without -# bin re-mapping. The original default of 15 degrees led to large -# beam-pair differences. Based on diagnostic plots it appears that -# only tilt angles smaller than 9 degrees or so are satisfactory. -# In case of the IWISE data set, such a tight constraint causes -# too many data gaps. The compromise of 12 degrees seems to work -# quite well, based on the p0 vs epsilon correlation across 5 -# data sets. +# Instrument Tilt # -# NB: if this default is changed, the usage message in [LADCP_w] -# needs to be updated as well. +# It is not fully clear what tilt angles are acceptable for +# obtaining good vertical velocities. Up to 2016 (V1.3) the +# default limit was 12 degrees based on an analysis of +# the 2010 IWISE data with inaccurate 2-beam transformations. +# As re-processing with a limit of 20 degrees improves the +# agreement between DL and UL data (R = 0.77/0.67 => 0.79/0.74) +# the limit was changed to 22 degrees, the same used +# in Martin Visbeck's inversion code. -&antsFloatOpt(\$opt_t,12); +&antsFloatOpt(\$opt_t,22); # The following sets the default error velocity limit; measurements @@ -224,8 +244,10 @@ # Previous Ping Interference editing as described in [edit_data.pl] # - enabled by default for WH150 data -# - the variable defines a string with a perl expression, which is -# evaluated once the data are loaded +# - PPI_seabed_editing_required defines a string with a perl expression +# that is evaluated once the data are loaded; if true, seabed PPI +# editing is enabled +# - to enable PPI editing without condition, set $PPI_editing = 1; # - 2014 CLIVAR P16 #47 has a slight discontinuity at 4000m; this # discontinuity is there without PPI filtering but gets slightly # worse with PPI filtering. Setting $PPI_extend_upper_limit to @@ -236,8 +258,9 @@ # set by the shortest acoustic path between the ADCP and the # seabed. -$PPI_editing_required = '($LADCP{BEAM_FREQUENCY} < 300)'; +$PPI_seabed_editing_required = '($LADCP{BEAM_FREQUENCY} < 300)'; +#$PPI_editing = 1; # uncomment to enable PPI always #$PPI_extend_upper_limit = 1.03; # see comments above @@ -311,7 +334,9 @@ # bin is chosen to construct a reference profile for Sv. The bin number # is automatically increased if the selected bin does not contain valid # data, i.e. the default value of 1 ensures that the closest valid bin -# is used to construct the reference profile. +# is used to construct the reference profile. The empirical correction +# causes artifacts every 100m. To disable the empirical +# correction, undefine the following variable. $Sv_ref_bin = 1;