13e12c5d1SDavid du Colombier #include "map.h" 23e12c5d1SDavid du Colombier 33e12c5d1SDavid du Colombier struct coord center; 43e12c5d1SDavid du Colombier 53e12c5d1SDavid du Colombier static int 6*219b2ee8SDavid du Colombier Xelliptic(struct place *place, double *x, double *y) 73e12c5d1SDavid du Colombier { 8*219b2ee8SDavid du Colombier double r1,r2; 93e12c5d1SDavid du Colombier r1 = acos(place->nlat.c*(place->wlon.c*center.c 103e12c5d1SDavid du Colombier - place->wlon.s*center.s)); 113e12c5d1SDavid du Colombier r2 = acos(place->nlat.c*(place->wlon.c*center.c 123e12c5d1SDavid du Colombier + place->wlon.s*center.s)); 133e12c5d1SDavid du Colombier *x = -(r1*r1 - r2*r2)/(4*center.l); 143e12c5d1SDavid du Colombier *y = (r1*r1+r2*r2)/2 - (center.l*center.l+*x**x); 153e12c5d1SDavid du Colombier if(*y < 0) 163e12c5d1SDavid du Colombier *y = 0; 173e12c5d1SDavid du Colombier *y = sqrt(*y); 183e12c5d1SDavid du Colombier if(place->nlat.l<0) 193e12c5d1SDavid du Colombier *y = -*y; 203e12c5d1SDavid du Colombier return(1); 213e12c5d1SDavid du Colombier } 223e12c5d1SDavid du Colombier 233e12c5d1SDavid du Colombier proj 24*219b2ee8SDavid du Colombier elliptic(double l) 253e12c5d1SDavid du Colombier { 263e12c5d1SDavid du Colombier l = fabs(l); 273e12c5d1SDavid du Colombier if(l>89) 283e12c5d1SDavid du Colombier return(0); 293e12c5d1SDavid du Colombier if(l<1) 303e12c5d1SDavid du Colombier return(Xazequidistant); 313e12c5d1SDavid du Colombier deg2rad(l,¢er); 323e12c5d1SDavid du Colombier return(Xelliptic); 333e12c5d1SDavid du Colombier } 34