#!/bin/bash

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

gnuplot<<EOF
    set terminal pngcairo enhanced color
    set xlabel 'r (mm)'
    set xrange [0:25]
    set key above

    FILEREF='reference/line.csv'
    FILE='../postProcessing/graph/1/line.csv'

    set ylabel 'alpha (-)'
    set output 'alpha.png'
    plot FILEREF u (\$1*1000):(column('alpha.air1')) t 'reference air1' w l lw 3 ,\
         FILEREF u (\$1*1000):(column('alpha.air2')) t 'reference air2' w l lw 3, \
         FILE u (\$1*1000):(column('alpha.air1')) t 'air1' w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('alpha.air2')) t 'air2' w l dt 2 lc black lw 3

    set ylabel 'w (m/s)'
    set output 'velocity.png'
    plot FILEREF u (\$1*1000):(column('U.air1_z')) t 'reference air1' w l lw 3 ,\
         FILEREF u (\$1*1000):(column('U.air2_z')) t 'reference air2' w l lw 3 ,\
         FILE u (\$1*1000):(column('U.air1_z')) t 'air1' w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('U.air2_z')) t 'air2' w l dt 2 lc black lw 3

    set ylabel 'k (m^2/s^2)'
    set output 'turbulence.png'
    plot FILEREF u (\$1*1000):(column('k.water')) t 'reference k.water' w l lw 3 ,\
         FILE u (\$1*1000):(column('k.water')) t 'k.water' w l dt 2 lc black lw 3

    set ylabel 'F_z (N/m^3)'

    set output 'forcesAxial.air.png'
    plot FILEREF u (\$1*1000):(column('dragForce.air1_z'))        t 'reference drag air1'         w l lw 3 ,\
         FILEREF u (\$1*1000):(column('dragForce.air2_z'))        t 'reference drag air2'         w l lw 3 ,\
         FILEREF u (\$1*1000):(column('virtualMassForce.air1_z')) t 'reference virtual mass air1' w l lw 3 ,\
         FILEREF u (\$1*1000):(column('virtualMassForce.air2_z')) t 'reference virtual mass air2' w l lw 3 , \
         FILE u (\$1*1000):(column('dragForce.air1_z'))        t 'drag air1'         w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('dragForce.air2_z'))        t 'drag air2'         w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('virtualMassForce.air1_z')) t 'virtual mass air1' w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('virtualMassForce.air2_z')) t 'virtual mass air2' w l dt 2 lc black lw 3

    set ylabel 'F_r (N/m^3)'

    set output 'forcesLateral.air.png'
    plot FILEREF u (\$1*1000):(column('liftForce.air1_x'))                t 'reference lift air1'                 w l lw 3 ,\
         FILEREF u (\$1*1000):(column('liftForce.air2_x'))                t 'reference lift air2'                 w l lw 3 ,\
         FILEREF u (\$1*1000):(column('wallLubricationForce.air1_x'))     t 'reference wall lubrication air1'     w l lw 3 ,\
         FILEREF u (\$1*1000):(column('wallLubricationForce.air2_x'))     t 'reference wall lubrication air2'     w l lw 3 ,\
         FILEREF u (\$1*1000):(column('turbulentDispersionForce.air1_x')) t 'reference turbulent dispersion air1' w l lw 3 ,\
         FILEREF u (\$1*1000):(column('turbulentDispersionForce.air2_x')) t 'reference turbulent dispersion air2' w l lw 3 ,\
         FILE u (\$1*1000):(column('liftForce.air1_x'))                t 'lift air1'                 w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('liftForce.air2_x'))                t 'lift air2'                 w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('wallLubricationForce.air1_x'))     t 'wall lubrication air1'     w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('wallLubricationForce.air2_x'))     t 'wall lubrication air2'     w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('turbulentDispersionForce.air1_x')) t 'turbulent dispersion air1' w l dt 2 lc black lw 3 ,\
         FILE u (\$1*1000):(column('turbulentDispersionForce.air2_x')) t 'turbulent dispersion air2' w l dt 2 lc black lw 3

EOF

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