135 void palPolmo (
double elongm,
double phim,
double xp,
double yp,
136 double *elong,
double *phi,
double *daz ) {
138 double sel,cel,sph,cph,xm,ym,zm,xnm,ynm,znm,
139 sxp,cxp,syp,cyp,zw,xt,yt,zt,xnt,ynt;
168 cph=sqrt(xt*xt+yt*yt);
169 if (cph == 0.0) xt=1.0;
174 if (xt != 0.0 || yt != 0.0) {
182 xnt=(xnm*cel+ynm*sel)*zt-znm*cph;
183 ynt=-xnm*sel+ynm*cel;
184 if (xnt != 0.0 || ynt != 0.0) {
185 *daz=atan2(-ynt,-xnt);
void palPolmo(double elongm, double phim, double xp, double yp, double *elong, double *phi, double *daz)