#!/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 -s air   setBubble -dict setBubbleAirDict
runApplication -s water setBubble -dict setBubbleWaterDict
runApplication foamRun

# Compare to reference solution
latestTime=$(foamListTimes -latestTime)
case=$(basename "$PWD")
cp validation/alpha.air.gz "${latestTime}/alpha.air.ref.gz"
cp validation/shearUnderResolution.air_water.gz \
    "${latestTime}/shearUnderResolution.air_water.ref.gz"
runApplication -s "foamRun.${case}.alpha.air" foamPostProcess \
    -func "deviation(field=alpha.air, reference=alpha.air.ref)" \
    -latestTime
runApplication -s "foamRun.${case}.shearUnderResolution.air_water" \
    foamPostProcess \
    -func \
     "deviation \
     ( \
         field=tmp<shearUnderResolution.air_water>, \
         reference=shearUnderResolution.air_water.ref \
     )" \
    -latestTime

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