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

# Clean everything
./Allclean

# Create mesh
m4 < m4/blockMeshDict.m4 > constant/polyMesh/blockMeshDict
m4 < m4/mirrorMeshDict.m4 > system/mirrorMeshDict
runApplication blockMesh

# Create second half of mesh
runApplication mirrorMesh
rm -fvf constant/polyMesh/*
mv -fv 0/polyMesh/* constant/polyMesh/.
rm -rvf 0/polyMesh

# Rotate mesh
transformPoints -rollPitchYaw "(90 0 0)" > log.transformPoints1
transformPoints -rollPitchYaw "(0 0 -90)" > log.transformPoints2

# Create missing patches
runApplication createPatch -overwrite

# Renumber mesh
runApplication renumberMesh -overwrite

# Run check mesh
runApplication checkMesh

# Run potentialFoam
runApplication potentialFoam -writep

# Get application name
application=$(getApplication)

# Run case - single
runApplication $application

## Run case - parallel
#runApplication decomposePar
#runParallel $application
