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

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

# Set application name
application="interFoam"

setDamBreakFine ()
{
    blockMeshDict="constant/polyMesh/blockMeshDict"
    controlDict="system/controlDict"
    sed \
    -e s/"23 8"/"46 10"/g \
    -e s/"19 8"/"40 10"/g \
    -e s/"23 42\(.*\) 1 1)"/"46 76\1 2 1)"/g \
    -e s/"4 42\(.*\) 1 1)"/"4 76\1 2 1)"/g \
    -e s/"19 42\(.*\) 1 1)"/"40 76\1 2 1)"/g \
    $blockMeshDict > temp.$$
    mv temp.$$ $blockMeshDict
    sed \
    -e s/"\(deltaT[ \t]*\) 0.001;"/"\1 5e-04;"/g \
    -e s/"\(endTime[ \t]*\) 1;"/"\1 0.4;"/g \
    $controlDict > temp.$$
    mv temp.$$ $controlDict
}

# Do damBreak
(
    cd damBreak || exit

    runApplication blockMesh
    runApplication setFields
    runApplication $application
)

# Clone case
cloneCase damBreak damBreakFine

(
    cd damBreakFine || exit

    # Modify case
    setDamBreakFine
    cp ../damBreak/0/alpha1.org 0/alpha1

    # And execute
    runApplication blockMesh
    runApplication setFields
    runApplication decomposePar
    runParallel $application 4
    runApplication reconstructPar
)

# Do capillaryRise
(
    cd capillaryRise || exit

    runApplication blockMesh
    runApplication setFields
    runApplication $application
)

# ----------------------------------------------------------------- end-of-file