#!/bin/bash

#Creates a directory for the outputs
DIR=output_1 #First test of units conversion
if [ -d "$DIR" ];
then
    echo "$DIR directory exists. Its content will be removed."
    rm $DIR/output_*
else
    echo "$DIR directory does not exists. It will be created."
    mkdir $DIR
fi

DIR=output_2 #Second test of units conversion
if [ -d "$DIR" ];
then
    echo "$DIR directory exists. Its content will be removed."
    rm $DIR/output_*
else
    echo "$DIR directory does not exists. It will be created."
    mkdir $DIR
fi

#Clears the previous figures
echo "Clearing existing figures."
if [ -f "circular_orbits_simulation_kpc.png" ];
then
    rm circular_orbits_simulation_kpc.png
fi

if [ -f "circular_orbits_simulation_Mpc.png" ];
then
    rm circular_orbits_simulation_Mpc.png
fi

if [ -f "deviation_simulation_kpc.png" ];
then
    rm deviation_simulation_kpc.png
fi
if [ -f "deviation_simulation_Mpc.png" ];
then
    rm deviation_simulation_Mpc.png
fi

if [ -f "deviation_from_original_data_simulation_kpc.png" ];
then
    rm deviation_from_original_data_simulation_kpc.png
fi

if [ -f "deviation_from_original_data_simulation_Mpc.png" ];
then
    rm deviation_from_original_data_simulation_Mpc.png
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 the circular orbits and of the errors"
if command -v python3 &>/dev/null; then
    python3 makePlots.py
else 
    python3 makePlots.py
fi