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

Definition at line 80 of file palEqgal.c.

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

Referenced by t_eqgal().

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( dr, dd, v1 );
107 
108  /* Equatorial to Galactic */
109  eraRxp( rmat, v1, v2 );
110 
111  /* Cartesian to spherical */
112  eraC2s( v2, dl, db );
113 
114  /* Express in conventional ranges */
115  *dl = eraAnp( *dl );
116  *db = eraAnpm( *db );
117 
118 }
void eraS2c(double theta, double phi, double c[3])
Definition: s2c.c:3
void eraC2s(double p[3], double *theta, double *phi)
Definition: c2s.c:3
void eraRxp(double r[3][3], double p[3], double rp[3])
Definition: rxp.c:3
double eraAnpm(double a)
Definition: anpm.c:3
double eraAnp(double a)
Definition: anp.c:3

+ Here is the call graph for this function:

+ Here is the caller graph for this function: