view release on metacpan or  search on metacpan

lib/Astro/PAL.pm  view on Meta::CPAN

themselves.

=head3 Anomalies

=over 4

=item palObs

palObs is special in that it returns an empty list if the return
status is bad. Additionally, palObs is called with a single
argument and the behaviour depends on whether the argument looks
like an integer or a string.

 ($ident, $name, $w, $p, $h) = palObs( 27 );
 ($ident, $name, $w, $p, $h) = palObs( "JCMT" );

=cut

# We need a local version of palObs that converts the single
# argument to the right form for the internal XS implementation

palsrc/palGeoc.c  view on Meta::CPAN

*-
*/

#include "pal.h"
#include "pal1sofa.h"

void palGeoc ( double p, double h, double *r, double *z ) {
  double xyz[3];
  const double elong = 0.0;   /* Use zero longitude */
  const double AU = 1.49597870E11;
  /* WGS84 looks to be the closest match */
  eraGd2gc( ERFA_WGS84, elong, p, h, xyz );
  *r = xyz[0] / (AU * cos(elong) );
  *z = xyz[2] / AU;
}