#!/bin/bash
cd "${0%/*}" || exit                                # Run from this directory
. ${WM_PROJECT_DIR:?}/bin/tools/RunFunctions        # Tutorial run functions
#------------------------------------------------------------------------------

plot_results() {

    image="plots/surface_compare.png"
    gnuplot<<EOF
    reset

    set terminal pngcairo font "helvetica,20" size 60.0in,70.0in
    set grid
    set xrange [0:1]
    set yrange [0:0.25]
    set xlabel "x [m]"
    set ylabel "y [m]"
    set output "$image"

    set style rect fc lt -1 fs solid 0.15 noborder
    set object  1 rect from 0.29, graph 0 to 0.59, graph 1

    array exp_data[10] = \
        ['4', '6', '8', '10', '12', '14', '16', '18', '20', '22']
    array time[10] = \
        ['0.4', '0.6', '0.8', '1', '1.2', '1.4', '1.6', '1.8', '2', '2.2']
    set multiplot \
        layout 5,2 \
        title "Free surface elevation y vs horizontal position x" \
        font "helvetica,42"

    do for[k=1:10] {
    set title "t = ".time[k]." [s]"
    plot \
        "postProcessing/freeSurface/".time[k]."/alpha.water_freeSurface.raw" \
            w points pt 7 ps 2 lc rgb "red"  t "OpenFOAM", \
        "experimentalData/data_".exp_data[k].".dat" \
            w points pt 7 ps 3 lc rgb "black"  t "Liu et al."
    }
    unset multiplot
EOF
}


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

# Requires gnuplot
command -v gnuplot >/dev/null || {
    echo "FOAM FATAL ERROR: gnuplot not found - skipping graph creation" 1>&2
    exit 1
}


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

dirPlots="plots"
[ -d "$dirPlots" ] || mkdir -p "$dirPlots"

plot_results


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