Merge pull request #1324 from kbevers/molodensky-eccentricity

Use 1st eccentricity instead of 2nd eccentricity
OSGeo · Mar 16, 2019 · eab2279 · eab2279
2 parents a9886a5 + 644a592
commit eab2279
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/transformations/molodensky.cpp b/src/transformations/molodensky.cpp
@@ -140,7 +140,7 @@ static PJ_LPZ calc_standard_params(PJ_LPZ lpz, PJ *P) {
 
     /* ellipsoid radii of curvature */
     double rho = RM(a, P->es, lpz.phi);
-    double nu  = RN(a, P->e2s, lpz.phi);
+    double nu  = RN(a, P->es, lpz.phi);
 
     /* delta phi */
     dphi  = (-dx*sphi*clam) - (dy*sphi*slam) + (dz*cphi)
@@ -193,7 +193,7 @@ static PJ_LPZ calc_abridged_params(PJ_LPZ lpz, PJ *P) {
 
     /* delta lambda */
     dlam = -dx*slam + dy*clam;
-    const double dlam_denom = RN(P->a, P->e2s, lpz.phi)*cphi;
+    const double dlam_denom = RN(P->a, P->es, lpz.phi)*cphi;
     if( dlam_denom == 0.0 ) {
         lpz.lam = HUGE_VAL;
         return lpz;