author | A.M. Thurnherr <athurnherr@yahoo.com> |
Sun, 13 Mar 2016 12:33:20 +0000 | |
changeset 39 | 91458506d56f |
parent 34 | e550db661c17 |
child 42 | f7690c7b92e0 |
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 |
39 | 4 |
# dlm: Sun Mar 13 10:50:09 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 |
39 | 6 |
# uE-Info: 130 13 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 |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
71 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
72 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
73 |
# Data Input |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
74 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
75 |
|
34 | 76 |
# Set $opt_4 to 1 (or use the -4 option) to suppress 3-beam LADCP |
77 |
# solutions |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
78 |
|
13 | 79 |
#$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
|
80 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
81 |
|
34 | 82 |
# The following variables allow bias-correcting the attiude |
83 |
# sensors. |
|
24 | 84 |
# NB: heading is not used for vertical-velocity processing! |
6 | 85 |
|
34 | 86 |
$pitch_bias = 0; |
87 |
$roll_bias = 0; |
|
88 |
$heading_bias = 0; |
|
6 | 89 |
|
90 |
||
34 | 91 |
# The following variable defines the minimum valid velocities |
92 |
# required in a LADCP file. If there are fewer data, an |
|
93 |
# error is produced |
|
29 | 94 |
|
95 |
$min_valid_vels = 50; |
|
96 |
||
97 |
||
34 | 98 |
# The -b option defines the range of bins to use in w calculations. |
99 |
# The '*' indicates the last bin in the ADCP file. For data |
|
100 |
# collected with non-zero blanking distance, -b '1,*' should |
|
101 |
# 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
|
102 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
103 |
$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
|
104 |
|
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 |
# Logging and Output |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
107 |
#====================================================================== |
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 |
# - there are 4 verbosity levels, selected by -v |
29 | 110 |
# 0 : errors |
111 |
# 1*: UNIX-like (warnings and info messages that are not produced for every cast; *DEFAULT) |
|
112 |
# 2 : progress messages and useful information |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
113 |
# >2: debug messges |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
114 |
# - the most useful ones of these are 1 & 2 |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
115 |
|
33 | 116 |
&antsCardOpt(\$opt_v,$ENV{VERB}); |
117 |
$opt_v = 1 unless numberp($opt_v); |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
118 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
119 |
|
34 | 120 |
# The "base name" of all output files (usually 0-padded 3-digits) |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
121 |
|
29 | 122 |
$out_basename = sprintf('%03d',$PROF); |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
123 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
124 |
|
34 | 125 |
# Output subdirectories |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
126 |
|
39 | 127 |
unless (-d $RUN) { |
128 |
warning(0,"Creating output directory ./$RUN\n"); |
|
129 |
mkdir($RUN); |
|
130 |
error("./$RUN: no such directory\n") unless (-d $RUN); |
|
131 |
} |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
132 |
$data_subdir = $plot_subdir = $log_subdir = $RUN; |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
133 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
134 |
|
34 | 135 |
# The -k option defines the minimum number of w samples required in each |
136 |
# vertical-velocity bin. The following sets the default value. |
|
30 | 137 |
|
138 |
&antsCardOpt(\$opt_k,20); |
|
139 |
||
140 |
||
34 | 141 |
# The -o option sets the output grid resolution in meters. The following |
142 |
# sets the default value. |
|
30 | 143 |
|
33 | 144 |
&antsFloatOpt(\$opt_o,20); |
30 | 145 |
|
146 |
||
34 | 147 |
# The following variables limit the bins used to grid w_oean |
30 | 148 |
# - in contrast to -b, the other bins are still used e.g. for BT |
149 |
# - values recorded in %outgrid_firstbin, %outgrid_lastbin |
|
31 | 150 |
# - values beyond range are: |
151 |
# - greyed out in *_mean_residuals.ps |
|
152 |
# - not used in *_w.ps, *_residuals.ps |
|
30 | 153 |
|
34 | 154 |
$outGrid_firstBin = '*'; # use $LADCP_firstBin (-b) |
155 |
$outGrid_lastBin = '*'; # use $LADCP_lastBin (-b) |
|
30 | 156 |
|
157 |
||
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
158 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
159 |
# Data Editing |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
160 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
161 |
|
34 | 162 |
# The following sets the default correlation limit; measurements with |
163 |
# 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
|
164 |
|
13 | 165 |
&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
|
166 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
167 |
|
34 | 168 |
# The following sets the default limit for instrument attitude |
169 |
# (pitch/roll). |
|
17
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
170 |
# |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
171 |
# The default value was established with IWISE profiles 004, 005, 045 |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
172 |
# and 049, which all show considerabe tilt-related discrepancies between |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
173 |
# the corresponding 2-beam solutions. The first and second pair of |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
174 |
# profiles were collected with 8 and 6m bins, respectively, without |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
175 |
# bin re-mapping. The original default of 15 degrees led to large |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
176 |
# beam-pair differences. Based on diagnostic plots it appears that |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
177 |
# only tilt angles smaller than 9 degrees or so are satisfactory. |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
178 |
# In case of the IWISE data set, such a tight constraint causes |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
179 |
# too many data gaps. The compromise of 12 degrees seems to work |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
180 |
# quite well, based on the p0 vs epsilon correlation across 5 |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
181 |
# data sets. |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
182 |
# |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
183 |
# NB: if this default is changed, the usage message in [LADCP_w] |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
184 |
# needs to be updated as well. |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
185 |
|
17
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
186 |
&antsFloatOpt(\$opt_t,12); |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
187 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
188 |
|
34 | 189 |
# The following sets the default error velocity limit; measurements |
190 |
# 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
|
191 |
|
13 | 192 |
&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
|
193 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
194 |
|
34 | 195 |
# The following variable allows editing the velocities farthest |
196 |
# from the transducer. It defines how many velocities are to be |
|
197 |
# 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
|
198 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
199 |
$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
|
200 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
201 |
|
34 | 202 |
# The following variable allows editing all data from a given |
203 |
# 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
|
204 |
|
13 | 205 |
$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
|
206 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
207 |
|
34 | 208 |
# The following sets the maximum gap length in the w time series that |
209 |
# 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
|
210 |
|
13 | 211 |
&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
|
212 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
213 |
|
34 | 214 |
# The following variable sets the max allowed vertical ocean velocity |
215 |
# 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
|
216 |
|
13 | 217 |
$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
|
218 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
219 |
|
34 | 220 |
# In each ensemble, vertical velocities differing more than this |
221 |
# parameter times the mean absolute deviation from the median, are |
|
222 |
# 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
|
223 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
224 |
$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
|
225 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
226 |
|
34 | 227 |
# Data from bins with less valid velocities than the following parameter |
228 |
# are considered bad and removed. It is not clear whether this really |
|
229 |
# makes sense, but this editing is likely safe because it only affects |
|
230 |
# 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
|
231 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
232 |
$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
|
233 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
234 |
|
34 | 235 |
# All ensembles recorded when the CTD is shallower than |
236 |
# the following parameter (depth in meters) are discarded. |
|
237 |
# Possible contamination: ship's hull, thrusters, bubble clouds |
|
238 |
# 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
|
239 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
240 |
$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
|
241 |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
242 |
|
34 | 243 |
# 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
|
244 |
# - enabled by default for WH150 data |
34 | 245 |
# - the variable defines a string with a perl expression, which is |
246 |
# evaluated once the data are loaded |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
247 |
# - 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
|
248 |
# 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
|
249 |
# 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
|
250 |
# 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
|
251 |
# 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
|
252 |
# 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
|
253 |
# 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
|
254 |
# 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
|
255 |
# seabed. |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
256 |
|
18 | 257 |
$PPI_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
|
258 |
|
34 | 259 |
#$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
|
260 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
261 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
262 |
# 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
|
263 |
# 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
|
264 |
# 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
|
265 |
# it is not necessary to edit DL data for surface contamination. However, |
34 | 266 |
# 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
|
267 |
# seabed should clearly be removed. |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
268 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
269 |
$sidelobe_editing_DL_surface = 0; |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
270 |
$sidelobe_editing_UL_seabed = 1; |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
271 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
272 |
# The following variable sets the depth for sidelobe contamination |
34 | 273 |
# from the surface. |
32
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
274 |
|
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
275 |
$vessel_draft = 6; # in meters |
6041a20feb39
V1.0 release candidate
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
31
diff
changeset
|
276 |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
277 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
278 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
279 |
# Time Lagging |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
280 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
281 |
|
34 | 282 |
# The -i option allows defining an initial guess for the time lag between |
283 |
# the LADCP and the CTD data. |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
284 |
|
13 | 285 |
# $opt_i = 567; |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
286 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
287 |
|
34 | 288 |
# The following variables define the bins used to calculate the reference- |
289 |
# 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
|
290 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
291 |
($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
|
292 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
293 |
|
34 | 294 |
# The -n option defines the number of windows used to calculate |
295 |
# 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
|
296 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
297 |
$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
|
298 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
299 |
|
34 | 300 |
# The -w option defines the width of the window (in seconds) used |
301 |
# to calculate the optimal time lag. There's one value for each |
|
302 |
# 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
|
303 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
304 |
$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
|
305 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
306 |
|
34 | 307 |
# The following parameters control whether the top three time lags |
308 |
# are accepted or not. |
|
309 |
# 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
|
310 |
# (in CTD scans) they must account for at least $TL_required_timelag_top_three_fraction |
34 | 311 |
# or an error is generated. |
312 |
# Notes: |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
313 |
# - $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
|
314 |
# 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
|
315 |
# - large spread may indicate dropped CTD scans |
34 | 316 |
# - the optimum value of $TL_max_allowed_three_lag_spread may be |
317 |
# 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
|
318 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
319 |
$TL_max_allowed_three_lag_spread = 3; |
13 | 320 |
&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
|
321 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
322 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
323 |
#====================================================================== |
29 | 324 |
# 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
|
325 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
326 |
|
29 | 327 |
# After applying the method of Deines (1999), an empirical correction |
328 |
# for Sv is applied to the data. The following variable determines which |
|
329 |
# bin is chosen to construct a reference profile for Sv. The bin number |
|
330 |
# is automatically increased if the selected bin does not contain valid |
|
331 |
# data, i.e. the default value of 1 ensures that the closest valid bin |
|
332 |
# is used to construct the reference profile. |
|
333 |
||
334 |
$Sv_ref_bin = 1; |
|
335 |
||
336 |
||
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
337 |
# 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
|
338 |
# 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
|
339 |
# which code is better. |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
340 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
341 |
$SS_use_BT = 0; |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
342 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
343 |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
344 |
# The following variable defines the minimum Sv signal in a bin (max - min) |
34 | 345 |
# 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
|
346 |
# 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
|
347 |
# 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
|
348 |
# 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
|
349 |
# 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
|
350 |
# 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
|
351 |
# (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
|
352 |
# 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
|
353 |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
354 |
$SS_min_signal = 30; |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
355 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
356 |
|
34 | 357 |
# Require at minimum nubmer of valid samples for seabed detection FROM ECHO |
358 |
# AMPLITUDES. Each sample is a bin with a clear seabed maximum. With a proper |
|
359 |
# 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
|
360 |
# 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
|
361 |
# 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
|
362 |
# 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
|
363 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
364 |
$SS_min_samp = 3; |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
365 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
366 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
367 |
# The following numbers define the valid range of height-above bottom |
34 | 368 |
# for seabed detection FROM ECHO AMPLITUDE. For data collected with WH300 |
369 |
# instruments and 8m bins, the maximum range needs to be greater than 250m |
|
370 |
# (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
|
371 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
372 |
$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
|
373 |
$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
|
374 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
375 |
|
34 | 376 |
# Number of ensembles around bottom to search sabed IN BT DATA. |
377 |
# 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
|
378 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
379 |
$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
|
380 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
381 |
|
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
382 |
# 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
|
383 |
# 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
|
384 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
385 |
$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
|
386 |
|
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 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
389 |
# Bottom Tracking |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
390 |
# - 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
|
391 |
# 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
|
392 |
#====================================================================== |
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 |
# 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
|
395 |
# 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
|
396 |
# 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
|
397 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
398 |
$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
|
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 |
# 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
|
402 |
# 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
|
403 |
# 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
|
404 |
# 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
|
405 |
# 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
|
406 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
407 |
$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
|
408 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
409 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
410 |
# 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
|
411 |
# 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
|
412 |
# 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
|
413 |
# 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
|
414 |
# 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
|
415 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
416 |
$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
|
417 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
418 |
#====================================================================== |
18 | 419 |
# ProcessingParams file selection |
420 |
#====================================================================== |
|
421 |
||
422 |
if (-r "ProcessingParams.$RUN") { |
|
423 |
$processing_param_file = "ProcessingParams.$RUN"; |
|
424 |
} elsif (-r "ProcessingParams.default") { |
|
425 |
$processing_param_file = "ProcessingParams.default"; |
|
426 |
} elsif (-r "ProcessingParams") { |
|
427 |
$processing_param_file = "ProcessingParams"; |
|
428 |
} else { |
|
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
429 |
error("$0: cannot find either <ProcessingParams.$RUN> or <ProcessingParams[.default]>\n"); |
18 | 430 |
} |
431 |
||
432 |
#---------------------------------------------------------------------- |
|
433 |
# Processing log (diagnostic messages) output |
|
434 |
#---------------------------------------------------------------------- |
|
435 |
||
436 |
$out_log = "$log_subdir/$out_basename.log"; |
|
437 |
||
438 |
||
439 |
#---------------------------------------------------------------------- |
|
440 |
# Vertical-velocity profile output and plots: |
|
441 |
# Data: |
|
29 | 442 |
# *.wprof vertical velocity profiles |
18 | 443 |
# Standard Plots: |
29 | 444 |
# *_wprof.ps vertical velocity profiles (main output plot) |
18 | 445 |
#---------------------------------------------------------------------- |
446 |
||
29 | 447 |
@out_profile = ("plot_wprof($plot_subdir/${out_basename}_wprof.ps)", |
448 |
"$data_subdir/$out_basename.wprof"); |
|
18 | 449 |
|
450 |
||
451 |
#---------------------------------------------------------------------- |
|
452 |
# Vertical-velocity sample data output and plots: |
|
453 |
# Data: |
|
31 | 454 |
# *.wsamp w sample data |
18 | 455 |
# Standard Plots: |
31 | 456 |
# *_wsamp.ps vertical velocity time-depth plot |
30 | 457 |
# *_residuals.ps residual vertical velocity time-depth plot |
458 |
# *_backscatter.ps volume scattering coefficient time-depth plot |
|
18 | 459 |
# Optional Plots: |
30 | 460 |
# *_correlation.ps correlation time-depth plot |
18 | 461 |
#---------------------------------------------------------------------- |
462 |
||
34 | 463 |
push(@out_wsamp,"$data_subdir/$out_basename.wsamp"); |
464 |
||
465 |
push(@out_wsamp,"plot_residuals($plot_subdir/${out_basename}_residuals.ps)"); |
|
466 |
push(@out_wsamp,"plot_backscatter($plot_subdir/${out_basename}_backscatter.ps)"); |
|
467 |
push(@out_wsamp,"plot_wsamp($plot_subdir/${out_basename}_wsamp.ps)"); |
|
468 |
||
469 |
#push(@out_wsamp,"plot_correlation($plot_subdir/${out_basename}_correlation.ps)"); |
|
18 | 470 |
|
471 |
||
472 |
#---------------------------------------------------------------------- |
|
473 |
# Time-series output |
|
474 |
# Data: |
|
475 |
# *.tis combined CTD/LADCP time-series data, including |
|
476 |
# package- and LADCP reference layer w |
|
477 |
#---------------------------------------------------------------------- |
|
478 |
||
30 | 479 |
@out_timeseries = ("$data_subdir/$out_basename.tis"); |
18 | 480 |
|
481 |
#---------------------------------------------------------------------- |
|
482 |
# Per-bin vertical-velocity residuals (plot only) |
|
483 |
#---------------------------------------------------------------------- |
|
484 |
||
30 | 485 |
@out_BR = ("plot_mean_residuals($plot_subdir/${out_basename}_mean_residuals.ps)"); |
18 | 486 |
|
487 |
||
488 |
#---------------------------------------------------------------------- |
|
489 |
# Time-lagging correlation statistics (plot only) |
|
490 |
#---------------------------------------------------------------------- |
|
491 |
||
30 | 492 |
@out_TL = ("plot_time_lags($plot_subdir/${out_basename}_time_lags.ps)"); |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
493 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
494 |
1; # return true |