author | A.M. Thurnherr <athurnherr@yahoo.com> |
Sat, 24 Jul 2021 10:35:41 -0400 | |
changeset 56 | 8f120b9f795a |
parent 54 | 828e5466391b |
child 57 | 69e39fcb7f41 |
permissions | -rw-r--r-- |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
1 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
2 |
# D E F A U L T S . P L |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
3 |
# doc: Tue Oct 11 17:11:21 2011 |
56 | 4 |
# dlm: Fri Jul 9 13:30:48 2021 |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
5 |
# (c) 2011 A.M. Thurnherr |
56 | 6 |
# uE-Info: 360 33 NIL 0 0 72 0 2 4 NIL ofnI |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
7 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
8 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
9 |
# HISTORY: |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
10 |
# Oct 11, 2011: - created |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
11 |
# Oct 12, 2011: - added $processing_param_file |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
12 |
# Oct 13, 2011: - added $out_basename, $TL_out, $TL_hist_out |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
13 |
# Oct 14, 2011: - added $w_out, $profile_out |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
14 |
# - renamed _out to out_; out_basename to out_basename |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
15 |
# Oct 17, 2011: - added {data,plot,log}_subdir |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
16 |
# - added $out_BR |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
17 |
# - adapted to new filter-plot scripts |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
18 |
# - changed -v default to 2 |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
19 |
# Oct 19, 2011: - added SS_max_allowed_range & renamed _min_ |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
20 |
# Oct 20, 2011: - added $out_timeseries default |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
21 |
# - added $per_bin_valid_frac_lim |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
22 |
# Oct 26, 2011: - added $first_guess_timelag |
6 | 23 |
# Oct 27, 2011: - modified ProcessingParam file loading |
24 |
# - added ${pitch,roll,heading}_bias |
|
10
7e76fdef8953
version prior to carrying out separate dc-/uc- time lagging
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
6
diff
changeset
|
25 |
# Oct 11, 2012: - added .TL output to defaults |
11 | 26 |
# Oct 15, 2012: - removed support for TLhist |
13 | 27 |
# Apr 22, 2013: - removed option variable aliases |
28 |
# May 14, 2013: - opt_m => w_max_lim |
|
29 |
# - BUG: default processing params file was still .default |
|
30 |
# May 15, 2013: - changed default profile plot to show 2-beam solutions |
|
31 |
# - BUG: out_TLhist entry was still there |
|
32 |
# May 16, 2013: - -a => -d |
|
33 |
# Jun 5, 2013: - renamed $discard_data_from_beam to $bad_beam |
|
34 |
# Jun 5, 2013: - made ProcessingParams (without .profiles) default file |
|
35 |
# Sep 5, 2013: - also allow ProcessingParams.default |
|
36 |
# - added LWplot_spec to default output |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
37 |
# May 20, 2014: - added support for $PPI_editing |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
38 |
# May 21, 2014: - added $PPI_extend_upper_limit |
17
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
39 |
# Oct 15, 2014: - investigated, modified and documented -t default |
18 | 40 |
# Oct 27, 2014: - removed CTD-acceleration-effects and w spectral plots |
41 |
# - removed time-lagging stats output file by default |
|
42 |
# Oct 31, 2014: - re-arranged order of things |
|
43 |
# - .w => .samp output |
|
44 |
# Nov 4, 2014: - BUG: PPI_editing did not work as advertised |
|
25
bd38f8bfb8e6
first step at re-designing plotting system
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
24
diff
changeset
|
45 |
# Apr 16, 2015: - turned output specifies into lists (re-design of |
bd38f8bfb8e6
first step at re-designing plotting system
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
24
diff
changeset
|
46 |
# plotting sub-system) |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
47 |
# - croak -> error |
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
48 |
# - added $SS_use_BT, $SS_min_signal, $SS_min_samp |
29 | 49 |
# Apr 20, 2015: - reduced value of $SS_min_allowed_range |
50 |
# - added $Sv_ref_bin |
|
51 |
# Apr 21: 2015: - BUG: typo in $Sv_ref_bin |
|
52 |
# - decreased default verbosity |
|
53 |
# May 15, 2015: - added $min_valid_vels |
|
54 |
# May 20, 2015: - STN -> PROF |
|
30 | 55 |
# Jul 26, 2015: - began adaptation to libGMT.pl |
29 | 56 |
# - changed .prof output .wprof |
57 |
# - -v docu was wrong |
|
30 | 58 |
# - added $outGrid_firstBin, $outGrid_lastBin |
59 |
# Jul 28, 2015: - implemented new plotting system |
|
60 |
# Jul 29, 2015: - implemented new plotting system |
|
31 | 61 |
# Sep 3, 2015: - renamed wsamp output and plot |
62 |
# - changed out_w to out_wsamp |
|
32
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
63 |
# Sep 26, 2015: - added sidelobe editing params |
33 | 64 |
# Oct 13, 2015: - addded support for $ENV{VERB} |
65 |
# Jan 4, 2016: - decreased default vertical resolution to 20m |
|
34 | 66 |
# Jan 22, 2016: - changed outGrid_firstBin default to 1 |
67 |
# Jan 26, 2016: - removed -d |
|
68 |
# - changed outGrid_firstBin default to '*', also lastBin |
|
69 |
# Jan 27, 2016: - added documentation |
|
39 | 70 |
# Mar 16, 2016: - added auto creation of output directory |
42 | 71 |
# Mar 18, 2016: - added comments about -l |
72 |
# Mar 19, 2016: - improved docu |
|
73 |
# Mar 29, 2016: - moved out dir creation to [LADCP_w_ocean] |
|
74 |
# - added opt_r support |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
75 |
# May 26, 2016: - added RDI_Coords::binMapping |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
76 |
# May 28, 2016: - added delta-residual filter (-r) |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
77 |
# Jun 1, 2016: - added $plotting_level |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
78 |
# Jun 2, 2016: - added $tilt_correction_* |
48 | 79 |
# Jun 6, 2016: - removed $tilt_correction_* |
80 |
# Jul 12, 2016: - updated docu on acoustic backscatter |
|
81 |
# Jul 31, 2016: - BUG: -d to disable bin mapping did not work, because |
|
82 |
# this file is read before the options are parsed |
|
83 |
# Aug 5, 2016: - updated header |
|
84 |
# Dec 22, 2016: - added $opt_p |
|
49 | 85 |
# Nov 27, 2017: - added @valid_ensemble_range |
86 |
# Nov 29, 2017: - replaced opt_i by initial_time_lag |
|
87 |
# Apr 24, 2018: - added $water_depth_db_cmd |
|
88 |
# May 2, 2018: - added max_hspeed |
|
89 |
# - replaced $PPI_seabed_editing_required by &PPI_seabed_editing_required |
|
54 | 90 |
# May 16, 2020: - updated doc |
56 | 91 |
# Jun 30, 2021: - ditto |
92 |
# Jul 9, 2021: - added $layer_residuals_rms_max |
|
93 |
# HISTORY END |
|
42 | 94 |
|
95 |
#====================================================================== |
|
96 |
# Output Log Files |
|
97 |
# - there are 4 verbosity levels, selected by -v |
|
98 |
# 0 : errors |
|
99 |
# 1*: UNIX-like (warnings and info messages that are not produced for every cast; *DEFAULT) |
|
100 |
# 2 : progress messages and useful information |
|
101 |
# >2: debug messges |
|
102 |
# - the most useful ones of these are 1 & 2 |
|
103 |
# - verbosity level can be set with the VERB shell variable |
|
104 |
#====================================================================== |
|
105 |
||
106 |
&antsCardOpt(\$opt_v,$ENV{VERB}); |
|
107 |
$opt_v = 1 unless numberp($opt_v); |
|
108 |
||
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
109 |
#====================================================================== |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
110 |
# Output Plots |
54 | 111 |
# - there are 4 plotting levels |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
112 |
# 0 : suppress all plots |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
113 |
# 1* : produce default plots; *DEFAULT |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
114 |
# 2 : produce default and diagnostic plots |
54 | 115 |
# 3 : produce debugging diagnostic plots |
116 |
# >3: produce all plots, including useless ones |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
117 |
#====================================================================== |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
118 |
|
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
119 |
$plotting_level = 1; |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
120 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
121 |
#====================================================================== |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
122 |
# Input Data |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
123 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
124 |
|
49 | 125 |
# The two elements in the @valid_ensemble_range array limit the minimum |
126 |
# and maximum ensemble numbers considered during processing. This is |
|
127 |
# useful primarily for files with lots of on-deck data. |
|
128 |
||
129 |
# @valid_ensemble_range = (3000,10000) |
|
130 |
||
34 | 131 |
# Set $opt_4 to 1 (or use the -4 option) to suppress 3-beam LADCP |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
132 |
# solutions. This only has an effect for beam-coordinate data. |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
133 |
|
13 | 134 |
#$opt_4 = 1; |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
135 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
136 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
137 |
# Set $RDI_Coords::binMapping to 'none' or use -d to disable |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
138 |
# linear bin interpolation (which is better than bin mapping). |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
139 |
# This only has an effect for beam-coordinate data. |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
140 |
|
48 | 141 |
# $RDI_Coords::binMapping = 'none'; |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
142 |
|
34 | 143 |
# The following variables allow bias-correcting the attiude |
144 |
# sensors. |
|
24 | 145 |
# NB: heading is not used for vertical-velocity processing! |
6 | 146 |
|
34 | 147 |
$pitch_bias = 0; |
148 |
$roll_bias = 0; |
|
149 |
$heading_bias = 0; |
|
6 | 150 |
|
151 |
||
34 | 152 |
# The following variable defines the minimum valid velocities |
153 |
# required in a LADCP file. If there are fewer data, an |
|
154 |
# error is produced |
|
29 | 155 |
|
156 |
$min_valid_vels = 50; |
|
157 |
||
158 |
||
34 | 159 |
# The -b option defines the range of bins to use in w calculations. |
160 |
# The '*' indicates the last bin in the ADCP file. For data |
|
161 |
# collected with non-zero blanking distance, -b '1,*' should |
|
162 |
# likely be used. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
163 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
164 |
$opt_b = '2,*' unless defined($opt_b); |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
165 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
166 |
#====================================================================== |
42 | 167 |
# Data Editing |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
168 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
169 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
170 |
# The following sets the max allowable rms residual w per ensemble, as |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
171 |
# well as the max allowable difference between the two beam-pair |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
172 |
# residuals. Measurements that fail either of these tests are are |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
173 |
# discarded. The limiting values were chosed by inspection of |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
174 |
# log files and diagnostic plots of a few example profiles. In case of |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
175 |
# the delta-residual limit, histograms of this parameter show a very |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
176 |
# steep cutoff at 0.05 cm/s for both IWISE and 2016_I08S data. |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
177 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
178 |
&antsFloatOpt(\$opt_r,'0.06,0.06'); |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
179 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
180 |
|
42 | 181 |
# By default, ensembles with uncertain time-lagging are discarded. |
182 |
# This allows profiles with dropped CTD scans to be processed without |
|
183 |
# manual intervention. For profiles collected in very calm conditions |
|
184 |
# (e.g. near the ice off Antarctica) time lagging is highly uncertain |
|
185 |
# most of the time --- setting $opt_l = 1 disables the lime-lagging |
|
186 |
# filter for those cases. |
|
30 | 187 |
|
42 | 188 |
# $opt_l = 1; |
30 | 189 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
190 |
|
34 | 191 |
# The following sets the default correlation limit; measurements with |
192 |
# correlations below this limit are discarded. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
193 |
|
13 | 194 |
&antsFloatOpt(\$opt_c,70); |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
195 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
196 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
197 |
# Instrument Tilt |
17
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
198 |
# |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
199 |
# It is not fully clear what tilt angles are acceptable for |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
200 |
# obtaining good vertical velocities. Up to 2016 (V1.3) the |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
201 |
# default limit was 12 degrees based on an analysis of |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
202 |
# the 2010 IWISE data with inaccurate 2-beam transformations. |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
203 |
# As re-processing with a limit of 20 degrees improves the |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
204 |
# agreement between DL and UL data (R = 0.77/0.67 => 0.79/0.74) |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
205 |
# the limit was changed to 22 degrees, the same used |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
206 |
# in Martin Visbeck's inversion code. |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
207 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
208 |
&antsFloatOpt(\$opt_t,22); |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
209 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
210 |
|
34 | 211 |
# The following sets the default error velocity limit; measurements |
212 |
# with error velocities below this limit are discarded. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
213 |
|
13 | 214 |
&antsFloatOpt(\$opt_e,0.1); |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
215 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
216 |
|
34 | 217 |
# The following variable allows editing the velocities farthest |
218 |
# from the transducer. It defines how many velocities are to be |
|
219 |
# removed from each ensemble. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
220 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
221 |
$truncate_farthest_valid_bins = 0; |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
222 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
223 |
|
34 | 224 |
# The following variable allows editing all data from a given |
225 |
# beam. Set to 1-4 to enable. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
226 |
|
13 | 227 |
$bad_beam = 0; |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
228 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
229 |
|
34 | 230 |
# The following sets the maximum gap length in the w time series that |
231 |
# is simply ignored. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
232 |
|
13 | 233 |
&antsFloatOpt(\$opt_g,60); |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
234 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
235 |
|
34 | 236 |
# The following variable sets the max allowed vertical ocean velocity |
237 |
# in m/s. Measurements with |w| this limit are discarded. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
238 |
|
13 | 239 |
$w_max_lim = 1; |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
240 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
241 |
|
34 | 242 |
# In each ensemble, vertical velocities differing more than this |
243 |
# parameter times the mean absolute deviation from the median, are |
|
244 |
# considered outliers and removed. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
245 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
246 |
$per_ens_outliers_mad_limit = 2; |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
247 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
248 |
|
34 | 249 |
# Data from bins with less valid velocities than the following parameter |
250 |
# are considered bad and removed. It is not clear whether this really |
|
251 |
# makes sense, but this editing is likely safe because it only affects |
|
252 |
# ensebles with the largest ranges. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
253 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
254 |
$per_bin_valid_frac_lim = 0.15; |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
255 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
256 |
|
34 | 257 |
# All ensembles recorded when the CTD is shallower than |
258 |
# the following parameter (depth in meters) are discarded. |
|
259 |
# Possible contamination: ship's hull, thrusters, bubble clouds |
|
260 |
# Inspired by 2011_IWISE station 8. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
261 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
262 |
$surface_layer_depth = 25; |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
263 |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
264 |
|
49 | 265 |
# Water depth is important for precious ping interference editing |
266 |
# (see below) and for setting the height-above bottom field. |
|
267 |
# - by default, water depth for dowwnward-facing ADCPs is determined |
|
268 |
# from the seabed echo return |
|
269 |
# - when water depth is set explicitly either via the $water_depth or |
|
270 |
# the $opt_h variable no search for the seabed is done |
|
271 |
# - the variable $water_depth_db_cmd can be set to the name of an |
|
272 |
# external command, which is used to get nominal water depth |
|
273 |
# from a data base if there is no other water depth information |
|
274 |
# (from echo return or supplied by user). The command will be |
|
275 |
# called with longitude and latitude as the only arguments and |
|
276 |
# is expected to return the water depth in meters. |
|
277 |
||
278 |
#$water_depth = 2048; # uncomment to set water depth to 2048m |
|
279 |
#$opt_h = 2048; # uncomment to set water depth to 2048m |
|
280 |
#$water_depth_db_cmd = 'waterdepth'; # uncomment to use 'waterdepth' command to get water depth |
|
281 |
||
282 |
||
34 | 283 |
# Previous Ping Interference editing as described in [edit_data.pl] |
49 | 284 |
# - enabled by default seabed editing of WH150 data but nothing else |
285 |
# - PPI_seabed_editing_required is a function that is called |
|
286 |
# once the data are loaded for the downlooker only; if it |
|
287 |
# returns true, seabed PPI editing is enabled |
|
288 |
# - PPI_surface_editing_required is a function that is called |
|
289 |
# once the data are loaded for the uplooker only; if it |
|
290 |
# returns true, sea surface PPI editing is enabled |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
291 |
# - 2014 CLIVAR P16 #47 has a slight discontinuity at 4000m; this |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
292 |
# discontinuity is there without PPI filtering but gets slightly |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
293 |
# worse with PPI filtering. Setting $PPI_extend_upper_limit to |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
294 |
# 1.03-1.05 partially removes the discontinuity but the profile |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
295 |
# never gets better than the profile wihtout PPI editing. Note |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
296 |
# the only reason why the upper PPI should be extended is if the |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
297 |
# recorded ping intervals are inaccurate as the upper limit is |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
298 |
# set by the shortest acoustic path between the ADCP and the |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
299 |
# seabed. |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
300 |
|
49 | 301 |
sub PPI_seabed_editing_required() |
302 |
{ |
|
303 |
# return 1; # uncomment to enable unconditional PPI editing |
|
304 |
return ($LADCP{BEAM_FREQUENCY} < 300); # low-frequency instruments require PPI editing |
|
305 |
} |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
306 |
|
49 | 307 |
sub PPI_surface_editing_required() |
308 |
{ |
|
309 |
return 0; # no sea surface PPI editing by default |
|
310 |
} |
|
311 |
||
312 |
#$PPI_extend_upper_limit = 1.03; # see comments above |
|
32
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
313 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
314 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
315 |
# The following variables control the "non-obvious" sidelobe editing for |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
316 |
# contamination from the seabed for the UL and from the sea surface for the |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
317 |
# DL. Tests with DoMORE-2 data (WH150 DL, WH300 UL) strongly suggest that |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
318 |
# it is not necessary to edit DL data for surface contamination. However, |
34 | 319 |
# at least for that instrument combination, UL (WH300) contamination from the |
32
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
320 |
# seabed should clearly be removed. |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
321 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
322 |
$sidelobe_editing_DL_surface = 0; |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
323 |
$sidelobe_editing_UL_seabed = 1; |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
324 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
325 |
# The following variable sets the depth for sidelobe contamination |
34 | 326 |
# from the surface. |
32
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
327 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
328 |
$vessel_draft = 6; # in meters |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
329 |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
330 |
|
49 | 331 |
# The following function, which is called after the LADCP data have been |
332 |
# read, must return the maximum horizontal reference-layer |
|
333 |
# speed that is allowed. The following values are based on 2018 GO-SHIP |
|
334 |
# S4P profile #106 where the CTD rosette was dragged quickly during |
|
335 |
# the latter part of the upcast. Of course, it is possible that the |
|
336 |
# differnces between the UL and DL data could be due to tilt- |
|
337 |
# sensor differences, rather than due to instrument type. |
|
338 |
||
51
0f6d9e64cc4f
before adding dropped-scan correction to LADCP_w_CTD
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
49
diff
changeset
|
339 |
$max_hspeed_300kHz = 0.55; # m/s |
0f6d9e64cc4f
before adding dropped-scan correction to LADCP_w_CTD
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
49
diff
changeset
|
340 |
$max_hspeed_150kHz = 0.35; # m/s |
0f6d9e64cc4f
before adding dropped-scan correction to LADCP_w_CTD
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
49
diff
changeset
|
341 |
|
49 | 342 |
sub max_hspeed() |
343 |
{ |
|
344 |
if (abs($LADCP{BEAM_FREQUENCY}-300) <= 25) { # 300kHz Workhorse |
|
51
0f6d9e64cc4f
before adding dropped-scan correction to LADCP_w_CTD
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
49
diff
changeset
|
345 |
$max_hspeed = $max_hspeed_300kHz; |
49 | 346 |
} elsif (abs($LADCP{BEAM_FREQUENCY}-150) <= 25) { # 150kHz Workhorse |
51
0f6d9e64cc4f
before adding dropped-scan correction to LADCP_w_CTD
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
49
diff
changeset
|
347 |
$max_hspeed = $max_hspeed_150kHz; |
49 | 348 |
} else { |
349 |
warning(2,"unknown horizontal speed limit for this instrument frequency ($LADCP{BEAM_FREQUENCY} kHz)\n"); |
|
350 |
$max_hspeed = 9e99; |
|
351 |
} |
|
352 |
} |
|
353 |
||
56 | 354 |
|
355 |
# After the beam-pair residuals have been binned into profiles, rms |
|
356 |
# values in 5-output-bin-thick layers (200m by default) are calculated. |
|
357 |
# Vertical velocities from bins with values greater than the following |
|
358 |
# limit are set to nan in the the output. |
|
359 |
||
360 |
$layer_residuals_rms_max = 0.003; |
|
361 |
||
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
362 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
363 |
# Time Lagging |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
364 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
365 |
|
49 | 366 |
# The following variable allows specifying an initial guess for the time |
367 |
# lag between the LADCP and the CTD data. The -i option overrides any value |
|
368 |
# set in the [ProcessingParams] file. |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
369 |
|
49 | 370 |
# $initial_time_lag = 567; |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
371 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
372 |
|
34 | 373 |
# The following variables define the bins used to calculate the reference- |
374 |
# layer velocities used for time lagging. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
375 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
376 |
($refLr_firstBin,$refLr_lastBin) = (2,6); |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
377 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
378 |
|
34 | 379 |
# The -n option defines the number of windows used to calculate |
380 |
# the optimal time lag. There's one value for each time-lagging step. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
381 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
382 |
$opt_n = '10,100' unless defined($opt_n); |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
383 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
384 |
|
34 | 385 |
# The -w option defines the width of the window (in seconds) used |
386 |
# to calculate the optimal time lag. There's one value for each |
|
387 |
# time-lagging step. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
388 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
389 |
$opt_w = '240,20' unless defined($opt_w); |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
390 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
391 |
|
34 | 392 |
# The following parameters control whether the top three time lags |
393 |
# are accepted or not. |
|
394 |
# If the top 3 lags have spread greater than $TL_max_allowed_three_lag_spread |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
395 |
# (in CTD scans) they must account for at least $TL_required_timelag_top_three_fraction |
34 | 396 |
# or an error is generated. |
397 |
# Notes: |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
398 |
# - $TL_max_allowed_three_lag_spread default was initially set to 2 but found to be |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
399 |
# violated quite often during 2011_IWISE |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
400 |
# - large spread may indicate dropped CTD scans |
34 | 401 |
# - the optimum value of $TL_max_allowed_three_lag_spread may be |
402 |
# cast-duration dependent |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
403 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
404 |
$TL_max_allowed_three_lag_spread = 3; |
13 | 405 |
&antsFloatOpt(\$opt_3,0.6); |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
406 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
407 |
|
48 | 408 |
# Time lagging is carried out in "partial cast pieces". By default |
409 |
# the down- and upcasts are treated separately, which is accomplished |
|
410 |
# by setting $opt_p to '+'. In general, the $opt_p variable contains |
|
411 |
# a comma-separated list of elapsed splitting times (in minutes), |
|
412 |
# with '+' denoting the time of maximum depth. |
|
413 |
||
414 |
$opt_p = '+' unless defined($opt_p); |
|
415 |
||
416 |
||
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
417 |
#====================================================================== |
29 | 418 |
# Acoustic Backscatter and Seabed Search |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
419 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
420 |
|
29 | 421 |
# After applying the method of Deines (1999), an empirical correction |
422 |
# for Sv is applied to the data. The following variable determines which |
|
423 |
# bin is chosen to construct a reference profile for Sv. The bin number |
|
48 | 424 |
# is automatically increased if its value is less than LADCP_firstBin (-b), |
425 |
# and also if the selected bin does not contain valid |
|
29 | 426 |
# data, i.e. the default value of 1 ensures that the closest valid bin |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
427 |
# is used to construct the reference profile. The empirical correction |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
428 |
# causes artifacts every 100m. To disable the empirical |
48 | 429 |
# correction, undefine the following variable ($Sv_ref_bin = undef;) |
430 |
# |
|
431 |
# NOTE: Accoustic backscatter data in the reference bin are not |
|
432 |
# corrected beyond the method of Deines (1999). |
|
29 | 433 |
|
434 |
$Sv_ref_bin = 1; |
|
435 |
||
436 |
||
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
437 |
# Set to folloing variable to 1 to use ADCP BT data to detect seabed |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
438 |
# instead of default code based on Sv (echo amplitude). I do not know |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
439 |
# which code is better. |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
440 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
441 |
$SS_use_BT = 0; |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
442 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
443 |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
444 |
# The following variable defines the minimum Sv signal in a bin (max - min) |
34 | 445 |
# required for reliable seabed detection FROM ECHO AMPLITUDES. A limit of 40dB is |
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
446 |
# indicated based on GoM#13, where the seabed is only visible in the last |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
447 |
# bin (#25). 30dB is chosen as the default to allow for variability. |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
448 |
# This value may need to be changed for data not collected with WH300 |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
449 |
# instruments with 8m bins, and perhaps also for different types of |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
450 |
# seafloor (soft sediments). To do this, set $SS_min_signal to a small value |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
451 |
# (e.g. 10) and inspect the \@SV_rng values reported in the log files. |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
452 |
# This parameter is only used when $SS_use_BT == 0. |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
453 |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
454 |
$SS_min_signal = 30; |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
455 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
456 |
|
34 | 457 |
# Require at minimum nubmer of valid samples for seabed detection FROM ECHO |
458 |
# AMPLITUDES. Each sample is a bin with a clear seabed maximum. With a proper |
|
459 |
# setting of $SS_min_signal, the algorithm is stable even with only a single |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
460 |
# sample (GoM#13). However, a default of 3 required samples is chosen |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
461 |
# to make seabed detection less sensitive to $SS_min_signal. |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
462 |
# This parameter is only used when $SS_use_BT == 0. |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
463 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
464 |
$SS_min_samp = 3; |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
465 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
466 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
467 |
# The following numbers define the valid range of height-above bottom |
34 | 468 |
# for seabed detection FROM ECHO AMPLITUDE. For data collected with WH300 |
469 |
# instruments and 8m bins, the maximum range needs to be greater than 250m |
|
470 |
# (based on # GoM#13). |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
471 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
472 |
$SS_min_allowed_range = 0; |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
473 |
$SS_max_allowed_range = 350; |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
474 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
475 |
|
34 | 476 |
# Number of ensembles around bottom to search sabed IN BT DATA. |
477 |
# Only used with $SS_use_BT == 1. |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
478 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
479 |
$SS_search_window_halfwidth = 200; |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
480 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
481 |
|
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
482 |
# Maximum allowed distance of seabed from mode of distribution. |
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
483 |
# Only used with $SS_use_BT == 1. |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
484 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
485 |
$SS_max_allowed_depth_range = 10; |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
486 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
487 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
488 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
489 |
# Bottom Tracking |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
490 |
# - at present, the ADCP BT data are ignored, i.e. "post-processed" |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
491 |
# BT data are used. |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
492 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
493 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
494 |
# Don't look for BT-referenced velocities if package is more than $BT_max_range |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
495 |
# above seabed. This parameter is frequency dependent and the current value is |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
496 |
# appropriate (if rather high) for 300kHz Workhorse intruments. |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
497 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
498 |
$BT_max_range = 300; |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
499 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
500 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
501 |
# The code only tries to bin BT-referenced velocities if a consistent bottom |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
502 |
# is available in all 4 beams. Ensembles where the range of bin numbers where |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
503 |
# the maximum echo is found is greater than $max_BIT_bin_range_diff are rejected. |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
504 |
# In addition to flukes this also rejects ensembles collected with large |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
505 |
# instrument tilts. The value of 3 is a first guess that has not been explored. |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
506 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
507 |
$BT_max_bin_range_diff = 3; |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
508 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
509 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
510 |
# If the difference between measured vertical velocity of the seabed (i.e. |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
511 |
# the package vertical velocity referenced by the seabed) and the vertical |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
512 |
# velocity of the CTD (from dp/dt) si greater than $BT_max_w_error the current |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
513 |
# ensemble is ignored and $nBTwFlag is increased. The value of |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
514 |
# 3cm/s is taken from listBT developed on A0304 cruise. |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
515 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
516 |
$BT_max_w_error = 0.03; |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
517 |
|
42 | 518 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
519 |
#====================================================================== |
42 | 520 |
# Gridded Velocity Profile Output |
18 | 521 |
#====================================================================== |
522 |
||
523 |
||
42 | 524 |
# The -k option defines the minimum number of w samples required in each |
525 |
# vertical-velocity bin. The following sets the default value. |
|
18 | 526 |
|
42 | 527 |
&antsCardOpt(\$opt_k,20); |
18 | 528 |
|
529 |
||
42 | 530 |
# The -o option sets the output grid resolution in meters. The following |
48 | 531 |
# sets the default value. It was increased from 20m to 40m in Feb 2017 |
532 |
# for V1.3 because this value was required for the DoMORE-2 data and |
|
533 |
# also improves the profiles from 2017 P18, as well as from a recent |
|
534 |
# GoM data set provided by J. Ochoa. |
|
18 | 535 |
|
48 | 536 |
&antsFloatOpt(\$opt_o,40); |
18 | 537 |
|
538 |
||
42 | 539 |
# The following variables limit the bins used to grid w_oean |
540 |
# - in contrast to -b, the other bins are still used e.g. for BT |
|
541 |
# - values recorded in %outgrid_firstbin, %outgrid_lastbin |
|
542 |
# - values beyond range are: |
|
56 | 543 |
# - greyed out in *_bin_residuals.ps |
42 | 544 |
# - not used in *_w.ps, *_residuals.ps |
18 | 545 |
|
42 | 546 |
$outGrid_firstBin = '*'; # use $LADCP_firstBin (-b) |
547 |
$outGrid_lastBin = '*'; # use $LADCP_lastBin (-b) |
|
18 | 548 |
|
549 |
||
550 |
||
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
551 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
552 |
1; # return true |
42 | 553 |