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 r0, a;
 6 28994509 2004-04-21 devnull 
 7 28994509 2004-04-21 devnull static int
 8 28994509 2004-04-21 devnull Xsimpleconic(struct place *place, double *x, double *y)
 9 28994509 2004-04-21 devnull {
10 28994509 2004-04-21 devnull 	double r = r0 - place->nlat.l;
11 28994509 2004-04-21 devnull 	double t = a*place->wlon.l;
12 28994509 2004-04-21 devnull 	*x = -r*sin(t);
13 28994509 2004-04-21 devnull 	*y = -r*cos(t);
14 28994509 2004-04-21 devnull 	return 1;
15 28994509 2004-04-21 devnull }
16 28994509 2004-04-21 devnull 
17 28994509 2004-04-21 devnull proj
18 28994509 2004-04-21 devnull simpleconic(double par0, double par1)
19 28994509 2004-04-21 devnull {
20 28994509 2004-04-21 devnull 	struct coord lat0;
21 28994509 2004-04-21 devnull 	struct coord lat1;
22 28994509 2004-04-21 devnull 	deg2rad(par0,&lat0);
23 28994509 2004-04-21 devnull 	deg2rad(par1,&lat1);
24 28994509 2004-04-21 devnull 	if(fabs(lat0.l+lat1.l)<.01)
25 28994509 2004-04-21 devnull 		return rectangular(par0);
26 28994509 2004-04-21 devnull 	if(fabs(lat0.l-lat1.l)<.01) {
27 28994509 2004-04-21 devnull 		a = lat0.s/lat0.l;
28 28994509 2004-04-21 devnull 		r0 = lat0.c/lat0.s + lat0.l;
29 28994509 2004-04-21 devnull 	} else {
30 28994509 2004-04-21 devnull 		a = (lat1.c-lat0.c)/(lat0.l-lat1.l);
31 28994509 2004-04-21 devnull 		r0 = ((lat0.c+lat1.c)/a + lat1.l + lat0.l)/2;
32 28994509 2004-04-21 devnull 	}
33 28994509 2004-04-21 devnull 	return Xsimpleconic;
34 28994509 2004-04-21 devnull }