#!/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 setFields
runApplication foamRun

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

# Compare to reference solution
latestTime=$(foamListTimes -latestTime)
case=$(basename "$PWD")
cp validation/alpha.air.gz "${latestTime}"/alpha.air.ref.gz
cp validation/k.water.gz   "${latestTime}"/k.water.ref.gz
runApplication -s "foamRun.${case}.alpha.air" foamPostProcess \
    -func "deviation \
    ( \
        field=alpha.air, \
        reference=alpha.air.ref, \
        maxAbsDev=1e-6 \
    )" \
    -latestTime
runApplication -s "foamRun.${case}.k.water" foamPostProcess \
    -func "deviation \
    ( \
        field=k.water, \
        reference=k.water.ref, \
        maxAbsDev=1e-8 \
    )" \
    -latestTime

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