83 double palPa(
double ha,
double dec,
double phi ) {
84 double cp, sqsz, cqsz;
88 cqsz = sin(phi) * cos(dec) - cp * sin(dec) * cos(ha);
89 if (sqsz == 0.0 && cqsz == 0.0) cqsz = 1.0;
90 return atan2( sqsz, cqsz );
double palPa(double ha, double dec, double phi)