Return polar radius based on projection

The polar radius is different if we are using a spherical projection
so return the value based on the projection being used.
Merge pull request #41 from cjh1/elliptical-vs-spherical-fix

Return polar radius based on projection