#!/bin/bash

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

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

runApplication blockMesh
runApplication createZones
runApplication createPatch
runApplication setFields
runApplication foamRun

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

# Compare to reference solution
latestTime=$(foamListTimes -latestTime)
case=$(basename "$PWD")
cp "validation/alpha.airD.gz" "${latestTime}/alpha.airD.ref.gz"
cp "validation/f2.airD.gz" "${latestTime}/f2.airD.ref.gz"
runApplication -s "${case}.alpha.airD" foamPostProcess \
    -func \
     "deviation \
     ( \
        field=alpha.airD, \
        reference=alpha.airD.ref, \
        maxAbsDev=1e-7  \
     )" \
    -latestTime
runApplication -s "${case}.f2.airD" foamPostProcess \
    -func \
     "deviation\
     ( \
        field=f2.airD, \
        reference=f2.airD.ref, \
        maxAbsDev=0.01 \
     )" \
    -latestTime

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