temp_ion.plt 272 Bytes
Newer Older
Francois Gygi committed
1
#!/bin/bash
2 3 4 5 6
if [ $1 == "-range" ]
then
  range=$2
  shift 2
fi
Francois Gygi committed
7
gnuplot -persist <<EOF
8 9 10 11
fit $range a*x+b "<grep -h temp_ion $*" u 0:2 via a,b
fit $range c "<grep -h temp_ion $*" u 0:2 via c
p $range "<grep -h temp_ion $*" u 2 w l, a*x+b, c
print "Tavg=",c,"   dT/dt=",a
Francois Gygi committed
12
EOF