#!/bin/sh

# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions

application="lesCavitatingFoam"

refineMeshByCellSet()
{
    echo "creating cell set for primary zone - $1"
    cp system/cellSetDict.$1 system/cellSetDict
    cellSet >& log.cellSet.$1

    echo "refining primary zone - $1"
    refineMesh -dict -overwrite >& log.refineMesh.$1
}

cd throttle
    runApplication blockMesh

    refineMeshByCellSet 1
    refineMeshByCellSet 2
    refineMeshByCellSet 3

    runApplication $application
cd ..

cd throttle3D
    cp -r 0.org 0

    runApplication blockMesh

    refineMeshByCellSet 1
    refineMeshByCellSet 2
    refineMeshByCellSet 3

    echo "mapping fields from 2D throttle case"
    mapFields ../throttle -sourceTime latestTime >& log.mapFields

    runApplication decomposePar
    hostname > system/machines
    runParallel $application 4 system/machines
    runApplication reconstructPar
cd ..