id	summary	reporter	owner	description	type	status	priority	milestone	component	version	resolution	keywords	cc
323	Rewrite RoC calculation	Ian Culverwell	Ian Culverwell	"The calculation of radius of curvature in '''ropp_utils/coordinates/curvature.f90''' is a bit dodgy: Ak2 is indeterminate (0/0) at lat = pi/2, and Ak1 includes such huge quantities (like ((SQRT(Re**4))**3) that it causes a floating point overflow if evaluated at single precision!

Suggest rewriting in terms of the principal radii R_NS and R_EW.

(Prompted by discussions with Josep Aparicio.)


"	defect	closed	normal	9.0	ropp_utils	6.1	fixed	RoC