Improvements to parallel inversion code

Can now cope with open field-lines, where zero-gradient BCs are currently used.

parderiv can now include additional terms for Z derivatives, mixed Y-Z derivatives,
and first-order Y derivs

SurfaceIter extended to test whether processor is on a boundary