author | A.M. Thurnherr <athurnherr@yahoo.com> |
Fri, 04 Sep 2015 06:38:44 +0000 | |
changeset 31 | d0ae3cb99021 |
parent 30 | 7fb67e771d85 |
child 32 | 6041a20feb39 |
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 |
31 | 4 |
# dlm: Thu Sep 3 11:59:09 2015 |
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 |
31 | 6 |
# uE-Info: 407 7 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 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
63 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
64 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
65 |
# Data Input |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
66 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
67 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
68 |
# CTD depth adjustment |
13 | 69 |
# - set with -d (-a up to 2013/05/16) |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
70 |
# - value is added to CTD pressure |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
71 |
# - use when CTD has -ve pressures |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
72 |
|
13 | 73 |
&antsFloatOpt(\$opt_d,0); |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
74 |
|
23 | 75 |
|
76 |
# set $opt_4 to 1 to suppress 3-beam LADCP solutions |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
77 |
|
13 | 78 |
#$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
|
79 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
80 |
|
6 | 81 |
# correct attiude sensors |
24 | 82 |
# NB: heading is not used for vertical-velocity processing! |
6 | 83 |
|
84 |
$pitch_bias = $roll_bias = $heading_bias = 0; |
|
85 |
||
86 |
||
29 | 87 |
# minimum valid velocities to require in a LADCP file |
88 |
# throws an error if not enough valid vels in a file |
|
89 |
||
90 |
$min_valid_vels = 50; |
|
91 |
||
92 |
||
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
93 |
# bins to use in w calculations |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
94 |
# - set with -b |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
95 |
# - defaults to 2-last |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
96 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
97 |
$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
|
98 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
99 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
100 |
# 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
|
101 |
#====================================================================== |
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 |
# - there are 4 verbosity levels, selected by -v |
29 | 104 |
# 0 : errors |
105 |
# 1*: UNIX-like (warnings and info messages that are not produced for every cast; *DEFAULT) |
|
106 |
# 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
|
107 |
# >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
|
108 |
# - 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
|
109 |
|
29 | 110 |
&antsCardOpt(\$opt_v,1); |
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
111 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
112 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
113 |
# output base name |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
114 |
|
29 | 115 |
$out_basename = sprintf('%03d',$PROF); |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
116 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
117 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
118 |
# output subdirectories |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
119 |
|
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
120 |
error("$RUN: no such directory\n") unless (-d $RUN); |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
121 |
$data_subdir = $plot_subdir = $log_subdir = $RUN; |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
122 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
123 |
|
30 | 124 |
# min w samples required for each vertical-velocity bin |
125 |
# - value recorded in %outgrid_minsamp and used by [LADCP_w_regrid] |
|
126 |
||
127 |
&antsCardOpt(\$opt_k,20); |
|
128 |
||
129 |
||
130 |
# output grid resolution in meters |
|
131 |
# - value recorded in %outgrid_dz and used by [LADCP_w_regrid] |
|
132 |
||
133 |
&antsFloatOpt(\$opt_o,10); |
|
134 |
||
135 |
||
136 |
# the following variables limit the bins used to grid w_oean |
|
137 |
# - in contrast to -b, the other bins are still used e.g. for BT |
|
138 |
# - values recorded in %outgrid_firstbin, %outgrid_lastbin |
|
31 | 139 |
# - values beyond range are: |
140 |
# - greyed out in *_mean_residuals.ps |
|
141 |
# - not used in *_w.ps, *_residuals.ps |
|
30 | 142 |
|
143 |
$outGrid_firstBin = 0; |
|
144 |
$outGrid_lastBin = 999; |
|
145 |
||
146 |
||
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
147 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
148 |
# Data Editing |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
149 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
150 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
151 |
# min correlation |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
152 |
|
13 | 153 |
&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
|
154 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
155 |
|
17
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
156 |
# max tilt (pitch/roll) |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
157 |
# |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
158 |
# 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
|
159 |
# 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
|
160 |
# 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
|
161 |
# 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
|
162 |
# 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
|
163 |
# 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
|
164 |
# 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
|
165 |
# 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
|
166 |
# 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
|
167 |
# 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
|
168 |
# data sets. |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
169 |
# |
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
170 |
# 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
|
171 |
# 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
|
172 |
|
17
fc83e436a800
just before removing code to calculated dc/uc average w
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
16
diff
changeset
|
173 |
&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
|
174 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
175 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
176 |
# max err vel |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
177 |
|
13 | 178 |
&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
|
179 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
180 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
181 |
# truncate farthest valid velocities |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
182 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
183 |
$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
|
184 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
185 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
186 |
# discard velocities from chosen beam (1-4) |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
187 |
|
13 | 188 |
$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
|
189 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
190 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
191 |
# max LADCP gap length in seconds |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
192 |
|
13 | 193 |
&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
|
194 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
195 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
196 |
# max allowed vertical velocity in m/s |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
197 |
|
13 | 198 |
$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
|
199 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
200 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
201 |
# in each ensemble, vertical velocities differing more than this |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
202 |
# parameter times mean absolute deviation from median, are considered |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
203 |
# 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
|
204 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
205 |
$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
|
206 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
207 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
208 |
# data from bins with less valid velocities than the following parameter |
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
209 |
# are considered bad 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
|
210 |
|
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
211 |
$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
|
212 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
213 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
214 |
# ensembles when instrument is shallower than |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
215 |
# $surface_layer_depth in meters are removed. |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
216 |
# possible contamination: ship's hull, thrusters, bubble clouds |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
217 |
# Inspired by 2011_IWISE station 8 |
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 |
$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
|
220 |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
221 |
|
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
222 |
# PPI editing as described in [edit_data.pl] |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
223 |
# - enabled by default for WH150 data |
18 | 224 |
# - variable sets an expression, to be evaluated once the data |
225 |
# have been loaded |
|
16
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
226 |
# - 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
|
227 |
# 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
|
228 |
# 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
|
229 |
# 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
|
230 |
# 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
|
231 |
# 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
|
232 |
# 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
|
233 |
# 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
|
234 |
# seabed. |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
235 |
|
18 | 236 |
$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
|
237 |
|
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
238 |
#$PPI_extend_upper_limit = 1.03; # arbitrarily increase calculated max dist from seabed by 3% |
29e867b3e070
whoosher version at beginning of FZ1
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
13
diff
changeset
|
239 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
240 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
241 |
# Time Lagging |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
242 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
243 |
|
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
244 |
# externally supplied lag initial guess |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
245 |
|
13 | 246 |
# $opt_i = 567; |
5
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
247 |
|
509cc9966b68
======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
4
diff
changeset
|
248 |
|
4
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
249 |
# reference layer bins for w for time matching |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
250 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
251 |
($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
|
252 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
253 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
254 |
# number of time lags during each of 2 lagging steps |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
255 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
256 |
$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
|
257 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
258 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
259 |
# time lag search window widths for each of 2 lagging steps |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
260 |
# - full width in seconds |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
261 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
262 |
$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
|
263 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
264 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
265 |
# if top 3 lags have spread greater than $TL_max_allowed_three_lag_spread |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
266 |
# (in CTD scans) they must account for at least $TL_required_timelag_top_three_fraction |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
267 |
# or there is an error |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
268 |
# - $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
|
269 |
# 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
|
270 |
# - large spread may indicate dropped CTD scans |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
271 |
# - the optimum value of $TL_max_allowed_three_lag_spread may be cast-duration dependent |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
272 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
273 |
$TL_max_allowed_three_lag_spread = 3; |
13 | 274 |
&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
|
275 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
276 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
277 |
#====================================================================== |
29 | 278 |
# 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
|
279 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
280 |
|
29 | 281 |
# After applying the method of Deines (1999), an empirical correction |
282 |
# for Sv is applied to the data. The following variable determines which |
|
283 |
# bin is chosen to construct a reference profile for Sv. The bin number |
|
284 |
# is automatically increased if the selected bin does not contain valid |
|
285 |
# data, i.e. the default value of 1 ensures that the closest valid bin |
|
286 |
# is used to construct the reference profile. |
|
287 |
||
288 |
$Sv_ref_bin = 1; |
|
289 |
||
290 |
||
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
291 |
# 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
|
292 |
# 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
|
293 |
# which code is better. |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
294 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
295 |
$SS_use_BT = 0; |
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
296 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
297 |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
298 |
# The following variable defines the minimum Sv signal in a bin (max - min) |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
299 |
# required for reliable seabed detection. A limiting value of 40dB is |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
300 |
# 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
|
301 |
# 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
|
302 |
# 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
|
303 |
# 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
|
304 |
# 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
|
305 |
# (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
|
306 |
# 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
|
307 |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
308 |
$SS_min_signal = 30; |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
309 |
|
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
310 |
|
28
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
311 |
# Require at minimum nubmer of valid samples for seabed detection. Each |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
312 |
# sample is a bin with a clear seabed maximum. With a proper setting |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
313 |
# of $SS_min_signal, the algorithm is stable even with only a single |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
314 |
# 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
|
315 |
# 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
|
316 |
# 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
|
317 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
318 |
$SS_min_samp = 3; |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
319 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
320 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
321 |
# The following numbers define the valid range of height-above bottom |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
322 |
# for seabed detection. For data collected with WH300 instruments and |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
323 |
# 8m bins, the maximum range needs to be greater than 250m (based on |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
324 |
# GoM#13). |
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
325 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
326 |
$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
|
327 |
$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
|
328 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
329 |
|
b07b23485336
after folding in EGU flight mods (empirical Sv correction)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
27
diff
changeset
|
330 |
# Number of ensembles around bottom to search. 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
|
331 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
332 |
$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
|
333 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
334 |
|
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
335 |
# 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
|
336 |
# 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
|
337 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
338 |
$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
|
339 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
340 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
341 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
342 |
# Bottom Tracking |
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
343 |
# - 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
|
344 |
# 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
|
345 |
#====================================================================== |
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
346 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
347 |
# 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
|
348 |
# 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
|
349 |
# 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
|
350 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
351 |
$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
|
352 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
353 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
354 |
# 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
|
355 |
# 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
|
356 |
# 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
|
357 |
# 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
|
358 |
# 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
|
359 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
360 |
$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
|
361 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
362 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
363 |
# 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
|
364 |
# 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
|
365 |
# 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
|
366 |
# 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
|
367 |
# 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
|
368 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
369 |
$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
|
370 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
371 |
#====================================================================== |
18 | 372 |
# ProcessingParams file selection |
373 |
#====================================================================== |
|
374 |
||
375 |
if (-r "ProcessingParams.$RUN") { |
|
376 |
$processing_param_file = "ProcessingParams.$RUN"; |
|
377 |
} elsif (-r "ProcessingParams.default") { |
|
378 |
$processing_param_file = "ProcessingParams.default"; |
|
379 |
} elsif (-r "ProcessingParams") { |
|
380 |
$processing_param_file = "ProcessingParams"; |
|
381 |
} else { |
|
27
2053d8de8d6b
updated plots & seabed code; GoM data process correctly
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
25
diff
changeset
|
382 |
error("$0: cannot find either <ProcessingParams.$RUN> or <ProcessingParams[.default]>\n"); |
18 | 383 |
} |
384 |
||
385 |
#---------------------------------------------------------------------- |
|
386 |
# Processing log (diagnostic messages) output |
|
387 |
#---------------------------------------------------------------------- |
|
388 |
||
389 |
$out_log = "$log_subdir/$out_basename.log"; |
|
390 |
||
391 |
||
392 |
#---------------------------------------------------------------------- |
|
393 |
# Vertical-velocity profile output and plots: |
|
394 |
# Data: |
|
29 | 395 |
# *.wprof vertical velocity profiles |
18 | 396 |
# Standard Plots: |
29 | 397 |
# *_wprof.ps vertical velocity profiles (main output plot) |
18 | 398 |
#---------------------------------------------------------------------- |
399 |
||
29 | 400 |
@out_profile = ("plot_wprof($plot_subdir/${out_basename}_wprof.ps)", |
401 |
"$data_subdir/$out_basename.wprof"); |
|
18 | 402 |
|
403 |
||
404 |
#---------------------------------------------------------------------- |
|
405 |
# Vertical-velocity sample data output and plots: |
|
406 |
# Data: |
|
31 | 407 |
# *.wsamp w sample data |
18 | 408 |
# Standard Plots: |
31 | 409 |
# *_wsamp.ps vertical velocity time-depth plot |
30 | 410 |
# *_residuals.ps residual vertical velocity time-depth plot |
411 |
# *_backscatter.ps volume scattering coefficient time-depth plot |
|
18 | 412 |
# Optional Plots: |
30 | 413 |
# *_correlation.ps correlation time-depth plot |
18 | 414 |
#---------------------------------------------------------------------- |
415 |
||
31 | 416 |
@out_wsamp = ("plot_residuals($plot_subdir/${out_basename}_residuals.ps)", |
417 |
"plot_backscatter($plot_subdir/${out_basename}_backscatter.ps)", |
|
418 |
# "plot_correlation($plot_subdir/${out_basename}_correlation.ps)", |
|
419 |
"plot_wsamp($plot_subdir/${out_basename}_wsamp.ps)", |
|
420 |
"$data_subdir/$out_basename.wsamp"); |
|
18 | 421 |
|
422 |
||
423 |
#---------------------------------------------------------------------- |
|
424 |
# Time-series output |
|
425 |
# Data: |
|
426 |
# *.tis combined CTD/LADCP time-series data, including |
|
427 |
# package- and LADCP reference layer w |
|
428 |
#---------------------------------------------------------------------- |
|
429 |
||
30 | 430 |
@out_timeseries = ("$data_subdir/$out_basename.tis"); |
18 | 431 |
|
432 |
#---------------------------------------------------------------------- |
|
433 |
# Per-bin vertical-velocity residuals (plot only) |
|
434 |
#---------------------------------------------------------------------- |
|
435 |
||
30 | 436 |
@out_BR = ("plot_mean_residuals($plot_subdir/${out_basename}_mean_residuals.ps)"); |
18 | 437 |
|
438 |
||
439 |
#---------------------------------------------------------------------- |
|
440 |
# Time-lagging correlation statistics (plot only) |
|
441 |
#---------------------------------------------------------------------- |
|
442 |
||
30 | 443 |
@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
|
444 |
|
e681262d0fd7
yay! first version I am really happy with. With plots 'n all.
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff
changeset
|
445 |
1; # return true |