static void pal__Geoc ( double  p, double  h, double *  r, double *  z  )

static

Definition at line 228 of file palAoppa.c.

Referenced by palAoppa().

                                                                   {
  /*  earth equatorial radius (metres) */
  const double A0=6378140.0;

  /*  reference spheroid flattening factor and useful function */
  const double f = 1.0/298.257;
  double b;

  /*  astronomical unit in metres */
  const double AU = 1.49597870e11;

  double sp,cp,c,s;

  b = pow( 1.0-f, 2.0 );

  /*  geodetic to geocentric conversion */
  sp = sin(p);
  cp = cos(p);
  c = 1.0/sqrt(cp*cp+b*sp*sp);
  s = b*c;
  *r = (A0*c+h)*cp/AU;
  *z = (A0*s+h)*sp/AU;

}

Here is the caller graph for this function: