Index: trunk/psphot/src/fit_galaxies.c
===================================================================
--- trunk/psphot/src/fit_galaxies.c	(revision 4574)
+++ trunk/psphot/src/fit_galaxies.c	(revision 4582)
@@ -14,5 +14,5 @@
 
     float MOMENT_R = psMetadataLookupF32 (&status, config, "GAL_MOMENTS_RADIUS");
-    float RADIUS   = psMetadataLookupF32 (&status, config, "FIT_RADIUS");
+    // float RADIUS   = psMetadataLookupF32 (&status, config, "FIT_RADIUS");
     float snFaint  = psMetadataLookupF32 (&status, config, "FAINT_SN_LIM");
     float OUTER    = psMetadataLookupF32 (&status, config, "OUTER_RADIUS");
@@ -22,5 +22,5 @@
     float FLUX_LIMIT  = FIT_NSIGMA * skyStats->sampleStdev;
 
-    psModelType   modelType   = psModelSetType ("PS_MODEL_RGAUSS");
+    psModelType   modelType   = psModelSetType ("PS_MODEL_SGAUSS");
     psModelRadius modelRadius = psModelRadius_GetFunction (modelType);
 
@@ -34,5 +34,9 @@
 
 	// recalculate the source moments using the galaxy radius (larger)
-	status = pmSourceMoments (source, MOMENT_R);
+	status = pmSourceMoments_EAM (source, MOMENT_R);
+	if (!status) {
+	  fprintf (stderr, "invalid moments, skipping\n");
+	  continue;
+	}
 
 	// use the source moments, etc to guess basic model parameters