void palPrebn	(	double	bep0,
		double	bep1,
		double	rmatp[3][3]
	)

Definition at line 76 of file palPrebn.c.

References PAL__DAS2R, palDeuler(), t, and telData::w.

Referenced by palPreces(), and t_prebn().

                                                               {

  double t,bigt, zeta, theta, z, tas2r, w;

  /* Interval between basic epoch B1850.0 and beginning epoch in TC */
  bigt = (bep0-1850)/100.;

  /*  Interval over which precession required, in tropical centuries */
  t = (bep1-bep0)/100.;

  /* Euler angles */
  tas2r = t * PAL__DAS2R;
  w = 2303.5548 + ( 1.39720 + 0.000059 * bigt) * bigt;

  zeta = ( w + ( 0.30242 - 0.000269 * bigt + 0.017996 * t ) * t ) * tas2r;
  z = ( w + ( 1.09478 + 0.000387 * bigt + 0.018324 * t ) * t ) * tas2r;
  theta = ( 2005.1125 + ( -0.85294 - 0.000365 * bigt ) * bigt +
            (-0.42647 - 0.000365 * bigt - 0.041802 * t ) * t ) * tas2r;

  /*  Rotation matrix */
  palDeuler("ZYZ", -zeta, theta, -z, rmatp);

}

Here is the call graph for this function:

Here is the caller graph for this function: