# -*- gnuplot -*-

fnH = "<grep -a ^H ".ARG1
fnP = "<./pressure.awk ".ARG1

set y2tics
set ytics nomirror
set ylab "Temperature [°C]"
set y2lab "Pressure [mbar]"
set xlab "Time [h]"

R1   = 3.3e3
R25  = 3.3e3
B25  = 4500.0
R(a) = R1*a/(4096-a)
degC(a) = B25/(log(R(a)/R25)+B25/298.0) - 273

t0=`date +%s -d 0`
t(x) = (x-t0)/3600.0

plot \
 fnH u (t($2)):(degC($7))     axis x1y1 tit "T_{FPGA}" w l, \
 fnH u (t($2)):($8<4090?degC($8):NaN)     axis x1y1 tit "T_{OUT}"  w l, \
 fnH u (t($2)):($9<4090?degC($9):NaN)     axis x1y1 tit "T_{DET}"  w l, \
 fnP."|grep ^PP0" u (t($2)):3 axis x1y1 tit "T_{BATE0}"  w l, \
 fnP."|grep ^PP1" u (t($2)):3 axis x1y1 tit "T_{BATE1}"  w l, \
 fnP."|grep ^PP0" u (t($2)):4 axis x1y2 tit "Pressure 0"  w l, \
 fnP."|grep ^PP1" u (t($2)):4 axis x1y2 tit "Pressure 1"  w l