From 5d39ae0bb01d142a69552ce3553dd61b75db4ae8 Mon Sep 17 00:00:00 2001 From: Phil Elson Date: Fri, 21 Dec 2018 14:04:04 +0000 Subject: Fixed rHealpix projection, closing #1206. Numerical test results based on function output, rather than mathematical derivation. Verified global coverage with graphical eyeballing through cartopy. --- src/PJ_healpix.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'src') diff --git a/src/PJ_healpix.c b/src/PJ_healpix.c index ee575abf..3a88695f 100644 --- a/src/PJ_healpix.c +++ b/src/PJ_healpix.c @@ -499,8 +499,8 @@ static XY combine_caps(double x, double y, int north_square, int south_square, double a[2]; /* Workaround cppcheck git issue */ double* pa = a; - pa[0] = -3*M_FORTPI + ((inverse == 0) ? 0 : capmap.cn) *M_HALFPI; - pa[1] = M_HALFPI; + pa[0] = -3*M_FORTPI + ((inverse == 0) ? pole : capmap.cn) *M_HALFPI; + pa[1] = ((capmap.region == north) ? 1 : -1) *M_HALFPI; vector_add(ret_dot, a, vector); } -- cgit v1.2.3