1 28994509 2004-04-21 devnull #include <u.h>
 2 28994509 2004-04-21 devnull #include <libc.h>
 3 28994509 2004-04-21 devnull #include "map.h"
 4 28994509 2004-04-21 devnull 
 5 28994509 2004-04-21 devnull static double a;
 6 28994509 2004-04-21 devnull 
 7 28994509 2004-04-21 devnull static int
 8 28994509 2004-04-21 devnull Xnewyorker(struct place *place, double *x, double *y)
 9 28994509 2004-04-21 devnull {
10 28994509 2004-04-21 devnull 	double r = PI/2 - place->nlat.l;
11 28994509 2004-04-21 devnull 	double s;
12 28994509 2004-04-21 devnull 	if(r<.001)	/* cheat to plot center */
13 28994509 2004-04-21 devnull 		s = 0;
14 28994509 2004-04-21 devnull 	else if(r<a)
15 28994509 2004-04-21 devnull 		return -1;
16 28994509 2004-04-21 devnull 	else
17 28994509 2004-04-21 devnull 		s = log(r/a);
18 28994509 2004-04-21 devnull 	*x = -s * place->wlon.s;
19 28994509 2004-04-21 devnull 	*y = -s * place->wlon.c;
20 28994509 2004-04-21 devnull 	return(1);
21 28994509 2004-04-21 devnull }
22 28994509 2004-04-21 devnull 
23 28994509 2004-04-21 devnull proj
24 28994509 2004-04-21 devnull newyorker(double a0)
25 28994509 2004-04-21 devnull {
26 28994509 2004-04-21 devnull 	a = a0*RAD;
27 28994509 2004-04-21 devnull 	return(Xnewyorker);
28 28994509 2004-04-21 devnull }