Index: branches/eam_branches/ipp-20130904/psphot/src/psphotGalaxyShape.c
===================================================================
--- branches/eam_branches/ipp-20130904/psphot/src/psphotGalaxyShape.c	(revision 36209)
+++ branches/eam_branches/ipp-20130904/psphot/src/psphotGalaxyShape.c	(revision 36210)
@@ -6,5 +6,5 @@
 
     fprintf (stdout, "\n");
-    psLogMsg ("psphot", PS_LOG_INFO, "--- psphot Kron Fluxes ---");
+    psLogMsg ("psphot", PS_LOG_INFO, "--- psphot Galaxy Shapes ---");
 
     // select the appropriate recipe information
@@ -58,9 +58,11 @@
 
     if (!sources->n) {
-        psLogMsg ("psphot", PS_LOG_INFO, "no sources, skipping kron fluxes");
+        psLogMsg ("psphot", PS_LOG_INFO, "no sources, skipping galaxy shapes");
         return true;
     }
 
     psTimerStart ("psphot.galaxy");
+
+    psphotInitRadiusEXT (recipe, readout);
 
     // determine the number of allowed threads
@@ -218,4 +220,5 @@
     // we are using fitOptions->mode : be sure this makes sense
     pmPCMdata *pcm = pmPCMinit (source, fitOptions, model, maskVal, psfSize);
+    if (!pcm) return false;
 
     // we are fitting only PM_PAR_I0; the shape elements are generated from a grid
@@ -302,5 +305,5 @@
     float Chisq = YY - 2 * Io * YM + Io * Io * MM;
 
-    fprintf (stderr, "Io: %f, Chisq: %f\n", Io, Chisq);
+    // fprintf (stderr, "Io: %f, Chisq: %f\n", Io, Chisq);
 
     psVectorAppend (source->galaxyFits->Io, Io);