Scaled atomistic VTK magnetic moments data in units of Bohr magnetons

Author: davidcorteso

fidimag/atomistic/llg.py

@@ -326,7 +326,18 @@ def add_monitor_at(self, i, j, k, name='p'):
         self.saver.update_entity_order()
 
     def save_vtk(self):
-        self.vtk.save_vtk(self.spin.reshape(-1, 3), self._mu_s, step=self.step)
+        """
+        Save a VTK file with the magnetisation vector field and magnetic
+        moments as cell data. Magnetic moments are saved in units of
+        Bohr magnetons
+
+        NOTE: It is recommended to use a *cell to point data* filter in
+        Paraview or Mayavi to plot the vector field
+        """
+        self.vtk.save_vtk(self.spin.reshape(-1, 3),
+                          self._mu_s / const.mu_B,
+                          step=self.step
+                          )
 
     def save_m(self):
         if not os.path.exists('%s_npys' % self.name):