bfast plot gnuplot

#!/bin/bash
#
set -e
#set -x

help()
{
  [ ".$1" != "." ] && echo "ERROR: $1"
  cat<<EOF
Usage:
$0 <data_f>
EOF
  exit 1
}

data_f=$1
[ ".$data_f"  == "." ] && help

plot()
{
  root=$1
  col=$2
  cfg_f="./.tmp.$$.txt"
  (
  cat <<-EOF
set datafile separator ","
set nokey
set terminal png nocrop size 1000,600
set xdata time
set timefmt "%Y-%m-%d %H:%M"
set t png
#set yrange [0:25]
set output "output/${root}.png"
set ylabel "Time computing: ${root}"
set title "bfast processing across SEAs (${root})"
plot "${data_f}" using 1:${col}
EOF
) > $cfg_f

  cat $cfg_f | gnuplot
  rm -f $cfg_f
}

i=2
mkdir -p plots
for r in "load_ref" "load_indexes" "searching" "merging" "elapsed"
do
  plot $r $i
  i=$[$i+1]
done