add cartesian to poloidal comparison subroutine
add subroutine to compute both main and aux fields
Implement field line following, still have to test for accuracy