/Proj4/PJ_goode.c
http://github.com/route-me/route-me · C · 49 lines · 49 code · 0 blank · 0 comment · 11 complexity · 723ddbf7acdfde8ff6bd99524c7cce07 MD5 · raw file
- #ifndef lint
- static const char SCCSID[]="@(#)PJ_goode.c 4.1 94/02/15 GIE REL";
- #endif
- #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)