xref: /freebsd-src/contrib/ntp/scripts/stats/ensemble.awk (revision cfe30d02adda7c3b5c76156ac52d50d8cab325d9)
1*c0b746e5SOllivier Robert# program to produce loran ensemble statistics from clockstats files
2*c0b746e5SOllivier Robert#
3*c0b746e5SOllivier Robert# usage: awk -f ensemble.awk clockstats
4*c0b746e5SOllivier Robert#
5*c0b746e5SOllivier Robert# format of input record (time values in seconds)
6*c0b746e5SOllivier Robert# 49165 8.628 127.127.10.1 93:178:00:00:07.241 LORAN ENSEMBLE
7*c0b746e5SOllivier Robert# -6.43E-08 +5.02E-08 .091 +5.98E-08 +1.59E-08 .909 +4.85E-08 +3.52E-08
8*c0b746e5SOllivier Robert#
9*c0b746e5SOllivier Robert# format of output record (time values in nanoseconds)
10*c0b746e5SOllivier Robert#  MJD       sec     GPS    wgt    LORAN   wgt      avg   sigma
11*c0b746e5SOllivier Robert# 49165     8.628   -64.3  0.091    59.8  0.909    48.5    35.2
12*c0b746e5SOllivier Robert#
13*c0b746e5SOllivier Robert# select LORAN ENSEMBLE records with valid format and weights
14*c0b746e5SOllivier Robert{
15*c0b746e5SOllivier Robert	if (NF >= 14 && $6 == "ENSEMBLE" && $9 > 0 && $12 > 0)
16*c0b746e5SOllivier Robert		printf "%5s %9.3f %7.1f %6.3f %7.1f %6.3f %7.1f %7.1f\n", $1, $2, $7*1e9, $9, $10*1e9, $12, $13*1e9, $14*1e9
17*c0b746e5SOllivier Robert}
18