Index: /branches/eam_02_branch/psphot/src/psphotOutput.c
===================================================================
--- /branches/eam_02_branch/psphot/src/psphotOutput.c	(revision 12980)
+++ /branches/eam_02_branch/psphot/src/psphotOutput.c	(revision 12981)
@@ -151,68 +151,2 @@
     return header;
 }
-
-# if (0)
-// output functions: we have several fixed modes (sx, obj, cmp)
-void psphotOutput (pmReadout *readout, psMetadata *arguments) {
-
-    bool status;
-    char *outputFile = NULL;
-
-    psMetadata *header = pmReadoutGetHeader (readout);
-
-    pmPSF *psf = psMetadataLookupPtr (&status, readout->analysis, "PSPHOT.PSF");
-    // sample PSF images??
-    if (psfSample != NULL) psphotSamplePSFs (psf, readout->image, psfSample);
-    if (psfSample != NULL) psphotSamplePSFs (psf, readout->image, psfSample);
-
-    if (psfFile != NULL) {
-	psMetadata *psfData = pmPSFtoMD (NULL, psf);
-	psMetadataConfigWrite (psfData, psfFile);
-	psFree (psfData);
-    }
-}
-
-psImage *pmModelPSFatXY (psImage *image, pmModel *modelEXT, pmPSF *psf, int x, int y, int dx, int dy) {
-
-    psRegion region = {x - dx, x + dx, y - dy, y + dy};
-    psImage *view = psImageSubset (image, region);
-    psImage *sample = psImageCopy (NULL, view, PS_TYPE_F32);
-    psImageInit (sample, 0);
-    modelEXT->params->data.F32[2] = x;
-    modelEXT->params->data.F32[3] = y;
-    pmModel *modelPSF = pmModelFromPSF (modelEXT, psf);
-    pmModelAdd (sample, NULL, modelPSF, PM_MODEL_OP_FULL);
-    psFree (modelPSF);
-    return (sample);
-}
-
-bool psphotSamplePSFs (pmPSF *psf, psImage *image, char *output) {
-
-    // make sample PSFs for 4 corners and the center
-    psImage *sample;
-
-    // optional dump of all rough source data
-    if (output[0] == 0) return false;
-
-    pmModel *modelEXT = pmModelAlloc (psf->type);
-    modelEXT->params->data.F32[0] = 0;
-    modelEXT->params->data.F32[1] = 1;
-
-    psFits *fits = psFitsOpen (output, "w");
-
-    // the centers are in parent coordinates; they do not need to correspond to valid pixels...
-    sample = pmModelPSFatXY (image, modelEXT, psf, 25, 25, 25, 25);
-    psFitsWriteImage (fits, NULL, sample, 0);
-    sample = pmModelPSFatXY (image, modelEXT, psf, image->numCols - 25, image->numRows - 25, 25, 25);
-    psFitsWriteImage (fits, NULL, sample, 0);
-    sample = pmModelPSFatXY (image, modelEXT, psf, image->numCols - 25, 25, 25, 25);
-    psFitsWriteImage (fits, NULL, sample, 0);
-    sample = pmModelPSFatXY (image, modelEXT, psf, 25, image->numRows - 25, 25, 25);
-    psFitsWriteImage (fits, NULL, sample, 0);
-    sample = pmModelPSFatXY (image, modelEXT, psf, image->numCols / 2, image->numRows / 2, 25, 25);
-    psFitsWriteImage (fits, NULL, sample, 0);
-
-    psFitsClose (fits);
-    return (TRUE);
-}
-# endif
