Revision Log

nortek.py

Revision 168 (checked in by cbc, 16 years ago)	Adding diwasp customizations.

Line
1	"""This code takes the raw data from the Nortek AWAC and prepares it
2	for matlab You input a raw binary *.WPR file from the AWAC and it outputs three Currents files with all of the currents
3	data, as well as a waves orbital velocity, range, pressure and sysinfo file for each wave burst"""
4
5
6	import sys, struct, math, shutil, os
7
8	try:
9	infileName = sys.argv[1]
10
11	except:
12	print 'error'
13	sys.exit(1)
14
15	ifile = open(infileName, 'rb')
16	testOfile=open('test','w')
17	currentsOfile = open('currents', 'w')
18	currentsHeaderOfile = open('currentsHeader','w')
19	currentsTiltOfile = open('currentsTilt','w')
20
21	currentsOfile.write('%this is the Nortek AWAC Currents data \n')
22	currentsHeaderOfile.write('%this is the Header data for the Nortek AWAC Currents data \n %it includes Header Info, Date info, Pressure and Temp \n')
23	currentsTiltOfile.write('%this is the Tilt readings (heading pitch and roll) for the Nortek AWAC Currents data \n')
24
25	# HARDWARE CONFIGURATION
26
27	#read header and make sure it is valid
28	HWconfID= 0x05a5
29
30	readHeaderID=ifile.read(2)
31	headerID=struct.unpack('H',readHeaderID)[0]
32	if headerID != HWconfID:
33	print 'error with Hardware Config Header'
34	sys.exit(1)
35
36	#size of structure in words (1 word = 2 bytes)
37	readSize=ifile.read(2)
38	size=struct.unpack('H',readSize)[0]
39	#compute the length in bytes - the checksum
40	length=(size*2)-2
41
42	#Serial Number (not sure what to do with this)
43	readSerial=ifile.read(14)
44
45	#Board Configuration
46	readBoardC=ifile.read(2)
47
48	#Board Frequency (kHz)
49	readBoardF=ifile.read(2)
50
51	#PIC Code Version #
52	readPIC=ifile.read(2)
53
54	#Hardware Revision
55	readHardRev=ifile.read(2)
56
57	#Recorder Size (*65536 bytes)
58	readRecSize=ifile.read(2)
59
60	#Status
61	readStatus=ifile.read(2)
62
63	#Spare
64	ifile.read(12)
65
66	#Firmware Version
67	readFirmVersion=ifile.read(4)
68	firmVersion=struct.unpack('I',readFirmVersion)[0]
69	currentsHeaderOfile.write(str(firmVersion))
70	currentsHeaderOfile.write(' %Firmware Version \n')
71
72	#Checksum
73	readCheckSum=ifile.read(2)
74
75
76
77	# HEAD CONFIGURATION
78
79	#read header and make sure it is valid
80	HeadConfID= 0x04a5
81
82	readHeaderID=ifile.read(2)
83	headerID=struct.unpack('H',readHeaderID)[0]
84	if headerID != HeadConfID:
85	print 'error with Head Config Header'
86	sys.exit(1)
87
88	#size of structure in words (1 word = 2 bytes)
89	readSize=ifile.read(2)
90	size=struct.unpack('H',readSize)[0]
91	#compute the length in bytes - the checksum
92	length=(size*2)-2
93
94	#Head Configuration
95	readHeadConfig=ifile.read(2)
96
97	#Head Frequency (kHz)
98	readHeadFreq=ifile.read(2)
99	headFreq=struct.unpack('H',readHeadFreq)[0]
100	currentsHeaderOfile.write(str(headFreq))
101	currentsHeaderOfile.write(' %head frequency in khz \n')
102
103	#Head Type
104	readHeadType=ifile.read(2)
105
106	#Head Serial Number
107	readHeadSerial=ifile.read(12)
108
109	#System Data
110	ifile.read(176)
111
112	#Spare
113	ifile.read(22)
114
115	#Number of beams
116	readNbeams=ifile.read(2)
117	Nbeams=struct.unpack('H',readNbeams)[0]
118	currentsHeaderOfile.write(str(Nbeams))
119	currentsHeaderOfile.write(' %number of beams \n')
120
121	#Checksum
122	readChecksum=ifile.read(2)
123
124
125
126	# USER CONFIGURATION
127
128	#read header and make sure it is valid
129	UserConfID= 0x00a5
130
131	readHeaderID=ifile.read(2)
132	headerID=struct.unpack('H',readHeaderID)[0]
133	if headerID != UserConfID:
134	print 'error with User Config Header'
135	sys.exit(1)
136
137	#size of structure in words (1 word = 2 bytes)
138	readSize=ifile.read(2)
139	size=struct.unpack('H',readSize)[0]
140	#compute the length in bytes - the checksum
141	length=(size*2)-2
142
143	#Transmit Pulse Length (counts)
144	readPulse=ifile.read(2)
145	pulse=struct.unpack('H',readPulse)[0]
146	currentsHeaderOfile.write(str(pulse))
147	currentsHeaderOfile.write(' %transmit pulse length (counts) \n')
148
149	#Blanking Distance (counts)
150	readBlanking=ifile.read(2)
151	blanking=struct.unpack('H',readBlanking)[0]
152
153	#Receive Length (counts)
154	readReceive=ifile.read(2)
155	receive=struct.unpack('H',readReceive)[0]
156	currentsHeaderOfile.write(str(receive))
157	currentsHeaderOfile.write(' %recieve length (counts) \n')
158
159	#Time Between Pings (counts)
160	readTimePing=ifile.read(2)
161	timePing=struct.unpack('H',readTimePing)[0]
162	currentsHeaderOfile.write(str(timePing))
163	currentsHeaderOfile.write(' %time between pings (counts) \n')
164
165	#Time Between Burst Sequence (counts)
166	readTimeBurst=ifile.read(2)
167	timeBurst=struct.unpack('H',readTimeBurst)[0]
168	currentsHeaderOfile.write(str(timeBurst))
169	currentsHeaderOfile.write(' %time between burst sequence (counts) \n')
170
171	#Number of Beam Sequences per Burst
172	readBeamSeq=ifile.read(2)
173	beamSeq=struct.unpack('H',readBeamSeq)[0]
174	currentsHeaderOfile.write(str(beamSeq))
175	currentsHeaderOfile.write(' %number of beam sequences per burst \n')
176
177	#Average interval (seconds)
178	readAvgInt=ifile.read(2)
179	avgInt=struct.unpack('H',readAvgInt)[0]
180	currentsHeaderOfile.write(str(avgInt))
181	currentsHeaderOfile.write(' %average interval (seconds) \n')
182
183
184	#Number of Beams (alread have this from before)
185	readNbeams=ifile.read(2)
186
187	#Timing Controller Mode
188	readTimeCont=ifile.read(2)
189
190	#Power Control Register
191	readPowerCont=ifile.read(2)
192
193	#Spare
194	ifile.read(6)
195
196	#Compass Update Rate
197	readCompRate=ifile.read(2)
198
199	#Coordinate System (0=ENU, 1=XYZ, 2=Beam)
200	readCoord=ifile.read(2)
201	coord=struct.unpack('H',readCoord)[0]
202	currentsHeaderOfile.write(str(coord))
203	currentsHeaderOfile.write(' %coord system, 0=ENU, 1=XYZ, 2=Beam \n')
204
205	#Number of Cells or Bins
206	readNcells=ifile.read(2)
207	Ncells=struct.unpack('H',readNcells)[0]
208	currentsHeaderOfile.write(str(Ncells))
209	currentsHeaderOfile.write(' %number of bins or cells \n')
210
211
212	#Cell Size
213	readCellSize=ifile.read(2)
214	cCellSize=struct.unpack('H',readCellSize)[0]
215
216	#compute the currents cell size
217	cell=cCellSize100.0(750/((headFreq1000.04)/255.0))
218	cellint = math.floor(cell)
219	rawCellSize=0.01*cellint/256.0
220	cellSize=math.cos(25(math.pi/180))rawCellSize
221	currentsHeaderOfile.write(str('%3.2f' %cellSize))
222	currentsHeaderOfile.write(' %bin cell size (m) \n')
223
224	#compute the blanking distance
225	pLength=rawCellSize*100.0
226	cellBlanking=blanking100.0(750.0/32768.0)-pLength
227	cellBlanking=0.01*math.floor(cellBlanking)
228	cellBlanking=math.cos(25(math.pi/180))cellBlanking
229	currentsHeaderOfile.write(str('%3.2f' %cellBlanking))
230	currentsHeaderOfile.write(' %blanking distance (m) \n')
231
232	#compute and write the cell depths
233	currentsHeaderOfile.write(' % cell heights from bottom (meters) \n')
234	depths=[]
235	for i in range(1,Ncells+1):
236	depths.append(cellBlanking+cellSize*i)
237	currentsHeaderOfile.write(str('%3.2f' %depths[i-1]))
238	currentsHeaderOfile.write('\n')
239
240	#Measurement Interval
241	readMeasInt=ifile.read(2)
242	measInt=struct.unpack('H',readMeasInt)[0]
243	currentsHeaderOfile.write(str(measInt))
244	currentsHeaderOfile.write(' %measurement interval (seconds) \n')
245
246	#Recorder Deployment Name
247	readDepName=ifile.read(6)
248
249	#Recorder Wrap Mode
250	readRecWrap=ifile.read(2)
251
252	#Deployment Start Time (we can use the start time of the samples instead)
253	readStarTime=ifile.read(6)
254
255	#Time Between Diagnostic Measurements (seconds)
256	readDiagTime=ifile.read(4)
257	diagTime=struct.unpack('I',readDiagTime)[0]
258
259	#Mode
260	readMode=ifile.read(2)
261
262	#User Input Sound Speed Adjustment Factor
263	readAdjSound=ifile.read(2)
264
265	# Number of Samples in Diagnostic Mode
266	readNsamples=ifile.read(2)
267	Nsamples=struct.unpack('H',readNsamples)[0]
268
269	# Beam and Cell Number to Measure in Diagnostic Mode
270	readNbeamCell=ifile.read(2)
271
272	#Number Pings in Diagnostic Mode
273	readNpingsDiag=ifile.read(2)
274
275	#Mode Test
276	readModeTest=ifile.read(2)
277
278	#Analog Input Address
279	readAnalAddr=ifile.read(2)
280
281	#Software Version
282	readSoftVers=ifile.read(2)
283
284	#Spare
285	ifile.read(2)
286
287	#Velocity Adjustment Table
288	readVelAdjTable=ifile.read(180)
289
290	#File Comments
291	readFileComments=ifile.read(180)
292
293	#Mode for data rate and wave cell position
294	readMode=ifile.read(2)
295	byteList=struct.unpack('2B', readMode)
296
297	def int2msbits(n, count=8):
298	"""returns the most-significant binary of integer n, using count number of digits"""
299	return "".join([str((n >> y) & 1) for y in range(0, count)])
300
301	bytemap = ''
302	bitmap = ''
303	for ind in range(len(byteList)):
304	byte = struct.pack('B', byteList[ind])
305	bytemap += byte
306	bits = int2msbits(byteList[ind])
307	bitmap += bits
308
309	#Percentage for Wave Cell Positioning
310	readPercWave=ifile.read(2)
311	PercWave=struct.unpack('H',readPercWave)[0]
312	PercWave=(PercWave/32767)*100
313
314	#Wave Transmit Pulse
315	readWavePulse=ifile.read(2)
316	wavePulse=struct.unpack('H',readWavePulse)[0]
317
318	#Fixed Wave Blanking Distance (counts)
319	readFixedBlanking=ifile.read(2)
320	fixedBlanking=struct.unpack('H',readFixedBlanking)[0]
321
322	#Wave Measurement Cell Size
323	readWaveCellSize=ifile.read(2)
324	waveCellSize=struct.unpack('H',readWaveCellSize)[0]
325
326	#Number of Diagnostic/Wave Samples
327	readWaveSamples=ifile.read(2)
328	waveSamples=struct.unpack('H',readWaveSamples)[0]
329
330	#Spare
331	ifile.read(8)
332
333	#Analog Output Scale Factor
334	readAnalogScale=ifile.read(2)
335
336	#Correlation Threshold
337	readCorrThresh=ifile.read(2)
338
339	#Spare
340	ifile.read(2)
341
342	#Transmit Pulse Length Second Lag (counts)
343	readPulseLag=ifile.read(2)
344
345	#Spare
346	ifile.read(30)
347
348	#Stage Match Filter Constants
349	readStageMatch=ifile.read(16)
350
351	#Checksum
352	readChecksum=ifile.read(2)
353
354
355
356	# VELOCITY PROFILE DATA
357
358	#first create the lists that we can write the data to and then write to the file
359	currentsTime=[]
360	currentsPress=[]
361	currentsTemp=[]
362	currentsHeading=[]
363	currentsPitch=[]
364	currentsRoll=[]
365	allVel1=[]
366	allVel2=[]
367	allVel3=[]
368	allAmp1=[]
369	allAmp2=[]
370	allAmp3=[]
371
372	VelID= 0x20a5
373
374	readHeaderID=ifile.read(2)
375	headerID=struct.unpack('H',readHeaderID)[0]
376	if headerID != VelID:
377	print 'error with Velocity Profile Header'
378	sys.exit(1)
379
380	while headerID == VelID:
381
382
383	#size in words
384	readSize=ifile.read(2)
385	if len(readSize) == 0:
386	break
387	size=struct.unpack('H',readSize)[0]
388	#compute the length in bytes - the checksum
389	length=(size*2)-2
390
391	# READ AND COMPUTE THE START TIME FOR THE CURRENT VEL>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
392	#minute
393	readMinute=ifile.read(1)
394	minute=struct.unpack('B',readMinute)[0]
395	minute=str(hex(minute))[2:4]
396	#Second
397	readSecond=ifile.read(1)
398	second=struct.unpack('B',readSecond)[0]
399	second=str(hex(second))[2:4]
400	#day
401	readDay=ifile.read(1)
402	day=struct.unpack('B',readDay)[0]
403	day=str(hex(day))[2:4]
404	#hour
405	readHour=ifile.read(1)
406	hour=struct.unpack('B',readHour)[0]
407	hour=str(hex(hour))[2:4]
408	#year
409	readYear=ifile.read(1)
410	year=struct.unpack('B',readYear)[0]
411	year=str(hex(year))[2:4]
412	#month
413	readMonth=ifile.read(1)
414	month=struct.unpack('B',readMonth)[0]
415	month=str(hex(month))[2:4]
416	#century
417	century = str(20)
418
419	startTime=century+(year.zfill(2))+(month.zfill(2))+(day.zfill(2))+(hour.zfill(2))+(minute.zfill(2))+(second.zfill(2))
420
421	#create a time stamp for the file names that only includes from year up to minutes
422	FileStamp=(year.zfill(2))+(month.zfill(2))+(day.zfill(2))+(hour.zfill(2))+(minute.zfill(2))
423	currentsTime.append(FileStamp)
424
425
426	#error code
427	readError=ifile.read(2)
428	error=struct.unpack('H',readError)[0]
429
430	#analog input 1
431	readAnalog1=ifile.read(2)
432	analog1=struct.unpack('H',readAnalog1)[0]
433
434	#Battery Voltage (.1 V)
435	readBattery=ifile.read(2)
436	battery=struct.unpack('H',readBattery)[0]
437
438	#Speed of Sound or Analog 2 (.1 m/s)
439	readSound=ifile.read(2)
440	sound=struct.unpack('H',readSound)[0]
441
442	#Heading (.1 deg)
443	readHeading=ifile.read(2)
444	heading=struct.unpack('H',readHeading)[0]
445	currentsHeading.append(heading)
446
447	#Pitch (.1 deg)
448	readPitch=ifile.read(2)
449	pitch=struct.unpack('h',readPitch)[0]
450	currentsPitch.append(pitch)
451
452	#Roll (.1 deg)
453	readRoll=ifile.read(2)
454	roll=struct.unpack('h',readRoll)[0]
455	currentsRoll.append(roll)
456
457	#pressure MSB in mm
458	readPressMSB=ifile.read(1)
459	pressureMSB=struct.unpack('B',readPressMSB)[0]
460
461	#status code
462	readStatus=ifile.read(1)
463
464	#pressure LSW in mm
465	readPressLSW=ifile.read(2)
466	pressureLSW=struct.unpack('H',readPressLSW)[0]
467
468	#compute the Pressure
469	pressure= 65536*pressureMSB + pressureLSW
470	currentsPress.append(pressure)
471
472	#Temperature (0.01 deg C)
473	readTemp=ifile.read(2)
474	temp=struct.unpack('h',readTemp)[0]
475	currentsTemp.append(temp)
476
477	#Spare
478	readSpare=ifile.read(88)
479
480
481	#Read in the velocity data for beam 1-3 (mm/s) from cell 1(bottom) to n
482	numCells=3*Ncells
483	vel1=[]
484	vel2=[]
485	vel3=[]
486	while numCells > 0:
487	readVel=ifile.read(2)
488	vel=struct.unpack('h',readVel)[0]
489	if numCells > 2*Ncells:
490	vel1.append(vel)
491	elif numCells > Ncells:
492	vel2.append(vel)
493	else:
494	vel3.append(vel)
495	numCells -=1
496	allVel1.append(vel1)
497	allVel2.append(vel2)
498	allVel3.append(vel3)
499
500	#Read in the amplitude data for beam 1-3 (counts) from cell 1(bottom) to n
501
502	numCells=3*Ncells
503	amp1=[]
504	amp2=[]
505	amp3=[]
506	while numCells > 0:
507	readAmp=ifile.read(1)
508	amp=struct.unpack('B',readAmp)[0]
509	if numCells > 2*Ncells:
510	amp1.append(amp)
511	elif numCells > Ncells:
512	amp2.append(amp)
513	else:
514	amp3.append(amp)
515	numCells -=1
516	allAmp1.append(amp1)
517	allAmp2.append(amp2)
518	allAmp3.append(amp3)
519
520	#Add a fill byte if there are an odd number of cells
521	if Ncells/2 != 0:
522	ifile.read(1)
523
524	else:
525	ifile.read(0)
526
527	#Compute the Checksum
528	ifile.seek(-(length+1),1)
529	dataStruct=ifile.read(length)
530
531	total =0
532	for byte in dataStruct:
533	value = struct.unpack('B', byte)[0]
534	total += value
535	#add 0xb58c in hex to put the cs in the right format
536	cs=total
537
538	#Read in the Checksum
539	readChecksum=ifile.read(2)
540	checksum=struct.unpack('H',readChecksum)[0]
541
542
543
544
545	# WAVE DATA HEADER
546
547	WdataID=0x31a5
548
549	#Read the Sync and ID as a 2 byte hex instead of seperately
550	readHeaderID=ifile.read(2)
551	if len(readHeaderID) == 0:
552	break
553	headerID=struct.unpack('H',readHeaderID)[0]
554
555
556
557	while headerID == WdataID:
558
559	if headerID != WdataID:
560	break
561
562	pressureOfile = open('pressure', 'w')
563	rangeOfile = open('range', 'w')
564	orbitOfile = open('orbit', 'w')
565	sysinfoOfile = open('sysinfo', 'w')
566
567
568	#write the data from the hardware/user config to the sysinfoOfile
569	sysinfoOfile.write(str(coord))
570	sysinfoOfile.write(' %coord system, 0=ENU, 1=XYZ, 2=Beam \n')
571
572	sysinfoOfile.write(str(bitmap))
573	sysinfoOfile.write(' % bit1:data rate 0=1Hz,1=2Hz bit2:wave cell position 0=fixed,1=dynamic bit3=type of pos 0=% of mean press, 1=% of min press \n')
574
575	sysinfoOfile.write(str(PercWave))
576	sysinfoOfile.write(' %percentage for wave cell position \n')
577
578
579	#size in words
580	readSize=ifile.read(2)
581	size=struct.unpack('H',readSize)[0]
582	#compute the length in bytes - the checksum
583	length=(size*2)-2
584
585	# READ AND COMPUTE THE START TIME FOR THE WAVE HEADER>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
586	#minute
587	readMinute=ifile.read(1)
588	minute=struct.unpack('B',readMinute)[0]
589	minute=str(hex(minute))[2:4]
590	#Second
591	readSecond=ifile.read(1)
592	second=struct.unpack('B',readSecond)[0]
593	second=str(hex(second))[2:4]
594	#day
595	readDay=ifile.read(1)
596	day=struct.unpack('B',readDay)[0]
597	day=str(hex(day))[2:4]
598	#hour
599	readHour=ifile.read(1)
600	hour=struct.unpack('B',readHour)[0]
601	hour=str(hex(hour))[2:4]
602	#year
603	readYear=ifile.read(1)
604	year=struct.unpack('B',readYear)[0]
605	year=str(hex(year))[2:4]
606	#month
607	readMonth=ifile.read(1)
608	month=struct.unpack('B',readMonth)[0]
609	month=str(hex(month))[2:4]
610	#century
611	century = str(20)
612
613	startTime=century+(year.zfill(2))+(month.zfill(2))+(day.zfill(2))+(hour.zfill(2))+(minute.zfill(2))+(second.zfill(2))
614
615	#create a time stamp for the file names that only includes from year up to minutes
616	FileStamp=(year.zfill(2))+(month.zfill(2))+(day.zfill(2))+(hour.zfill(2))+(minute.zfill(2))
617
618	sysinfoOfile.write(startTime)
619	sysinfoOfile.write(' % this is the start time \n')
620
621
622
623	#Number of Wave Data records to follow
624	readWaveRecords=ifile.read(2)
625	waveRecords=struct.unpack('H',readWaveRecords)[0]
626
627	#Blanking Distance (counts)
628	readBlanking=ifile.read(2)
629	blanking=struct.unpack('H',readBlanking)[0]
630
631	#Battery Voltage
632	readBattery=ifile.read(2)
633
634	#Sound Speed
635	readSound=ifile.read(2)
636	sound=struct.unpack('H',readSound)[0]
637	sysinfoOfile.write(str(sound))
638	sysinfoOfile.write(' %speed of sound (0.1m/s) \n')
639
640	#Heading
641	readHeading=ifile.read(2)
642	heading=struct.unpack('H',readHeading)[0]
643	sysinfoOfile.write(str(heading))
644	sysinfoOfile.write(' %heading \n')
645
646	#Pitch
647	readPitch=ifile.read(2)
648	pitch=struct.unpack('h',readPitch)[0]
649	sysinfoOfile.write(str(pitch))
650	sysinfoOfile.write(' %pitch \n')
651
652	#Roll
653	readRoll=ifile.read(2)
654	roll=struct.unpack('h',readRoll)[0]
655	sysinfoOfile.write(str(roll))
656	sysinfoOfile.write(' %roll \n')
657
658	#Min Pressure Value (mm)
659	readMinPress=ifile.read(2)
660	minPress=struct.unpack('H',readMinPress)[0]
661	sysinfoOfile.write(str(minPress))
662	sysinfoOfile.write(' %min press (mm) \n')
663
664	#Max Pressure Value (mm)
665	readMaxPress=ifile.read(2)
666
667	#Temperature
668	readTemp=ifile.read(2)
669
670	#hCell Size (counts)
671	readCellSize=ifile.read(2)
672	hCellSize=struct.unpack('H',readCellSize)[0]
673
674
675	#compute the wave cell size
676	cell=hCellSize100.0(750/((headFreq1000.04)/255.0))
677	cellint = math.floor(cell)
678	cellSize=0.01*cellint/256.0
679	sysinfoOfile.write(str('%3.2f' %cellSize))
680	sysinfoOfile.write(' %wave cell size (m) \n')
681
682	#compute the wave cell position
683	pLength=cellSize*100.0
684	cellDist=blanking100.0(750.0/32768.0)-pLength
685	cellDist=0.01*math.floor(cellDist)
686	beamCellPos=cellDist+cellSize
687	cellPos=math.cos(25(math.pi/180))beamCellPos
688	sysinfoOfile.write(str('%3.2f' %cellPos))
689	sysinfoOfile.write(' %wave cell position (m) \n')
690	testOfile.write(str('%3.2f' %cellPos))
691	testOfile.write(' \n')
692
693	#Noise Amplitude Beam1
694	readAmpBeam1=ifile.read(1)
695
696	#Noise Amplitude Beam2
697	readAmpBeam2=ifile.read(1)
698
699	#Noise Amplitude Beam3
700	readAmpBeam3=ifile.read(1)
701
702	#Noise Amplitude Beam4
703	readAmpBeam4=ifile.read(1)
704
705	#Processing Mag Beam1
706	readMagBeam1=ifile.read(2)
707
708	#Processing Mag Beam2
709	readMagBeam2=ifile.read(2)
710
711	#Processing Mag Beam3
712	readMagBeam3=ifile.read(2)
713
714	#Processing Mag Beam4
715	readMagBeam4=ifile.read(2)
716
717	#Spare
718	ifile.read(14)
719
720	#Checksum
721	readChecksum=ifile.read(2)
722
723
724
725
726	#WAVE DATA
727
728
729
730	#what is the wave ID or the extra data structure ID
731	waveID= 0x30a5
732	extraID= 0x42a5
733
734
735	headerID=waveID
736
737	while (headerID == waveID) or (headerID == extraID):
738
739	#Read the Sync and ID as a 2 byte hex instead of seperately
740	readHeaderID=ifile.read(2)
741	if len(readHeaderID) == 0:
742	break
743	headerID=struct.unpack('H',readHeaderID)[0]
744
745	if (headerID != waveID) and (headerID != extraID):
746	break
747
748
749	#if the extra data structure exists, read through it
750	if headerID == extraID:
751	readSize=ifile.read(2)
752	size=struct.unpack('H',readSize)[0]
753	ifile.read(size*2-2)
754
755
756	#size in words
757	readSize=ifile.read(2)
758	size=struct.unpack('H',readSize)[0]
759	#compute the length in bytes - the checksum
760	length=(size*2)-2
761
762	#Pressure
763	readPressure=ifile.read(2)
764	pressure=struct.unpack('H',readPressure)[0]
765	pressureOfile.write(str(pressure))
766	pressureOfile.write('\n')
767
768	#AST 1
769	readAST1=ifile.read(2)
770	AST1=struct.unpack('H',readAST1)[0]
771	rangeOfile.write(str('%06d' % AST1))
772	rangeOfile.write(' \n')
773
774	#Analog input
775	aInput=ifile.read(2)
776
777	#Velocity Beam 1 (mm/s)
778	readVel1=ifile.read(2)
779	vel1=struct.unpack('h',readVel1)[0]
780	orbitOfile.write(str('%06d' % vel1))
781	orbitOfile.write(' ')
782
783	#Velocity Beam 2 (mm/s)
784	readVel2=ifile.read(2)
785	vel2=struct.unpack('h',readVel2)[0]
786	orbitOfile.write(str('%06d' % vel2))
787	orbitOfile.write(' ')
788
789	#Velocity Beam 3 (mm/s)
790	readVel3=ifile.read(2)
791	vel3=struct.unpack('h',readVel3)[0]
792	orbitOfile.write(str('%06d' % vel3))
793	orbitOfile.write('\n')
794
795	#AST 2
796	readAST2=ifile.read(2)
797	AST2=struct.unpack('H',readAST2)[0]
798	rangeOfile.write(str('%06d' % AST2))
799	rangeOfile.write(' \n')
800
801	#Noise Amplitude for Beam 1
802	readAmp1=ifile.read(1)
803	amp1=struct.unpack('B',readAmp1)[0]
804
805	#Noise Amplitude for Beam 2
806	readAmp2=ifile.read(1)
807	amp2=struct.unpack('B',readAmp2)[0]
808
809	#Noise Amplitude for Beam 3
810	readAmp3=ifile.read(1)
811	amp3=struct.unpack('B',readAmp3)[0]
812
813	#AST quality
814	readASTQuality=ifile.read(1)
815	ASTQuality=struct.unpack('B',readASTQuality)[0]
816
817	#CheckSum
818	readChecksum=ifile.read(2)
819
820
821
822	pressureOfile.close()
823	rangeOfile.close()
824	orbitOfile.close()
825	sysinfoOfile.close()
826
827	shutil.copy('pressure','pressure_'+FileStamp+'.txt')
828	shutil.copy('orbit','orbit_'+FileStamp+'.txt')
829	shutil.copy('range','range_'+FileStamp+'.txt')
830	shutil.copy('sysinfo','sysinfo_'+FileStamp+'.txt')
831
832	os.remove('pressure')
833	os.remove('orbit')
834	os.remove('range')
835	os.remove('sysinfo')
836
837
838	#write the rest of the data to the currents file and rename
839
840	currentsHeaderOfile.write('%List of time records for this data set, YYMMDDHHMM \n')
841
842	for date in currentsTime:
843	currentsHeaderOfile.write(str(date))
844	currentsHeaderOfile.write('\n')
845
846	currentsHeaderOfile.write('%List of pressure readings for each current burst \n')
847
848	for press in currentsPress:
849	currentsHeaderOfile.write(str(press/1000.000))
850	currentsHeaderOfile.write('\n')
851
852	currentsHeaderOfile.write('%List of temp readings in degrees C for each current burst \n')
853
854	for temp in currentsTemp:
855	currentsHeaderOfile.write(str(temp/100.00))
856	currentsHeaderOfile.write('\n')
857
858	currentsTiltOfile.write('%List heading,pitch and roll values in degrees for each current burst \n')
859
860	for value in range(len(currentsHeading)):
861	currentsTiltOfile.write(str('%- 8.1f' %(currentsHeading[value]/10.0)))
862	currentsTiltOfile.write(str('%- 8.1f' %(currentsPitch[value]/10.0)))
863	currentsTiltOfile.write(str('%- 8.1f' %(currentsRoll[value]/10.0)))
864	currentsTiltOfile.write('\n')
865
866	#write the vel and amplitude data to the file
867
868	currentsOfile.write('% velocities and amplitudes with lines of Ncells for each burst \n')
869	currentsOfile.write('% velocities from lowest cell up, in m/s ---- amplitudes in counts \n')
870	currentsOfile.write('% vel1(beam1, x or u) vel2(beam2, y or v) vel3(beam3, z or w) ---- amp1(beam1) amp2(beam2) amp3(beam3) \n')
871
872
873	for list in range(len(allVel1)):
874	for i in range(0,Ncells):
875	currentsOfile.write(str('%- 8.3f' % (allVel1[list][i]/1000.000)))
876	currentsOfile.write(str('%- 8.3f' % (allVel2[list][i]/1000.000)))
877	currentsOfile.write(str('%- 8.3f' % (allVel3[list][i]/1000.000)))
878	currentsOfile.write(str('%-8d' % allAmp1[list][i]))
879	currentsOfile.write(str('%-8d' % allAmp2[list][i]))
880	currentsOfile.write(str('%-8d' % allAmp3[list][i]))
881	currentsOfile.write( '\n')
882
883
884	currentsOfile.close()
885	currentsHeaderOfile.close()
886	currentsTiltOfile.close()
887	shutil.copy('currents','currents_'+currentsTime[0]+'.txt')
888	shutil.copy('currentsHeader','currentsHeader_'+currentsTime[0]+'.txt')
889	shutil.copy('currentsTilt','currentsTilt_'+currentsTime[0]+'.txt')
890
891	print 'End Of File, Processing Complete'
892
893	currentsOfile.close()
894	currentsHeaderOfile.close()
895	currentsTiltOfile.close()
896	pressureOfile.close()
897	rangeOfile.close()
898	orbitOfile.close()
899	sysinfoOfile.close()
900
901	os.remove('currents')
902	os.remove('currentsHeader')
903	os.remove('currentsTilt')
904
905	ifile.close()
906	testOfile.close()

Note: See TracBrowser for help on using the browser.

root/DPWP/trunk/DPWavesProc/nortek/nortek.py

Download in other formats: