#!/bin/sh
#------------------------------------------------------------------------------
# =========                 |
# \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
#  \\    /   O peration     |
#   \\  /    A nd           | Copyright (C) 2009-2009 OpenCFD Ltd.
#    \\/     M anipulation  |
#-------------------------------------------------------------------------------
# License
#     This file is part of OpenFOAM.
#
#     OpenFOAM is free software; you can redistribute it and/or modify it
#     under the terms of the GNU General Public License as published by the
#     Free Software Foundation; either version 2 of the License, or (at your
#     option) any later version.
#
#     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
#     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
#     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
#     for more details.
#
#     You should have received a copy of the GNU General Public License
#     along with OpenFOAM; if not, write to the Free Software Foundation,
#     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
#
# Script
#     createGraphs
#
# Description
#     Creates .eps graphs of OpenFOAM results vs experiment for the buoyant
#     cavity case
#
#------------------------------------------------------------------------------

createEpsT() {
    index=$1
    OF=$2
    EXPT=$3

    gnuplot<<EOF
    set terminal postscript eps color enhanced
    set output "OF_vs_EXPT_T$i.eps"
    set xlabel "Channel width, x / [m]"
    set ylabel "Temperature / [K]"
    set grid
    set key left top
    set size 0.6, 0.6
    set xrange [0:0.08]
    set yrange [285:310]
    plot \
        "$EXPT" u (\$1/1000):(\$2+273.15) title "Expt 0.$index" with points lt 1 pt 6, \
        "$OF" title "OpenFOAM 0.$index" with lines linetype -1
EOF
}


createEpsU() {
    index=$1
    OF=$2
    EXPT=$3

    gnuplot<<EOF
    set terminal postscript eps color enhanced
    set output "OF_vs_EXPT_U$i.eps"
    set xlabel "Channel width, x / [m]"
    set ylabel "Vertical velocity component, Uy / [m/s]"
    set grid
    set key left top
    set size 0.6, 0.6
    set xrange [0:0.08]
    set yrange [-0.2:0.2]
    plot \
        "$EXPT" u (\$1/1000):(\$2) title "Expt 0.$index" with points lt 1 pt 6, \
        "$OF" u 1:3 title "OpenFOAM 0.$index" with lines linetype -1
EOF
}


# test if gnuplot exists on the system
type -P gnuplot &>/dev/null || {
    echo "gnuplot not found - skipping graph creation" >&2; exit 1;
}

# paths to data
LATESTTIME=`ls ../sets`
OFDATAROOT=../sets/$LATESTTIME

EXPTDATAROOT=./exptData

# generate temperature profiles
TSets="1 3 4 5 6 7 9"
for i in $TSets; do
    echo "    processing temperature profile at y/yMax of 0.$i"

    OF="$OFDATAROOT/y0.${i}_T.xy"
    EXPT="$EXPTDATAROOT/mt_z0_${i}0_lo.dat"

    createEpsT $i $OF $EXPT
done


# generate velocity profiles
USets="1 3 4 5 6 7 9"
for i in $USets; do
    echo "    processing velocity profile at y/yMax of 0.$i"

    OF="$OFDATAROOT/y0.${i}_U.xy"
    EXPT="$EXPTDATAROOT/mv_z0_${i}0_lo.dat"

    createEpsU $i $OF $EXPT
done

echo "done"

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