void palAoppa	(	double	date,
		double	dut,
		double	elongm,
		double	phim,
		double	hm,
		double	xp,
		double	yp,
		double	tdk,
		double	pmb,
		double	rh,
		double	wl,
		double	tlr,
		double	aoprms[14]
	)

Definition at line 172 of file palAoppa.c.

References C, PAL__D2PI, PAL__DS2R, pal__Eqeqx(), pal__Geoc(), palAoppat(), and palRefco().

Referenced by PointingModel::CalcPointingPos(), palAop(), palOap(), and t_aop().

                                                                      {

  /* Constants */
  const double C = 173.14463331; /* Speed of light (AU per day) */
  const double SOLSID = 1.0027379093; /* Ratio between solar and sidereal time */

  /* Local variables */
  double cphim,xt,yt,zt,xc,yc,zc,elong,phi,uau,vau;

  /*  Observer's location corrected for polar motion */
  cphim = cos(phim);
  xt = cos(elongm)*cphim;
  yt = sin(elongm)*cphim;
  zt = sin(phim);
  xc = xt-xp*zt;
  yc = yt+yp*zt;
  zc = xp*xt-yp*yt+zt;
  if (xc == 0.0 && yc == 0.0) {
    elong = 0.0;
  } else {
    elong = atan2(yc,xc);
  }
  phi = atan2(zc,sqrt(xc*xc+yc*yc));
  aoprms[0] = phi;
  aoprms[1] = sin(phi);
  aoprms[2] = cos(phi);

  /*  magnitude of the diurnal aberration vector */
  pal__Geoc(phi,hm,&uau,&vau);
  aoprms[3] = PAL__D2PI*uau*SOLSID/C;

  /*  copy the refraction parameters and compute the a & b constants */
  aoprms[4] = hm;
  aoprms[5] = tdk;
  aoprms[6] = pmb;
  aoprms[7] = rh;
  aoprms[8] = wl;
  aoprms[9] = tlr;
  palRefco(hm,tdk,pmb,rh,wl,phi,tlr,1e-10,
           &aoprms[10],&aoprms[11]);

  /*  longitude + equation of the equinoxes + sidereal equivalent of DUT
   *  (ignoring change in equation of the equinoxes between UTC and TDB) */
  aoprms[12] = elong+pal__Eqeqx(date)+dut*SOLSID*PAL__DS2R;

  /*  sidereal time */
  palAoppat(date,aoprms);

}

Here is the call graph for this function:

Here is the caller graph for this function: