From 28994509cc11ac6a5443054dfae1fedfb69039bc Mon Sep 17 00:00:00 2001 From: rsc Date: Wed, 21 Apr 2004 22:19:33 +0000 Subject: Why not? --- src/cmd/map/libmap/mollweide.c | 25 +++++++++++++++++++++++++ 1 file changed, 25 insertions(+) create mode 100644 src/cmd/map/libmap/mollweide.c (limited to 'src/cmd/map/libmap/mollweide.c') diff --git a/src/cmd/map/libmap/mollweide.c b/src/cmd/map/libmap/mollweide.c new file mode 100644 index 00000000..3284c495 --- /dev/null +++ b/src/cmd/map/libmap/mollweide.c @@ -0,0 +1,25 @@ +#include +#include +#include "map.h" + +static int +Xmollweide(struct place *place, double *x, double *y) +{ + double z; + double w; + z = place->nlat.l; + if(fabs(z)<89.9*RAD) + do { /*newton for 2z+sin2z=pi*sin(lat)*/ + w = (2*z+sin(2*z)-PI*place->nlat.s)/(2+2*cos(2*z)); + z -= w; + } while(fabs(w)>=.00001); + *y = sin(z); + *x = - (2/PI)*cos(z)*place->wlon.l; + return(1); +} + +proj +mollweide(void) +{ + return(Xmollweide); +} -- cgit v1.2.3