author | A.M. Thurnherr <athurnherr@yahoo.com> |
Thu, 29 May 2014 08:41:19 +0000 | |
changeset 18 | bb7bb9f83db9 |
parent 13 | b176da8559b3 |
child 21 | 0b5bbe60131c |
permissions | -rwxr-xr-x |
0 | 1 |
#!/usr/bin/perl |
2 |
#====================================================================== |
|
3 |
# L I S T E N S |
|
4 |
# doc: Sat Jan 18 18:41:49 2003 |
|
18 | 5 |
# dlm: Tue Mar 4 12:56:36 2014 |
0 | 6 |
# (c) 2003 A.M. Thurnherr |
18 | 7 |
# uE-Info: 46 64 NIL 0 0 72 2 2 4 NIL ofnI |
0 | 8 |
#====================================================================== |
9 |
||
10 |
# Print useful info from the ensemble list or dump ensembles to |
|
11 |
# separate files. |
|
12 |
||
13 |
# HISTORY: |
|
14 |
# Jan 18, 2003: - created |
|
15 |
# Mar 18, 2004: - updated |
|
16 |
# Sep 15, 2005: - made ESW optional (BB150) |
|
17 |
# - change RDI binread library name |
|
18 |
# Aug 25, 2006: - added -r)ange |
|
19 |
# - added write -E)nsembles |
|
20 |
# Aug 26, 2006: - added -M)agdecl |
|
21 |
# Sep 19, 2007: - adapted to new [RDI_BB_Read.pl] (not tested) |
|
22 |
# Jan 26, 2008: - BUG: diagnostic output had been written to STDOUT |
|
23 |
# Feb 1, 2008: - BUG: still more diagnostic output written to STDOUT |
|
24 |
# - BUG: -E/-A combo had ignored -E |
|
25 |
# - changed %-good fieldnames for earth coordinates |
|
26 |
# - allowed for 3-beam solutions |
|
27 |
# Feb 7, 2008: - added -f)ields |
|
28 |
# Apr 4, 2008: - made -f output nan on undefined values |
|
29 |
# - BUG: -f fields did not allow array indices |
|
30 |
# - added in-w)ater data only |
|
31 |
# - restructured for simplicity |
|
32 |
# Mar 2, 2009: - added # of valid bin-1 vels to non-ANTS output |
|
33 |
# Jul 30, 2009: - NaN => nan |
|
1 | 34 |
# Aug 15, 2010: - BUG: usage typo |
5 | 35 |
# Dec 10, 2010: - changed ANTS output to time/date instead of UNIX time |
36 |
# Jan 5, 2011: - added -b |
|
6 | 37 |
# May 12, 2011: - disabled error exit on built-in-test errors when ensembles are dumped to stdout |
10 | 38 |
# Mar 20, 2013: - removed DATA_FORMAT stuff |
11 | 39 |
# Apr 12, 2013: - added bin number to -E output |
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
40 |
# Aug 7, 2013: - added -d, -4 |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
41 |
# - added w12, w34 to -E output |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
42 |
# - added sounspeed correction warning |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
43 |
# - changed -E from prefix to suffix |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
44 |
# - added active header line to -E output |
18 | 45 |
# Mar 4, 2014: - added partial support for DATA_SOURCE_ID |
46 |
# - added support for missing PITCH/ROLL/HEADING |
|
0 | 47 |
|
48 |
# Notes: |
|
5 | 49 |
# - -E outputs data in earth coordinates, unless -b is set also |
0 | 50 |
# - -E output is always in ANTS format, ignoring -A |
51 |
# - no soundspeed correction |
|
52 |
||
53 |
require "getopts.pl"; |
|
54 |
$0 =~ m{(.*/)[^/]+}; |
|
55 |
require "$1RDI_BB_Read.pl"; |
|
56 |
require "$1RDI_Coords.pl"; |
|
57 |
||
58 |
die("Usage: $0 [-A)nts] [-Q)uiet (errcheck only)] " . |
|
59 |
"[-f)ields <[name=]FIELD[,...]>] " . |
|
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
60 |
"[require -4)-beam solutions] [-d)iscard <beam#>] " . |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
61 |
"[write -E)nsemples <.suff> [-M)agnetic <declination>] [min -p)ercent-good <#>] [keep -b)eam coords]] " . |
1 | 62 |
"[-r)ange <first_ens,last_ens>] [in-w)ater ensembles only] " . |
0 | 63 |
"<RDI file...>\n") |
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
64 |
unless (&Getopts("4Abd:E:f:M:p:Qr:w") && $#ARGV >= 0); |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
65 |
|
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
66 |
print(STDERR "WARNING: no soundspeed correction applied!\n"); |
0 | 67 |
|
68 |
print(STDERR "WARNING: magnetic declination not set!\n") |
|
69 |
if defined($opt_E) && !defined($opt_M); |
|
70 |
||
71 |
die("$0: illegal option combination\n") |
|
5 | 72 |
if ($opt_Q && $opt_A) || ((defined($opt_M) || defined($opt_p) || defined($opt_b)) && !defined($opt_E)); |
0 | 73 |
|
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
74 |
die("$0: -4 and -d are mutually exclusive\n") |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
75 |
if ($opt_4 && defined($opt_d)); |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
76 |
|
0 | 77 |
($first_ens,$last_ens) = split(',',$opt_r) |
78 |
if defined($opt_r); |
|
79 |
||
80 |
undef($opt_A) if defined($opt_E); |
|
81 |
||
82 |
$opt_p = 0 unless defined($opt_p); |
|
83 |
||
84 |
if ($opt_f) { # additional fields |
|
85 |
@addFields = split(',',$opt_f); |
|
86 |
foreach my $f (@addFields) { |
|
87 |
$f =~ s/\s//g; # remove spaces |
|
88 |
@def = split('=',$f); |
|
89 |
if (@def == 2) { # name=field |
|
90 |
$addLayout .= $opt_A ? " {$def[0]}" : " $def[0]"; |
|
91 |
$f = $def[1]; |
|
92 |
} else { # field |
|
93 |
$addLayout .= " {$f}"; |
|
94 |
} |
|
95 |
} |
|
96 |
# print(STDERR "addLayout = $addLayout\n"); |
|
97 |
# print(STDERR "\@addFields = @addFields\n"); |
|
98 |
} |
|
99 |
||
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
100 |
$RDI_Coords::minValidVels = 4 if ($opt_4); # no 3-beam solutions |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
101 |
|
0 | 102 |
#---------------------------------------------------------------------- |
103 |
# MAIN |
|
104 |
#---------------------------------------------------------------------- |
|
105 |
||
106 |
while (-f $ARGV[0]) { |
|
107 |
if ($opt_A && !$opt_E) { |
|
108 |
print("#ANTS#PARAMS# RDI_file{$ARGV[0]}\n"); |
|
109 |
} elsif (!$opt_Q) { |
|
110 |
print(STDERR "$ARGV[0]: "); |
|
111 |
} |
|
112 |
readData(@ARGV,\%dta); |
|
113 |
printf(STDERR "%d complete ensembles...\n",scalar(@{$dta{ENSEMBLE}})) |
|
114 |
unless ($opt_Q); |
|
115 |
$dta{HEADING_BIAS} = -$opt_M; # magnetic declination |
|
116 |
shift; |
|
117 |
||
118 |
if ($dta{BEAM_COORDINATES}) { # coords used |
|
119 |
$beamCoords = 1; |
|
120 |
} elsif (!$dta{EARTH_COORDINATES}) { |
|
121 |
die("$ARGV[0]: only beam and earth coordinates implemented so far\n"); |
|
122 |
} |
|
123 |
||
5 | 124 |
die("$ARGV[0]: -b only makes sense for beam-coordinate data\n") |
125 |
if ($opt_b && !$beamCoords); |
|
126 |
||
0 | 127 |
if ($opt_A) { # select output fmt: ANTS |
128 |
unless ($opt_Q) { |
|
129 |
printf("#ANTS#PARAMS# N_ensembles{%d}\n",scalar(@{$dta{ENSEMBLE}})); |
|
130 |
print('#ANTS#FIELDS# {ens} {time} {xducer_up} {temp} {hdg} {pitch} {roll} {XMIT_VOLTAGE} {XMIT_CURRENT}'); |
|
10 | 131 |
print(' {ESW}') if ($dta{FIXED_LEADER_BYTES} >= 53); |
0 | 132 |
print("$addLayout\n"); |
133 |
} |
|
134 |
||
135 |
$dumpEns = sub ($) |
|
136 |
{ |
|
137 |
my($e) = @_; |
|
138 |
||
18 | 139 |
printf('%d %lf %d %g', |
0 | 140 |
$dta{ENSEMBLE}[$e]->{NUMBER}, |
141 |
$dta{ENSEMBLE}[$e]->{UNIX_TIME}, |
|
142 |
$dta{ENSEMBLE}[$e]->{XDUCER_FACING_UP} ? 1 : 0, |
|
143 |
$dta{ENSEMBLE}[$e]->{TEMPERATURE}, |
|
18 | 144 |
); |
145 |
if (defined($dta{ENSEMBLE}[$e]->{HEADING})) { printf(' %g',$dta{ENSEMBLE}[$e]->{HEADING}); } |
|
146 |
else { printf(' nan'); } |
|
147 |
if (defined($dta{ENSEMBLE}[$e]->{PITCH})) { printf(' %g',$dta{ENSEMBLE}[$e]->{PITCH}); } |
|
148 |
else { printf(' nan'); } |
|
149 |
if (defined($dta{ENSEMBLE}[$e]->{ROLL})) { printf(' %g',$dta{ENSEMBLE}[$e]->{ROLL}); } |
|
150 |
else { printf(' nan'); } |
|
151 |
printf(' %g %g', |
|
0 | 152 |
$dta{ENSEMBLE}[$e]->{ADC_XMIT_VOLTAGE}, |
153 |
$dta{ENSEMBLE}[$e]->{ADC_XMIT_CURRENT}, |
|
154 |
); |
|
155 |
printf(' %08X',$dta{ENSEMBLE}[$e]->{ERROR_STATUS_WORD}) |
|
10 | 156 |
if ($dta{FIXED_LEADER_BYTES} >= 53); |
0 | 157 |
foreach my $f (@addFields) { |
158 |
my($fn,$fi) = ($f =~ m{([^[]*)(\[.*)}); |
|
159 |
$fn = $f unless defined($fn); |
|
160 |
my($v) = eval("\$dta{ENSEMBLE}[$e]->{$fn}$fi"); |
|
161 |
print(defined($v) ? " $v" : " nan"); |
|
162 |
} |
|
163 |
print("\n"); |
|
164 |
} |
|
165 |
||
166 |
} elsif ($opt_E) { # one file per ens |
|
167 |
||
168 |
$dumpEns = sub ($) |
|
169 |
{ |
|
170 |
my($e) = @_; |
|
171 |
my($b,$i); |
|
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
172 |
my($file) = "$dta{ENSEMBLE}[$e]->{NUMBER}$opt_E"; |
0 | 173 |
|
174 |
open(P,">$file") || die("$file: $!\n"); |
|
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
175 |
print(P "#!/usr/bin/perl -S list\n"); |
0 | 176 |
print(P "#ANTS#PARAMS# " . |
177 |
"BT_u{$dta{ENSEMBLE}[$e]->{BT_VELOCITY}[0]} " . |
|
178 |
"BT_v{$dta{ENSEMBLE}[$e]->{BT_VELOCITY}[1]} " . |
|
179 |
"BT_w{$dta{ENSEMBLE}[$e]->{BT_VELOCITY}[2]} " . |
|
180 |
"BT_e{$dta{ENSEMBLE}[$e]->{BT_VELOCITY}[3]} " . |
|
181 |
"BT_cor1{$dta{ENSEMBLE}[$e]->{BT_CORRELATION}[0]} " . |
|
182 |
"BT_cor2{$dta{ENSEMBLE}[$e]->{BT_CORRELATION}[1]} " . |
|
183 |
"BT_cor3{$dta{ENSEMBLE}[$e]->{BT_CORRELATION}[2]} " . |
|
184 |
"BT_cor4{$dta{ENSEMBLE}[$e]->{BT_CORRELATION}[3]} " . |
|
185 |
"BT_amp1{$dta{ENSEMBLE}[$e]->{BT_AMPLITUDE}[0]} " . |
|
186 |
"BT_amp2{$dta{ENSEMBLE}[$e]->{BT_AMPLITUDE}[1]} " . |
|
187 |
"BT_amp3{$dta{ENSEMBLE}[$e]->{BT_AMPLITUDE}[2]} " . |
|
188 |
"BT_amp4{$dta{ENSEMBLE}[$e]->{BT_AMPLITUDE}[3]} " . |
|
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
189 |
"soundspeed_correction{NONE!} " . |
0 | 190 |
"\n" |
191 |
); |
|
192 |
print(P "#ANTS#FIELDS# " . |
|
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
193 |
"{bin} {dz} {u} {v} {w} {e} {w12} {w34} {cor1} {cor2} {cor3} {cor4} " . |
0 | 194 |
"{amp1} {amp2} {amp3} {amp4} " |
195 |
); |
|
196 |
if ($beamCoords) { |
|
197 |
print(P "{pcg1} {pcg2} {pcg3} {pcg4}"); |
|
198 |
} else { |
|
199 |
print(P "{pc3beam} {pcBadErrVel} {pc1or2beam} {pc4beam}"); |
|
200 |
} |
|
201 |
print(P "$addLayout\n"); |
|
202 |
||
203 |
for (my($b)=0; $b<$dta{N_BINS}; $b++) { |
|
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
204 |
my(@v,$w12,$w34); |
0 | 205 |
my($dz) = $dta{DISTANCE_TO_BIN1_CENTER} + $b*$dta{BIN_LENGTH}; |
206 |
||
207 |
if ($beamCoords) { |
|
208 |
undef($dta{ENSEMBLE}[$e]->{VELOCITY}[$b][0]) |
|
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
209 |
if (($dta{ENSEMBLE}[$e]->{PERCENT_GOOD}[$b][0] < $opt_p) || ($opt_d == 1)); |
0 | 210 |
undef($dta{ENSEMBLE}[$e]->{VELOCITY}[$b][1]) |
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
211 |
if (($dta{ENSEMBLE}[$e]->{PERCENT_GOOD}[$b][1] < $opt_p) || ($opt_d == 2)); |
0 | 212 |
undef($dta{ENSEMBLE}[$e]->{VELOCITY}[$b][2]) |
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
213 |
if (($dta{ENSEMBLE}[$e]->{PERCENT_GOOD}[$b][2] < $opt_p) || ($opt_d == 3)); |
0 | 214 |
undef($dta{ENSEMBLE}[$e]->{VELOCITY}[$b][3]) |
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
215 |
if (($dta{ENSEMBLE}[$e]->{PERCENT_GOOD}[$b][3] < $opt_p) || ($opt_d == 4)); |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
216 |
($dummy,$w12,$dummy,$w34) = |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
217 |
velBeamToBPEarth(\%dta,$e,@{$dta{ENSEMBLE}[$e]->{VELOCITY}[$b]}); |
5 | 218 |
@v = $opt_b ? @{$dta{ENSEMBLE}[$e]->{VELOCITY}[$b]} : |
219 |
velInstrumentToEarth(\%dta,$e, |
|
220 |
velBeamToInstrument(\%dta, |
|
221 |
@{$dta{ENSEMBLE}[$e]->{VELOCITY}[$b]})); |
|
0 | 222 |
} else { |
223 |
@v = velApplyHdgBias(\%dta,$e,@{$dta{ENSEMBLE}[$e]->{VELOCITY}[$b]}); |
|
224 |
} |
|
225 |
||
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
226 |
$v[0] = nan unless defined($v[0]); # u |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
227 |
$v[1] = nan unless defined($v[1]); # v |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
228 |
$v[2] = nan unless defined($v[2]); # w |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
229 |
$v[3] = nan unless defined($v[3]); # err_vel |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
230 |
$w12 = nan unless defined($w12); # w from beams 1&2 |
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
231 |
$w34 = nan unless defined($w34); # w from beams 3&4 |
0 | 232 |
my(@out) = ( |
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
233 |
$b,$dz,$v[0],$v[1],$v[2],$v[3],$w12,$w34, |
0 | 234 |
@{$dta{ENSEMBLE}[$e]->{CORRELATION}[$b]}, |
235 |
@{$dta{ENSEMBLE}[$e]->{ECHO_AMPLITUDE}[$b]}, |
|
236 |
@{$dta{ENSEMBLE}[$e]->{PERCENT_GOOD}[$b]} |
|
237 |
); |
|
238 |
foreach my $f (@addFields) { |
|
239 |
my($fn,$fi) = ($f =~ m{([^[]*)(\[.*)}); |
|
240 |
$fn = $f unless defined($fn); |
|
241 |
push(@out,eval("\$dta{ENSEMBLE}[$e]->{$fn}$fi")); |
|
242 |
} |
|
13
b176da8559b3
before implementing WBWens (PD0 writing)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
11
diff
changeset
|
243 |
for ($i=0; $i<19+@addFields; $i++) { |
0 | 244 |
$out[$i] = nan unless defined($out[$i]); |
245 |
} |
|
246 |
print(P "@out\n"); |
|
247 |
} |
|
248 |
close(P); |
|
249 |
} |
|
250 |
||
251 |
} else { # neither ANTS nor ens files |
|
252 |
unless ($opt_Q) { |
|
10 | 253 |
if ($dta{FIXED_LEADER_BYTES} >= 53) { |
18 | 254 |
printf(" # Date Time XD Temp Headng Pitch Roll #vv DSID ESW$addLayout\n"); |
0 | 255 |
printf("-----------------------------------------------------------------------\n"); |
256 |
} else { |
|
18 | 257 |
printf(" # Date Time XD Temp Headng Pitch Roll #vv DSID$addLayout\n"); |
0 | 258 |
printf("-------------------------------------------------------------------\n"); |
259 |
} |
|
260 |
} |
|
261 |
||
262 |
$dumpEns = sub ($) |
|
263 |
{ |
|
264 |
my($e) = @_; |
|
265 |
||
18 | 266 |
printf('%5d %s %s %s %5.1f', |
0 | 267 |
$dta{ENSEMBLE}[$e]->{NUMBER}, |
268 |
$dta{ENSEMBLE}[$e]->{DATE}, |
|
269 |
$dta{ENSEMBLE}[$e]->{TIME}, |
|
270 |
$dta{ENSEMBLE}[$e]->{XDUCER_FACING_UP} ? "UP" : "DN", |
|
271 |
$dta{ENSEMBLE}[$e]->{TEMPERATURE}, |
|
18 | 272 |
); |
273 |
if (defined($dta{ENSEMBLE}[$e]->{HEADING})) { printf(' %6.1f',$dta{ENSEMBLE}[$e]->{HEADING}); } |
|
274 |
else { printf(' nan'); } |
|
275 |
if (defined($dta{ENSEMBLE}[$e]->{PITCH})) { printf(' %5.1f',$dta{ENSEMBLE}[$e]->{PITCH}); } |
|
276 |
else { printf(' nan'); } |
|
277 |
if (defined($dta{ENSEMBLE}[$e]->{ROLL})) { printf(' %5.1f',$dta{ENSEMBLE}[$e]->{ROLL}); } |
|
278 |
else { printf(' nan'); } |
|
279 |
printf(' %3d 0x%02X', |
|
0 | 280 |
$dta{ENSEMBLE}[$e]->{BIN1VELS}, |
18 | 281 |
$dta{ENSEMBLE}[$e]->{DATA_SOURCE_ID}, |
0 | 282 |
); |
283 |
printf(' 0x%08X',$dta{ENSEMBLE}[$e]->{ERROR_STATUS_WORD}) |
|
10 | 284 |
if ($dta{FIXED_LEADER_BYTES} >= 53); |
0 | 285 |
foreach my $f (@addFields) { |
286 |
my($fn,$fi) = ($f =~ m{([^[]*)(\[.*)}); |
|
287 |
$fn = $f unless defined($fn); |
|
288 |
my($v) = eval("\$dta{ENSEMBLE}[$e]->{$fn}$fi"); |
|
289 |
print(defined($v) ? " $v" : " nan"); |
|
290 |
} |
|
6 | 291 |
print(" BUILT-IN-TEST ERROR") |
292 |
if defined($dta{ENSEMBLE}[$e]->{BUILT_IN_TEST_ERROR}); |
|
0 | 293 |
print("\n"); |
294 |
} |
|
295 |
||
296 |
} |
|
297 |
||
298 |
for ($e=0; $e<=$#{$dta{ENSEMBLE}}; $e++) { |
|
299 |
next if (defined($first_ens) && |
|
300 |
$dta{ENSEMBLE}[$e]->{NUMBER} < $first_ens); |
|
301 |
last if (defined($last_ens) && |
|
302 |
$dta{ENSEMBLE}[$e]->{NUMBER} > $last_ens); |
|
303 |
$dta{ENSEMBLE}[$e]->{BIN1VELS} = |
|
304 |
defined($dta{ENSEMBLE}[$e]->{VELOCITY}[1][0]) + |
|
305 |
defined($dta{ENSEMBLE}[$e]->{VELOCITY}[1][1]) + |
|
306 |
defined($dta{ENSEMBLE}[$e]->{VELOCITY}[1][2]) + |
|
307 |
defined($dta{ENSEMBLE}[$e]->{VELOCITY}[1][3]); |
|
308 |
next if ($opt_w && $dta{ENSEMBLE}[$e]->{BIN1VELS}<3); |
|
309 |
||
310 |
die("3-beams used in ensemble #$dta{ENSEMBLE}[$e]->{NUMBER}\n") |
|
311 |
if ($dta{ENSEMBLE}[$e]->{N_BEAMS_USED} < 4); |
|
312 |
die("BIT error in ensemble $dta{ENSEMBLE}[$e]->{NUMBER}\n") |
|
6 | 313 |
if ($opt_Q || $opt_A || $opt_E) && defined($dta{ENSEMBLE}[$e]->{BUILT_IN_TEST_ERROR}); |
0 | 314 |
die("Low gain in ensemble #$dta{ENSEMBLE}[$e]->{NUMBER}\n") |
315 |
if ($dta{ENSEMBLE}[$e]->{LOW_GAIN}); |
|
316 |
||
317 |
&$dumpEns($e) |
|
318 |
unless ($opt_Q); |
|
319 |
} |
|
320 |
} |
|
321 |
||
322 |
exit(0); |