5 static struct coord stdp0, stdp1;
9 Xlambert(struct place *place, double *x, double *y)
12 if(place->nlat.l < -80.*RAD)
14 if(place->nlat.l > 89.*RAD)
17 r = stdp0.c*exp(0.5*k*log(
18 (1+stdp0.s)*(1-place->nlat.s)/((1-stdp0.s)*(1+place->nlat.s))));
21 *x = - r*sin(k * place->wlon.l);
22 *y = - r*cos(k * place->wlon.l);
27 lambert(double par0, double par1)
30 if(fabs(par0)>fabs(par1)){
35 deg2rad(par0, &stdp0);
36 deg2rad(par1, &stdp1);
37 if(fabs(par1+par0)<.1)
39 if(fabs(par1-par0)<.1)
40 return(perspective(-1.));
41 if(fabs(par0)>89.5||fabs(par1)>89.5)
43 k = 2*log(stdp1.c/stdp0.c)/log(
44 (1+stdp0.s)*(1-stdp1.s)/((1-stdp0.s)*(1+stdp1.s)));