Help with magnetostatic parameters

[Peela baradwaj naidu]

Hello VAMPIRE Users, 
I need some help with magnetosatic parameters in my simulation. I am studying ultrafast demagnetisation under applied field and temperature pulses and for some reason when i switch off dipole solvers and apply a field strength of 3.5 T, I egt a demagnetisation behaviour that is fairly expected based on my experiments. However the real problem arises when I switch on the dipole solvers (commented out in the code). Firstly my magnetisation before starting applied field doesn't look equilibrated compared to my result when dipole parameters are off. I also looked at the magnetostatic energy output to see if it could provide any additional info. My goal is to reduce the applied required for similar demagnetisation to closer to <0.9 T. For this reason I simply cannot disregard the contribution of magnetostaic energy, however including magnetosatatics is causing further problems. Can anyone please advise how can I improve my dipole parameters?
Thank You
I cannot upload files for some reason so here is a drive link to the files necessary link