If at least one of the input arguments x and y is of type Matrix, then all output arguments are of type Matrix, with the dimensions of the first input argument of type Matrix. If one input argument is of type Number and the other of type Matrix, then the Number type argument is "expanded" to become a Matrix filled uniformly with the value of the Number.
The coordinates system expected by MAGNETIC_FIELD is the internal frame. If you describe your cell in Cartesian coordinates, then the internal frame is simply the Cartesian coordinate system. If use polar coordinates however, then you should supply the following parameters:
for x: log(r) for y: \φ, in radians
These mappings can be performed with the POLAR_TO_INTERNAL procedure.
The internal coordinates system is Cartesian when the cell is described in Cartesian coordinates, but when you describe the cell in polar coordinates, you've to apply the following mapping:
Br = Bx/r Bphi = By/r Bz = Bz
The internal units for the magnetic field are 100\ G or 0.01\ T.
Formatted on 21/01/18 at 16:55.