diff options
Diffstat (limited to 'src/PJ_lask.c')
| -rw-r--r-- | src/PJ_lask.c | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/src/PJ_lask.c b/src/PJ_lask.c new file mode 100644 index 00000000..34a61d94 --- /dev/null +++ b/src/PJ_lask.c @@ -0,0 +1,29 @@ +#ifndef lint +static const char SCCSID[]="@(#)PJ_lask.c 4.1 94/02/15 GIE REL"; +#endif +#define PROJ_PARMS__ +#define PJ_LIB__ +#include <projects.h> +PROJ_HEAD(lask, "Laskowski") "\n\tMisc Sph, no inv."; +#define a10 0.975534 +#define a12 -0.119161 +#define a32 -0.0143059 +#define a14 -0.0547009 +#define b01 1.00384 +#define b21 0.0802894 +#define b03 0.0998909 +#define b41 0.000199025 +#define b23 -0.0285500 +#define b05 -0.0491032 +FORWARD(s_forward); /* sphere */ + double l2, p2; + + l2 = lp.lam * lp.lam; + p2 = lp.phi * lp.phi; + xy.x = lp.lam * (a10 + p2 * (a12 + l2 * a32 + p2 * a14)); + xy.y = lp.phi * (b01 + l2 * (b21 + p2 * b23 + l2 * b41) + + p2 * (b03 + p2 * b05)); + return (xy); +} +FREEUP; if (P) pj_dalloc(P); } +ENTRY0(lask) P->fwd = s_forward; P->inv = 0; P->es = 0.; ENDENTRY(P) |