48 lines
836 B
Awk
48 lines
836 B
Awk
|
#!/usr/bin/awk -f
|
||
|
|
|
||
|
|
BEGIN{ if (!U) U="THREED" }
|
||
|
|
@include "RPi.awk"
|
||
|
|
|
||
|
|
BEGIN {
|
||
|
|
T = 500 # integration time
|
||
|
|
dT0 = 0
|
||
|
|
tmax = 180 # time between P
|
||
|
|
min_dtime = 10000
|
||
|
|
pmin = 500
|
||
|
|
pmax = 1100
|
||
|
|
}
|
||
|
|
|
||
|
|
isGPS() {}
|
||
|
|
|
||
|
|
isH() { iTime = Time - pTime }
|
||
|
|
|
||
|
|
isP(1) >= 2 {
|
||
|
|
if (p && iTime>0) {
|
||
|
|
dt = iTime
|
||
|
|
if (dt>tmax) dt = tmax
|
||
|
|
tt += dt
|
||
|
|
tp += dt*p
|
||
|
|
}
|
||
|
|
if (iTime > tmax || tt+dT0 >= T) {
|
||
|
|
if (tt>0) {
|
||
|
|
printf "%.0f %.1f", T0 + tt/2 + file_epoch, tp/tt
|
||
|
|
for (j=0; j<8; j++) printf " %g %g", R[j]/tt, sqrt(R[j]+1)/tt
|
||
|
|
printf "\n"
|
||
|
|
}
|
||
|
|
if (iTime > tmax) dT0=0
|
||
|
|
else dT0 = tt-T
|
||
|
|
tp = 0
|
||
|
|
tt = 0
|
||
|
|
delete R
|
||
|
|
T0 = Time
|
||
|
|
}
|
||
|
|
if (!p) T0=Time
|
||
|
|
p = Pressure
|
||
|
|
pTime = Time
|
||
|
|
iTime = 0
|
||
|
|
}
|
||
|
|
|
||
|
|
p && isEE() && iTime>=0 && iTime<=tmax && $2==0 {
|
||
|
|
l2 = strtonum($4)
|
||
|
|
for (j=0; j<8; j++) if (and(l2,lshift(1,j))) R[j]++
|
||
|
|
}
|