Index: /trunk/psphot/src/psphotSourceFits.c
===================================================================
--- /trunk/psphot/src/psphotSourceFits.c	(revision 10267)
+++ /trunk/psphot/src/psphotSourceFits.c	(revision 10268)
@@ -291,6 +291,7 @@
     moments.y2 = source->moments->Sy;
     moments.xy = source->moments->Sxy;
-    axes = psEllipseMomentsToAxes (moments);
-
+    axes = psEllipseMomentsToAxes (moments, 20.0);
+
+    // XXX this is really arbitrary: 4 pixel separation?
     dx = 2 * cos (axes.theta);
     dy = 2 * sin (axes.theta);
@@ -305,12 +306,12 @@
     DBL = pmModelCopy (PSF);
     DBL->params->data.F32[1] *= 0.5;
-    DBL->params->data.F32[2] = source->moments->x + dx;
-    DBL->params->data.F32[3] = source->moments->y + dy;
+    DBL->params->data.F32[2] = source->peak->xf + dx;
+    DBL->params->data.F32[3] = source->peak->yf + dy;
     modelSet->data[0] = DBL;
 
     DBL = pmModelCopy (PSF);
     DBL->params->data.F32[1] *= 0.5;
-    DBL->params->data.F32[2] = source->moments->x - dx;
-    DBL->params->data.F32[3] = source->moments->y - dy;
+    DBL->params->data.F32[2] = source->peak->xf - dx;
+    DBL->params->data.F32[3] = source->peak->yf - dy;
     modelSet->data[1] = DBL;