#!/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.airD.gz  "${latestTime}"/alpha.airD.ref.gz
cp validation/alpha.water.gz "${latestTime}"/alpha.water.ref.gz
cp validation/U.airD.gz      "${latestTime}"/U.airD.ref.gz
runApplication -s "${case}.alpha.airD" foamPostProcess \
    -func "deviation(field=alpha.airD, reference=alpha.airD.ref)" -latestTime
runApplication -s "${case}.alpha.water" foamPostProcess \
    -func "deviation(field=alpha.water, reference=alpha.water.ref)" -latestTime
runApplication -s "${case}.U.airD" foamPostProcess \
    -func "deviation(field=U.airD, reference=U.airD.ref)" -latestTime

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