#!/bin/bash

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

gnuplot<<'EOF'

set linetype 1 lc rgb "blue"   lw 3 dt 1 pt 7 ps 0.5
set linetype 2 lc rgb "cyan"   lw 3 dt 1 pt 5 ps 0.5
set linetype 3 lc rgb "green"  lw 3 dt 1
set linetype 2 lc rgb "orange" lw 3 dt 1
set linetype 5 lc rgb "pink"   lw 3 dt 1
set linetype 6 lc rgb "black"  lw 3 dt 1
set linetype 7 lc rgb "black"  lw 2 dt 2

set terminal pngcairo enhanced color

PP  = "../postProcessing"
REF = "./reference"

set output "volumePerPhase.png"
set xlabel "t [s]"
set ylabel "Volume per phase [m^3]"
plot PP."/phasesVolume/0/volFieldValue.dat"  u 1:2 w l ls 1 t "airC", \
     PP."/phasesVolume/0/volFieldValue.dat"  u 1:3 w l ls 2 t "airD", \
     PP."/phasesVolume/0/volFieldValue.dat"  u 1:2 w l ls 3 t "water", \
     REF."/phasesVolume/0/volFieldValue.dat" u 1:2 w l ls 7 t "airC (ref)", \
     REF."/phasesVolume/0/volFieldValue.dat" u 1:3 w l ls 7 t "airD (ref)", \
     REF."/phasesVolume/0/volFieldValue.dat" u 1:2 w l ls 7 t "water (ref)"

set output "volumePerSizeClass.png"
set xlabel "t [s]"
set ylabel "Volume per size class [m^3]"
plot PP."/volumeClass0/0/volFieldValue.dat"  u 1:2 w l ls 1 t "class0", \
     PP."/volumeClass1/0/volFieldValue.dat"  u 1:2 w l ls 2 t "class1", \
     PP."/volumeClass2/0/volFieldValue.dat"  u 1:2 w l ls 3 t "class2", \
     PP."/volumeClass3/0/volFieldValue.dat"  u 1:2 w l ls 2 t "class3", \
     PP."/volumeClass2/0/volFieldValue.dat"  u 1:2 w l ls 5 t "class2", \
     REF."/volumeClass0/0/volFieldValue.dat" u 1:2 w l ls 7 t "class0 (ref)", \
     REF."/volumeClass1/0/volFieldValue.dat" u 1:2 w l ls 7 t "class1 (ref)", \
     REF."/volumeClass2/0/volFieldValue.dat" u 1:2 w l ls 7 t "class2 (ref)", \
     REF."/volumeClass3/0/volFieldValue.dat" u 1:2 w l ls 7 t "class3 (ref)", \
     REF."/volumeClass2/0/volFieldValue.dat" u 1:2 w l ls 7 t "class2 (ref)"

set output 'numberDensity.png'
set grid
set xlabel 'diameter [m]'
set ylabel 'number density [m^{-2}]'
set boxwidth 0.5
plot PP."/numberDensity/2/numberDensity.xy" \
     u 1:2 w boxes fs solid 1 lw 0 t 'Result (t=2 s)', \
     REF."/numberDensity/2/numberDensity.xy" \
     u 1:2 w boxes fill empty lw 2 dt 2 lc rgb 'black' t 'Reference'

set output "totalNumberOfBubbles.png"
set xlabel "t [s]"
set ylabel "Total number of bubbles [-]"
plot PP."/totalNumberOfBubbles/0/volFieldValue.dat"  u 1:2 w l ls 1 t "Total", \
     REF."/totalNumberOfBubbles/0/volFieldValue.dat" u 1:2 w l ls 7 t "Ref."

EOF
