rpirena/l2_rates.awk
stephan 14f8f107d3 cospi: analysis work, epoch
git-svn-id: svn+ssh://asterix.ieap.uni-kiel.de/home/subversion/stephan/solo/eda/cospi/host@8187 bc5caf13-1734-44f8-af43-603852e9ee25
2021-03-12 20:38:26 +00:00

48 lines
836 B
Awk
Executable file

#!/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]++
}