Index: /branches/eam_branches/ipp-20130904/psphot/src/psphotGalaxyShape.c
===================================================================
--- /branches/eam_branches/ipp-20130904/psphot/src/psphotGalaxyShape.c	(revision 36246)
+++ /branches/eam_branches/ipp-20130904/psphot/src/psphotGalaxyShape.c	(revision 36247)
@@ -326,4 +326,21 @@
     float dflux = flux * (dIo / Io);
 
+// Set this to be the model of choice.
+// XXX: Maybe: figure out which is the best and save that one.
+{
+    psFree (source->modelEXT);
+
+    source->modelEXT = psMemIncrRefCounter (pcm->modelConv);
+    source->type = PM_SOURCE_TYPE_EXTENDED;
+    source->mode |= PM_SOURCE_MODE_EXTMODEL;
+    source->mode |= PM_SOURCE_MODE_NONLINEAR_FIT;
+
+    // adjust the window so the subtraction covers the faint wings
+    // psphotSetRadiusMoments(&fitRadius, &windowRadius, readout, source, markVal);
+
+    // cache the model flux
+    pmPCMCacheModel (source, maskVal, psfSize);
+}
+
     psVectorAppend (source->galaxyFits->Flux, flux);
     psVectorAppend (source->galaxyFits->dFlux, dflux);