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 struct coord stdpar0, stdpar1;
 6 28994509 2004-04-21 devnull static double k;
 7 28994509 2004-04-21 devnull static double yeq;
 8 28994509 2004-04-21 devnull 
 9 28994509 2004-04-21 devnull static int
10 28994509 2004-04-21 devnull Xtrapezoidal(struct place *place, double *x, double *y)
11 28994509 2004-04-21 devnull {
12 28994509 2004-04-21 devnull 	*y = yeq + place->nlat.l;
13 28994509 2004-04-21 devnull 	*x = *y*k*place->wlon.l;
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 trapezoidal(double par0, double par1)
19 28994509 2004-04-21 devnull {
20 28994509 2004-04-21 devnull 	if(fabs(fabs(par0)-fabs(par1))<.1)
21 28994509 2004-04-21 devnull 		return rectangular(par0);
22 28994509 2004-04-21 devnull 	deg2rad(par0,&stdpar0);
23 28994509 2004-04-21 devnull 	deg2rad(par1,&stdpar1);
24 28994509 2004-04-21 devnull 	if(fabs(par1-par0) < .1)
25 28994509 2004-04-21 devnull 		k = stdpar1.s;
26 28994509 2004-04-21 devnull 	else
27 28994509 2004-04-21 devnull 		k = (stdpar1.c-stdpar0.c)/(stdpar0.l-stdpar1.l);
28 28994509 2004-04-21 devnull 	yeq = -stdpar1.l - stdpar1.c/k;
29 28994509 2004-04-21 devnull 	return Xtrapezoidal;
30 28994509 2004-04-21 devnull }