#!/bin/sh

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

time=$(foamListTimes -case .. -latestTime)

gnuplot<<EOF
    set terminal postscript eps color enhanced font "Helvetica,20"
    set out "graphs.eps"

    set style line 1 lt 1 linecolor rgb "blue" linewidth 1.5
    set style line 2 lt 2 linecolor rgb "black" linewidth 1.5

    set xlabel "x"
    set ylabel "u'"
    set title "T3A - Flat Plate - turbulent intensity"
    plot [:1.5][:0.05] \
        "../postProcessing/kGraph/$time/line.xy" \
        u (\$1-0.04):(1./5.4*sqrt(2./3.*\$2))title "kOmegaSSTLM" w l ls 1, \
        "exptData/T3A.dat" u (\$1/1000):(\$3/100) title "Exp T3A" w p ls 2

    set xlabel "Re_x"
    set ylabel "c_f"
    set title "T3A - Flat Plate - C_f"
    plot [:6e+5][0:0.01] \
        "../postProcessing/wallShearStressGraph/$time/line.xy" \
        u ((\$1-0.04)*5.4/1.5e-05):(-\$2/0.5/5.4**2) title "kOmegaSSTLM" w l ls 1, \
        "exptData/T3A.dat" u (\$1/1000*5.4/1.51e-05):2 title "Exp" w p ls 2
EOF

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