To run the program simply type:
remanence_anisotropy_magic.py 

