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

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

# Run cases
cipsaMultiphaseEuler/Allrun
addonMultiphaseEuler/Allrun

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

# Compare to reference solution
latestTime=$(foamListTimes -case cipsaMultiphaseEuler/ -latestTime)
case=$(basename "$PWD")
cp "addonMultiphaseEuler/${latestTime}/U.gas1" \
    "cipsaMultiphaseEuler/${latestTime}/U.gas1.ref"
cp "addonMultiphaseEuler/${latestTime}/T.gas1" \
    "cipsaMultiphaseEuler/${latestTime}/T.gas1.ref"
runApplication -s "${case}.U" foamPostProcess -case cipsaMultiphaseEuler/ \
    -func "deviation(field=U.gas1, reference=U.gas1.ref, maxAbsDev=2e-3)" \
    -latestTime
runApplication -s "${case}.T" foamPostProcess -case cipsaMultiphaseEuler/ \
    -func "deviation(field=T.gas1, reference=T.gas1.ref, maxAbsDev=1e-2)" \
    -latestTime

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