Generalised Spherical coordinates into own function; test routines.
Fixed theta output routine (was relying on old temporary variables)