1# awk program to scan loopstats files and report errors/statistics 2# 3# usage: awk -f loop.awk loopstats 4# 5# format of loopstats record 6# MJD sec time (s) freq (ppm) poll 7# 49235 3.943 0.000016 22.4716 6 8# 9# format of output dataset (time values in milliseconds, freq in ppm) 10# loopstats.19960706 11# loop 1180, 0+/-11.0, rms 2.3, freq -24.45+/-0.045, var 0.019 12# 13BEGIN { 14 loop_tmax = loop_fmax = -1e9 15 loop_tmin = loop_fmin = 1e9 16} 17# 18# scan all records in file 19# 20{ 21 if (NF >= 5) { 22 loop_count++ 23 if ($3 > loop_tmax) 24 loop_tmax = $3 25 if ($3 < loop_tmin) 26 loop_tmin = $3 27 if ($4 > loop_fmax) 28 loop_fmax = $4 29 if ($4 < loop_fmin) 30 loop_fmin = $4 31 loop_time += $3 32 loop_time_rms += $3 * $3 33 loop_freq += $4 34 loop_freq_rms += $4 * $4 35 } 36} END { 37 if (loop_count > 0) { 38 loop_time /= loop_count 39 loop_time_rms = sqrt(loop_time_rms / loop_count - loop_time * loop_time) 40 loop_freq /= loop_count 41 loop_freq_rms = sqrt(loop_freq_rms / loop_count - loop_freq * loop_freq) 42 printf "loop %d, %.0f+/-%.1f, rms %.1f, freq %.2f+/-%0.3f, var %.3f\n", loop_count, (loop_tmax + loop_tmin) / 2 * 1e6, (loop_tmax - loop_tmin) / 2 * 1e6, loop_time_rms * 1e6, (loop_fmax + loop_fmin) / 2, (loop_fmax - loop_fmin) / 2, loop_freq_rms 43 } 44} 45 46