mpi_rmdir
store_step
mpi_wrap_master_thread
sorts.simulation.
Wrap function to only execute on thread rank=0