- declination = geoDec - mpar * rho * sin (gclat) * sin (g - geoDec) / sin(g);
-
- /* FG_LOG( FG_GENERAL, FG_INFO,
+ if (fabs(lat) > 0) {
+ declination
+ = geoDec - mpar * rho * sin (gclat) * sin (g - geoDec) / sin(g);
+ } else {
+ declination = geoDec;
+ // cerr << "Geocentric vs. Topocentric position" << endl;
+ // cerr << "RA (difference) : "
+ // << SGD_RADIANS_TO_DEGREES * (geoRa - rightAscension) << endl;
+ // cerr << "Dec (difference) : "
+ // << SGD_RADIANS_TO_DEGREES * (geoDec - declination) << endl;
+ }
+
+ /* SG_LOG( SG_GENERAL, SG_INFO,