Index: /trunk/psastro/src/psastroOneChip.c
===================================================================
--- /trunk/psastro/src/psastroOneChip.c	(revision 13164)
+++ /trunk/psastro/src/psastroOneChip.c	(revision 13165)
@@ -88,4 +88,5 @@
     // XXX allow statitic to be set by the user
     psStats *fitStats = psStatsAlloc (PS_STAT_CLIPPED_MEAN | PS_STAT_CLIPPED_STDEV);
+    // psStats *fitStats = psStatsAlloc (PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
     fitStats->clipSigma = psMetadataLookupF32 (&status, recipe, "PSASTRO.CHIP.NSIGMA");
     fitStats->clipIter = psMetadataLookupS32 (&status, recipe, "PSASTRO.CHIP.NITER");
@@ -107,7 +108,9 @@
     // pixError is the average 1D scatter in pixels ('results' are in FPA units = microns)
     float pixError = 0.5*(results->xStats->clippedStdev + results->yStats->clippedStdev) / pixelScale;
+    // float pixError = 0.5*(results->xStats->robustStdev + results->yStats->robustStdev) / pixelScale;
 
     // astError is the average 1D scatter in arcsec ('results' are in FPA units = microns)
     float astError = 0.5*(results->xStats->clippedStdev + results->yStats->clippedStdev) * plateScale;
+    // float astError = 0.5*(results->xStats->robustStdev + results->yStats->robustStdev) * plateScale;
     int astNstar = results->yStats->clippedNvalues;
 
@@ -144,4 +147,7 @@
 	psastroDumpRawstars (rawstars, fpa, chip);
     }
+    if (psTraceGetLevel("psastro.dump") > 0) {
+	psastroDumpStars (refstars, "refstars.cal.dat");
+    }
 
     if (psTraceGetLevel("psastro.plot") > 0) {