#!/bin/bash

cd zevsm
> /root/zevs.dat

for T in *
do
   awk '{
      printf( "%d %.6f %.6e\n", T, $1, $2);
   }' T=$T < $T >> /root/zevs.dat

   echo >> /root/zevs.dat
done

T1=`vtdate -in '2017-05-28 06:00'`
T2=`vtdate -in '2017-05-28 21:30'`

echo "
   set terminal png small size 2000,400
   set ylabel 'Frequency, Hz'
   set xlabel 'UTC HH:MM'
   set palette color maxcolors 4096
   set pm3d at b interpolate 450,136e-6
   set view map
   set style data pm3d
   set style func pm3d
   set nokey
   set xdata time
   set timefmt '%s'
   set format x '%H:%M'
   set xtics border mirror out 1800
   set xrange ['$T1':'$T2']
   set yrange [80.5:83.5]
   set xtics 1800
   set mxtics 6
   set ytics border mirror out 0.2
   set cblabel 'Relative power'
        
   set title 'Todmorden H-field 2017-05-26 82Hz. Bearing 34 deg. Relative power in 16.6mHz.  Sample width 60 secs, stepping 60 secs'
   splot '/root/zevs.dat' using 1:2:(\$3 * \$3 * 1e9)
         " | gnuplot > /root/zevs.png