Implemented periodic boundary conditions for parallel simulation

Modified generation of halo atoms to include replication of system generating periodic boundaries. Initialisation of total number of super cells has been necessarily moved to master create function. Neighbourlist calculation has been modified to accept negative supercell indices. Tidied up tabs.
Implemented macrocell data output including position, magnetisation and demagnetisation field data.