Skip to content
Fetching contributors…
Cannot retrieve contributors at this time
47 lines (46 sloc) 1.03 KB
#define PROJ_PARMS__ \
struct PJconsts *sinu; \
struct PJconsts *moll;
#define PJ_LIB__
#include <projects.h>
PROJ_HEAD(goode, "Goode Homolosine") "\n\tPCyl, Sph.";
C_NAMESPACE PJ
*pj_sinu(PJ *), *pj_moll(PJ *);
#define Y_COR 0.05280
#define PHI_LIM .71093078197902358062
FORWARD(s_forward); /* spheroid */
if (fabs(lp.phi) <= PHI_LIM)
xy = P->sinu->fwd(lp, P->sinu);
else {
xy = P->moll->fwd(lp, P->moll);
xy.y -= lp.phi >= 0.0 ? Y_COR : -Y_COR;
}
return (xy);
}
INVERSE(s_inverse); /* spheroid */
if (fabs(xy.y) <= PHI_LIM)
lp = P->sinu->inv(xy, P->sinu);
else {
xy.y += xy.y >= 0.0 ? Y_COR : -Y_COR;
lp = P->moll->inv(xy, P->moll);
}
return (lp);
}
FREEUP;
if (P) {
if (P->sinu)
(*(P->sinu->pfree))(P->sinu);
if (P->moll)
(*(P->moll->pfree))(P->moll);
pj_dalloc(P);
}
}
ENTRY2(goode, sinu, moll)
P->es = 0.;
if (!(P->sinu = pj_sinu(0)) || !(P->moll = pj_moll(0)))
E_ERROR_0;
if (!(P->sinu = pj_sinu(P->sinu)) || !(P->moll = pj_moll(P->moll)))
E_ERROR_0;
P->fwd = s_forward;
P->inv = s_inverse;
ENDENTRY(P)
Jump to Line
Something went wrong with that request. Please try again.