void eraPvtob	(	double	elong,
		double	phi,
		double	hm,
		double	xp,
		double	yp,
		double	sp,
		double	theta,
		double	pv[2][3]
	)

Definition at line 3 of file pvtob.c.

References eraGd2gc(), eraPom00(), eraTrxp(), ERFA_D2PI, and ERFA_DAYSEC.

Referenced by eraApco(), eraApio(), and t_pvtob().

{
/* Earth rotation rate in radians per UT1 second */
   const double OM = 1.00273781191135448 * ERFA_D2PI / ERFA_DAYSEC;

   double xyzm[3], rpm[3][3], xyz[3], x, y, z, s, c;

/* Geodetic to geocentric transformation (ERFA_WGS84). */
   (void) eraGd2gc(1, elong, phi, hm, xyzm);

/* Polar motion and TIO position. */
   eraPom00(xp, yp, sp, rpm);
   eraTrxp(rpm, xyzm, xyz);
   x = xyz[0];
   y = xyz[1];
   z = xyz[2];

/* Functions of ERA. */
   s = sin(theta);
   c = cos(theta);

/* Position. */
   pv[0][0] = c*x - s*y;
   pv[0][1] = s*x + c*y;
   pv[0][2] = z;

/* Velocity. */
   pv[1][0] = OM * ( -s*x - c*y );
   pv[1][1] = OM * (  c*x - s*y );
   pv[1][2] = 0.0;

/* Finished. */

}

Here is the call graph for this function:

Here is the caller graph for this function: