/contrib/ntp/scripts/stats/clock.awk

  1# awk program to scan clockstat files and report errors/statistics
  2#
  3# usage: awk -f check.awk clockstats
  4#
  5# This program works for the following radios:
  6# PST/Traconex 1020 WWV reciever
  7# Arbiter 1088 GPS receiver
  8# Spectracom 8170/Netclock-2 WWVB receiver
  9# IRIG audio decoder
 10# Austron 2200A/2201A GPS receiver (see README.austron file)
 11#
 12BEGIN {
 13	etf_min = osc_vmin = osc_tmin = 1e9
 14	etf_max = osc_vmax = osc_tmax = -1e9
 15}
 16#
 17# scan all records in file
 18#
 19{
 20	#
 21	# select PST/Traconex WWV records
 22	# 00:00:37.234  96/07/08/190 O6@0:5281825C07510394
 23	#
 24	if (NF >= 4 && $3 == "127.127.3.1") {
 25		if (substr($6, 14, 4) > "0010")
 26			wwv_sync++
 27		if (substr($6, 13, 1) == "C")
 28			wwv_wwv++
 29		if (substr($6, 13, 1) == "H")
 30			wwv_wwvh++
 31		x = substr($6, 12, 1)
 32		if (x == "1")
 33			wwv_2.5++
 34		else if (x == "2")
 35			wwv_5++
 36		else if (x == "3")
 37			wwv_10++
 38		else if (x == "4")
 39			wwv_15++
 40		else if (x == "5")
 41			wwv_20++
 42		continue
 43	}
 44	#
 45	# select Arbiter GPS records
 46	# 96 190 00:00:37.000 0 V=08 S=44 T=3 P=10.6 E=00
 47	# N39:42:00.951 W075:46:54.880 210.55      2.50 0.00
 48	#
 49	if (NF >= 4 && $3 == "127.127.11.1") {
 50		if (NF > 8) {
 51			arb_count++
 52			if ($7 != 0)
 53				arb_sync++
 54			x = substr($10, 3, 1)
 55			if (x == "0")
 56				arb_0++
 57			else if (x == "1")
 58				arb_1++
 59			else if (x == "2")
 60				arb_2++
 61			else if (x == "3")
 62				arb_3++
 63			else if (x == "4")
 64				arb_4++
 65			else if (x == "5")
 66				arb_5++
 67			else if (x == "6")
 68			arb_6++
 69		} else if (NF == 8) {
 70			arbn++
 71			arb_mean += $7
 72			arb_rms += $7 * $7
 73			if (arbn > 0) {
 74				x = $7 - arb_val
 75				arb_var += x * x
 76			}
 77			arb_val = $7
 78		}
 79		continue
 80	}
 81	#
 82	# select Spectracom WWVB records
 83	# see summary for decode
 84	#   96 189 23:59:32.248  D
 85	#
 86	if (NF >= 4 && $3 == "127.127.4.1") {
 87		if ($4 == "SIGNAL" || NF > 7)
 88			printf "%s\n", $0
 89		else {
 90			wwvb_count++
 91			if ($4 ~ /\?/)
 92				wwvb_x++
 93			else if ($4 ~ /A/)
 94				wwvb_a++
 95			else if ($4 ~ /B/)
 96				wwvb_b++
 97			else if ($4 ~ /C/)
 98				wwvb_c++
 99			else if ($4 ~ /D/)
100				wwvb_d++
101		}
102		continue
103	}
104	#
105	# select IRIG audio decoder records
106	# see summary for decode
107	#
108	if (NF >= 4 && $3 == "127.127.6.0") {
109		irig_count++
110		if ($5 ~ /\?/)
111			irig_error++
112		continue
113	}
114	#
115	# select Austron GPS LORAN ENSEMBLE records
116	# see summary for decode
117	#
118	else if (NF >= 13 && $6 == "ENSEMBLE") {
119		ensemble_count++
120		if ($9 <= 0)
121			ensemble_badgps++
122		else if ($12 <= 0)
123			ensemble_badloran++
124		else {
125			if ($13 > 200e-9 || $13 < -200e-9)
126				ensemble_200++
127			else if ($13 > 100e-9 || $13 < -100e-9)
128				ensemble_100++
129			ensemble_mean += $13
130			ensemble_rms += $13 * $13
131		}
132		continue
133	}
134	#
135	# select Austron LORAN TDATA records
136	# see summary for decode; note that signal quality log is simply
137	# copied to output
138	#
139	else if (NF >= 7 && $6 == "TDATA") {
140                tdata_count++
141                for (i = 7; i < NF; i++) {
142                        if ($i == "M" && $(i+1) == "OK") {
143                                i += 5
144                                m += $i
145                		tdata_m++
146		        }
147                        else if ($i == "W" && $(i+1) == "OK") {
148                                i += 5
149                                w += $i
150                        	tdata_w++
151			}
152                        else if ($i == "X" && $(i+1) == "OK") {
153                                i += 5
154                                x += $i
155                        	tdata_x++
156			}
157                        else if ($i == "Y" && $(i+1) == "OK") {
158                                i += 5
159                                y += $i
160                        	tdata_y++
161			}
162                        else if ($i == "Z" && $(i+1) == "OK") {
163                                i += 5
164                                z += $i
165                        	tdata_z++
166			}
167		}	
168		continue
169	}
170	#
171	# select Austron ITF records
172	# see summary for decode
173	#
174	else if (NF >= 13 && $5 == "ITF" && $12 >= 100) {
175		itf_count++
176		if ($9 > 200e-9 || $9 < -200e-9)
177			itf_200++
178		else if ($9 > 100e-9 || $9 < -100e-9)
179			itf_100++
180		itf_mean += $9
181		itf_rms += $9 * $9
182		itf_var += $10 * $10
183		continue
184	}
185	#
186	# select Austron ETF records
187	# see summary for decode
188	#
189	else if (NF >= 13 && $5 == "ETF" && $13 >= 100) {
190		etf_count++
191		if ($6 > etf_max)
192			etf_max = $6
193		else if ($6 < etf_min)
194			etf_min = $6
195		etf_mean += $6
196		etf_rms += $6 * $6
197		etf_var += $9 * $9
198		continue
199	}
200	#
201	# select Austron TRSTAT records
202	# see summary for decode
203	#
204	else if (NF >= 5 && $5 == "TRSTAT") {
205		trstat_count++
206		j = 0
207		for (i = 6; i <= NF; i++)
208			if ($i == "T")
209				j++
210		trstat_sat[j]++
211		continue
212	}
213	#
214	# select Austron ID;OPT;VER records
215	#
216	# config GPS 2201A TTY1 TC1 LORAN IN OUT1 B.00 B.00 28-Apr-93
217	#
218	# GPS 2201A	receiver model
219	# TTY1		rs232 moduel
220	# TC1		IRIG module
221	# LORAN		LORAN assist module
222	# IN		input module
223	# OUT1		output module
224	# B.00 B.00	firmware revision
225	# 28-Apr-9	firmware date3
226        #
227	else if (NF >= 5 && $5 == "ID;OPT;VER") {
228		id_count++
229		id_temp = ""
230		for (i = 6; i <= NF; i++)
231			id_temp = id_temp " " $i
232		if (id_string != id_temp)
233			printf "config%s\n", id_temp
234		id_string = id_temp
235		continue	
236	}
237	#
238	# select Austron POS;PPS;PPSOFF records
239	#
240	# position +39:40:48.425 -075:45:02.392 +74.09 Stored UTC 0 200 0
241	#
242	# +39:40:48.425	position north latitude
243	# -075:45:02.392 position east longitude
244	# +74.09	elevation (meters)
245	# Stored	position is stored
246	# UTC		time is relative to UTC
247	# 0 200 0	PPS offsets
248	#
249	else if (NF >= 5 && $5 == "POS;PPS;PPSOFF") {
250		pos_count++
251		pos_temp = ""
252		for (i = 6; i <= NF; i++)
253			pos_temp = pos_temp " " $i
254		if (pos_string != pos_temp)
255			printf "position%s\n", pos_temp
256		pos_string = pos_temp
257	continue
258	}
259	#
260	# select Austron OSC;ET;TEMP records
261	#
262	# loop 1121 Software Control Locked
263	#
264	# 1121		oscillator type
265	# Software Control loop is under software control
266	# Locked	loop is locked
267	#
268	else if (NF >= 5 && $5 == "OSC;ET;TEMP") {
269		osc_count++
270		osc_temp = $6 " " $7 " " $8 " " $9
271		if (osc_status != osc_temp)
272			printf "loop %s\n", osc_temp
273		osc_status = osc_temp
274		if ($10 > osc_vmax)
275			osc_vmax = $10
276		if ($10 < osc_vmin)
277			osc_vmin = $10
278		if ($11 > osc_tmax)
279			osc_tmax = $11
280		if ($11 < osc_tmin)
281			osc_tmin = $11
282	continue
283	}
284	#
285	# select Austron UTC records
286	# these ain't ready yet
287	#
288	else if (NF >= 5 && $5 == "UTC") {
289		utc_count++
290		utc_temp = ""
291		for (i = 6; i <= NF; i++)
292			utc_temp = utc_temp " " $i
293		if (utc_string != utc_temp)
294#			printf "utc%s\n", utc_temp
295                utc_string = utc_temp
296	continue
297	}
298} END {
299#
300# PST/Traconex WWV summary data
301#
302	if (wwv_wwv + wwv_wwvh > 0)
303		printf "wwv %d, wwvh %d, err %d, MHz (2.5) %d, (5) %d, (10) %d, (15) %d, (20) %d\n", wwv_wwv, wwv_wwvh, wwv_sync, wwv_2.5, wwv_5, wwv_10, wwv_15, wwv_20
304#
305# Arbiter 1088 summary data
306#
307# gps		record count
308# err		error count
309# sats(0-6)	satellites tracked
310# mean		1 PPS mean (us)
311# rms		1 PPS rms error (us)
312# var		1 PPS Allan variance
313#
314	if (arb_count > 0) {
315		printf "gps %d, err %d, sats(0-6) %d %d %d %d %d %d %d", arb_count, arb_sync, arb_0, arb_1, arb_2, arb_3, arb_4, arb_5, arb_6
316		if (arbn > 1) {
317			arb_mean /= arbn
318			arb_rms = sqrt(arb_rms / arbn - arb_mean * arb_mean)
319			arb_var = sqrt(arb_var / (2 * (arbn - 1)))
320			printf ", mean %.2f, rms %.2f, var %.2e\n", arb_mean, arb_rms, arb_var * 1e-6
321		} else {
322			printf "\n"
323		}
324	}
325#
326# ensemble summary data
327#
328# ensemble	record count
329# badgps	gps data unavailable
330# badloran	loran data unavailable
331# rms		ensemble rms error (ns)
332# >200		ensemble error >200 ns
333# >100		100 ns < ensemble error < 200 ns
334#
335	if (ensemble_count > 0) {
336		ensemble_mean /= ensemble_count
337		ensemble_rms = sqrt(ensemble_rms / ensemble_count - ensemble_mean * ensemble_mean) * 1e9 
338		printf "ensemble %d, badgps %d, badloran %d, rms %.1f, >200 %d, >100 %d\n", ensemble_count, ensemble_badgps, ensemble_badloran, ensemble_rms, ensemble_200, ensemble_100
339	}
340#
341# wwvb summary data
342#
343# wwvb		record count
344# ?		unsynchronized
345# >1		error > 1 ms
346# >10		error > 10 ms
347# >100		error > 100 ms
348# >500		error > 500 ms
349#
350	if (wwvb_count > 0)
351		printf "wwvb %d, ? %d, >1 %d, >10 %d, >100 %d, >500 %d\n", wwvb_count, wwvb_x, wwvb_a, wwvb_b, wwvb_c, wwvb_d
352#
353# irig summary data
354#
355# irig		record count
356# err		error count
357#
358	if (irig_count > 0)
359		printf "irig %d, err %d\n", irig_count, irig_error
360#
361# tdata summary data
362#
363# tdata		record count
364# m		M master OK-count, mean level (dB)
365# w		W slave OK-count, mean level (dB)
366# x		X slave OK-count, mean level (dB)
367# y		Y slave OK-count, mean level (dB)
368# z		Z slave OK-count, mean level (dB)
369#
370	if (tdata_count > 0 ) {
371		if (tdata_m > 0)
372			m /= tdata_count
373		if (tdata_x > 0)
374			w /= tdata_count
375		if (tdata_x > 0)
376			x /= tdata_count
377		if (tdata_y > 0)
378			y /= tdata_count
379		if (tdata_z > 0)
380			z /= tdata_count
381		printf "tdata %d, m %d %.1f, w %d %.1f, x %d %.1f, y %d %.1f, z %d %.1f\n", tdata_count, tdata_m, m, tdata_w, w, tdata_x, x, tdata_y, y, tdata_z, z
382	}
383#
384# itf summary data
385#
386# itf		record count
387# rms		itf rms error (ns)
388# >200		itf error > 200 ns
389# >100		itf error > 100 ns
390# var		Allan variance
391#
392	if (itf_count > 1) { 
393		itf_mean /= itf_count
394		itf_rms = sqrt(itf_rms / itf_count - itf_mean * itf_mean) * 1e9
395		itf_var = sqrt(itf_var / (2 * (itf_count - 1)))
396		printf "itf %d, rms %.1f, >200 %d, >100 %d, var %.2e\n", itf_count, itf_rms, itf_200, itf_100, itf_var
397	}
398#
399# etf summary data
400#
401# etf		record count
402# mean		etf mean (ns)
403# rms		etf rms error (ns)
404# max		etf maximum (ns)
405# min		etf minimum (ns)
406# var		Allan variance
407#
408	if (etf_count > 0) {
409                etf_mean /= etf_count
410		etf_rms = sqrt(etf_rms / etf_count - etf_mean * etf_mean)
411		etf_var = sqrt(etf_var / (2 * (etf_count - 1)))
412		printf "etf %d, mean %.1f, rms %.1f, max %d, min %d, var %.2e\n", etf_count, etf_mean, etf_rms, etf_max, etf_min, etf_var
413	}
414#
415# trstat summary data
416#
417# trstat	record count
418# sat		histogram of tracked satellites (0 - 7)
419#
420	if (trstat_count > 0)
421		printf "trstat %d, sat %d %d %d %d %d %d %d %d\n", trstat_count, trstat_sat[0], trstat_sat[1], trstat_sat[2], trstat_sat[2], trstat_sat[3], trstat_sat[4], trstat_sat[5], trstat_sat[6], trstat_sat[7]
422#
423# osc summary data
424#
425# osc		record count
426# control	control midrange (V) +/- deviation (mV)
427# temp		oven temperature midrange +/- deviation (deg C)
428#
429	if (osc_count > 0)
430		printf "osc %d, control %.3f+/-%.3f, temp %.1f+/-%.2f\n", osc_count, (osc_vmax + osc_vmin) / 2, (osc_vmax - osc_vmin) / 2 * 1e3, (osc_tmax + osc_tmin) / 2, (osc_tmax - osc_tmin) / 2
431}