void palMappa	(	double	eq,
		double	date,
		double	amprms[21]
	)

Definition at line 94 of file palMappa.c.

References eraEpj(), eraPn(), i, PAL__CR, PAL__MJD0, palEvp(), and palPrenut().

Referenced by PointingModel::CalcPointingPos(), palAmp(), palMap(), t_ampqk(), t_mappa(), and t_mapqkz().

                                                          {

/* Local constants */

/*  Gravitational radius of the Sun x 2 (2*mu/c**2, AU) */
  const double GR2 = 2.0 * 9.87063e-9;

/* Local Variables; */
   int i;
   double ebd[ 3 ], ehd[ 3 ], eh[ 3 ], e, vn[ 3 ], vm;

/* Initialise so that unsused values are returned holding zero */
   memset( amprms, 0, 21*sizeof( *amprms ) );

/* Time interval for proper motion correction. */
   amprms[ 0 ] = eraEpj( PAL__MJD0, date ) - eq;

/* Get Earth barycentric and heliocentric position and velocity. */
   palEvp( date, eq, ebd, &amprms[ 1 ], ehd, eh );

/* Heliocentric direction of Earth (normalized) and modulus. */
   eraPn( eh, &e, &amprms[ 4 ] );

/* Light deflection parameter */
   amprms[7] = GR2 / e;

/* Aberration parameters. */
   for( i = 0; i < 3; i++ ) {
      amprms[ i + 8 ] = ebd[ i ]*PAL__CR;
   }
   eraPn( &amprms[8], &vm, vn );
   amprms[ 11 ] = sqrt( 1.0 - vm*vm );

/* NPB matrix. */
   palPrenut( eq, date, (double(*)[ 3 ]) &amprms[ 12 ] );
}

Here is the call graph for this function:

Here is the caller graph for this function: