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

runApplication foamPostProcess -latestTime -func "
    graphCell
    (
        name=graph,
        start=( 0     0.0  0.05),
        end=  ( 0.025 0.0  0.05),
        axis=x,
        fields=( tmp<((3*sqr(min((max((yWall-zeroWallDist),0)|(Cd*d.air)),1)))\
-(2*pow3(min((max((yWall-zeroWallDist),0)|(Cd*d.air)),1))))> ),
        setFormat=csv,
        separator=\" \"
    )"

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

# Compare to reference solution
case=$(basename "$PWD")
runCommand --suffix "${case}" -- diff -sb \
    validation/line.csv \
    postProcessing/graph/0.2/line.csv
#------------------------------------------------------------------------------
