Index: trunk/psLib/src/astronomy/psCoord.c
===================================================================
--- trunk/psLib/src/astronomy/psCoord.c	(revision 2583)
+++ trunk/psLib/src/astronomy/psCoord.c	(revision 2600)
@@ -8,8 +8,8 @@
 *  @ingroup CoordinateTransform
 *
-*  @author George Gusciora, MHPCC
+*  @author GLG, MHPCC
 *
-*  @version $Revision: 1.36 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-12-01 19:56:05 $
+*  @version $Revision: 1.37 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2004-12-02 21:12:51 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -164,14 +164,14 @@
 This function prototype has been modified since the SDRS.
  *****************************************************************************/
-psSphereTransform* psSphereTransformAlloc(double NPlat,
-        double Xo,
-        double xo)
+psSphereTransform* psSphereTransformAlloc(double alphaP,
+        double deltaP,
+        double phiP)
 {
     psSphereTransform* tmp = (psSphereTransform* ) psAlloc(sizeof(psSphereTransform));
 
-    tmp->sinPhi = sin(NPlat);
-    tmp->cosPhi = cos(NPlat);
-    tmp->Xo = Xo;
-    tmp->xo = xo;
+    tmp->cosDeltaP = cos(deltaP);
+    tmp->sinDeltaP = cos(deltaP);
+    tmp->alphaP = alphaP;
+    tmp->phiP = phiP;
 
     return (tmp);
@@ -207,11 +207,11 @@
     x = coord->r;
     y = coord->d;
-    dx = x - transform->xo;
-    sinY = cos(y) * sin(dx) * transform->sinPhi + sin(y) * transform->cosPhi;
+    dx = x - transform->phiP;
+    sinY = cos(y) * sin(dx) * transform->sinDeltaP + sin(y) * transform->cosDeltaP;
     cosY = sqrt(1.0 - sinY * sinY);
-    sinX = (cos(y) * sin(dx) * transform->cosPhi - sin(y) * transform->sinPhi) / cos(y);
+    sinX = (cos(y) * sin(dx) * transform->cosDeltaP - sin(y) * transform->sinDeltaP) / cos(y);
     cosX = cos(y) * cos(dx) / cos(y);
 
-    out->r = atan2(sinX, cosX) + transform->Xo;
+    out->r = atan2(sinX, cosX) + transform->alphaP;
     out->d = atan2(sinY, cosY);
 
@@ -227,9 +227,9 @@
     double T = year / 100.0;
     double phi = -23.452294 + 0.013013 * T + 0.000001639 * T * T - 0.000000503 * T * T * T;
-    double Xo = 0.0;
-    double xo = 0.0;
+    double alphaP = 0.0;
+    double phiP = 0.0;
 
     psFree(tmTime);
-    return (psSphereTransformAlloc(phi, Xo, xo));
+    return (psSphereTransformAlloc(phi, alphaP, phiP));
 }
 
@@ -242,9 +242,9 @@
     double T = year / 100.0;
     double phi = +23.452294 - 0.013013 * T - 0.000001639 * T * T + 0.000000503 * T * T * T;
-    double Xo = 0.0;
-    double xo = 0.0;
+    double alphaP = 0.0;
+    double phiP = 0.0;
 
     psFree(tmTime);
-    return (psSphereTransformAlloc(phi, Xo, xo));
+    return (psSphereTransformAlloc(phi, alphaP, phiP));
 }
 
@@ -536,10 +536,9 @@
                           const psTime *toTime)
 {
-    psF64 fromMJD = fromTime->sec/86400.0 + fromTime->usec/86400000000.0 + 40587.0;
-    psF64 toMJD = toTime->sec/86400.0 + toTime->usec/86400000000.0 + 40587.0;
-    psF64 = (toMJD - fromMJD) / 36525.0;
-
-
-
+    /*
+        psF64 fromMJD = fromTime->sec/86400.0 + fromTime->usec/86400000000.0 + 40587.0;
+        psF64 toMJD = toTime->sec/86400.0 + toTime->usec/86400000000.0 + 40587.0;
+        psF64 = (toMJD - fromMJD) / 36525.0;
+    */
     return(NULL);
 }