author | Andreas Thurnherr <ant@ldeo.columbia.edu> |
Mon, 08 Mar 2021 07:58:53 -0500 | |
changeset 57 | 5a59411306ba |
parent 56 | 82384f531dba |
child 58 | 78607e2e8add |
permissions | -rw-r--r-- |
14 | 1 |
#====================================================================== |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
2 |
# R D I _ P D 0 _ I O . P L |
14 | 3 |
# doc: Sat Jan 18 14:54:43 2003 |
57
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
4 |
# dlm: Sat Mar 6 14:13:23 2021 |
14 | 5 |
# (c) 2003 A.M. Thurnherr |
57
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
6 |
# uE-Info: 643 43 NIL 0 0 72 2 2 4 NIL ofnI |
14 | 7 |
#====================================================================== |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
8 |
|
37 | 9 |
# Read RDI PD0 binary data files (*.[0-9][0-9][0-9]) |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
10 |
|
14 | 11 |
# HISTORY: |
12 |
# Jan 18, 2003: - incepted aboard the Aurora Australis (KAOS) |
|
13 |
# Jan 19, 2003: - continued |
|
14 |
# Jan 20, 2003: - replaced INTENSITY by AMPLITUDE |
|
15 |
# Jan 21, 2003: - changed heading-correction field names |
|
16 |
# Jan 27, 2003: - cosmetics |
|
17 |
# Feb 14, 2003: - moved BT setup into header |
|
18 |
# Mar 15, 2003: - moved 10th xmit voltage into header as BATTERY field |
|
19 |
# - removed again, because values are not meaningful |
|
20 |
# Feb 24, 2004: - continued aboard Nathaniel B. Palmer (Anslope II) |
|
21 |
# - BUG: ensemble # was wrong on error messages |
|
22 |
# Feb 26, 2004: - removed ESW_ERROR and bitmasking of ERROR_STATUS_WORD |
|
23 |
# Feb 27, 2004: - removed some unused (commented-out) baggage |
|
24 |
# Mar 11, 2004: - BUG: renamed ACD -> ADC |
|
25 |
# Mar 18, 2004: - cosmetics |
|
26 |
# Mar 30, 2004: - rewrote to speed up reading; new version takes |
|
27 |
# ~40% less time |
|
28 |
# Sep 14, 2005: - made BT data optional (NUMBER_OF_DATA_TYPES) |
|
29 |
# - added DATA_FORMAT |
|
30 |
# Sep 15, 2005: - debugged |
|
31 |
# - implement checksum to robustly find EOF |
|
32 |
# - renamed to RDI_BB_Read.pl |
|
33 |
# - BUG: had used POSIX::mktime with wrong month def! |
|
34 |
# Oct 30, 2005: - added WH300 FW16.27 file format |
|
35 |
# - added DATA_FORMAT_VARIANT |
|
36 |
# - changed semantics so that first valid ensemble is |
|
39 | 37 |
# in E[0] (instead of E[$ensNo-1]) |
14 | 38 |
# Nov 8, 2005: - replaced UNIXTIME by UNIX_TIME |
39 |
# - added SECNO |
|
40 |
# Aug 31: 2006: - added DAYNO |
|
41 |
# Aug 1, 2007: - BUG: typo in monthLength() |
|
42 |
# Sep 18, 2007: - modified readHeader() readDta() WBRhdr() WBRens() to |
|
43 |
# conserve memory (no large arrays as retvals) |
|
44 |
# Jun 4, 2008: - BUG: BB150 code was not considered on Sep 18, 2007 |
|
45 |
# Aug 15, 2010: - downgraded "unexpected number of data types" from error to warning |
|
46 |
# - BUG: WBRcfn had not been set correctly |
|
47 |
# - modified to allow processing files without time info |
|
48 |
# May 12, 2011: - added code to report built-in-test errors |
|
49 |
# Mar 19, 2013: - added support for WH600 data file (58 fixed leader bytes) |
|
50 |
# Mar 20, 2013: - removed DATA_FORMAT stuff |
|
51 |
# - added support for BT data in subset of ensembles |
|
52 |
# Apr 29, 2013: - changed semantics to assume EOF when unexpected number of data types |
|
53 |
# are present in an ensemble |
|
54 |
# Nov 25, 2013: - renamed from [RDI_BB_Read.pl] |
|
55 |
# - begin implementing WBWens() |
|
56 |
# - checkEnsemble() expunged |
|
18 | 57 |
# Mar 3, 2014: - BUG: WBPens() did not handle incomple ensembles at EOF correctly |
58 |
# Mar 4, 2014: - added support for DATA_SOURCE_ID |
|
59 |
# Apr 24, 2014: - added debug statements to log %-GOOD values |
|
60 |
# May 6, 2014: - loosened input format checks |
|
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
61 |
# May 7, 2014: - removed BT_PRESENT flag |
21 | 62 |
# Sep 6, 2014: - adapted WBRhdr to >7 data types |
63 |
# Oct 15, 2014: - implemented work-around for readData() not recognizing |
|
64 |
# incomplete ensemble at the end, which seems to imply that there is |
|
39 | 65 |
# a garbage final ensemble that passes the checksum test??? |
26
02776c168135
made RDI_BB_IO work with files containing leading garbage
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
24
diff
changeset
|
66 |
# Oct 2, 2015: - added &skip_initial_trash() |
28 | 67 |
# Dec 18, 2015: - added most data types to WBPofs() |
68 |
# - BUG: WBPens() requires round() for scaled values |
|
31 | 69 |
# Jan 9, 2016: - removed system() from writeData() |
70 |
# - BUG: WBRhdr() did not set DATA_SOURCE_ID |
|
29 | 71 |
# - added PRODUCER |
31 | 72 |
# - BUG: writeData() did not work correctly for ECOGIG OC26 moored data (spaces in filename?) |
73 |
# - added support for patching coordinate system |
|
32 | 74 |
# Feb 16, 2016: - added transducer orientation to WBPens() |
75 |
# - BUG: most WBPens() error messages used wrong file name |
|
76 |
# Feb 23, 2016: - changed WBRhdr() to use 2nd ensemble (with correct data-source id) |
|
77 |
# Feb 26, 2016: - added basic BT data to WBPens(); not BT_RL_* and BT_SIGNAL_STRENGTH |
|
36
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
78 |
# Feb 29, 2016: - LEAP DAY: actually got BT data patching to work |
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
79 |
# Jul 30, 2016: - BUG: incomplete last ensemble with garbage content was returned on reading |
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
80 |
# WH300 data |
37 | 81 |
# Aug 5, 2016: - cosmetics |
82 |
# Aug 23, 2016: - added &clearEns() |
|
83 |
# Nov 9, 2016: - made WBRhdr() return undef on "empty" files |
|
84 |
# Nov 18, 2016: - BUG: ensNo was not reported correctly in format errors |
|
85 |
# Nov 23, 2016: - no longer set pitch/roll/heading to undef in clearEns() |
|
86 |
# Mar 7, 2016: - renamed round() to stop clashing with ANTSLIB |
|
39 | 87 |
# May 18, 2017: - maded readHeader() more permissive (checkFmt flag) |
88 |
# - added partial support for Ocean Surveyor data |
|
89 |
# Aug 1, 2017: BUG: minor typo in ocean surveyor code (err check not done?) |
|
90 |
# Aug 7, 2017: - added LAG_LENGTH |
|
91 |
# - added SPEED_OF_SOUND to header |
|
92 |
# Aug 8, 2017: - replaced croak() by die() |
|
93 |
# - added actual transducer frequencies |
|
40
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
94 |
# Nov 22, 2017: - BUG: dayNo() and monthLength() clashed with [libconv.pl] |
42 | 95 |
# - added support for RDI_PD0_IO::IGNORE_Y2K_CLOCK |
96 |
# Dec 7, 2017: - added suppress_error to readHeader() |
|
97 |
# Dec 23, 2017: - BUG: could no longer read Anslope II raw files |
|
98 |
# - added support for patching ADCP time data |
|
99 |
# - added support for RDI_PD0_IO::OVERRIDE_Y2K_CLOCK |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
100 |
# Feb 6, 2018: - added supprort for first_ens, last_ens in readData() |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
101 |
# Feb 7, 2018: - BUG: new params were not optional |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
102 |
# - made read routines more permissive by allowing |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
103 |
# garbage between ensembles |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
104 |
# Mar 15, 2018: - BUG: WBPens() did not work for files with garbage |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
105 |
# Mar 16, 2018: - added skipped garbage warning |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
106 |
# - BUG: garbage skipping did not work correctly for files w/o BT |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
107 |
# Mar 20, 2018: - BUG: garbage skipping STILL did not work correctly for files w/o BT |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
108 |
# Apr 9, 2018: - added last_bin argument to readData() |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
109 |
# Apr 10, 2018: - slight improvement (parameter range check) |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
110 |
# Apr 23, 2018: - make WBRens() work (again?) with BB150 data froim 1996 |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
111 |
# Apr 24, 2018: - BUG: undefined lat_bin argument to readData() did not work |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
112 |
# Apr 30, 2018: - added support for repeated ensembles |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
113 |
# - added warning on wrong ensemble length |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
114 |
# Jun 9, 2018: - removed double \n from warnings |
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
115 |
# Jun 12, 2018: - BUG: IMPed files did not pass the garbage detection |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
116 |
# Jun 13, 2019: - adapted reading routines to RTI files (free order of data types) |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
117 |
# - removed old BT_PRESENT code |
50 | 118 |
# Jun 28, 2019: - renamed SECONDS to SECOND for consistency |
119 |
# Jun 30, 2019: - added dirty flag to prevent bad PD0 patching |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
120 |
# - added $show_progress |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
121 |
# Feb 13, 2020: - added support for $readDataProgress (to Jun 13 version) |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
122 |
# Apr 12, 2020: - merged laptop and whoosher versions which both implemented |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
123 |
# progreass; not sure whether this merge is successful |
52
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
124 |
# - disabled duplicate ens detection, which turns reading into |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
125 |
# O(N^2) process |
55 | 126 |
# Apr 14, 2020: - BUG: WBPens did not work for ens# > 65535 (high byte) |
127 |
# Mar 3, 2021: - adapted to Nortek PD0 files |
|
128 |
# END OF HISTORY |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
129 |
|
14 | 130 |
# FIRMWARE VERSIONS: |
131 |
# It appears that different firmware versions generate different file |
|
132 |
# structures. Currently (Sep 2005) these routines have been tested |
|
133 |
# with the following firmware versions (as reported by [listHdr]): |
|
134 |
# |
|
39 | 135 |
# Firmw. DATA_FORMAT(_VARIANT) Owner Cruise FIXED_LEADER_LENGTH |
14 | 136 |
#------------------------------------------------------------ |
39 | 137 |
# 05.52 BB150 (1) UH CLIVAR/P16S 42 |
14 | 138 |
# 16.12 WH300 (1) FSU A0304 53 |
39 | 139 |
# 16.21 WH300 (1) LDEO NBP0402 53 |
140 |
# 16.27 WH300 (2) Nash ? 59 |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
141 |
|
29 | 142 |
# PD0 FILE FORMAT EXTENSIONS: |
143 |
# |
|
32 | 144 |
# - file creator encoded in DATA_SOURCE_ID |
145 |
# |
|
146 |
# - first ensemble uses default RDI DATA_SOURCE_ID because the LDEO_IX |
|
147 |
# software assumes this |
|
148 |
# |
|
39 | 149 |
# - DATA_SOURCE_ID = 0x7F original TRDI PD0 file |
29 | 150 |
# |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
151 |
# - DATA_SOURCE_ID = 0xA0 | PATCHED_MASK produced by IMP+LADCP, KVH+LADCP |
29 | 152 |
# PATCHED_MASK & 0x04: pitch value has been patched |
153 |
# PATCHED_MASK & 0x02: roll value has been patched |
|
154 |
# PATCHED_MASK & 0x01: heading value has been patched |
|
155 |
# - PITCH & ROLL can be missing (0x8000 badval as in velocities) |
|
156 |
# - HEADING can be missing (0xF000 badval, as 0x8000 is valid 327.68 heading) |
|
157 |
# |
|
31 | 158 |
# - DATA_SOURCE_ID = 0xE0 produced by editPD0 |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
159 |
|
14 | 160 |
# NOTES: |
161 |
# - RDI stores data in VAX/Intel byte order (little endian) |
|
162 |
# - the output data structure does not exactly mirror the file data |
|
163 |
# structure; the header is not stored at all and the fixed leader |
|
164 |
# data are not duplicated in every ensemble |
|
165 |
# - in the RDI files some fields that logically belong into the header |
|
166 |
# or the fixed leader (e.g. BT_MIN_CORRELATION) appear in the |
|
167 |
# ensemble data --- these are not read on input |
|
168 |
# - the field names are generally unabbreviated except for |
|
169 |
# BT (= Bottom Track), RL (= Reference Layer), MIN and MAX |
|
170 |
# - all arrays are 0-referenced, but the ensemble number is not! |
|
171 |
# - a list of filenames can be be passed to readData() so that |
|
172 |
# files split onto several memory cards (typically .000 .001 &c) |
|
173 |
# can be read --- not sure if this works, actually |
|
174 |
# - the RDI manuals are not entirely clear everywhere; I have made |
|
175 |
# guesses in some cases, but they should not affect the main |
|
176 |
# fields of interest |
|
177 |
# - some fields in the fixed leader are not really fixed in a LADCP |
|
178 |
# setting (e.g. xducer orientation); I'v made an educated guess |
|
179 |
# as to which fields to move to the ENS array |
|
180 |
# - all units except pressure are SI, i.e. in m and m/s |
|
181 |
# - I don't understand the ERROR_STATUS_WORD; here's what 3 different |
|
182 |
# instruments returned: |
|
39 | 183 |
# 0x88000100 FSU instrument during A0304 (Firmware 16.12) |
14 | 184 |
# 0x88008180 LDEO uplooker (slave) during NBP0402 (Firmware 16.21) |
185 |
# 0x00008180 LDEO downlooker (master) during NBP0402 (Firmware 16.21) |
|
186 |
# According to the manual (January 2001 version) this would, for example, |
|
187 |
# indicate power failures on both FSU and LDEO slave instruments... |
|
40
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
188 |
# - defining the variable "$RDI_PD0_IO::IGNORE_Y2K_CLOCK" before calling &readData() |
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
189 |
# makes the code ignore the Y2K clock and use the old clock instead; this |
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
190 |
# is used for Dan Torres' KVH system |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
191 |
|
14 | 192 |
# &readData() returns perl obj (ref to anonymous hash) with the following |
193 |
# structure: |
|
194 |
# |
|
18 | 195 |
# DATA_SOURCE_ID scalar 0x7f (Workhorse, also DVL) |
14 | 196 |
# NUMBER_OF_DATA_TYPES scalar 6 (no BT) or 7 |
197 |
# ENSEMBLE_BYTES scalar ?, number of bytes w/o checksum |
|
198 |
# HEADER_BYTES scalar ? |
|
199 |
# FIXED_LEADER_BYTES scalar 42 for BB150; 53 for WH300, 58 for WH600, 59 for WH300(Nash) |
|
200 |
# VARIABLE_LEADER_BYTES scalar ? |
|
39 | 201 |
# VELOCITY_DATA_BYTES scalar ? |
14 | 202 |
# CORRELATION_DATA_BYTES scalar ? |
203 |
# ECHO_INTENSITY_DATA_BYTES scalar ? |
|
39 | 204 |
# PERCENT_GOOD_DATA_BYTES scalar ? |
14 | 205 |
# CPU_FW_VER scalar 0--255 |
206 |
# CPU_FW_REV scalar 0--255 |
|
207 |
# BEAM_FREQUENCY scalar 75, 150, 300, 600, 1200, 2400 [kHz] |
|
39 | 208 |
# LAG_LENGTH scalar ??? |
209 |
# CONVEX_BEAM_PATTERN bool undefined, 1 |
|
14 | 210 |
# CONCAVE_BEAM_PATTERN bool undefined, 1 |
211 |
# SENSOR_CONFIG scalar 1--3 |
|
212 |
# XDUCER_HEAD_ATTACHED bool undefined, 1 |
|
213 |
# BEAM_ANGLE scalar 15,20,30,undefined=other [deg] |
|
39 | 214 |
# N_BEAMS scalar 4--5 |
14 | 215 |
# N_DEMODS scalar 2--3(???),undefined=n/a |
216 |
# N_BINS scalar 1--128 |
|
217 |
# PINGS_PER_ENSEMBLE scalar 0--16384 |
|
218 |
# BIN_LENGTH scalar 0.01--64 [m] |
|
219 |
# BLANKING_DISTANCE scalar 0-99.99 [m] |
|
39 | 220 |
# MIN_CORRELATION scalar 0--255 |
14 | 221 |
# N_CODE_REPETITIONS scalar 0--255 |
222 |
# MIN_PERCENT_GOOD scalar 1--100 [%] |
|
223 |
# MAX_ERROR_VELOCITY scalar 0--5 [m/s] |
|
224 |
# TIME_BETWEEN_PINGS scalar 0--? [s] |
|
225 |
# BEAM_COORDINATES bool undefined,1 |
|
226 |
# INSTRUMENT_COORDINATES bool undefined,1 |
|
227 |
# SHIP_COORDINATES bool undefined,1 |
|
228 |
# EARTH_COORDINATES bool undefined,1 |
|
39 | 229 |
# PITCH_AND_ROLL_USED bool undefined,1 |
14 | 230 |
# USE_3_BEAM_ON_LOW_CORR bool undefined,1 |
39 | 231 |
# BIN_MAPPING_ALLOWED bool undefined,1 |
232 |
# HEADING_ALIGNMENT scalar -179.99..180 [deg] |
|
233 |
# HEADING_BIAS scalar -179.99..180 [deg] |
|
57
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
234 |
# SSPEED_CALCULATED bool undefined,1 |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
235 |
# SSPEED_USING_PRESS bool undefined,1 |
39 | 236 |
# USE_COMPASS bool undefined,1 |
14 | 237 |
# USE_PITCH_SENSOR bool undefined,1 |
39 | 238 |
# USE_ROLL_SENSOR bool undefined,1 |
239 |
# USE_CONDUCTIVITY_SENSOR bool undefined,1 |
|
14 | 240 |
# USE_TEMPERATURE_SENSOR bool undefined,1 |
241 |
# SPEED_OF_SOUND_CALCULATED bool undefined,1 |
|
242 |
# PRESSURE_SENSOR_AVAILABLE bool undefined,1 |
|
243 |
# COMPASS_AVAILABLE bool undefined,1 |
|
244 |
# PITCH_SENSOR_AVAILABLE bool undefined,1 |
|
245 |
# ROLL_SENSOR_AVAILABLE bool undefined,1 |
|
246 |
# CONDUCTIVITY_SENSOR_AVAILABLE bool undefined,1 |
|
247 |
# TEMPERATURE_SENSOR_AVAILABLE bool undefined,1 |
|
39 | 248 |
# DISTANCE_TO_BIN1_CENTER scalar 0--655.35 [m] |
14 | 249 |
# TRANSMITTED_PULSE_LENGTH scalar 0--655.35 [m] |
250 |
# RL_FIRST_BIN scalar 1--128 |
|
39 | 251 |
# RL_LAST_BIN scalar 1--128 |
14 | 252 |
# FALSE_TARGET_THRESHOLD scalar 0--254, undefined=disabled |
39 | 253 |
# LOW_LATENCY_SETTING scalar 0--5(???) |
14 | 254 |
# TRANSMIT_LAG_DISTANCE scalar 0--655.35 [m] |
255 |
# CPU_SERIAL_NUMBER scalar undefined, 0--65535 if WH300 |
|
256 |
# NARROW_BANDWIDTH bool undefined,1 (only set if WH300) |
|
257 |
# WIDE_BANDWIDTH bool undefined,1 (only set if WH300) |
|
258 |
# TRANSMIT_POWER scalar undefined, 0--255(high) if WH300 |
|
39 | 259 |
# TRANSMIT_POWER_HIGH bool undefined,1 (only set if WH300) |
14 | 260 |
# BT_PINGS_PER_ENSEMBLE scalar 0--999 |
261 |
# BT_DELAY_BEFORE_REACQUIRE scalar 0--999 |
|
262 |
# BT_MIN_CORRELATION scalar 0--255 |
|
263 |
# BT_MIN_EVAL_AMPLITUDE scalar 0--255 |
|
39 | 264 |
# BT_MIN_PERCENT_GOOD scalar 0--100 [%] |
265 |
# BT_MODE scalar 4,5,6(?) |
|
14 | 266 |
# BT_MAX_ERROR_VELOCITY scalar 0--5 [m/s], undef=not screened |
267 |
# BT_RL_MIN_SIZE scalar 0--99.9 [m] |
|
268 |
# BT_RL_NEAR scalar 0--999.9 [m] |
|
269 |
# BT_RL_FAR scalar 0--999.9 [m] |
|
270 |
# BT_MAX_TRACKING_DEPTH scalar 8--999.9 [m] |
|
39 | 271 |
# ENSEMBLE[ensemble_no-1] array ensemble info |
14 | 272 |
# XDUCER_FACING_UP bool undefined, 1 |
273 |
# XDUCER_FACING_DOWN bool undefined, 1 |
|
274 |
# N_BEAMS_USED scalar 3,4,5(?) |
|
275 |
# NUMBER scalar 1--16777215 |
|
39 | 276 |
# BUILT_IN_TEST_ERROR scalar ?,undefined=none |
14 | 277 |
# SPEED_OF_SOUND scalar 1400--1600 [m/s] |
278 |
# XDUCER_DEPTH scalar 0.1--999.9 [m] |
|
39 | 279 |
# HEADING scalar 0--359.99 [deg] --- IMP EXTENSION: undef |
18 | 280 |
# PITCH scalar -20.00-20.00 [deg] --- IMP EXTENSION: undef |
281 |
# ROLL scalar -20.00-20.00 [deg] --- IMP EXTENSION: undef |
|
14 | 282 |
# SALINITY scalar 0-40 [psu] |
39 | 283 |
# TEMPERATURE scalar -5.00--40.00 [deg] |
14 | 284 |
# MIN_PRE_PING_WAIT_TIME scalar ? [s] |
285 |
# HEADING_STDDEV scalar 0-180 [deg] |
|
286 |
# PITCH_STDDEV scalar 0.0-20.0 [deg] |
|
39 | 287 |
# ROLL_STDDEV scalar 0.0-20.0 [deg] |
14 | 288 |
# ADC_XMIT_CURRENT scalar 0--255 |
289 |
# ADC_XMIT_VOLTAGE scalar 0--255 |
|
39 | 290 |
# ADC_AMBIENT_TEMPERATURE scalar 0--255 |
14 | 291 |
# ADC_PRESSURE_PLUS scalar 0--255 |
292 |
# ADC_PRESSURE_MINUS scalar 0--255 |
|
293 |
# ADC_ATTITUDE_TEMPERATURE scalar 0--255 |
|
294 |
# ADC_ATTITUDE scalar 0--255 |
|
295 |
# ADC_CONTAMINATION scalar 0--255 |
|
296 |
# ERROR_STATUS_WORD scalar undefined, ? (only set if WH300) |
|
297 |
# PRESSURE scalar undefined, ?-? [dbar] (only set if WH300) |
|
39 | 298 |
# PRESSURE_STDDEV scalar undefined, ?-? [dbar] (only set if WH300) |
14 | 299 |
# DATE string MM/DD/YYYY |
300 |
# YEAR scalar ? |
|
301 |
# MONTH scalar 1--12 |
|
39 | 302 |
# DAY scalar 1--31 |
14 | 303 |
# TIME string HH:MM:SS.hh |
304 |
# HOUR scalar 0--23 |
|
305 |
# MINUTE scalar 0--59 |
|
50 | 306 |
# SECOND scalar 0--59.99 |
14 | 307 |
# UNIX_TIME scalar 0--? |
308 |
# SECNO scalar 0--? (number of seconds since daystart) |
|
309 |
# DAYNO double fractional day number since start of current year (1.0 is midnight Jan 1st) |
|
39 | 310 |
# VELOCITY[bin][beam] scalars -32.767--32.768 [m/s], undef=bad |
311 |
# CORRELATION[bin][beam] scalars 1--255, undefined=bad |
|
312 |
# ECHO_AMPLITUDE[bin][beam] scalars 0--255 |
|
313 |
# PERCENT_GOOD[bin][beam] scalars 0--255 |
|
314 |
# BT_RANGE[beam] scalars tons [m] |
|
315 |
# BT_VELOCITY[beam] scalars see VELOCITY |
|
316 |
# BT_CORRELATION[beam] scalars see CORRELATION |
|
317 |
# BT_EVAL_AMPLITUDE[beam] scalars 0--255 |
|
318 |
# BT_PERCENT_GOOD[beam] scalars see PERCENT_GOOD |
|
319 |
# BT_RL_VELOCITY[beam] scalars see VELOCITY |
|
320 |
# BT_RL_CORRELATION[beam] scalars see CORRELATION |
|
321 |
# BT_RL_ECHO_AMPLITUDE[beam] scalars see ECHO_AMPLITUDE |
|
322 |
# BT_RL_PERCENT_GOOD[beam] scalars see PERCENT_GOOD |
|
323 |
# BT_SIGNAL_STRENGTH[beam] scalars 0--255 |
|
14 | 324 |
# HIGH_GAIN bool 1, undefined |
325 |
# LOW_GAIN bool 1, undefined |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
326 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
327 |
use strict; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
328 |
use Time::Local; # timegm() |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
329 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
330 |
$RDI_PD0_IO::show_progress = 0; # when set, print . every 1000 ens read |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
331 |
|
14 | 332 |
#---------------------------------------------------------------------- |
333 |
# Time Conversion Subroutines |
|
40
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
334 |
# - prepended with _ to avoid conflicts with [libconv.pl] |
14 | 335 |
#---------------------------------------------------------------------- |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
336 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
337 |
sub _monthLength($$) # of days in month |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
338 |
{ |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
339 |
my($Y,$M) = @_; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
340 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
341 |
return 31 if ($M==1 || $M==3 || $M==5 || $M==7 || |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
342 |
$M==8 || $M==10 || $M==12); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
343 |
return 30 if ($M==4 || $M==6 || $M==9 || $M==11); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
344 |
return 28 if ($Y%4 != 0); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
345 |
return 29 if ($Y%100 != 0); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
346 |
return 28 if ($Y%400 > 0); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
347 |
return 29; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
348 |
} |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
349 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
350 |
{ my($epoch,$lM,$lD,$lY,$ldn); # static scope |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
351 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
352 |
sub _dayNo($$$$$$) |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
353 |
{ |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
354 |
my($Y,$M,$D,$h,$m,$s) = @_; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
355 |
my($dn); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
356 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
357 |
if ($Y==$lY && $M==$lM && $D==$lD) { # same day as last samp |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
358 |
$dn = $ldn; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
359 |
} else { # new day |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
360 |
$epoch = $Y unless defined($epoch); # 1st call |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
361 |
$lY = $Y; $lM = $M; $lD = $D; # store |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
362 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
363 |
for ($dn=0,my($cY)=$epoch; $cY<$Y; $cY++) { # multiple years |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
364 |
$dn += 337 + &_monthLength($Y,$M); |
14 | 365 |
} |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
366 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
367 |
$dn += $D; # day in month |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
368 |
while (--$M > 0) { # preceding months |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
369 |
$dn += &_monthLength($Y,$M); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
370 |
} |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
371 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
372 |
$ldn = $dn; # store |
14 | 373 |
} |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
374 |
return $dn + $h/24 + $m/24/60 + $s/24/3600; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
375 |
} |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
376 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
377 |
} # static scope |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
378 |
|
14 | 379 |
#---------------------------------------------------------------------- |
380 |
# Read Data |
|
381 |
#---------------------------------------------------------------------- |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
382 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
383 |
my($WBRcfn,$WBPcfn); # current file names for reading/patching |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
384 |
my($BIT_errors) = 0; # built-in-test errors |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
385 |
|
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
386 |
my($FmtErr) = "%s: illegal %s Id 0x%04x at ensemble %d"; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
387 |
|
14 | 388 |
#---------------------------------------------------------------------- |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
389 |
# skip to next valid start of ensemble (skip over garbage) |
26
02776c168135
made RDI_BB_IO work with files containing leading garbage
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
24
diff
changeset
|
390 |
#---------------------------------------------------------------------- |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
391 |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
392 |
sub goto_next_ens(@) |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
393 |
{ |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
394 |
my($fh,$return_skipped) = @_; # if return_skipped not set, return file pos |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
395 |
my($buf,$dta); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
396 |
|
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
397 |
my($found) = 0; # zero consecutive 0x7f found |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
398 |
my($skipped) = 0; my($garbage_start); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
399 |
while ($found < 2) { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
400 |
sysread($fh,$buf,1) == 1 || last; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
401 |
($dta) = unpack('C',$buf); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
402 |
if ($dta == 0x7f) { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
403 |
$found++; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
404 |
} elsif ($found==1 && |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
405 |
($dta==0xE0 || # from editPD0 |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
406 |
(($dta&0xF0)==0xA0 && ($dta&0x0F)<8))) { # from IMP+LADCP or KVH+LADCP |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
407 |
$found++; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
408 |
} elsif ($found == 0) { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
409 |
$garbage_start = sysseek($fh,0,1)-1 unless defined($garbage_start); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
410 |
$skipped++; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
411 |
} else { # here, found == 1 but 2nd byte was not found |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
412 |
$garbage_start = sysseek($fh,0,1)-$found unless defined($garbage_start); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
413 |
$skipped += $found; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
414 |
$found = 0; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
415 |
} |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
416 |
} |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
417 |
my($fpos) = ($found < 2) ? undef : sysseek($fh,-2,1); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
418 |
return $skipped if ($return_skipped); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
419 |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
420 |
if ($skipped) { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
421 |
if (eof($fh)) { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
422 |
# 04/18/18: disabled the following line of code because it is very common at |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
423 |
# least with the older RDI instruments I am looking at in the context |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
424 |
# of the SR1b repeat section analysis |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
425 |
# print(STDERR "WARNING (RDI_PD0_IO): PD0 file ends with $skipped garbage bytes\n"); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
426 |
} elsif ($garbage_start == 0) { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
427 |
print(STDERR "WARNING (RDI_PD0_IO): PD0 file starts with $skipped garbage bytes\n"); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
428 |
} else { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
429 |
print(STDERR "WARNING (RDI_PD0_IO): $skipped garbage bytes in PD0 file beginning at byte $garbage_start\n"); |
26
02776c168135
made RDI_BB_IO work with files containing leading garbage
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
24
diff
changeset
|
430 |
} |
02776c168135
made RDI_BB_IO work with files containing leading garbage
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
24
diff
changeset
|
431 |
} |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
432 |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
433 |
return $fpos; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
434 |
} |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
435 |
|
26
02776c168135
made RDI_BB_IO work with files containing leading garbage
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
24
diff
changeset
|
436 |
#---------------------------------------------------------------------- |
14 | 437 |
# readHeader(file_name,^dta) WBRhdr(^data) |
438 |
# - read header data |
|
439 |
# - also includes some data from 1st ens |
|
440 |
#---------------------------------------------------------------------- |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
441 |
|
14 | 442 |
sub readHeader(@) |
443 |
{ |
|
42 | 444 |
my($fn,$dta,$suppress_error) = @_; |
14 | 445 |
$WBRcfn = $fn; |
39 | 446 |
open(WBRF,$WBRcfn) || die("$WBRcfn: $!"); |
42 | 447 |
if (WBRhdr($dta)) { |
448 |
return 1; |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
449 |
} elsif ($suppress_error) { |
42 | 450 |
return undef; |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
451 |
} else { |
42 | 452 |
die("$WBRcfn: Insufficient data\n"); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
453 |
} |
14 | 454 |
} |
455 |
||
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
456 |
sub WBRhdr($) |
14 | 457 |
{ |
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
458 |
my($dta) = @_; |
56
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
459 |
my($hid,$did,$B,$W,$i,$dummy); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
460 |
local our($ndt,$buf,$id,$start_ens,@WBRofs); |
14 | 461 |
my($B1,$B2,$B3,$B4,$B5,$B6,$B7,$W1,$W2,$W3,$W4,$W5); |
39 | 462 |
|
14 | 463 |
#-------------------- |
464 |
# HEADER |
|
465 |
#-------------------- |
|
466 |
||
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
467 |
my($skipped) = goto_next_ens(\*WBRF,1); |
50 | 468 |
if ($skipped > 0) { |
469 |
$RDI_PD0_IO::File_Dirty = 1; |
|
470 |
printf(STDERR "WARNING: %d bytes of initial garbage\n",$skipped); |
|
471 |
} |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
472 |
|
37 | 473 |
sysread(WBRF,$buf,6) == 6 || return undef; |
32 | 474 |
($hid,$did,$dta->{ENSEMBLE_BYTES},$dummy,$dta->{NUMBER_OF_DATA_TYPES}) |
475 |
= unpack('CCvCC',$buf); |
|
36
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
476 |
$hid == 0x7f || die(sprintf($FmtErr,$WBRcfn,"Header (hid)",$hid,0)); |
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
477 |
$did == 0x7f || die(sprintf($FmtErr,$WBRcfn,"Header (did)",$did,0)); |
56
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
478 |
$ndt = $dta->{NUMBER_OF_DATA_TYPES}; |
32 | 479 |
|
37 | 480 |
$start_ens = sysseek(WBRF,$dta->{ENSEMBLE_BYTES}-6+2,1) || return undef; |
481 |
sysread(WBRF,$buf,6) == 6 || return undef; |
|
14 | 482 |
($hid,$did,$dta->{ENSEMBLE_BYTES},$dummy,$dta->{NUMBER_OF_DATA_TYPES}) |
483 |
= unpack('CCvCC',$buf); |
|
36
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
484 |
$hid == 0x7f || die(sprintf($FmtErr,$WBRcfn,"Header (hid2)",$hid,0)); |
29 | 485 |
$dta->{DATA_SOURCE_ID} = $did; |
486 |
if ($did == 0x7f) { |
|
57
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
487 |
$dta->{PRODUCER} = 'TRDI or Nortek ADCP'; |
31 | 488 |
} elsif (($did&0xF0) == 0xA0) { |
489 |
$dta->{PRODUCER} = 'IMP+LADCP (Thurnherr software)'; |
|
490 |
} elsif (($did&0xF0) == 0xE0) { |
|
491 |
$dta->{PRODUCER} = 'editPD0 (Thurnherr software)'; |
|
29 | 492 |
} else { |
31 | 493 |
$dta->{PRODUCER} = sprintf('unknown (0x%02X)'); |
29 | 494 |
} |
495 |
||
14 | 496 |
sysread(WBRF,$buf,2*$dta->{NUMBER_OF_DATA_TYPES}) |
497 |
== 2*$dta->{NUMBER_OF_DATA_TYPES} |
|
498 |
|| die("$WBRcfn: $!"); |
|
499 |
@WBRofs = unpack("v$dta->{NUMBER_OF_DATA_TYPES}",$buf); |
|
56
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
500 |
|
14 | 501 |
# for ($i=0; $i<$dta->{NUMBER_OF_DATA_TYPES}; $i++) { |
56
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
502 |
# printf(STDERR "WBRofs[$i] = %d\n",$WBRofs[$i]); |
14 | 503 |
# } |
504 |
||
39 | 505 |
$dta->{HEADER_BYTES} = $WBRofs[0]; |
506 |
$dta->{FIXED_LEADER_BYTES} = $WBRofs[1] - $WBRofs[0]; |
|
507 |
$dta->{VARIABLE_LEADER_BYTES} = $WBRofs[2] - $WBRofs[1]; |
|
508 |
||
509 |
if ($dta->{FIXED_LEADER_BYTES} == 42) { # Eric Firing's old instrument I used in 2004 |
|
510 |
$dta->{INSTRUMENT_TYPE} = 'BB150'; |
|
511 |
} elsif ($dta->{FIXED_LEADER_BYTES} == 53) { # old firmware: no serial numbers |
|
512 |
$dta->{INSTRUMENT_TYPE} = 'Workhorse'; |
|
513 |
} elsif ($dta->{FIXED_LEADER_BYTES} == 59) { # new firmware: with serial numbers |
|
57
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
514 |
$dta->{INSTRUMENT_TYPE} = 'Workhorse or Signature'; |
39 | 515 |
} elsif ($dta->{FIXED_LEADER_BYTES} == 58) { # DVL |
516 |
$dta->{INSTRUMENT_TYPE} = 'Explorer'; |
|
517 |
} elsif ($dta->{FIXED_LEADER_BYTES} == 60) { # OS75 |
|
518 |
$dta->{INSTRUMENT_TYPE} = 'Ocean Surveyor'; |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
519 |
} else { |
39 | 520 |
$dta->{INSTRUMENT_TYPE} = 'unknown'; |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
521 |
} |
39 | 522 |
|
56
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
523 |
my($BT_dt) = WBRdtaIndex(0x0600); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
524 |
$dta->{BT_PRESENT} = defined($BT_dt); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
525 |
@{$dta->{DATA_TYPES}} = WBRdtaTypes(); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
526 |
|
39 | 527 |
#-------------------------------- |
528 |
# Variable Leader: SPEED_OF_SOUND |
|
529 |
#-------------------------------- |
|
14 | 530 |
|
26
02776c168135
made RDI_BB_IO work with files containing leading garbage
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
24
diff
changeset
|
531 |
sysseek(WBRF,$start_ens+$WBRofs[1],0) || die("$WBRcfn: $!"); |
14 | 532 |
sysread(WBRF,$buf,2) == 2 || die("$WBRcfn: $!"); |
533 |
$id = unpack('v',$buf); |
|
39 | 534 |
if ($dta->{INSTRUMENT_TYPE} eq 'Ocean Surveyor') { |
535 |
$id == 0x0081 || printf(STDERR $FmtErr."\n",$WBRcfn,"Variable Leader",$id,1); |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
536 |
} else { |
39 | 537 |
$id == 0x0080 || printf(STDERR $FmtErr."\n",$WBRcfn,"Variable Leader",$id,1); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
538 |
} |
39 | 539 |
sysseek(WBRF,12,1) || die("$WBRcfn: $!"); # skip up to speed of sound |
14 | 540 |
sysread(WBRF,$buf,2) == 2 || die("$WBRcfn: $!"); |
39 | 541 |
$dta->{SPEED_OF_SOUND} = unpack('v',$buf); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
542 |
|
14 | 543 |
#-------------------- |
544 |
# FIXED LEADER |
|
545 |
#-------------------- |
|
546 |
||
26
02776c168135
made RDI_BB_IO work with files containing leading garbage
A.M. Thurnherr <ant@ldeo.columbia.edu>
parents:
24
diff
changeset
|
547 |
sysseek(WBRF,$start_ens+$WBRofs[0],0) || die("$WBRcfn: $!"); |
14 | 548 |
sysread(WBRF,$buf,42) == 42 || die("$WBRcfn: $!"); |
39 | 549 |
($id,$dta->{CPU_FW_VER},$dta->{CPU_FW_REV},$B1,$B2,$dummy, |
550 |
$dta->{LAG_LENGTH},$dummy, |
|
14 | 551 |
$dta->{N_BINS},$dta->{PINGS_PER_ENSEMBLE},$dta->{BIN_LENGTH}, |
552 |
$dta->{BLANKING_DISTANCE},$dummy,$dta->{MIN_CORRELATION}, |
|
553 |
$dta->{N_CODE_REPETITIONS},$dta->{MIN_PERCENT_GOOD}, |
|
554 |
$dta->{MAX_ERROR_VELOCITY},$dta->{TIME_BETWEEN_PINGS},$B3,$B4,$B5, |
|
555 |
$dta->{HEADING_ALIGNMENT},$dta->{HEADING_BIAS},$B6,$B7, |
|
556 |
$dta->{DISTANCE_TO_BIN1_CENTER},$dta->{TRANSMITTED_PULSE_LENGTH}, |
|
557 |
$dta->{REF_LAYER_FIRST_BIN},$dta->{REF_LAYER_LAST_BIN}, |
|
558 |
$dta->{FALSE_TARGET_THRESHOLD},$dta->{LOW_LATENCY_SETTING}, |
|
559 |
$dta->{TRANSMIT_LAG_DISTANCE}) = |
|
560 |
unpack('vCCCCC3CvvvCCCCvCCCCvvCCvvCCCCv',$buf); |
|
561 |
||
39 | 562 |
if ($dta->{INSTRUMENT_TYPE} eq 'Ocean Surveyor') { |
563 |
$id == 0x0001 || printf(STDERR $FmtErr."\n",$WBRcfn,"Fixed Leader",$id,0); |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
564 |
} else { |
39 | 565 |
$id == 0x0000 || printf(STDERR $FmtErr."\n",$WBRcfn,"Fixed Leader",$id,0); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
566 |
} |
14 | 567 |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
568 |
# $dta->{BEAM_FREQUENCY} = 2**($B1 & 0x07) * 75; # nominal |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
569 |
if (($B1&0x07) == 0b000) { $dta->{BEAM_FREQUENCY} = 76.8; } # actual |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
570 |
elsif (($B1&0x07) == 0b001) { $dta->{BEAM_FREQUENCY} = 153.6; } |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
571 |
elsif (($B1&0x07) == 0b010) { $dta->{BEAM_FREQUENCY} = 307.2; } |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
572 |
elsif (($B1&0x07) == 0b011) { $dta->{BEAM_FREQUENCY} = 614.4; } |
39 | 573 |
elsif (($B1&0x07) == 0b100) { $dta->{BEAM_FREQUENCY} = 1228.8; } |
574 |
elsif (($B1&0x07) == 0b101) { $dta->{BEAM_FREQUENCY} = 2457.6; } |
|
575 |
else { die(sprintf("$WBRcfn: cannot decode BEAM_FREQUENCY (%03b)\n",$B1&0x07)); } |
|
576 |
||
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
577 |
$dta->{CONVEX_BEAM_PATTERN} = 1 if ($B1 & 0x08); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
578 |
$dta->{CONCAVE_BEAM_PATTERN} = 1 if (!($B1 & 0x08)); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
579 |
$dta->{SENSOR_CONFIG} = ($B1 & 0x30) >> 4; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
580 |
$dta->{XDUCER_HEAD_ATTACHED} = 1 if ($B1 & 0x40); |
14 | 581 |
|
582 |
if (($B2 & 0x03) == 0x00) { $dta->{BEAM_ANGLE} = 15; } |
|
583 |
elsif (($B2 & 0x03) == 0x01) { $dta->{BEAM_ANGLE} = 20; } |
|
584 |
elsif (($B2 & 0x03) == 0x02) { $dta->{BEAM_ANGLE} = 30; } |
|
585 |
if (($B2 & 0xF0) == 0x40) { $dta->{N_BEAMS} = 4; } |
|
586 |
elsif (($B2 & 0xF0) == 0x50) { $dta->{N_BEAMS} = 5; $dta->{N_DEMODS} = 3; } |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
587 |
elsif (($B2 & 0xF0) == 0xF0) { $dta->{N_BEAMS} = 5; $dta->{N_DEMODS} = 2; } |
14 | 588 |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
589 |
$dta->{BIN_LENGTH} /= 100; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
590 |
$dta->{BLANKING_DISTANCE} /= 100; |
14 | 591 |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
592 |
$dta->{MAX_ERROR_VELOCITY} /= 1000; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
593 |
$dta->{TIME_BETWEEN_PINGS} *= 60; |
14 | 594 |
$dta->{TIME_BETWEEN_PINGS} += $B3 + $B4/100; |
595 |
||
596 |
$dta->{BEAM_COORDINATES} = 1 if (($B5 & 0x18) == 0x00); |
|
597 |
$dta->{INSTRUMENT_COORDINATES} = 1 if (($B5 & 0x18) == 0x08); |
|
598 |
$dta->{SHIP_COORDINATES} = 1 if (($B5 & 0x18) == 0x10); |
|
599 |
$dta->{EARTH_COORDINATES} = 1 if (($B5 & 0x18) == 0x18); |
|
600 |
$dta->{PITCH_AND_ROLL_USED} = 1 if ($B5 & 0x04); |
|
601 |
$dta->{USE_3_BEAM_ON_LOW_CORR} = 1 if ($B5 & 0x02); |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
602 |
$dta->{BIN_MAPPING_ALLOWED} = 1 if ($B5 & 0x01); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
603 |
|
14 | 604 |
$dta->{HEADING_ALIGNMENT} = |
605 |
($dta->{EARTH_COORDINATES} || $dta->{SHIP_COORDINATES}) ? |
|
606 |
$dta->{HEADING_ALIGNMENT} / 100 : undef; |
|
607 |
$dta->{HEADING_BIAS} = |
|
608 |
($dta->{EARTH_COORDINATES} || $dta->{SHIP_COORDINATES}) ? |
|
609 |
$dta->{HEADING_BIAS} / 100 : undef; |
|
610 |
||
57
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
611 |
$dta->{SSPEED_FROM_EC_SETTING} = 1 unless ($B6 & 0x40); # these are the aspirations |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
612 |
$dta->{DEPTH_FROM_ED_SETTING} = 1 unless ($B6 & 0x20); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
613 |
$dta->{HEADING_FROM_EH_SETTING} = 1 unless ($B6 & 0x10); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
614 |
$dta->{PITCH_FROM_EP_SETTING} = 1 unless ($B6 & 0x08); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
615 |
$dta->{ROLL_FROM_ER_SETTING} = 1 unless ($B6 & 0x04); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
616 |
$dta->{SALIN_FROM_ES_SETTING} = 1 unless ($B6 & 0x02); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
617 |
$dta->{TEMP_FROM_ET_SETTING} = 1 unless ($B6 & 0x01); |
14 | 618 |
|
57
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
619 |
$dta->{SSPEED_CALCULATED} = 1 if ($B7 & 0x40); # and this is how it is in practice |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
620 |
$dta->{PRESSURE_SENSOR} = 1 if ($B7 & 0x20); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
621 |
$dta->{COMPASS} = 1 if ($B7 & 0x10); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
622 |
$dta->{PITCH_SENSOR} = 1 if ($B7 & 0x08); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
623 |
$dta->{ROLL_SENSOR} = 1 if ($B7 & 0x04); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
624 |
$dta->{CONDUCTIVITY_SENSOR} = 1 if ($B7 & 0x02); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
625 |
$dta->{TEMPERATURE_SENSOR} = 1 if ($B7 & 0x01); |
14 | 626 |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
627 |
$dta->{DISTANCE_TO_BIN1_CENTER} /= 100; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
628 |
$dta->{TRANSMITTED_PULSE_LENGTH} /= 100; |
14 | 629 |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
630 |
$dta->{FALSE_TARGET_THRESHOLD} = undef |
14 | 631 |
if ($dta->{FALSE_TARGET_THRESHOLD} == 255); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
632 |
$dta->{TRANSMIT_LAG_DISTANCE} /= 100; |
14 | 633 |
|
57
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
634 |
if ($dta->{INSTRUMENT_TYPE} eq 'Workhorse or Signature') { |
14 | 635 |
sysread(WBRF,$buf,11) == 11 || die("$WBRcfn: $!"); |
636 |
($W1,$W2,$W3,$W4,$W5,$dta->{TRANSMIT_POWER}) = |
|
637 |
unpack('vvvvvC',$buf); |
|
638 |
||
57
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
639 |
if ($W1 == 0x6953 && $W2 == 0x6E67 && |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
640 |
$W3 == 0x7461 && $W4 == 0x0075) { |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
641 |
$dta->{INSTRUMENT_TYPE} = 'Signature'; |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
642 |
$dta->{TRANSMIT_POWER_MAX} = ($dta->{TRANSMIT_POWER} == 122); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
643 |
$dta->{PRODUCER} = 'Nortek ADCP' if ($dta->{PRODUCER} eq 'TRDI or Nortek ADCP'); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
644 |
} else { |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
645 |
$dta->{INSTRUMENT_TYPE} = 'Workhorse'; |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
646 |
$dta->{CPU_SERIAL_NUMBER} = sprintf("%04X%04X%04X%04X",$W1,$W2,$W3,$W4); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
647 |
$dta->{TRANSMIT_POWER_MAX} = ($dta->{TRANSMIT_POWER} == 255); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
648 |
$dta->{PRODUCER} = 'TRDI ADCP' if ($dta->{PRODUCER} eq 'TRDI or Nortek ADCP'); |
5a59411306ba
updated for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
56
diff
changeset
|
649 |
} |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
650 |
|
14 | 651 |
$dta->{NARROW_BANDWIDTH} = ($W5 == 1); |
652 |
$dta->{WIDE_BANDWIDTH} = ($W5 == 0); |
|
653 |
||
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
654 |
if ($dta->{FIXED_LEADER_BYTES} == 59) { # new style with serial number |
14 | 655 |
sysread(WBRF,$buf,6) == 6 || die("$WBRcfn: $!"); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
656 |
($dummy,$dta->{SERIAL_NUMBER},$dummy) = # last bytes is beam angle, but that info has |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
657 |
unpack('CVC',$buf); # already been provided above |
14 | 658 |
} |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
659 |
} |
14 | 660 |
|
661 |
if ($dta->{INSTRUMENT_TYPE} eq 'Explorer') { |
|
662 |
sysread(WBRF,$buf,16) == 16 || die("$WBRcfn: $!"); |
|
663 |
($dummy,$dummy,$W5,$dummy,$dta->{SERIAL_NUMBER}) = |
|
664 |
unpack('VVvvV',$buf); |
|
665 |
$dta->{NARROW_BANDWIDTH} = ($W5 == 1); |
|
666 |
$dta->{WIDE_BANDWIDTH} = ($W5 == 0); |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
667 |
} |
14 | 668 |
|
669 |
#----------------------- |
|
670 |
# 1st ENSEMBLE, BT Setup |
|
671 |
#----------------------- |
|
672 |
||
56
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
673 |
if ($dta->{BT_PRESENT}) { |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
674 |
sysseek(WBRF,$start_ens+$WBRofs[$BT_dt],0) || die("$WBRcfn: $!"); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
675 |
sysread(WBRF,$buf,12) == 12 || die("$WBRcfn: $!"); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
676 |
($id,$dta->{BT_PINGS_PER_ENSEMBLE},$dta->{BT_DELAY_BEFORE_REACQUIRE}, |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
677 |
$dta->{BT_MIN_CORRELATION},$dta->{BT_MIN_EVAL_AMPLITUDE}, |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
678 |
$dta->{BT_MIN_PERCENT_GOOD},$dta->{BT_MODE}, |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
679 |
$dta->{BT_MAX_ERROR_VELOCITY}) = unpack('vvvCCCCv',$buf); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
680 |
|
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
681 |
$id == 0x0600 || |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
682 |
printf(STDERR $FmtErr."\n",$WBRcfn,"Bottom Track",$id,0,tell(WBRF)); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
683 |
|
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
684 |
$dta->{BT_MAX_ERROR_VELOCITY} = |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
685 |
$dta->{BT_MAX_ERROR_VELOCITY} ? $dta->{BT_MAX_ERROR_VELOCITY} / 1000 |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
686 |
: undef; |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
687 |
|
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
688 |
sysseek(WBRF,28,1) || die("$WBRcfn: $!"); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
689 |
sysread(WBRF,$buf,6) == 6 || die("$WBRcfn: $!"); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
690 |
($dta->{BT_RL_MIN_SIZE},$dta->{BT_RL_NEAR},$dta->{BT_RL_FAR}) |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
691 |
= unpack('vvv',$buf); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
692 |
|
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
693 |
$dta->{BT_RL_MIN_SIZE} /= 10; |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
694 |
$dta->{BT_RL_NEAR} /= 10; |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
695 |
$dta->{BT_RL_FAR} /= 10; |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
696 |
|
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
697 |
sysseek(WBRF,20,1) || die("$WBRcfn: $!"); # skip data |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
698 |
sysread(WBRF,$buf,2) == 2 || die("$WBRcfn: $!"); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
699 |
$dta->{BT_MAX_TRACKING_DEPTH} = unpack('v',$buf) / 10; |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
700 |
} |
14 | 701 |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
702 |
return $dta; |
14 | 703 |
} |
704 |
||
705 |
#---------------------------------------------------------------------- |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
706 |
# readData(file_name,^data[,first_ens,last_ens[,last_bin]]) |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
707 |
# - read ensembles |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
708 |
# - read all ensembles unless first_ens and last_ens are given |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
709 |
# - read all bins unless last_bin is given |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
710 |
# - NB: CODE GETS PROGRESSIVELY SLOWER => HUGE PAIN FOR LARGE FILES |
14 | 711 |
#---------------------------------------------------------------------- |
712 |
||
713 |
sub readData(@) |
|
714 |
{ |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
715 |
my($fn,$dta,$fe,$le,$lb) = @_; |
14 | 716 |
$WBRcfn = $fn; |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
717 |
open(WBRF,$WBRcfn) || die("$WBRcfn: $!\n"); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
718 |
WBRhdr($dta) || die("$WBRcfn: Insufficient Data\n"); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
719 |
$lb = $dta->{N_BINS} |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
720 |
unless (numberp($lb) && $lb>=1 && $lb<=$dta->{N_BINS}); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
721 |
WBRens($lb,$dta->{FIXED_LEADER_BYTES},\@{$dta->{ENSEMBLE}},$fe,$le); |
14 | 722 |
print(STDERR "$WBRcfn: $BIT_errors built-in-test errors\n") |
723 |
if ($BIT_errors); |
|
724 |
} |
|
725 |
||
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
726 |
sub WBRens(@) |
14 | 727 |
{ |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
728 |
my($nbins,$fixed_leader_bytes,$E,$fe,$le) = @_; |
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
729 |
my($B1,$B2,$B3,$B4,$I,$bin,$beam,$dummy,@dta,$i,$cs); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
730 |
my($ens,$ensNo,$dayStart,$ens_length,$hid,$did,$el); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
731 |
local our($ndt,$buf,$id,$start_ens,@WBRofs); |
14 | 732 |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
733 |
sysseek(WBRF,0,0) || die("$WBRcfn: $!"); |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
734 |
ENSEMBLE: |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
735 |
for ($ens=0; 1; $ens++) { |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
736 |
print(STDERR '.') if ($RDI_PD0_IO::show_progress && $ens % 1000 == 0); |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
737 |
$start_ens = goto_next_ens(\*WBRF); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
738 |
last unless defined($start_ens); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
739 |
|
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
740 |
#---------------------------------------- |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
741 |
# Handle first_ens and last_ens |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
742 |
#---------------------------------------- |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
743 |
|
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
744 |
if (defined($fe) && $ens>0 && ${$E}[$ens-1]->{NUMBER}<$fe) { # delete previous ensemble |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
745 |
pop(@{$E}); $ens--; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
746 |
} |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
747 |
|
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
748 |
if (defined($le) && $ens>0 && ${$E}[$ens-1]->{NUMBER}>$le) { # delete previous ensemble and finish |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
749 |
pop(@{$E}); $ens--; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
750 |
last; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
751 |
} |
14 | 752 |
|
753 |
#---------------------------------------- |
|
754 |
# Get ensemble length and # of data types |
|
755 |
#---------------------------------------- |
|
756 |
||
757 |
sysseek(WBRF,$start_ens,0) || die("$WBRcfn: $!"); |
|
758 |
sysread(WBRF,$buf,6) == 6 || last; |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
759 |
($hid,$did,$el,$dummy,$ndt) = unpack('CCvCC',$buf); |
37 | 760 |
$hid == 0x7f || die(sprintf($FmtErr,$WBRcfn,"Header",$hid,defined($ensNo)?$ensNo+1:0)); |
29 | 761 |
${$E}[$ens]->{DATA_SOURCE_ID} = $did; |
762 |
if ($did == 0x7f) { |
|
763 |
${$E}[$ens]->{PRODUCER} = 'TRDI ADCP'; |
|
764 |
} elsif ($did&0xF0 == 0xA0) { |
|
765 |
${$E}[$ens]->{PRODUCER} = 'IMP+LADCP (Thurnherr software)'; |
|
31 | 766 |
} elsif ($did&0xF0 == 0xE0) { |
29 | 767 |
${$E}[$ens]->{PRODUCER} = 'editPD0 (Thurnherr software)'; |
768 |
} else { |
|
769 |
${$E}[$ens]->{PRODUCER} = 'unknown'; |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
770 |
} |
29 | 771 |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
772 |
if (defined($ens_length) && ($el != $ens_length)) { |
50 | 773 |
$RDI_PD0_IO::File_Dirty = 1; |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
774 |
print(STDERR "WARNING (RDI_PD0_IO): ensemble ${$E}[$#{$E}]->{NUMBER} skipped (unexpected length)\n"); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
775 |
pop(@{$E}); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
776 |
$ens--; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
777 |
next; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
778 |
} |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
779 |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
780 |
$ens_length = $el; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
781 |
|
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
782 |
## printf(STDERR "$WBRcfn: WARNING: unexpected number of data types (%d, ens=$ens)\n",$ndt),last |
18 | 783 |
## unless ($ndt == 6 || $ndt == 7); |
50 | 784 |
my($nread) = sysread(WBRF,$buf,2*$ndt); # 2019 EPR test |
785 |
if ($nread != 2*$ndt) { |
|
786 |
printf(STDERR "$WBRcfn: WARNING: expected to read %d bytes, got only %d in ensemble %d\n", |
|
787 |
2*$ndt,$nread,${$E}[$ens]->{NUMBER}); |
|
788 |
last; |
|
789 |
} |
|
790 |
||
791 |
||
14 | 792 |
@WBRofs = unpack("v$ndt",$buf); |
793 |
$fixed_leader_bytes = $WBRofs[1] - $WBRofs[0]; |
|
794 |
# print(STDERR "@WBRofs\n"); |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
795 |
|
14 | 796 |
#------------------------------- |
797 |
# Make Sure Ensemble is Complete |
|
798 |
#------------------------------- |
|
799 |
||
800 |
# UH BB150 writes incomplete ensembles (i.e. short read |
|
801 |
# indicates EOF). FSU WH300 has bogus data in incomplete |
|
802 |
# final ensemble. |
|
803 |
||
804 |
sysseek(WBRF,$start_ens,0) || die("$WBRcfn: $!"); |
|
37 | 805 |
unless ((sysread(WBRF,$buf,$ens_length) == $ens_length) && # incomplete ensemble |
36
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
806 |
(sysread(WBRF,$cs,2) == 2)) { |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
807 |
# print(STDERR "INCOMPLETE ENSEMBLE\n"); |
36
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
808 |
pop(@{$E}); |
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
809 |
last; |
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
810 |
} |
14 | 811 |
|
37 | 812 |
unless (unpack('%16C*',$buf) == unpack('v',$cs)) { # bad checksum |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
813 |
# print(STDERR "BAD CHECKSUM\n"); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
814 |
pop(@{$E}); $ens--; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
815 |
next; |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
816 |
} |
14 | 817 |
|
818 |
#------------------------------ |
|
819 |
# Variable Leader |
|
820 |
#------------------------------ |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
821 |
|
37 | 822 |
my($lastEns) = $ensNo; |
14 | 823 |
sysseek(WBRF,$start_ens+$WBRofs[1],0) || die("$WBRcfn: $!"); |
824 |
sysread(WBRF,$buf,4) == 4 || die("$WBRcfn: $!"); |
|
37 | 825 |
($id,$ensNo) = unpack("vv",$buf); # only lower two bytes!!! |
14 | 826 |
|
39 | 827 |
if (${$E}[$ens]->{INSTRUMENT_TYPE} eq 'Ocean Surveyor') { |
828 |
$id == 0x0081 || |
|
829 |
die(sprintf($FmtErr,$WBRcfn,"Variable Leader",$id,$ensNo + ($lastEns - ($lastEns & 0xFFFF)))); |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
830 |
} else { |
39 | 831 |
$id == 0x0080 || |
832 |
die(sprintf($FmtErr,$WBRcfn,"Variable Leader",$id,$ensNo + ($lastEns - ($lastEns & 0xFFFF)))); |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
833 |
} |
14 | 834 |
|
42 | 835 |
# if ($fixed_leader_bytes==42 || $fixed_leader_bytes==58) { # BB150 & Explorer DVL (if DISABLED!) |
836 |
sysread(WBRF,$buf,7) == 7 || die("$WBRcfn: $!"); # always read pre-Y2K clock |
|
14 | 837 |
(${$E}[$ens]->{YEAR},${$E}[$ens]->{MONTH}, |
838 |
${$E}[$ens]->{DAY},${$E}[$ens]->{HOUR},${$E}[$ens]->{MINUTE}, |
|
50 | 839 |
${$E}[$ens]->{SECOND},$B4) = unpack('CCCCCCC',$buf); |
840 |
${$E}[$ens]->{SECOND} += $B4/100; |
|
14 | 841 |
${$E}[$ens]->{YEAR} += (${$E}[$ens]->{YEAR} > 80) ? 1900 : 2000; |
40
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
842 |
# } else { |
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
843 |
# sysseek(WBRF,7,1) || die("$WBRcfn: $!"); # use Y2K RTC instead |
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
844 |
# } |
14 | 845 |
|
846 |
sysread(WBRF,$buf,1) == 1 || die("$WBRcfn: $!"); |
|
847 |
$ensNo += unpack('C',$buf) << 16; |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
848 |
|
52
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
849 |
## CODE DISABLED BECAUSE IT TAKES TOO MUCH TIME ON LARGE FILES. |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
850 |
## NOT SURE HOW COMMON DUPLICATE ENSEMBLES ARE. |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
851 |
## |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
852 |
## for (my($i)=$ens; $i>0; $i--) { # check for duplicate ens; e.g. 2018 S4P 24UL |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
853 |
## if (${$E}[$i]->{NUMBER} == $ensNo) { |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
854 |
## $RDI_PD0_IO::File_Dirty = 1; |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
855 |
## print(STDERR "WARNING (RDI_PD0_IO): duplicate ensemble $ensNo skipped\n"); |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
856 |
## pop(@{$E}); |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
857 |
## $ens--; |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
858 |
## next ENSEMBLE; |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
859 |
## } |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
860 |
## } |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
861 |
|
14 | 862 |
${$E}[$ens]->{NUMBER} = $ensNo; |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
863 |
|
14 | 864 |
sysread(WBRF,$buf,30) == 30 || die("$WBRcfn: $!"); |
865 |
(${$E}[$ens]->{BUILT_IN_TEST_ERROR},${$E}[$ens]->{SPEED_OF_SOUND}, |
|
866 |
${$E}[$ens]->{XDUCER_DEPTH},${$E}[$ens]->{HEADING}, |
|
867 |
${$E}[$ens]->{PITCH},${$E}[$ens]->{ROLL}, |
|
868 |
${$E}[$ens]->{SALINITY},${$E}[$ens]->{TEMPERATURE}, |
|
869 |
${$E}[$ens]->{MIN_PRE_PING_WAIT_TIME},$B1,$B2, |
|
870 |
${$E}[$ens]->{HEADING_STDDEV},${$E}[$ens]->{PITCH_STDDEV}, |
|
871 |
${$E}[$ens]->{ROLL_STDDEV},${$E}[$ens]->{ADC_XMIT_CURRENT}, |
|
872 |
${$E}[$ens]->{ADC_XMIT_VOLTAGE},${$E}[$ens]->{ADC_AMBIENT_TEMPERATURE}, |
|
873 |
${$E}[$ens]->{ADC_PRESSURE_PLUS},${$E}[$ens]->{ADC_PRESSURE_MINUS}, |
|
874 |
${$E}[$ens]->{ADC_ATTITUDE_TEMPERATURE},${$E}[$ens]->{ADC_ATTITUDE}, |
|
875 |
${$E}[$ens]->{ADC_CONTAMINATION}) |
|
876 |
= unpack('vvvvvvvvCCCCCCCCCCCCCC',$buf); |
|
877 |
||
878 |
${$E}[$ens]->{BUILT_IN_TEST_ERROR} = undef |
|
879 |
unless (${$E}[$ens]->{BUILT_IN_TEST_ERROR}); |
|
880 |
$BIT_errors++ if (${$E}[$ens]->{BUILT_IN_TEST_ERROR}); |
|
881 |
||
882 |
${$E}[$ens]->{XDUCER_DEPTH} /= 10; |
|
18 | 883 |
|
884 |
#------------------------------------------------- |
|
885 |
# IMP EXTENSION: PITCH/ROLL/HEADING CAN BE MISSING |
|
886 |
#------------------------------------------------- |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
887 |
|
18 | 888 |
${$E}[$ens]->{HEADING} = (${$E}[$ens]->{HEADING} == 0xF000) |
889 |
? undef |
|
890 |
: ${$E}[$ens]->{HEADING} / 100; |
|
891 |
${$E}[$ens]->{PITCH} = (${$E}[$ens]->{PITCH} == 0x8000) |
|
892 |
? undef |
|
893 |
: unpack('s',pack('S',${$E}[$ens]->{PITCH})) / 100; |
|
894 |
${$E}[$ens]->{ROLL} = (${$E}[$ens]->{ROLL} == 0x8000) |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
895 |
? undef |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
896 |
: unpack('s',pack('S',${$E}[$ens]->{ROLL})) / 100; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
897 |
|
18 | 898 |
${$E}[$ens]->{TEMPERATURE} = unpack('s',pack('S',${$E}[$ens]->{TEMPERATURE})) / 100; |
14 | 899 |
${$E}[$ens]->{MIN_PRE_PING_WAIT_TIME} *= 60; |
900 |
${$E}[$ens]->{MIN_PRE_PING_WAIT_TIME} += $B1 + $B2/100; |
|
901 |
${$E}[$ens]->{PITCH_STDDEV} /= 10; |
|
902 |
${$E}[$ens]->{ROLL_STDDEV} /= 10; |
|
903 |
||
40
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
904 |
if (($fixed_leader_bytes==53 || $fixed_leader_bytes==59) && # Workhorse instruments |
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
905 |
!defined($RDI_PD0_IO::IGNORE_Y2K_CLOCK)) { |
14 | 906 |
sysread(WBRF,$buf,23) == 23 || die("$WBRcfn: $!"); |
907 |
(${$E}[$ens]->{ERROR_STATUS_WORD}, |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
908 |
$dummy,${$E}[$ens]->{PRESSURE},${$E}[$ens]->{PRESSURE_STDDEV}, |
14 | 909 |
$dummy,${$E}[$ens]->{YEAR},$B3,${$E}[$ens]->{MONTH}, |
910 |
${$E}[$ens]->{DAY},${$E}[$ens]->{HOUR},${$E}[$ens]->{MINUTE}, |
|
50 | 911 |
${$E}[$ens]->{SECOND},$B4) |
14 | 912 |
= unpack('VvVVCCCCCCCCC',$buf); |
913 |
||
914 |
${$E}[$ens]->{PRESSURE} /= 1000; |
|
915 |
${$E}[$ens]->{PRESSURE_STDDEV} /= 1000; |
|
916 |
${$E}[$ens]->{YEAR} *= 100; ${$E}[$ens]->{YEAR} += $B3; |
|
50 | 917 |
${$E}[$ens]->{SECOND} += $B4/100; |
14 | 918 |
} |
919 |
||
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
920 |
# THE FOLLOWING LINE OF CODE WAS REMOVED 7/30/2016 WHEN I ADDED A POP |
36
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
921 |
# TO THE last STATEMENT ABOVE (INCOMPLETE ENSEMBLE) |
42 | 922 |
# THE LINE WAS RE-ENABLED ON 12/23/2017 BECAUSE OTHERWISE THE |
923 |
# ANSLOPE II PROFILES IN THE HOWTO CANNOT BE READ. |
|
924 |
pop(@{$E}),last if (${$E}[$ens]->{MONTH}>12); # 10/15/2014; IWISE#145 UL ??? |
|
21 | 925 |
|
14 | 926 |
if ($fixed_leader_bytes == 58) { # Explorer DVL |
927 |
sysread(WBRF,$buf,14) == 14 || die("$WBRcfn: $!"); |
|
928 |
(${$E}[$ens]->{ERROR_STATUS_WORD}, |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
929 |
$dummy,${$E}[$ens]->{PRESSURE},${$E}[$ens]->{PRESSURE_STDDEV}) |
14 | 930 |
= unpack('VvVV',$buf); |
931 |
${$E}[$ens]->{PRESSURE} /= 1000; |
|
932 |
${$E}[$ens]->{PRESSURE_STDDEV} /= 1000; |
|
933 |
} |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
934 |
|
14 | 935 |
${$E}[$ens]->{DATE} |
936 |
= sprintf("%02d/%02d/%d",${$E}[$ens]->{MONTH}, |
|
937 |
${$E}[$ens]->{DAY}, |
|
938 |
${$E}[$ens]->{YEAR}); |
|
939 |
${$E}[$ens]->{TIME} |
|
940 |
= sprintf("%02d:%02d:%05.02f",${$E}[$ens]->{HOUR}, |
|
941 |
${$E}[$ens]->{MINUTE}, |
|
50 | 942 |
${$E}[$ens]->{SECOND}); |
14 | 943 |
${$E}[$ens]->{DAYNO} |
40
6a46e9d31106
after Hamburg visit
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
39
diff
changeset
|
944 |
= &_dayNo(${$E}[$ens]->{YEAR},${$E}[$ens]->{MONTH},${$E}[$ens]->{DAY}, |
50 | 945 |
${$E}[$ens]->{HOUR},${$E}[$ens]->{MINUTE},${$E}[$ens]->{SECOND}); |
14 | 946 |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
947 |
# when analyzing an STA file from an OS75 SADCP (Poseidon), |
14 | 948 |
# I noticed that there is no time information. This causes |
949 |
# timegm to bomb. |
|
950 |
if (${$E}[$ens]->{MONTH} == 0) { # no time info |
|
951 |
${$E}[$ens]->{UNIX_TIME} = 0; |
|
952 |
${$E}[$ens]->{SECNO} = 0; |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
953 |
} else { |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
954 |
# print(STDERR "[$ens]->${$E}[$ens]->{MINUTE}:${$E}[$ens]->{HOUR},${$E}[$ens]->{DAY},${$E}[$ens]->{MONTH},${$E}[$ens]->{YEAR}-<\n"); |
14 | 955 |
${$E}[$ens]->{UNIX_TIME} |
956 |
= timegm(0,${$E}[$ens]->{MINUTE}, |
|
957 |
${$E}[$ens]->{HOUR}, |
|
958 |
${$E}[$ens]->{DAY}, |
|
959 |
${$E}[$ens]->{MONTH}-1, # timegm jan==0!!! |
|
960 |
${$E}[$ens]->{YEAR}) |
|
50 | 961 |
+ ${$E}[$ens]->{SECOND}; |
14 | 962 |
|
963 |
$dayStart = timegm(0,0,0,${$E}[$ens]->{DAY}, |
|
964 |
${$E}[$ens]->{MONTH}-1, |
|
965 |
${$E}[$ens]->{YEAR}) |
|
966 |
unless defined($dayStart); |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
967 |
${$E}[$ens]->{SECNO} = ${$E}[$ens]->{UNIX_TIME} - $dayStart; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
968 |
} |
14 | 969 |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
970 |
sysseek(WBRF,$start_ens+$WBRofs[0]+4,0) # System Config / Fixed Leader |
14 | 971 |
|| die("$WBRcfn: $!"); |
972 |
||
973 |
sysread(WBRF,$buf,5) == 5 || die("$WBRcfn: $!"); |
|
974 |
($B1,$dummy,$dummy,$dummy,${$E}[$ens]->{N_BEAMS_USED}) |
|
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
975 |
= unpack('CCCCC',$buf); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
976 |
${$E}[$ens]->{XDUCER_FACING_UP} = 1 if ($B1 & 0x80); |
14 | 977 |
${$E}[$ens]->{XDUCER_FACING_DOWN} = 1 unless ($B1 & 0x80); |
978 |
||
979 |
#-------------------- |
|
980 |
# Velocity Data |
|
981 |
#-------------------- |
|
982 |
||
983 |
my($ndata) = $nbins * 4; |
|
984 |
||
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
985 |
my($vel_di) = WBRdtaIndex(0x0100); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
986 |
die("no velocity data in ensemble #$ensNo\n") |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
987 |
unless defined($vel_di); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
988 |
|
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
989 |
sysseek(WBRF,$start_ens+$WBRofs[$vel_di],0) || die("$WBRcfn: $!"); |
14 | 990 |
sysread(WBRF,$buf,2+$ndata*2) == 2+$ndata*2 || die("$WBRcfn: $!"); |
991 |
($id,@dta) = unpack("vv$ndata",$buf); |
|
992 |
||
993 |
for ($i=0,$bin=0; $bin<$nbins; $bin++) { |
|
994 |
for ($beam=0; $beam<4; $beam++,$i++) { |
|
995 |
${$E}[$ens]->{VELOCITY}[$bin][$beam] = |
|
996 |
unpack('s',pack('S',$dta[$i])) / 1000 |
|
997 |
if ($dta[$i] != 0x8000); |
|
998 |
} |
|
999 |
} |
|
1000 |
||
1001 |
#-------------------- |
|
1002 |
# Correlation Data |
|
1003 |
#-------------------- |
|
1004 |
||
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1005 |
my($corr_di) = WBRdtaIndex(0x0200); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1006 |
die("no correlation data in ensemble #$ensNo\n") |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1007 |
unless defined($corr_di); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
1008 |
|
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1009 |
sysseek(WBRF,$start_ens+$WBRofs[$corr_di],0) || die("$WBRcfn: $!"); |
14 | 1010 |
sysread(WBRF,$buf,2+$ndata) == 2+$ndata || die("$WBRcfn: $!"); |
1011 |
($id,@dta) = unpack("vC$ndata",$buf); |
|
1012 |
||
1013 |
for ($i=0,$bin=0; $bin<$nbins; $bin++) { |
|
1014 |
for ($beam=0; $beam<4; $beam++,$i++) { |
|
1015 |
${$E}[$ens]->{CORRELATION}[$bin][$beam] = $dta[$i] |
|
1016 |
if ($dta[$i]); |
|
1017 |
} |
|
1018 |
} |
|
1019 |
||
1020 |
#-------------------- |
|
1021 |
# Echo Intensity Data |
|
1022 |
#-------------------- |
|
1023 |
||
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1024 |
my($echo_di) = WBRdtaIndex(0x0300); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1025 |
die("no echo intensity data in ensemble #$ensNo\n") |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1026 |
unless defined($echo_di); |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
1027 |
|
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1028 |
sysseek(WBRF,$start_ens+$WBRofs[$echo_di],0) || die("$WBRcfn: $!"); |
14 | 1029 |
sysread(WBRF,$buf,2+$ndata) == 2+$ndata || die("$WBRcfn: $!"); |
1030 |
($id,@dta) = unpack("vC$ndata",$buf); |
|
1031 |
||
1032 |
$id == 0x0300 || |
|
37 | 1033 |
die(sprintf($FmtErr,$WBRcfn,"Echo Intensity",$id,$ensNo)); |
14 | 1034 |
|
1035 |
for ($i=0,$bin=0; $bin<$nbins; $bin++) { |
|
1036 |
for ($beam=0; $beam<4; $beam++,$i++) { |
|
1037 |
${$E}[$ens]->{ECHO_AMPLITUDE}[$bin][$beam] = $dta[$i]; |
|
1038 |
} |
|
1039 |
} |
|
1040 |
||
1041 |
#-------------------- |
|
1042 |
# Percent Good Data |
|
1043 |
#-------------------- |
|
1044 |
||
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1045 |
my($pctg_di) = WBRdtaIndex(0x0400); |
55 | 1046 |
if (defined($pctg_di)) { |
1047 |
sysseek(WBRF,$start_ens+$WBRofs[$pctg_di],0) || die("$WBRcfn: $!"); |
|
1048 |
sysread(WBRF,$buf,2+$ndata) == 2+$ndata || die("$WBRcfn: $!"); |
|
1049 |
($id,@dta) = unpack("vC$ndata",$buf); |
|
1050 |
$id == 0x0400 || |
|
1051 |
die(sprintf($FmtErr,$WBRcfn,"Percent-Good Data",$id,$ensNo)); |
|
1052 |
for ($i=0,$bin=0; $bin<$nbins; $bin++) { |
|
1053 |
for ($beam=0; $beam<4; $beam++,$i++) { |
|
1054 |
${$E}[$ens]->{PERCENT_GOOD}[$bin][$beam] = $dta[$i]; |
|
1055 |
} |
|
1056 |
} |
|
1057 |
} else { |
|
1058 |
printf(STDERR "\nWARNING (RDI_PD0_IO): Percent-Good data missing from PD0 file\n") |
|
1059 |
if ($ens == 0); |
|
1060 |
for ($i=0,$bin=0; $bin<$nbins; $bin++) { |
|
1061 |
for ($beam=0; $beam<4; $beam++,$i++) { |
|
1062 |
${$E}[$ens]->{PERCENT_GOOD}[$bin][$beam] = 100; # fake it |
|
1063 |
} |
|
1064 |
} |
|
1065 |
} |
|
14 | 1066 |
|
18 | 1067 |
#----------------------------------------- |
14 | 1068 |
# Bottom-Track Data |
18 | 1069 |
# - scan through remaining data types |
1070 |
#----------------------------------------- |
|
14 | 1071 |
|
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1072 |
my($bt_di) = WBRdtaIndex(0x0600); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1073 |
unless (defined($pctg_di)) { # no BT found => next ens |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1074 |
sysseek(WBRF,$start_ens+$ens_length+2,0) || die("$WBRcfn: $!"); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1075 |
next; |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
1076 |
} |
18 | 1077 |
|
32 | 1078 |
sysseek(WBRF,14,1) || die("$WBRcfn: $!"); # BT range, velocity, corr, %-good, ... |
18 | 1079 |
sysread(WBRF,$buf,28) == 28 || die("$WBRcfn: $!"); |
1080 |
@dta = unpack('v4v4C4C4C4',$buf); |
|
1081 |
for ($beam=0; $beam<4; $beam++) { |
|
32 | 1082 |
${$E}[$ens]->{BT_RANGE}[$beam] = $dta[$beam] / 100 # lower 2 bytes only! |
1083 |
if ($dta[$beam]); # (see below for high bytes) |
|
18 | 1084 |
} |
1085 |
for ($beam=0; $beam<4; $beam++) { |
|
1086 |
${$E}[$ens]->{BT_VELOCITY}[$beam] = |
|
1087 |
unpack('s',pack('S',$dta[4+$beam])) / 1000 |
|
1088 |
if ($dta[4+$beam] != 0x8000); |
|
1089 |
} |
|
1090 |
for ($beam=0; $beam<4; $beam++) { |
|
1091 |
${$E}[$ens]->{BT_CORRELATION}[$beam] = $dta[8+$beam] |
|
1092 |
if ($dta[8+$beam]); |
|
1093 |
} |
|
32 | 1094 |
for ($beam=0; $beam<4; $beam++) { # BT filter parameter |
18 | 1095 |
${$E}[$ens]->{BT_EVAL_AMPLITUDE}[$beam] = $dta[12+$beam]; |
1096 |
} |
|
1097 |
for ($beam=0; $beam<4; $beam++) { |
|
1098 |
${$E}[$ens]->{BT_PERCENT_GOOD}[$beam] = $dta[16+$beam]; |
|
1099 |
} |
|
1100 |
||
32 | 1101 |
sysseek(WBRF,6,1) || die("$WBRcfn: $!"); # BT ref level stuff |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
1102 |
if (sysread(WBRF,$buf,20) != 20) { # EN642/PITA1 |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
1103 |
pop(@{$E}); |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
1104 |
last; |
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
1105 |
} |
18 | 1106 |
@dta = unpack('v4C4C4C4',$buf); |
1107 |
for ($beam=0; $beam<4; $beam++) { |
|
1108 |
${$E}[$ens]->{BT_RL_VELOCITY}[$beam] = |
|
1109 |
unpack('s',pack('S',$dta[$beam])) / 1000 |
|
1110 |
if ($dta[$beam] != 0x8000); |
|
1111 |
} |
|
1112 |
for ($beam=0; $beam<4; $beam++) { |
|
1113 |
${$E}[$ens]->{BT_RL_CORRELATION}[$beam] = $dta[4+$beam] |
|
1114 |
if ($dta[4+$beam]); |
|
1115 |
} |
|
1116 |
for ($beam=0; $beam<4; $beam++) { |
|
1117 |
${$E}[$ens]->{BT_RL_ECHO_AMPLITUDE}[$beam] = $dta[8+$beam]; |
|
1118 |
} |
|
1119 |
for ($beam=0; $beam<4; $beam++) { |
|
1120 |
${$E}[$ens]->{BT_RL_PERCENT_GOOD}[$beam] = $dta[12+$beam]; |
|
1121 |
} |
|
1122 |
||
32 | 1123 |
sysseek(WBRF,2,1) || die("$WBRcfn: $!"); # BT signal strength & BT range high bytes |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1124 |
# sysread(WBRF,$buf,9) == 9 || die("$WBRcfn: $!"); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1125 |
if (sysread(WBRF,$buf,9) == 9) { # SR1b JR16 BB150 data files require this |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1126 |
@dta = unpack('C4CC4',$buf); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1127 |
for ($beam=0; $beam<4; $beam++) { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1128 |
${$E}[$ens]->{BT_SIGNAL_STRENGTH}[$beam] = $dta[$beam]; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1129 |
} |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1130 |
${$E}[$ens]->{HIGH_GAIN} if ($dta[4]); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1131 |
${$E}[$ens]->{LOW_GAIN} unless ($dta[4]); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1132 |
for ($beam=0; $beam<4; $beam++) { # high bytes (1 byte per beam) |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1133 |
${$E}[$ens]->{BT_RANGE}[$beam] += $dta[5+$beam] * 655.36 |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1134 |
if ($dta[5+$beam]); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1135 |
} |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1136 |
# sysseek(WBRF,8,1) || die("$WBRcfn: $!"); # remainder of ensemble |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1137 |
} |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1138 |
sysseek(WBRF,$start_ens+$ens_length+2,0) || die("$WBRcfn: $!"); |
14 | 1139 |
} # ens loop |
51
148c092b3a09
after somewhat suspicous RDI_PD0_IO.pl merge and just before adding
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
50
diff
changeset
|
1140 |
print(STDERR "\n") if ($RDI_PD0_IO::show_progress); |
14 | 1141 |
} |
1142 |
||
52
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
1143 |
sub WBRdtaIndex($) # return index of particular data type |
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
1144 |
{ # where index refers to data sections inside PD0 files |
48
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1145 |
my($trgid) = @_; |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1146 |
our($ndt,$buf,$id,$start_ens,@WBRofs); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1147 |
|
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1148 |
for (my($di)=2; $di<$ndt; $di++) { |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1149 |
sysseek(WBRF,$start_ens+$WBRofs[$di],0) || die("$WBRcfn: $!"); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1150 |
sysread(WBRF,$buf,2) == 2 || die("$WBRcfn: $!"); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1151 |
$id = unpack('v',$buf); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1152 |
return $di if ($id == $trgid); |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1153 |
} |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1154 |
return undef; |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1155 |
} |
cdc74ebada81
after partial adaptation to RTI files
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
45
diff
changeset
|
1156 |
|
56
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1157 |
sub WBRdtaTypes() # return list of data types |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1158 |
{ |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1159 |
my(@dt); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1160 |
our($ndt,$buf,$id,$start_ens,@WBRofs); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1161 |
|
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1162 |
for (my($di)=2; $di<$ndt; $di++) { |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1163 |
sysseek(WBRF,$start_ens+$WBRofs[$di],0) || die("$WBRcfn: $!"); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1164 |
sysread(WBRF,$buf,2) == 2 || die("$WBRcfn: $!"); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1165 |
$id = unpack('v',$buf); |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1166 |
if ($id == 0x0000) { push(@dt,'FixedLeader'); } |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1167 |
elsif ($id == 0x0080) { push(@dt,'VariableLeader'); } |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1168 |
elsif ($id == 0x0081) { push(@dt,'VariableLeader'); } # TRDI Ocean Surveyor |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1169 |
elsif ($id == 0x0100) { push(@dt,'VELOCITY'); } |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1170 |
elsif ($id == 0x0200) { push(@dt,'CORRELATION'); } |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1171 |
elsif ($id == 0x0300) { push(@dt,'ECHO_AMPLITUDE'); } |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1172 |
elsif ($id == 0x0400) { push(@dt,'PERCENT_GOOD'); } |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1173 |
elsif ($id == 0x0600) { push(@dt,'BOTTOM_TRACK'); } |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1174 |
else { push(@dt,'Unknown'); } |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1175 |
} |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1176 |
return @dt; |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1177 |
} |
82384f531dba
improvements for Nortek
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
55
diff
changeset
|
1178 |
|
14 | 1179 |
#---------------------------------------------------------------------- |
28 | 1180 |
# writeData(output_file_name,^data) WBPens(nbins,fixed_leader_bytes,^data) |
1181 |
# - writeData() copies file previously read with readData() to output_file_name |
|
1182 |
# - WBPens() patches new PD0 file with ^data |
|
1183 |
# - ^data is modified!!!! |
|
1184 |
# - output file must already exist and have correct structure |
|
1185 |
# - only subset of data structure is patched: |
|
1186 |
# - Header: Data Source Id |
|
1187 |
# - Var Ldr: Soundspeed, Depth, Heading, Pitch, Roll, Temp, Salin |
|
1188 |
# - Data: Velocity, Correlation, Echo Amp, %-Good, |
|
14 | 1189 |
#---------------------------------------------------------------------- |
1190 |
||
1191 |
sub writeData(@) |
|
1192 |
{ |
|
1193 |
my($fn,$dta) = @_; |
|
1194 |
||
50 | 1195 |
die("writeData() needs \$WBRcfn from previous readData()\n") |
31 | 1196 |
unless (length($WBRcfn) > 0); |
50 | 1197 |
die("writeData() only works with clean PD0 files\n") |
1198 |
if ($RDI_PD0_IO::File_Dirty); |
|
31 | 1199 |
|
1200 |
sysseek(WBRF,0,0) || die("$WBRcfn: $!"); # rewind input file |
|
1201 |
$WBPcfn = $fn; # set patch file name for error messages |
|
1202 |
open(WBPF,"+>$WBPcfn") || die("$WBPcfn: $!"); # open patch file for r/w |
|
14 | 1203 |
|
31 | 1204 |
while (1) { # copy input file to patch file |
1205 |
my($buf); |
|
1206 |
my($nread) = sysread(WBRF,$buf,100*1024); |
|
1207 |
die("$WBRcfn: $!\n") if ($nread < 0); |
|
1208 |
last if ($nread == 0); |
|
1209 |
my($nwritten) = syswrite(WBPF,$buf,100*1024); |
|
1210 |
die("$WBPcfn: $! ($nwritten of $nread written)\n") |
|
1211 |
unless ($nwritten = $nread); |
|
1212 |
} |
|
1213 |
sysseek(WBPF,0,0) || die("$WBPcfn: $!"); # rewind patch file |
|
1214 |
||
1215 |
WBPens($dta->{N_BINS},$dta->{FIXED_LEADER_BYTES},$dta); |
|
14 | 1216 |
} |
1217 |
||
37 | 1218 |
sub _round(@) |
28 | 1219 |
{ |
1220 |
return $_[0] >= 0 ? int($_[0] + 0.5) |
|
1221 |
: int($_[0] - 0.5); |
|
1222 |
} |
|
1223 |
||
1224 |
||
14 | 1225 |
sub WBPens($$$) |
1226 |
{ |
|
31 | 1227 |
my($nbins,$fixed_leader_bytes,$dta) = @_; |
14 | 1228 |
my($start_ens,$B1,$B2,$B3,$B4,$I,$id,$bin,$beam,$buf,$dummy,@dta,$i,$cs,@WBPofs); |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1229 |
my($ens,$dayStart,$ens_length,$hid,$ndt,$el); |
14 | 1230 |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1231 |
# for ($ens=$start_ens=0; $ens<=$#{$dta->{ENSEMBLE}}; $ens++,$start_ens+=$ens_length+2) { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1232 |
for ($ens=0; $ens<=$#{$dta->{ENSEMBLE}}; $ens++) { |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1233 |
$start_ens = goto_next_ens(\*WBPF); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1234 |
die("ens = $ens\n") unless defined($start_ens); |
14 | 1235 |
|
18 | 1236 |
#------------------------------ |
1237 |
# Patch Header (Data Source Id) |
|
1238 |
#------------------------------ |
|
14 | 1239 |
|
1240 |
sysseek(WBPF,$start_ens,0) || die("$WBPcfn: $!"); |
|
18 | 1241 |
sysread(WBPF,$buf,1) || die("$WBPcfn: unexpected EOF"); |
1242 |
($hid) = unpack('C',$buf); |
|
1243 |
$hid == 0x7f || die(sprintf($FmtErr,$WBPcfn,"Header",$hid,$ens)); |
|
1244 |
||
31 | 1245 |
$buf = pack('C',$dta->{ENSEMBLE}[$ens]->{DATA_SOURCE_ID}); |
18 | 1246 |
my($nw) = syswrite(WBPF,$buf,1); |
1247 |
$nw == 1 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1248 |
||
1249 |
sysread(WBPF,$buf,4) == 4 || die("$WBPcfn: unexpected EOF"); |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1250 |
($el,$dummy,$ndt) = unpack('vCC',$buf); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1251 |
$ens--,next if (defined($ens_length) && ($el != $ens_length)); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1252 |
$ens_length = $el; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1253 |
|
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1254 |
printf(STDERR "$WBPcfn: WARNING: unexpected number of data types (%d, ens=$ens)\n",$ndt),last |
14 | 1255 |
unless ($ndt == 6 || $ndt == 7); |
18 | 1256 |
|
14 | 1257 |
sysread(WBPF,$buf,2*$ndt) == 2*$ndt || die("$WBPcfn: $!"); |
1258 |
@WBPofs = unpack("v$ndt",$buf); |
|
1259 |
$fixed_leader_bytes = $WBPofs[1] - $WBPofs[0]; |
|
31 | 1260 |
|
1261 |
#-------------------- |
|
1262 |
# Fixed Leader |
|
1263 |
#-------------------- |
|
14 | 1264 |
|
32 | 1265 |
sysseek(WBPF,$start_ens+$WBPofs[0]+4,0) # system config (transducer orientation) |
1266 |
|| die("$WBPcfn: $!"); |
|
1267 |
sysread(WBPF,$buf,1) == 1 || die("$WBPcfn: $!"); |
|
1268 |
$B1 = unpack('C',$buf); |
|
1269 |
$B1 |= 0x80 if ($dta->{ENSEMBLE}[$ens]->{XDUCER_FACING_UP}); |
|
1270 |
$B1 &= 0x7F if ($dta->{ENSEMBLE}[$ens]->{XDUCER_FACING_DOWN}); |
|
1271 |
$buf = pack('C',$B1); |
|
1272 |
sysseek(WBPF,$start_ens+$WBPofs[0]+4,0) |
|
1273 |
|| die("$WBPcfn: $!"); |
|
1274 |
syswrite(WBPF,$buf,1) == 1 || die("$WBPcfn: $!"); |
|
1275 |
||
1276 |
sysseek(WBPF,$start_ens+$WBPofs[0]+25,0) || die("$WBPcfn: $!"); # EX / coord-transformation |
|
31 | 1277 |
sysread(WBPF,$buf,1) == 1 || die("$WBPcfn: $!"); |
1278 |
my($EX) = unpack('C',$buf); |
|
1279 |
if ($dta->{BEAM_COORDINATES}) { |
|
1280 |
$EX &= ~0x18; |
|
1281 |
} elsif ($dta->{EARTH_COORDINATES}) { |
|
1282 |
$EX |= 0x18; |
|
1283 |
} else { |
|
1284 |
die("$WBPcfn: only beam- and earth coordinates are supported (implementation restriction)\n"); |
|
1285 |
} |
|
1286 |
$buf = pack('C',$EX); |
|
1287 |
sysseek(WBPF,$start_ens+$WBPofs[0]+25,0) || die("$WBPcfn: $!"); |
|
1288 |
syswrite(WBPF,$buf,1) == 1 || die("$WBPcfn: $!"); |
|
1289 |
||
42 | 1290 |
#---------------------------------------------------------------------- |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1291 |
# Variable Leader #0 |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1292 |
# - read ensNo for debugging purposes |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1293 |
#---------------------------------------------------------------------- |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1294 |
|
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1295 |
sysseek(WBPF,$start_ens+$WBPofs[1]+2,0) || die("$WBPcfn: $!"); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1296 |
sysread(WBPF,$buf,2) == 2 || die("$WBPcfn: $!"); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1297 |
my($ensNo) = unpack("v",$buf); # only lower two bytes!!! |
52
5b07a9b89aee
- adapted to work with IMPed moored ADCP data
Andreas Thurnherr <ant@ldeo.columbia.edu>
parents:
51
diff
changeset
|
1298 |
sysseek(WBPF,$start_ens+$WBPofs[1]+11,0) || die("$WBPcfn: $!"); # jump to high byte |
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1299 |
sysread(WBPF,$buf,1) == 1 || die("$WBPcfn: $!"); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1300 |
$ensNo += unpack('C',$buf) << 16; |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1301 |
die("ensNo = $ensNo (should be $dta->{ENSEMBLE}[$ens]->{NUMBER})\n") |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1302 |
unless ($ensNo == $dta->{ENSEMBLE}[$ens]->{NUMBER}); |
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1303 |
|
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1304 |
#---------------------------------------------------------------------- |
42 | 1305 |
# Variable Leader #1 |
1306 |
# - if $RDI_PD0_IO::OVERRIDE_Y2K_CLOCK is set, the data from the pre-Y2K |
|
1307 |
# clock are used to override the ADCP clock values; this allows |
|
1308 |
# a better time to be recorded by the data acquisition system |
|
1309 |
# without overwriting the main instrument clock data |
|
1310 |
#---------------------------------------------------------------------- |
|
31 | 1311 |
|
42 | 1312 |
if ($RDI_PD0_IO::OVERRIDE_Y2K_CLOCK) { |
1313 |
sysseek(WBPF,$start_ens+$WBPofs[1]+4,0) || die("$WBPcfn: $!"); # jump to RTC_YEAR |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1314 |
sysread(WBPF,$buf,7) == 7 || die("$WBPcfn: $!"); # read pre-Y2K clock |
42 | 1315 |
($dta->{ENSEMBLE}[$ens]->{YEAR}, |
1316 |
$dta->{ENSEMBLE}[$ens]->{MONTH}, |
|
1317 |
$dta->{ENSEMBLE}[$ens]->{DAY}, |
|
1318 |
$dta->{ENSEMBLE}[$ens]->{HOUR}, |
|
1319 |
$dta->{ENSEMBLE}[$ens]->{MINUTE}, |
|
50 | 1320 |
$dta->{ENSEMBLE}[$ens]->{SECOND},$B4) = |
42 | 1321 |
unpack('CCCCCCC',$buf); |
50 | 1322 |
$dta->{ENSEMBLE}[$ens]->{SECOND} += $B4/100; |
42 | 1323 |
$dta->{ENSEMBLE}[$ens]->{YEAR} += ($dta->{ENSEMBLE}[$ens]->{YEAR} > 80) ? 1900 : 2000; |
1324 |
} |
|
1325 |
||
1326 |
#---------------------------------------------------------------------- |
|
1327 |
# Variable Leader #2 |
|
43
b63fa355644c
commit to merge with changes from EN620
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
42
diff
changeset
|
1328 |
# - read ensemble number for debugging purposes |
42 | 1329 |
# - patch everything from SPEED_OF_SOUND to TEMPERATURE |
1330 |
# - at one stage, IMP allowed for missing values in pitch/roll and heading; |
|
1331 |
# on 12/23/2017 the corresponding code was disabled (replaced by assertion) |
|
1332 |
#---------------------------------------------------------------------- |
|
1333 |
||
1334 |
sysseek(WBPF,$start_ens+$WBPofs[1]+14,0) || die("$WBPcfn: $!"); # jump to SPEED_OF_SOUND |
|
14 | 1335 |
|
37 | 1336 |
$dta->{ENSEMBLE}[$ens]->{XDUCER_DEPTH} = _round($dta->{ENSEMBLE}[$ens]->{XDUCER_DEPTH}*10); |
18 | 1337 |
|
42 | 1338 |
# $dta->{ENSEMBLE}[$ens]->{HEADING} = defined($dta->{ENSEMBLE}[$ens]->{HEADING}) |
1339 |
# ? _round($dta->{ENSEMBLE}[$ens]->{HEADING}*100) |
|
1340 |
# : 0xF000; |
|
1341 |
# $dta->{ENSEMBLE}[$ens]->{PITCH} = defined($dta->{ENSEMBLE}[$ens]->{PITCH}) |
|
1342 |
# ? unpack('S',pack('s',_round($dta->{ENSEMBLE}[$ens]->{PITCH}*100))) |
|
1343 |
# : 0x8000; |
|
1344 |
# $dta->{ENSEMBLE}[$ens]->{ROLL} = defined($dta->{ENSEMBLE}[$ens]->{ROLL}) |
|
1345 |
# ? unpack('S',pack('s',_round($dta->{ENSEMBLE}[$ens]->{ROLL}*100))) |
|
1346 |
# : 0x8000; |
|
28 | 1347 |
|
42 | 1348 |
croak("$0: assertion failed") unless defined($dta->{ENSEMBLE}[$ens]->{HEADING}) && |
1349 |
defined($dta->{ENSEMBLE}[$ens]->{PITCH}) && |
|
1350 |
defined($dta->{ENSEMBLE}[$ens]->{ROLL}); |
|
1351 |
$dta->{ENSEMBLE}[$ens]->{HEADING} = _round($dta->{ENSEMBLE}[$ens]->{HEADING}*100); |
|
1352 |
$dta->{ENSEMBLE}[$ens]->{PITCH} = unpack('S',pack('s',_round($dta->{ENSEMBLE}[$ens]->{PITCH}*100))); |
|
1353 |
$dta->{ENSEMBLE}[$ens]->{ROLL} = unpack('S',pack('s',_round($dta->{ENSEMBLE}[$ens]->{ROLL}*100))); |
|
1354 |
$dta->{ENSEMBLE}[$ens]->{TEMPERATURE} = unpack('S',pack('s',_round($dta->{ENSEMBLE}[$ens]->{TEMPERATURE}*100))); |
|
14 | 1355 |
|
1356 |
$buf = pack('vvvvvvv', |
|
31 | 1357 |
$dta->{ENSEMBLE}[$ens]->{SPEED_OF_SOUND}, |
1358 |
$dta->{ENSEMBLE}[$ens]->{XDUCER_DEPTH},$dta->{ENSEMBLE}[$ens]->{HEADING}, |
|
1359 |
$dta->{ENSEMBLE}[$ens]->{PITCH},$dta->{ENSEMBLE}[$ens]->{ROLL}, |
|
1360 |
$dta->{ENSEMBLE}[$ens]->{SALINITY},$dta->{ENSEMBLE}[$ens]->{TEMPERATURE}); |
|
14 | 1361 |
|
1362 |
my($nw) = syswrite(WBPF,$buf,14); |
|
1363 |
$nw == 14 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1364 |
||
28 | 1365 |
|
42 | 1366 |
#---------------------------------------------------------------------- |
1367 |
# Variable Leader #3 |
|
1368 |
# - patch Y2K RTC |
|
1369 |
#---------------------------------------------------------------------- |
|
1370 |
||
1371 |
sysseek(WBPF,$start_ens+$WBPofs[1]+57,0) || die("$WBPcfn: $!"); # jump to RTC_CENTURY |
|
1372 |
||
1373 |
my($century) = int($dta->{ENSEMBLE}[$ens]->{YEAR} / 100); |
|
1374 |
my($year) = $dta->{ENSEMBLE}[$ens]->{YEAR} % 100; |
|
50 | 1375 |
my($seconds) = int($dta->{ENSEMBLE}[$ens]->{SECOND}); |
1376 |
my($hundredths) = 100 * ($dta->{ENSEMBLE}[$ens]->{SECOND} - $seconds); |
|
42 | 1377 |
$buf = pack('CCCCCCCC',$century,$year,$dta->{ENSEMBLE}[$ens]->{MONTH}, |
1378 |
$dta->{ENSEMBLE}[$ens]->{DAY},$dta->{ENSEMBLE}[$ens]->{HOUR}, |
|
1379 |
$dta->{ENSEMBLE}[$ens]->{MINUTE},$seconds,$hundredths); |
|
1380 |
my($nw) = syswrite(WBPF,$buf,8); |
|
1381 |
$nw == 8 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1382 |
||
28 | 1383 |
#-------------------- |
1384 |
# Velocity Data |
|
1385 |
#-------------------- |
|
1386 |
||
32 | 1387 |
sysseek(WBPF,$start_ens+$WBPofs[2]+2,0) || die("$WBPcfn: $!"); # skip velocity data id (assume it is correct) |
28 | 1388 |
for ($bin=0; $bin<$nbins; $bin++) { |
1389 |
for ($beam=0; $beam<4; $beam++) { |
|
31 | 1390 |
$dta->{ENSEMBLE}[$ens]->{VELOCITY}[$bin][$beam] = defined($dta->{ENSEMBLE}[$ens]->{VELOCITY}[$bin][$beam]) |
37 | 1391 |
? _round($dta->{ENSEMBLE}[$ens]->{VELOCITY}[$bin][$beam]*1000) |
28 | 1392 |
: 0x8000; |
31 | 1393 |
$buf = pack('v',unpack('S',pack('s',$dta->{ENSEMBLE}[$ens]->{VELOCITY}[$bin][$beam]))); |
28 | 1394 |
my($nw) = syswrite(WBPF,$buf,2); |
1395 |
$nw == 2 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1396 |
} |
|
1397 |
} |
|
1398 |
||
1399 |
#-------------------- |
|
1400 |
# Correlation Data |
|
1401 |
#-------------------- |
|
1402 |
||
32 | 1403 |
sysseek(WBPF,$start_ens+$WBPofs[3]+2,0) || die("$WBPcfn: $!"); |
28 | 1404 |
for ($bin=0; $bin<$nbins; $bin++) { |
1405 |
for ($beam=0; $beam<4; $beam++) { |
|
31 | 1406 |
$buf = pack('C',$dta->{ENSEMBLE}[$ens]->{CORRELATION}[$bin][$beam]); |
28 | 1407 |
my($nw) = syswrite(WBPF,$buf,1); |
1408 |
$nw == 1 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1409 |
} |
|
1410 |
} |
|
1411 |
||
1412 |
#-------------------- |
|
1413 |
# Echo Intensity Data |
|
1414 |
#-------------------- |
|
1415 |
||
32 | 1416 |
sysseek(WBPF,$start_ens+$WBPofs[4]+2,0) || die("$WBPcfn: $!"); |
28 | 1417 |
|
1418 |
for ($bin=0; $bin<$nbins; $bin++) { |
|
1419 |
for ($beam=0; $beam<4; $beam++) { |
|
31 | 1420 |
$buf = pack('C',$dta->{ENSEMBLE}[$ens]->{ECHO_AMPLITUDE}[$bin][$beam]); |
28 | 1421 |
my($nw) = syswrite(WBPF,$buf,1); |
1422 |
$nw == 1 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1423 |
} |
|
1424 |
} |
|
1425 |
||
1426 |
#-------------------- |
|
1427 |
# Percent Good Data |
|
1428 |
#-------------------- |
|
1429 |
||
32 | 1430 |
sysseek(WBPF,$start_ens+$WBPofs[5]+2,0) || die("$WBPcfn: $!"); |
28 | 1431 |
|
1432 |
for ($i=0,$bin=0; $bin<$nbins; $bin++) { |
|
1433 |
for ($beam=0; $beam<4; $beam++,$i++) { |
|
31 | 1434 |
$buf = pack('C',$dta->{ENSEMBLE}[$ens]->{PERCENT_GOOD}[$bin][$beam]); |
28 | 1435 |
my($nw) = syswrite(WBPF,$buf,1); |
1436 |
$nw == 1 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1437 |
} |
|
1438 |
} |
|
1439 |
||
1440 |
#----------------------------------------- |
|
1441 |
# Bottom-Track Data |
|
1442 |
# - scan through remaining data types |
|
1443 |
#----------------------------------------- |
|
1444 |
||
1445 |
my($nxt); |
|
1446 |
for ($nxt=6; $nxt<$ndt; $nxt++) { # scan until BT found |
|
32 | 1447 |
sysseek(WBPF,$start_ens+$WBPofs[$nxt],0) || die("$WBPcfn: $!"); |
36
515b06dae59c
version at end of ECOGIG EN586 cruise
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
32
diff
changeset
|
1448 |
sysread(WBPF,$buf,2) == 2 || die("$WBPcfn: $! [ens=$ens]"); |
28 | 1449 |
$id = unpack('v',$buf); |
1450 |
last if ($id == 0x0600); |
|
1451 |
} |
|
1452 |
||
1453 |
unless ($nxt == $ndt) { # BT found |
|
32 | 1454 |
sysseek(WBPF,14,1) || die("$WBPcfn: $!"); # skip BT config |
1455 |
for ($beam=0; $beam<4; $beam++) { # BT range low bytes (2 per beam) |
|
37 | 1456 |
$buf = pack('v',_round($dta->{ENSEMBLE}[$ens]->{BT_RANGE}[$beam] * 100) & 0xFFFF); |
32 | 1457 |
my($nw) = syswrite(WBPF,$buf,2); |
1458 |
$nw == 2 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1459 |
} |
|
1460 |
||
1461 |
for ($beam=0; $beam<4; $beam++) { # BT velocities |
|
1462 |
$buf = pack('v',unpack('S',pack('s', |
|
1463 |
defined($dta->{ENSEMBLE}[$ens]->{BT_VELOCITY}[$beam]) |
|
37 | 1464 |
? _round($dta->{ENSEMBLE}[$ens]->{BT_VELOCITY}[$beam]*1000) |
32 | 1465 |
: 0x8000))); |
1466 |
my($nw) = syswrite(WBPF,$buf,2); |
|
1467 |
$nw == 2 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1468 |
} |
|
1469 |
||
1470 |
for ($beam=0; $beam<4; $beam++) { # BT correlation |
|
1471 |
$buf = pack('C',$dta->{ENSEMBLE}[$ens]->{BT_CORRELATION}[$beam]); |
|
1472 |
my($nw) = syswrite(WBPF,$buf,1); |
|
1473 |
$nw == 1 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1474 |
} |
|
1475 |
||
1476 |
for ($beam=0; $beam<4; $beam++) { # BT evaluation amp of matching filter |
|
1477 |
$buf = pack('C',$dta->{ENSEMBLE}[$ens]->{BT_EVAL_AMPLITUDE}[$beam]); |
|
1478 |
my($nw) = syswrite(WBPF,$buf,1); |
|
1479 |
$nw == 1 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1480 |
} |
|
1481 |
||
1482 |
for ($beam=0; $beam<4; $beam++) { # BT percent good |
|
1483 |
$buf = pack('C',$dta->{ENSEMBLE}[$ens]->{BT_PERCENT_GOOD}[$beam]); |
|
1484 |
my($nw) = syswrite(WBPF,$buf,1); |
|
1485 |
$nw == 1 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1486 |
} |
|
28 | 1487 |
|
32 | 1488 |
sysseek(WBPF,33,1) || die("$WBPcfn: $!"); # BT range high bytes (1 per beam) |
1489 |
for ($beam=0; $beam<4; $beam++) { |
|
37 | 1490 |
$buf = pack('C',(_round($dta->{ENSEMBLE}[$ens]->{BT_RANGE}[$beam]*100) & 0xFF0000) >> 16); |
32 | 1491 |
my($nw) = syswrite(WBPF,$buf,1); |
1492 |
$nw == 1 || die("$WBPcfn: $nw bytes written ($!)"); |
|
1493 |
} |
|
1494 |
||
1495 |
} # if BT found |
|
28 | 1496 |
|
14 | 1497 |
#---------------- |
1498 |
# Update Checksum |
|
1499 |
#---------------- |
|
1500 |
||
1501 |
sysseek(WBPF,$start_ens,0) || die("$WBPcfn: $!"); |
|
1502 |
sysread(WBPF,$buf,$ens_length) == $ens_length || die("$WBPcfn: $!"); |
|
1503 |
$cs = unpack('%16C*',$buf); |
|
1504 |
$buf = pack('v',$cs); |
|
1505 |
$nw = syswrite(WBPF,$buf,2); |
|
1506 |
$nw == 2 || die("$WBPcfn: $nw bytes written, ens=$ens ($!)"); |
|
1507 |
||
1508 |
} # ens loop |
|
1509 |
} |
|
1510 |
||
37 | 1511 |
#---------------------------------------------------------------------- |
1512 |
# &clearEns(^data,ens-index) |
|
1513 |
# - undefine all velocities in ensemble, including BT |
|
1514 |
# - this is required for the LDEO_IX software, |
|
1515 |
# which does not recognize missing attitude values |
|
1516 |
# - set percent good to zero |
|
1517 |
# - this is done for consistency |
|
1518 |
# - DO NOT undef heading, pitch and roll |
|
1519 |
# - the LDEO software does not recognize missing attitude vals |
|
1520 |
# and, therefore, misinterprets those |
|
1521 |
# - while this should not matter because all the velocities are |
|
1522 |
# also deleted, it was found that setting only the heading to |
|
1523 |
# undef'd and leaving pitch and roll unchanged causes |
|
1524 |
# significant errors in GPS velocity referencing! This |
|
1525 |
# must be a bug |
|
1526 |
# - also, if attitudes are undef'd the LDEO software |
|
1527 |
# cannto determine the instrument offset from pitch/roll |
|
1528 |
# and the pitch/roll DL vs UL plots are bogus |
|
1529 |
#---------------------------------------------------------------------- |
|
1530 |
||
1531 |
sub clearEns($$) |
|
1532 |
{ |
|
1533 |
my($dta,$ens) = @_; |
|
39 | 1534 |
die("clearEns: ens-index $ens out of range\n") |
37 | 1535 |
unless ($ens>=0 && $ens<=$#{$dta->{ENSEMBLE}}); |
1536 |
for (my($bin)=0; $bin<$dta->{N_BINS}; $bin++) { |
|
1537 |
undef(@{$dta->{ENSEMBLE}[$ens]->{VELOCITY}[$bin]}); |
|
1538 |
@{$dta->{ENSEMBLE}[$ens]->{PERCENT_GOOD}[$bin]} = (0,0,0,0); |
|
1539 |
} |
|
1540 |
undef(@{$dta->{ENSEMBLE}[$ens]->{BT_VELOCITY}}); |
|
1541 |
@{$dta->{ENSEMBLE}[$ens]->{BT_PERCENT_GOOD}} = (0,0,0,0); |
|
1542 |
# undef($dta->{ENSEMBLE}[$ens]->{HEADING}); |
|
1543 |
# undef($dta->{ENSEMBLE}[$ens]->{PITCH}); |
|
1544 |
# undef($dta->{ENSEMBLE}[$ens]->{ROLL}); |
|
1545 |
} |
|
1546 |
||
14 | 1547 |
1; # return true for all the world to see |