FACT++  1.0
void palGaleq ( double  dl,
double  db,
double *  dr,
double *  dd 
)

Definition at line 80 of file palGaleq.c.

References eraAnp(), eraAnpm(), eraC2s(), eraS2c(), and eraTrxp().

Referenced by t_galeq().

80  {
81 
82  double v1[3];
83  double v2[3];
84 
85 /*
86 * L2,B2 system of galactic coordinates
87 *
88 * P = 192.25 RA of galactic north pole (mean B1950.0)
89 * Q = 62.6 inclination of galactic to mean B1950.0 equator
90 * R = 33 longitude of ascending node
91 *
92 * P,Q,R are degrees
93 *
94 * Equatorial to galactic rotation matrix (J2000.0), obtained by
95 * applying the standard FK4 to FK5 transformation, for zero proper
96 * motion in FK5, to the columns of the B1950 equatorial to
97 * galactic rotation matrix:
98 */
99  double rmat[3][3] = {
100  { -0.054875539726,-0.873437108010,-0.483834985808 },
101  { +0.494109453312,-0.444829589425,+0.746982251810 },
102  { -0.867666135858,-0.198076386122,+0.455983795705 }
103  };
104 
105  /* Spherical to Cartesian */
106  eraS2c( dl, db, v1 );
107 
108  /* Galactic to equatorial */
109  eraTrxp( rmat, v1, v2 );
110 
111  /* Cartesian to spherical */
112  eraC2s( v2, dr, dd );
113 
114  /* Express in conventional ranges */
115  *dr = eraAnp( *dr );
116  *dd = eraAnpm( *dd );
117 
118 }
eraS2c
void eraS2c(double theta, double phi, double c[3])
Definition: s2c.c:3
eraC2s
void eraC2s(double p[3], double *theta, double *phi)
Definition: c2s.c:3
eraTrxp
void eraTrxp(double r[3][3], double p[3], double trp[3])
Definition: trxp.c:3
eraAnpm
double eraAnpm(double a)
Definition: anpm.c:3
eraAnp
double eraAnp(double a)
Definition: anp.c:3

+ Here is the call graph for this function:

+ Here is the caller graph for this function: