From 2cac138a22e8f01e34e9cfafb6b526694548b269 Mon Sep 17 00:00:00 2001
From: Kristian Evers <kristianevers@gmail.com>
Date: Wed, 19 Oct 2016 18:21:04 +0200
Subject: Incorrect formula used, changed to comply with Snyder and other
 sources. Fixes #331

---
 src/PJ_omerc.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'src')

diff --git a/src/PJ_omerc.c b/src/PJ_omerc.c
index f0a660bd..6abd2d33 100644
--- a/src/PJ_omerc.c
+++ b/src/PJ_omerc.c
@@ -219,7 +219,7 @@ PJ *PROJECTION(omerc) {
     if (no_off)
         Q->u_0 = 0;
     else {
-        Q->u_0 = fabs(Q->ArB * atan2(sqrt(D * D - 1.), cos(alpha_c)));
+        Q->u_0 = fabs(Q->ArB * atan(sqrt(D * D - 1.) / cos(alpha_c)));
         if (P->phi0 < 0.)
             Q->u_0 = - Q->u_0;
     }
-- 
cgit v1.2.3