  We present a method for performing atomistic spin dynamic simulations. A
comprehensive summary of all pertinent details for performing the simulations
such as equations of motions, models for including temperature, methods of
extracting data and numerical schemes for performing the simulations is given.
The method can be applied in a first principles mode, where all interatomic
exchange is calculated self-consistently, or it can be applied with frozen
parameters estimated from experiments or calculated for a fixed
spin-configuration. Areas of potential applications to different magnetic
questions are also discussed. The method is finally applied to one situation
where the macrospin model breaks down; magnetic switching in ultra strong
fields.
