#!/bin/bash

GEOM="LDH"
RE_NUMBERS=("20000" "30000" "35000" "40000" "50000" "60000" "65000" "70000" "80000")
CORES=6
SOLVER="buoyantBoussinesqSimpleFoam"

echo "==================================================="
echo "  Starting Execution for $GEOM Geometry Cases"
echo "==================================================="

for re in "${RE_NUMBERS[@]}"; do
    CASE_DIR="${GEOM}_${re}"
    
    if [ -d "$CASE_DIR" ]; then
        echo "--> Running Case: $CASE_DIR"
        cd "$CASE_DIR" || exit
        
        blockMesh > log.blockMesh 2>&1
        decomposePar -force > log.decomposePar 2>&1
        mpirun -np $CORES $SOLVER -parallel > log.$SOLVER 2>&1
        reconstructPar > log.reconstructPar 2>&1
        
        cd ../
    else
        echo "--> Error: Directory $CASE_DIR not found!"
    fi
done

echo ""
echo "==================================================="
echo "  All $GEOM Re cases executed successfully."
echo "==================================================="
