#!/bin/bash
cd "${0%/*}" || exit 1    # Run from this directory

# Source run functions
. "$WM_PROJECT_SITE/$WM_PROJECT_VERSION/bin/tools/RunFunctions"

runApplication blockMesh
runApplication foamRun

case=$(basename "$PWD")
latestTime=$(foamListTimes -latestTime)

# Test solver output
cp "validation/reference/alpha.water.gz" \
   "$latestTime/alpha.water.ref.gz"
runApplication -s "$case" foamPostProcess \
    -func "deviation \
    ( \
        field=alpha.water,  \
        reference=alpha.water.ref,  \
        maxAbsDev=3e-3 \
    )" \
    -latestTime

# Test function object output
actual=validation/reference/line.csv
desired=postProcessing/graph/1/line.csv
mpytest almostequal --decimal 6 "$actual" "$desired"

# Create plots
( cd validation && ./createGraphs )

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