proc_num=8

startTotal=$(date "+%s")

startMesh=$(date "+%s")

blockMesh
surfaceFeatureExtract -includedAngle 150 -case . ./constant/triSurface/entities.stl features
decomposePar
changeDictionary -enableFunctionEntries -dict system/changeDictionaryDict-ptscotch -instance system
mpiexec -n $proc_num snappyHexMesh -parallel -overwrite
changeDictionary -enableFunctionEntries -dict system/changeDictionaryDict-scotch -instance system
reconstructParMesh -mergeTol 1E-6 -constant

nowMesh=$(date "+%s")
timeMesh=$((nowMesh-startMesh))
echo "Meshing time used:$timeMesh seconds"


startCalc=$(date "+%s")

decomposePar -force
mpiexec -n $proc_num porousSimpleFoam -parallel
reconstructPar

echo "Meshing time used:$timeMesh seconds"

nowCalc=$(date "+%s")
timeCalc=$((nowCalc-startCalc))
echo "Calc time used:$timeCalc seconds"


nowTotal=$(date "+%s")
timeTotal=$((nowTotal-startTotal))
echo "Total time used:$timeTotal seconds"
