From 565a4bd035b9d4a83955808efef20f1d8dfa24cf Mon Sep 17 00:00:00 2001
From: Frank Warmerdam <warmerdam@pobox.com>
Date: Thu, 18 Mar 1999 16:34:52 +0000
Subject: New

git-svn-id: http://svn.osgeo.org/metacrs/proj/trunk@776 4e78687f-474d-0410-85f9-8d5e500ac6b2
---
 src/PJ_wag2.c | 24 ++++++++++++++++++++++++
 1 file changed, 24 insertions(+)
 create mode 100644 src/PJ_wag2.c

(limited to 'src/PJ_wag2.c')

diff --git a/src/PJ_wag2.c b/src/PJ_wag2.c
new file mode 100644
index 00000000..d6a0e3f2
--- /dev/null
+++ b/src/PJ_wag2.c
@@ -0,0 +1,24 @@
+#ifndef lint
+static const char SCCSID[]="@(#)PJ_wag2.c	4.1	94/02/15	GIE	REL";
+#endif
+#define PJ_LIB__
+# include	<projects.h>
+PROJ_HEAD(wag2, "Wagner II") "\n\tPCyl., Sph.";
+#define C_x 0.92483
+#define C_y 1.38725
+#define C_p1 0.88022
+#define C_p2 0.88550
+FORWARD(s_forward); /* spheroid */
+	lp.phi = aasin(C_p1 * sin(C_p2 * lp.phi));
+	xy.x = C_x * lp.lam * cos(lp.phi);
+	xy.y = C_y * lp.phi;
+	return (xy);
+}
+INVERSE(s_inverse); /* spheroid */
+	lp.phi = xy.y / C_y;
+	lp.lam = xy.x / (C_x * cos(lp.phi));
+	lp.phi = aasin(sin(lp.phi) / C_p1) / C_p2;
+	return (lp);
+}
+FREEUP; if (P) pj_dalloc(P); }
+ENTRY0(wag2) P->es = 0.; P->inv = s_inverse; P->fwd = s_forward; ENDENTRY(P)
-- 
cgit v1.2.3