Hi,
I am simulating a system with circular 
exchange bias by taking the field generated after relaxing the structure
 at some defined points of the mesh using:
  h=sim.get_subfield('H_total_Co', SI(1, 'A/m'))
After treating the numpy array "h" I plug it to the system using: 
  sim.set_H_ext(h, unit=SI('A/m'))
After
 doing so, the numpy "h" gets multiplied by 1e-6, why is that? It does 
not happen when I use a vector or a function instead of a numpy array...
thanks!!!