#!/bin/bash

if ! which gnuplot > /dev/null 2>&1
then
    echo 'gnuplot not found - skipping graph creation' >&2
    exit 1
fi

gnuplot<<EOF
    set grid
    set key outside below

    #---PNG---
    set term pngcairo size 800,600
    FType = '.png'

    #---Track-Bubbles---
    PATH =      "../postProcessing/trackPhaseStructures/0/"
    LIST =      system("cd ".PATH." && ls -1 [0-9]* | sort -g")
    NUM =       words(LIST)
    IDS(n) =    word(LIST,n)
    FILES(n) =  PATH.IDS(n)
    TIMES(n) =  FILES(n)
    do for [i=1:NUM] {
        TIMES(i) =  system("awk 'FNR == 1 {print \$2}' ".FILES(i))
    }

    RPATH =     "./"
    RLIST =     system("cd ".PATH." && ls -1 [0-9]* | sort -g")
    RNUM =      words(RLIST)
    RIDS(n) =   word(RLIST,n)
    RFILES(n) = RPATH.RIDS(n)
    RTIMES(n) =  RFILES(n)
    do for [i=1:RNUM] {
        RTIMES(i) =  system("awk 'FNR == 1 {print \$2}' ".RFILES(i))
    }

    set output "volume".FType
    set xlabel "t"
    set ylabel "V"
    set format y "%.1e"
    set yrange [0:]
    plot for [i=1:RNUM] RFILES(i) u 2:9 w l lt i dt 2 lw 2 notitle, \
         for [i=1:NUM]  FILES(i)  u 2:9 w l lt i notitle,           \
         1/0 lw 2 dt 2 lc 0 t "reference"
    set format y "%g"

    set output "position".FType
    set xlabel "t - t_0"
    set ylabel "y"
    plot for [i=1:RNUM] RFILES(i) u (\$2-(RTIMES(i)+0)):4   \
            w l lt i lw 2 dt 2 notitle,                     \
         for [i=1:NUM]  FILES(i)  u (\$2-(TIMES(i) +0)):4   \
            w l lt i           notitle,                     \
         1/0 lw 2 dt 2 lc 0 t "reference"

    set output "verticalVelocity".FType
    set ylabel "U.y"
    plot for [i=1:RNUM] RFILES(i) u (\$2-(RTIMES(i)+0)):7   \
            w l lt i lw 2 dt 2 notitle,                     \
         for [i=1:NUM]  FILES(i)  u (\$2-(TIMES(i) +0)):7   \
            w l lt i notitle,                               \
         1/0 lw 2 dt 2 lc 0 t "reference"


    #---PNG---
    set term pngcairo size 350,600
    set key outside below


    set output "path".FType
    set size ratio -1
    set xlabel "x"
    set ylabel "y"
    set xrange [0:0.05]
    set yrange [0:0.1]
    plot for [i=1:RNUM] RFILES(i) u 3:4 w l lt i lw 2 dt 2 notitle, \
         for [i=1:NUM]  FILES(i)  u 3:4 w l lt i notitle, \
         1/0 lw 2 dt 2 lc 0 t "reference"

EOF

#------------------------------------------------------------------------------
