#!/bin/bash # make run.sh fail if a subcommand fails set -e set -o pipefail if [ ! -e glassCube_64.hdf5 ] then echo "Fetching initial glass file for Strömgen Sphere 3D example ..." ./getGlass.sh fi if [ ! -f 'stromgrenSphere-3D-HHe.hdf5' ]; then echo "Generating ICs" python3 makeIC_HHe.py fi # Run SWIFT with RT ../../../swift \ --hydro --threads=4 --stars --external-gravity \ --feedback --radiation \ stromgrenSphere-3D-MFHHe.yml 2>&1 | tee output.log # option with mpi # mpirun -np 2 ../../../swift_mpi --hydro --threads=14 --stars --external-gravity --feedback --radiation stromgrenSphere-3D-MFHHe.yml 2>&1 | tee output.log # Plot the Stromgren 3D with hydrogen and helium checks. python3 ./plotStromgren3DMFHHeCheck.py 10