rpirena/hour_rate.awk

#!/usr/bin/awk -f

@include "RPi.awk"

BEGIN {
    T = 3600  # counting cadence
    tmax = 60 # max time between timestamps
}

T0 && isE() { incrRATE() }

function printRate() {
    if (tt) {
	printf "%.0f", T0+T/2
	for (i in RR) printf " %g %g", RR[i]/tt, sqrt(RR[i]+1)/tt
	printf "\n"
    }
}

END { printRate() }

isH() {
    if (!T0) T0 = Time
    dT = Time - tTime
    tTime = Time
    if (dT>0 && dT<tmax) {
	for (i in R) RR[i] += R[i]
	tt += dT
    }
    for (i in R) R[i] = 0
    if (Time > T0+T) {
	printRate()
	for (i in RR) RR[i] = 0
	T0 = Time
	tt = 0
    }
}