Index: lwspheroid.c
===================================================================
--- lwspheroid.c	(revision 5178)
+++ lwspheroid.c	(working copy)
@@ -315,6 +315,7 @@
 	double oneOver2e = 1.0 / (2.0 * e);
 	double t2p1 = oneOver2e * log((1.0 + e * sinPhi1) / (1.0 - e * sinPhi1));
 	double t2p2 = oneOver2e * log((1.0 + e * sinPhi2) / (1.0 - e * sinPhi2));
+        char s[20]; sprintf(s, "%.12g", e); 
 	return z0 * (t1p2 + t2p2) - z0 * (t1p1 + t2p1);
 }