author | A.M. Thurnherr <athurnherr@yahoo.com> |
Fri, 05 Aug 2016 11:02:51 -0400 | |
changeset 47 | 2ccb81b7cea5 |
parent 42 | f7690c7b92e0 |
child 48 | d9309804b6cf |
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 |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
4 |
# dlm: Mon Jun 6 22:07:19 2016 |
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 |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
6 |
# uE-Info: 179 18 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_* |
42 | 79 |
|
80 |
#====================================================================== |
|
81 |
# Output Log Files |
|
82 |
# - there are 4 verbosity levels, selected by -v |
|
83 |
# 0 : errors |
|
84 |
# 1*: UNIX-like (warnings and info messages that are not produced for every cast; *DEFAULT) |
|
85 |
# 2 : progress messages and useful information |
|
86 |
# >2: debug messges |
|
87 |
# - the most useful ones of these are 1 & 2 |
|
88 |
# - verbosity level can be set with the VERB shell variable |
|
89 |
#====================================================================== |
|
90 |
||
91 |
&antsCardOpt(\$opt_v,$ENV{VERB}); |
|
92 |
$opt_v = 1 unless numberp($opt_v); |
|
93 |
||
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
94 |
#====================================================================== |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
95 |
# Output Plots |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
96 |
# - there are 3 plotting levels |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
97 |
# 0 : suppress all plots |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
98 |
# 1* : produce default plots; *DEFAULT |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
99 |
# 2 : produce default and diagnostic plots |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
100 |
# >2: produce all plots, including useless ones |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
101 |
#====================================================================== |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
102 |
|
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
103 |
$plotting_level = 1; |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
104 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
105 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
106 |
#====================================================================== |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
107 |
# 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
|
108 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
109 |
|
34 | 110 |
# 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
|
111 |
# 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
|
112 |
|
13 | 113 |
#$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
|
114 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
115 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
116 |
# 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
|
117 |
# 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
|
118 |
# 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
|
119 |
|
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
120 |
$RDI_Coords::binMapping = $opt_d ? 'none' : 'linterp'; |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
121 |
|
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
122 |
|
34 | 123 |
# The following variables allow bias-correcting the attiude |
124 |
# sensors. |
|
24 | 125 |
# NB: heading is not used for vertical-velocity processing! |
6 | 126 |
|
34 | 127 |
$pitch_bias = 0; |
128 |
$roll_bias = 0; |
|
129 |
$heading_bias = 0; |
|
6 | 130 |
|
131 |
||
34 | 132 |
# The following variable defines the minimum valid velocities |
133 |
# required in a LADCP file. If there are fewer data, an |
|
134 |
# error is produced |
|
29 | 135 |
|
136 |
$min_valid_vels = 50; |
|
137 |
||
138 |
||
34 | 139 |
# The -b option defines the range of bins to use in w calculations. |
140 |
# The '*' indicates the last bin in the ADCP file. For data |
|
141 |
# collected with non-zero blanking distance, -b '1,*' should |
|
142 |
# 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
|
143 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
144 |
$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
|
145 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
146 |
#====================================================================== |
42 | 147 |
# 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
|
148 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
149 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
150 |
# 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
|
151 |
# 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
|
152 |
# 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
|
153 |
# 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
|
154 |
# 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
|
155 |
# 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
|
156 |
# 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
|
157 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
158 |
&antsFloatOpt(\$opt_r,'0.06,0.06'); |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
159 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
160 |
|
42 | 161 |
# By default, ensembles with uncertain time-lagging are discarded. |
162 |
# This allows profiles with dropped CTD scans to be processed without |
|
163 |
# manual intervention. For profiles collected in very calm conditions |
|
164 |
# (e.g. near the ice off Antarctica) time lagging is highly uncertain |
|
165 |
# most of the time --- setting $opt_l = 1 disables the lime-lagging |
|
166 |
# filter for those cases. |
|
30 | 167 |
|
42 | 168 |
# $opt_l = 1; |
30 | 169 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
170 |
|
34 | 171 |
# The following sets the default correlation limit; measurements with |
172 |
# 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
|
173 |
|
13 | 174 |
&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
|
175 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
176 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
177 |
# Instrument Tilt |
17
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
178 |
# |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
179 |
# 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
|
180 |
# 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
|
181 |
# 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
|
182 |
# 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
|
183 |
# 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
|
184 |
# 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
|
185 |
# 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
|
186 |
# 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
|
187 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
188 |
&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
|
189 |
|
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 error velocity limit; measurements |
192 |
# 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
|
193 |
|
13 | 194 |
&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
|
195 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
196 |
|
34 | 197 |
# The following variable allows editing the velocities farthest |
198 |
# from the transducer. It defines how many velocities are to be |
|
199 |
# 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
|
200 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
201 |
$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
|
202 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
203 |
|
34 | 204 |
# The following variable allows editing all data from a given |
205 |
# 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
|
206 |
|
13 | 207 |
$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
|
208 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
209 |
|
34 | 210 |
# The following sets the maximum gap length in the w time series that |
211 |
# 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
|
212 |
|
13 | 213 |
&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
|
214 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
215 |
|
34 | 216 |
# The following variable sets the max allowed vertical ocean velocity |
217 |
# 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
|
218 |
|
13 | 219 |
$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
|
220 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
221 |
|
34 | 222 |
# In each ensemble, vertical velocities differing more than this |
223 |
# parameter times the mean absolute deviation from the median, are |
|
224 |
# 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
|
225 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
226 |
$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
|
227 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
228 |
|
34 | 229 |
# Data from bins with less valid velocities than the following parameter |
230 |
# are considered bad and removed. It is not clear whether this really |
|
231 |
# makes sense, but this editing is likely safe because it only affects |
|
232 |
# 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
|
233 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
234 |
$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
|
235 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
236 |
|
34 | 237 |
# All ensembles recorded when the CTD is shallower than |
238 |
# the following parameter (depth in meters) are discarded. |
|
239 |
# Possible contamination: ship's hull, thrusters, bubble clouds |
|
240 |
# 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
|
241 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
242 |
$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
|
243 |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
244 |
|
34 | 245 |
# Previous Ping Interference editing as described in [edit_data.pl] |
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
246 |
# - enabled by default for WH150 data |
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
247 |
# - PPI_seabed_editing_required defines a string with a perl expression |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
248 |
# that is evaluated once the data are loaded; if true, seabed PPI |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
249 |
# editing is enabled |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
250 |
# - to enable PPI editing without condition, set $PPI_editing = 1; |
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
251 |
# - 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
|
252 |
# 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
|
253 |
# 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
|
254 |
# 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
|
255 |
# 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
|
256 |
# 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
|
257 |
# 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
|
258 |
# 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
|
259 |
# seabed. |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
260 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
261 |
$PPI_seabed_editing_required = '($LADCP{BEAM_FREQUENCY} < 300)'; |
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
262 |
|
47
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
263 |
#$PPI_editing = 1; # uncomment to enable PPI always |
34 | 264 |
#$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
|
265 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
266 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
267 |
# 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
|
268 |
# 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
|
269 |
# 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
|
270 |
# it is not necessary to edit DL data for surface contamination. However, |
34 | 271 |
# 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
|
272 |
# seabed should clearly be removed. |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
273 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
274 |
$sidelobe_editing_DL_surface = 0; |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
275 |
$sidelobe_editing_UL_seabed = 1; |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
276 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
277 |
# The following variable sets the depth for sidelobe contamination |
34 | 278 |
# from the surface. |
32
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
279 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
280 |
$vessel_draft = 6; # in meters |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
281 |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
282 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
283 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
284 |
# Time Lagging |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
285 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
286 |
|
34 | 287 |
# The -i option allows defining an initial guess for the time lag between |
288 |
# the LADCP and the CTD data. |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
289 |
|
13 | 290 |
# $opt_i = 567; |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
291 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
292 |
|
34 | 293 |
# The following variables define the bins used to calculate the reference- |
294 |
# 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
|
295 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
296 |
($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
|
297 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
298 |
|
34 | 299 |
# The -n option defines the number of windows used to calculate |
300 |
# 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
|
301 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
302 |
$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
|
303 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
304 |
|
34 | 305 |
# The -w option defines the width of the window (in seconds) used |
306 |
# to calculate the optimal time lag. There's one value for each |
|
307 |
# 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
|
308 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
309 |
$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
|
310 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
311 |
|
34 | 312 |
# The following parameters control whether the top three time lags |
313 |
# are accepted or not. |
|
314 |
# 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
|
315 |
# (in CTD scans) they must account for at least $TL_required_timelag_top_three_fraction |
34 | 316 |
# or an error is generated. |
317 |
# Notes: |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
318 |
# - $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
|
319 |
# 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
|
320 |
# - large spread may indicate dropped CTD scans |
34 | 321 |
# - the optimum value of $TL_max_allowed_three_lag_spread may be |
322 |
# 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
|
323 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
324 |
$TL_max_allowed_three_lag_spread = 3; |
13 | 325 |
&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
|
326 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
327 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
328 |
#====================================================================== |
29 | 329 |
# 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
|
330 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
331 |
|
29 | 332 |
# After applying the method of Deines (1999), an empirical correction |
333 |
# for Sv is applied to the data. The following variable determines which |
|
334 |
# bin is chosen to construct a reference profile for Sv. The bin number |
|
335 |
# is automatically increased if the selected bin does not contain valid |
|
336 |
# 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
|
337 |
# 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
|
338 |
# causes artifacts every 100m. To disable the empirical |
2ccb81b7cea5
version found on whoosher after repair
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
339 |
# correction, undefine the following variable. |
29 | 340 |
|
341 |
$Sv_ref_bin = 1; |
|
342 |
||
343 |
||
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
344 |
# 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
|
345 |
# 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
|
346 |
# which code is better. |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
347 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
348 |
$SS_use_BT = 0; |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
349 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
350 |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
351 |
# The following variable defines the minimum Sv signal in a bin (max - min) |
34 | 352 |
# 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
|
353 |
# 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
|
354 |
# 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
|
355 |
# 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
|
356 |
# 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
|
357 |
# 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
|
358 |
# (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
|
359 |
# 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
|
360 |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
361 |
$SS_min_signal = 30; |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
362 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
363 |
|
34 | 364 |
# Require at minimum nubmer of valid samples for seabed detection FROM ECHO |
365 |
# AMPLITUDES. Each sample is a bin with a clear seabed maximum. With a proper |
|
366 |
# 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
|
367 |
# 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
|
368 |
# 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
|
369 |
# 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
|
370 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
371 |
$SS_min_samp = 3; |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
372 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
373 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
374 |
# The following numbers define the valid range of height-above bottom |
34 | 375 |
# for seabed detection FROM ECHO AMPLITUDE. For data collected with WH300 |
376 |
# instruments and 8m bins, the maximum range needs to be greater than 250m |
|
377 |
# (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
|
378 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
379 |
$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
|
380 |
$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
|
381 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
382 |
|
34 | 383 |
# Number of ensembles around bottom to search sabed IN BT DATA. |
384 |
# 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
|
385 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
386 |
$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
|
387 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
388 |
|
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
389 |
# 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
|
390 |
# 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
|
391 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
392 |
$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
|
393 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
394 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
395 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
396 |
# Bottom Tracking |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
397 |
# - 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
|
398 |
# 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
|
399 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
400 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
401 |
# 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
|
402 |
# 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
|
403 |
# 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
|
404 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
405 |
$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
|
406 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
407 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
408 |
# 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
|
409 |
# 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
|
410 |
# 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
|
411 |
# 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
|
412 |
# 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
|
413 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
414 |
$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
|
415 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
416 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
417 |
# 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
|
418 |
# 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
|
419 |
# 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
|
420 |
# 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
|
421 |
# 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
|
422 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
423 |
$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
|
424 |
|
42 | 425 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
426 |
#====================================================================== |
42 | 427 |
# Gridded Velocity Profile Output |
18 | 428 |
#====================================================================== |
429 |
||
430 |
||
42 | 431 |
# The -k option defines the minimum number of w samples required in each |
432 |
# vertical-velocity bin. The following sets the default value. |
|
18 | 433 |
|
42 | 434 |
&antsCardOpt(\$opt_k,20); |
18 | 435 |
|
436 |
||
42 | 437 |
# The -o option sets the output grid resolution in meters. The following |
438 |
# sets the default value. |
|
18 | 439 |
|
42 | 440 |
&antsFloatOpt(\$opt_o,20); |
18 | 441 |
|
442 |
||
42 | 443 |
# The following variables limit the bins used to grid w_oean |
444 |
# - in contrast to -b, the other bins are still used e.g. for BT |
|
445 |
# - values recorded in %outgrid_firstbin, %outgrid_lastbin |
|
446 |
# - values beyond range are: |
|
447 |
# - greyed out in *_mean_residuals.ps |
|
448 |
# - not used in *_w.ps, *_residuals.ps |
|
18 | 449 |
|
42 | 450 |
$outGrid_firstBin = '*'; # use $LADCP_firstBin (-b) |
451 |
$outGrid_lastBin = '*'; # use $LADCP_lastBin (-b) |
|
18 | 452 |
|
453 |
||
454 |
||
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
455 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
456 |
1; # return true |
42 | 457 |