101 void palRdplan(
double date,
int np,
double elong,
double phi,
102 double * ra,
double * dec,
double * diam ) {
105 const double AUKM = 1.49597870e8;
108 const double EQRAU[] = {
130 double dx,dy,dz,r,tl;
133 if (np < 0 || np > 8 ) np=0;
144 eraRxp( rmat, &(v[3]), &(vgm[3]) );
150 for (i=0; i<6; i++) {
164 eraRxp( rmat, &(v[3]), &(vse[3]) );
167 for (i=0; i<6; i++) {
168 vsg[
i] = vse[
i] - 0.012150581 * vgm[
i];
175 for (i=0; i<6; i++) {
186 eraRxp( rmat, &(v[3]), &(vsp[3]) );
189 for (i=0; i<6; i++) {
190 v[
i] = vsp[
i] - vsg[
i];
199 for (i=0; i<6; i++) {
207 r = sqrt( dx*dx + dy*dy + dz*dz );
213 for (i=0; i<3; i++) {
222 *diam = 2.0 * asin( EQRAU[np] / (r * AUKM ) );
static const double PAL__CR
void palNut(double date, double rmatn[3][3])
void palPrenut(double epoch, double date, double rmatpn[3][3])
double palEpj(double date)
double palDt(double epoch)
void palDmoon(double date, double pv[6])
void eraC2s(double p[3], double *theta, double *phi)
void eraRxp(double r[3][3], double p[3], double rp[3])
void palPvobs(double p, double h, double stl, double pv[6])
void palPlanet(double date, int np, double pv[6], int *j)
void palRdplan(double date, int np, double elong, double phi, double *ra, double *dec, double *diam)
double palGmst(double ut1)