#!/bin/bash

# Generate the initial conditions if they are not present.
if [ ! -e homogeneousCube_32.hdf5 ]
then
    echo "Fetching initial conditions..."
    ./getICS.sh
fi


# Get the Grackle cooling table
if [ ! -e CloudyData_UVB=HM2012.h5 ]
then
    echo "Fetching the Cloudy tables required by Grackle..."
    ../getGrackleCoolingTable.sh
fi


# Create output directory
rm -rf snap
mkdir snap

# Run SWIFT
../../../swift --hydro --cooling --cosmology --threads=14  CoolingWithPrimordialElements.yml

# Check energy conservation and cooling rate
python3 plotAbundances.py snap/snapshot*.hdf5