Fix a mistake in the spherical polar line element formula
Set up toroidal field in spherical coordinates with axisymmetry