If you want it at T = 0 K then you could use sim:program = static-hysteresis-loop instead of sim:program = hysteresis-loop. If it's not at 0 K then you have to do the dynamic simulation and it might look like the material has much higher coersivity than it should simply because the sweep is too fast. One trick is to set critical damping, which you already have (alpha = 1.0), and another thing to do is to just increase the simulation time, i.e. the loop time steps.
I have tried with just Py but the simulation time became too long before I could get reasonable coercivity. There is an article about it where they use "kinetic Monte Carlo" instead of llg:
Kinetic Monte Carlo simulations of M-H Loops for HAMR recording media: Comparison with MOKE data, Plumer M., et. al.,
https://www.researchgate.net/publication/263022496_Kinetic_Monte_Carlo_simulations_of_M-H_Loops_for_HAMR_recording_media_Comparison_with_MOKE_data .
If you find a good solution in Vampire I would appreciate it if you shared it.