#!/bin/bash #Clears the previous figures echo "Clearing existing figures." if [ -f "orbits_simulation_kpc.png" ]; then rm orbits_simulation_kpc.png fi if [ -f "orbits_simulation_Mpc.png" ]; then rm orbits_simulation_Mpc.png fi if [ ! -e orbit.csv ] then echo "Fetching solution..." wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ReferenceSolutions/MWPotential_2014/orbit.csv fi #Clears the IC file if [ -f "circular_orbits_MW.hdf5" ]; then rm circular_orbits_MW.hdf5 fi #Generates the initial conditions echo "Generate initial conditions for circular orbits" if command -v python3 &>/dev/null; then python3 makeIC.py else python3 makeIC.py fi #Runs the simulation # self gravity G, external potential g, hydro s, threads t and high verbosity v echo "Starting the simulation with the first type of units (kpc)... Look at output_1.log for the simulation details." ../../../swift --external-gravity --threads=8 params_unit_1.yml 2>&1 > output_1.log echo "Simulation ended." echo "Starting the simulation with the second type of units (Mpc)... Look at output_2.log for the simulation details." ../../../swift --external-gravity --threads=8 params_unit_2.yml 2>&1 > output_2.log echo "Simulation ended." #Saves the plots echo "Save plots of orbits." if command -v python3 &>/dev/null; then python3 makePlots.py else python3 makePlots.py fi