author | A.M. Thurnherr <athurnherr@yahoo.com> |
Mon, 11 Nov 2013 11:20:45 -0400 | |
changeset 22 | f6635c0384b7 |
parent 21 | ec19ba9622f3 |
child 26 | d778b73f2a43 |
permissions | -rw-r--r-- |
0
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
1 |
#====================================================================== |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
2 |
# L A D C P P R O C . D E F A U L T S |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
3 |
# doc: Fri Sep 17 09:44:21 2010 |
22
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
4 |
# dlm: Thu Sep 19 13:34:34 2013 |
0
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
5 |
# (c) 2010 A.M. Thurnherr |
22
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
6 |
# uE-Info: 142 75 NIL 0 0 72 0 2 4 NIL ofnI |
0
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
7 |
#====================================================================== |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
8 |
|
9
48b2d27aaebf
after fixing bad magdec bug
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
9 |
# default parameters for [LADCPproc] |
0
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
10 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
11 |
# NOTES: |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
12 |
# - defaults are taken: |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
13 |
# 1) from my current UH processing merge control files |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
14 |
# 2) from the defaults set in Eric's [merge.c] |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
15 |
# - the default version in the ANTS bin dir is always loaded |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
16 |
# - if there is a version in the current processing directory it is loaded |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
17 |
# afterwards |
3 | 18 |
# - per-station parameters can be chosed based on $LADCP_file and $CTD_file |
19 |
# - for additional notes, see [LADCPproc] |
|
0
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
20 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
21 |
# HISTORY: |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
22 |
# Sep 17, 2010: - created |
1 | 23 |
# Dec 9, 2010: - added doc for ASCII CTD file support |
2 | 24 |
# Jun 15, 2011: - added Svbin_start, Svbin_end |
3 | 25 |
# Jul 7, 2011: - added BT processing parameters |
4 | 26 |
# Jul 12, 2011: - added $PPI_editing_enabled |
9
48b2d27aaebf
after fixing bad magdec bug
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
27 |
# Apr 11, 2012: - cosmetics |
12 | 28 |
# May 16, 2012: - BUG: comment var name typo |
19 | 29 |
# Jan 8, 2013: - added CTD_ASCII_header_lines |
20 | 30 |
# Apr 22, 2013: - BUG: BT_begin_search_above value of 300m did not make sense |
31 |
# Jun 5, 2013: - added $bad_beam |
|
21 | 32 |
# Sep 6, 2013: - BUG: BT_begin_search_above value of 300m was correct; |
33 |
# the original bug was in the documentation |
|
22
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
34 |
# Sep 19, 2013: - added support for $BT_range_method |
20 | 35 |
|
36 |
#---------------------------------------------------------------------- |
|
37 |
# Data editing |
|
38 |
#---------------------------------------------------------------------- |
|
39 |
||
40 |
# Sometimes, an ADCP beam degrades to the point where its data contaminate |
|
41 |
# the 4-beam solutions. Set the following variable to the number of the |
|
42 |
# bad beam to discard its data. This variable can also be used for crude |
|
43 |
# wake editing in cases where only a single beam is affected, which can |
|
44 |
# happen in particular with CTD platforms that tend to "weather vane". |
|
45 |
||
46 |
# $bad_beam = 1; |
|
1 | 47 |
|
48 |
#---------------------------------------------------------------------- |
|
49 |
# ASCII CTD file support |
|
50 |
#---------------------------------------------------------------------- |
|
51 |
||
52 |
# By default, [LADCPproc] reads a standard seabird .cnv file. Alternatively, |
|
53 |
# the CTD data can be supplied in a plain ASCII file. To do so, the following |
|
54 |
# variables must be defined in the setup file (-s). Each variable defines |
|
55 |
# the field (column) number of the corresponding data field. Fields are |
|
56 |
# numbered beginning with 1. |
|
57 |
||
58 |
# $CTD_ASCII_sampfreq = 1; ## in Hz |
|
59 |
||
60 |
# $CTD_ASCII_press_field = 1; |
|
61 |
# $CTD_ASCII_temp_field = 2; |
|
62 |
# $CTD_ASCII_salin_field = 3; |
|
63 |
# $CTD_ASCII_lat_field = 4; |
|
64 |
# $CTD_ASCII_lon_field = 5; |
|
65 |
||
19 | 66 |
# The following variables are optional: |
1 | 67 |
|
19 | 68 |
# $CTD_ASCII_header_lines = 4; # if this variable is defined, #-comments are disallowed |
1 | 69 |
# $CTD_ASCII_badval = 999; |
70 |
||
71 |
#---------------------------------------------------------------------- |
|
2 | 72 |
# Backscatter Profile Parameters |
73 |
#---------------------------------------------------------------------- |
|
74 |
||
75 |
# The output profile of volume-scattering coefficient is calculated |
|
76 |
# from a subset of the bins only. This is based on the observation, |
|
77 |
# from on a single P403 profile, that the volume-scattering coefficient |
|
78 |
# correction of Deines (IEEE 1999) yield similar results only for |
|
79 |
# bins 3-9. |
|
80 |
||
81 |
# The seabed-finding routine, on the other hand, uses acoustic |
|
82 |
# backscatter from all bins, as it is possible that the seabed is only |
|
83 |
# seen in the farthest bins. |
|
84 |
||
85 |
$Svbin_start = 3; |
|
86 |
$Svbin_end = 9; |
|
87 |
||
88 |
#---------------------------------------------------------------------- |
|
3 | 89 |
# Bottom Track Parameters |
90 |
#---------------------------------------------------------------------- |
|
91 |
||
92 |
# First bin to consider when looking for seabed to calculated post- |
|
93 |
# processed BT. NB: For consistency with EF's code, bin 1 is the |
|
94 |
# first bin. |
|
95 |
||
96 |
$BT_bin_start = 1; |
|
97 |
||
21 | 98 |
# Apparently valid BT data collected at distances greater than 300m |
99 |
# from the seabed are discarded. |
|
3 | 100 |
|
21 | 101 |
$BT_begin_search_above = 300; |
3 | 102 |
|
103 |
# Maximum allowed spread of bin number where max echo is found. Large values |
|
104 |
# imply sloping seabed and/or large instrument tilt. The default value of |
|
105 |
# 3 is probably too tight. The range spread calculation should be modified |
|
106 |
# to take instrument tilt into consideration. |
|
107 |
||
108 |
$BT_max_bin_spread = 3; |
|
109 |
||
110 |
# For tricky BT cases, the acoustic backscatter from the seabed can be plotted |
|
111 |
# and the bottom range set manually. This has been found to be very powerful |
|
112 |
# in case of 2011_IWISE yoyo profile 160, where the bottom was around 150m |
|
113 |
# from the lower turning points and where the method of Deines [1999] clearly |
|
114 |
# does not work. |
|
115 |
||
116 |
#$BT_min_depth = xxx; |
|
12 | 117 |
#$BT_max_depth = xxx; |
3 | 118 |
|
119 |
# Maximum difference between water depth and average distance of echo max., |
|
120 |
# if $BT_min_depth and $BT_max_depth are not set. The stddev of the detected |
|
22
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
121 |
# water depth is added to this number. |
3 | 122 |
|
123 |
$BT_max_depth_error = 20; |
|
124 |
||
125 |
# Maximum allowed difference between reference-layer w and BT w. Note |
|
126 |
# that this must be relaxed in regions of significant near-bottom |
|
127 |
# vertical velocity, e.g. due to cross-slope currents. |
|
128 |
||
129 |
$BT_max_w_difference = 0.03; |
|
130 |
||
22
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
131 |
# When the BT velocity is taken from the wrong bin, the measurement is not |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
132 |
# taken from the main acoustic beam but, rather, from a side lobe. As a |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
133 |
# result, the BT measurement is biased, causing biased near-bottom LADCP |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
134 |
# velocities when the instrument is moving horizontally. |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
135 |
# Methods: |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
136 |
# 0: Take BT velocity from bin with max(Sv). Not recommended. |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
137 |
# 1: Chose either bin at max(Sv) or one of its neighbors, |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
138 |
# depending which shows the smallest discrepancy between w_BT and w_reflr. |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
139 |
# Causes underestimation of instrument speed in P403 tow-yos. |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
140 |
# 2: Visbeck (2002): use median from 3 bins. Requires $BT_range_Visbeck_center |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
141 |
# to define where the 3 bins should be centered wrt. the max(Sv) bin. |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
142 |
# Visbeck suggests +1. Based on P403 tow-yo 030 it looks more like -2. |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
143 |
|
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
144 |
$BT_range_method = 1; |
f6635c0384b7
beginning of DIMES US5 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
21
diff
changeset
|
145 |
|
3 | 146 |
#---------------------------------------------------------------------- |
1 | 147 |
# Shear Processing Parameters |
148 |
#---------------------------------------------------------------------- |
|
0
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
149 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
150 |
## u_bin0, u_bin1, w_bin0, w_bin1: These set the first and |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
151 |
## last bin indices used when integrating horizontal and |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
152 |
## vertical velocity, respectively. The indices start |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
153 |
## from 1 and are inclusive. These parameters do not |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
154 |
## affect the calculation of shear. They should specify |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
155 |
## a good reference layer (for the calculation of the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
156 |
## barotropic component) such as bins 1 to 4. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
157 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
158 |
$wbin_start = 1; ## These parameters start from 1, not zero |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
159 |
$wbin_end = 5; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
160 |
$ubin_start = 1; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
161 |
$ubin_end = 5; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
162 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
163 |
## sh_bin0, sh_bin1: These are like the above except that |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
164 |
## they affect the shear calculation only. They can |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
165 |
## normally be left at the default range of 1 to 128 |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
166 |
## (disabled). |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
167 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
168 |
$shbin_start = 1; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
169 |
$shbin_end = $LADCP{N_BINS}; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
170 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
171 |
## w_ref_bin, w_dif: These control one of the editing |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
172 |
## criteria recommended by Fischer and Visbeck, or rather |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
173 |
## a modification of it. All velocity data are rejected |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
174 |
## below the point at which the vertical velocity |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
175 |
## estimate is w_dif larger or smaller than the estimate |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
176 |
## in w_ref_bin. I have not found this criterion |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
177 |
## helpful. 94/12/22: changed, so that only those points |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
178 |
## where w actually deviates from the mean from 0 to w_ref_bin |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
179 |
## are flagged. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
180 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
181 |
$w_ref_bin = 10; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
182 |
$w_dif = 0.05; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
183 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
184 |
## wake_hd_dif, wake_ang_min, min_wake_w, n_wake_bins: These control |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
185 |
## editing based upon the direction and inclination of the package wake, |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
186 |
## calculated from reference layer U,V, and W. wake_hd_dif |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
187 |
## sets how close to the heading of the wake must be to |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
188 |
## the heading of any beam for interference (in degrees). |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
189 |
## wake_ang_min sets the minimum wake angle from the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
190 |
## vertical for interference (in degrees). min_wake_w sets |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
191 |
## the minimum package speed required for wake interference |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
192 |
## (although upward speed is measured negative, this parameter |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
193 |
## is input as a positive value). All three criteria |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
194 |
## must be satisfied for the wake flag to be set. n_wake_bins |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
195 |
## determines how many bins from the top are removed from a |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
196 |
## flagged profile, the default is 1 (top bin only). |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
197 |
## Additionally, if the previous ensemble met all |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
198 |
## interference criteria, the present ensemble |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
199 |
## will be flagged even if the criteria are not met. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
200 |
## Default values are wake_hd_dif=0.0, wake_ang_min=90.0, |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
201 |
## and min_wake_w=0.1, which results in no wake editing. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
202 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
203 |
$wake_hd_dif = 0.0; ## set wake editing defaults so that no wake editing |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
204 |
$wake_ang_min = 90.0; ## occurs |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
205 |
$min_wake_w = 0.1; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
206 |
$n_wake_bins = 1; ## wake editing default - top bin only |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
207 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
208 |
## e_max is the maximum error velocity. An error velocity |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
209 |
## greater than e_max will flag the other velocity |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
210 |
## components. For the BB with 2-ping ensembles, 0.01 |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
211 |
## m/s looks about right for this parameter. It |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
212 |
## knocks out the incorrect ambiguity glitches as well |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
213 |
## as some smaller but still significant glitches. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
214 |
## Note that for the BB, the PG criteria operate on |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
215 |
## percent 4-beam solutions, so PG combined with the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
216 |
## e_max provides a reliable filter against big |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
217 |
## glitches. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
218 |
## For the NB with 18-ping ensembles, it is not yet |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
219 |
## clear whether there is any benefit in using e_max |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
220 |
## at all. To be effective, the value will have to be |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
221 |
## small, something like 0.015 or 0.02. Note also |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
222 |
## that for the NB, the pg array holds pg counting |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
223 |
## both 4-beam and 3-beam solutions, so 3-beam |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
224 |
## glitches will slip through the e_max net. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
225 |
## Default is e_max= 10.0, which effectively disables |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
226 |
## it. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
227 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
228 |
$e_max = 0.1; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
229 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
230 |
## min_correlation (BB only) is the minimum correlation, in |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
231 |
## counts, for each beam in a given bin. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
232 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
233 |
$min_cor = 70; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
234 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
235 |
## Shear editing: This requires 2 passes through the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
236 |
## database, one using the option |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
237 |
## binned_shear_time_range: |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
238 |
## to set the time range for a pass during which the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
239 |
## shear statistics will be calculated on a relatively |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
240 |
## coarsely grid by binning rather than interpolation, |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
241 |
## then using the usual option |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
242 |
## time_range: |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
243 |
## for the normal pass in which the results of the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
244 |
## previous pass can be used to flag bad velocity values |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
245 |
## based on anomalous shears. The first pass can include |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
246 |
## both up and down casts if desired, in which case it |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
247 |
## would be followed by a second pass for the up and |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
248 |
## downcasts separately. The binned shear statistics are |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
249 |
## saved until they are explicitly recalculated with the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
250 |
## "binned_shear_time_range:" option. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
251 |
## The shear editing is controlled by these parameters: |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
252 |
## shear_dev_max= x.x, where x.x is a floating point |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
253 |
## number giving a threshold in standard |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
254 |
## deviations. Any shear component deviating from |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
255 |
## the binned mean by more than this times the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
256 |
## local standard deviation will raise a flag. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
257 |
## Suitable values for this parameter are probably |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
258 |
## in the range 3-5. Less than 3 is likely to |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
259 |
## start trimming too many valid samples, more than |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
260 |
## 4 or 5 is likely to do nothing at all. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
261 |
## shear_sum_dev_max= x.x gives a threshold used in |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
262 |
## detecting isolated bad velocity points. If the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
263 |
## sum of two successive bad shears (based on |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
264 |
## shear_dev_max), divided by the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
265 |
## standard deviation, is LESS THAN this amount, |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
266 |
## then the common velocity point is considered to |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
267 |
## be an isolated glitch, and it alone is flagged. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
268 |
## Otherwise, both velocity samples contributing to |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
269 |
## a flagged shear will be flagged. A reasonable |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
270 |
## value for this parameter is probably around 1-2, |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
271 |
## but I don't yet have enough experience to be |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
272 |
## sure. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
273 |
## Warning: don't forget to reset the first_x=, |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
274 |
## first_y=, and first_z= parameters after pass 1 |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
275 |
## (with binned_shear_time_range:). If you don't, and |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
276 |
## the depths are not in the database, then the depths |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
277 |
## will be completely wrong in pass 2 (with |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
278 |
## time_range:). You may also want to use a different |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
279 |
## output file name for pass 1, so that the binned |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
280 |
## shear statistics file will not be overwritten |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
281 |
## during pass 2. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
282 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
283 |
$max_shdev = 3.5; ## to disable, set to nan |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
284 |
$max_shdev_sum = 1.5; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
285 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
286 |
## previous ping bottom bounce interference editing is controlled by: |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
287 |
## clip_margin 0.0 turns off this editing; otherwise, it is the margin in |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
288 |
## meters (on each side) by which the calculated range of the |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
289 |
## interference is expanded before clipping. A reasonable value would |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
290 |
## be something like 32 (2 depth bins). |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
291 |
|
4 | 292 |
undef($PPI_editing_enabled); # set to 1 to enable PPI editing |
293 |
||
0
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
294 |
# NOTES: |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
295 |
# - default value (90m) taken from comment in merge control files |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
296 |
# - clipping is disabled if water_depth cannot be determined |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
297 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
298 |
$clip_margin = 90; ## default: no previous ping bottom bounce editing |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
299 |
$first_clip_bin = 1; ## default: apply previous ping clipping to all bins */ |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
300 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
301 |
# tilt editing as in Visbeck's code |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
302 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
303 |
$max_tilt = 22; # max allowed angle from vertical |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
304 |
$max_delta_tilt = 4; # max allowed ping-to-ping tilt difference |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
305 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
306 |
# On DIMES US2 stations in Drake Passage it was found that the ambiguity velocity had been |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
307 |
# set too low. It appears, however, that aliased velocities are easy to detect because aliasing |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
308 |
# causes a large positive velocity to appear as a large negative one, and vice versa. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
309 |
# The following defines the maximum allowed discrepancy between the vertical velocity from CTD |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
310 |
# pressure and the LADCP reference-layer w. This parameter should probably be set to something |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
311 |
# similar as the ambiguity velocity. |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
312 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
313 |
$w_max_err = 2.5; |
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
314 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
315 |
|
de00d0f32431
initial version to process P403 yoyo
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
316 |
1; |