#!/bin/sh
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
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
-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
application=`getApplication`
runApplication blockMesh
runApplication setFields
runApplication $application
cd ..
# Clone case
cloneCase damBreak damBreakFine
cd damBreakFine
# Modify case
setDamBreakFine
cp ../damBreak/0/alpha1.org 0/alpha1
# And execute
runApplication decomposePar
runParallel $application 4
runApplication reconstructPar