Index: /trunk/psLib/src/astro/psCoord.c
===================================================================
--- /trunk/psLib/src/astro/psCoord.c	(revision 2930)
+++ /trunk/psLib/src/astro/psCoord.c	(revision 2931)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.41 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-12-21 20:42:07 $
+*  @version $Revision: 1.42 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-01-07 23:52:15 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -31,11 +31,5 @@
 /******************************************************************************/
 
-#define PS_COT(X) (1.0 / atan(X))
-#define DEG_TO_RAD(DEGREES) ((DEGREES) * M_PI / 180.0)
-#define MIN_TO_RAD(MINUTES) ((MINUTES) * M_PI / (180.0 * 60.0))
-#define SEC_TO_RAD(SECONDS) ((SECONDS) * M_PI / (180.0 * 60.0 * 60.0))
-#define RAD_TO_DEG(RADIANS) ((RADIANS) * 180 / M_PI)
-#define RAD_TO_MIN(RADIANS) ((RADIANS) * 180 * 60.0 / M_PI)
-#define RAD_TO_SEC(RADIANS) ((RADIANS) * 180 * 60.0 * 60.0 / M_PI)
+// None
 
 /******************************************************************************/
@@ -181,5 +175,5 @@
 
     tmp->cosDeltaP = cos(deltaP);
-    tmp->sinDeltaP = cos(deltaP);
+    tmp->sinDeltaP = sin(deltaP);
     tmp->alphaP = alphaP;
     tmp->phiP = phiP;
@@ -317,4 +311,5 @@
         tmpF64 = 1.0 / tmpF64;
         tmp->x = 2.0 * tmpF64 * cos(coord->r) * sin(0.5 * coord->d);
+        tmp->y = sin(coord->r) * tmpF64;
 
     } else if (projection->type == PS_PROJ_PAR) {
@@ -553,6 +548,4 @@
 psSpherePrecess(coords, fromTime, toTime):
  
-XXX: What algorithm should be used?
-XXX: Verify that this is correct.
 XXX: Use static memory for tmpST.
  *****************************************************************************/
Index: /trunk/psLib/src/astronomy/psAstrometry.c
===================================================================
--- /trunk/psLib/src/astronomy/psAstrometry.c	(revision 2930)
+++ /trunk/psLib/src/astronomy/psAstrometry.c	(revision 2931)
@@ -8,6 +8,6 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.53 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2004-12-30 20:18:36 $
+ *  @version $Revision: 1.54 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2005-01-07 23:52:15 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -120,4 +120,5 @@
     double F = transform->y->coeff[0][1];
 
+    // XXX: Use the constructor here.
     psPlaneTransform *out = psAlloc(sizeof(psPlaneTransform));
 
Index: /trunk/psLib/src/astronomy/psCoord.c
===================================================================
--- /trunk/psLib/src/astronomy/psCoord.c	(revision 2930)
+++ /trunk/psLib/src/astronomy/psCoord.c	(revision 2931)
@@ -10,6 +10,6 @@
 *  @author GLG, MHPCC
 *
-*  @version $Revision: 1.41 $ $Name: not supported by cvs2svn $
-*  @date $Date: 2004-12-21 20:42:07 $
+*  @version $Revision: 1.42 $ $Name: not supported by cvs2svn $
+*  @date $Date: 2005-01-07 23:52:15 $
 *
 *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -31,11 +31,5 @@
 /******************************************************************************/
 
-#define PS_COT(X) (1.0 / atan(X))
-#define DEG_TO_RAD(DEGREES) ((DEGREES) * M_PI / 180.0)
-#define MIN_TO_RAD(MINUTES) ((MINUTES) * M_PI / (180.0 * 60.0))
-#define SEC_TO_RAD(SECONDS) ((SECONDS) * M_PI / (180.0 * 60.0 * 60.0))
-#define RAD_TO_DEG(RADIANS) ((RADIANS) * 180 / M_PI)
-#define RAD_TO_MIN(RADIANS) ((RADIANS) * 180 * 60.0 / M_PI)
-#define RAD_TO_SEC(RADIANS) ((RADIANS) * 180 * 60.0 * 60.0 / M_PI)
+// None
 
 /******************************************************************************/
@@ -181,5 +175,5 @@
 
     tmp->cosDeltaP = cos(deltaP);
-    tmp->sinDeltaP = cos(deltaP);
+    tmp->sinDeltaP = sin(deltaP);
     tmp->alphaP = alphaP;
     tmp->phiP = phiP;
@@ -317,4 +311,5 @@
         tmpF64 = 1.0 / tmpF64;
         tmp->x = 2.0 * tmpF64 * cos(coord->r) * sin(0.5 * coord->d);
+        tmp->y = sin(coord->r) * tmpF64;
 
     } else if (projection->type == PS_PROJ_PAR) {
@@ -553,6 +548,4 @@
 psSpherePrecess(coords, fromTime, toTime):
  
-XXX: What algorithm should be used?
-XXX: Verify that this is correct.
 XXX: Use static memory for tmpST.
  *****************************************************************************/