Index: /trunk/psphot/src/psphotFullForceSummary.c
===================================================================
--- /trunk/psphot/src/psphotFullForceSummary.c	(revision 40808)
+++ /trunk/psphot/src/psphotFullForceSummary.c	(revision 40809)
@@ -301,9 +301,9 @@
         // there is now only a single chip (multiple readouts?). loop over it and process
         while ((cell = pmFPAviewNextCell (view, input->fpa, 1)) != NULL) {
- //           psLogMsg ("psphotFullForceSummary", 5, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
+	  // psLogMsg ("psphotFullForceSummary", 5, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
 
             // process each of the readouts
             while ((readout = pmFPAviewNextReadout (view, input->fpa, 1)) != NULL) {
-//                psLogMsg ("psphotFullForceSummary", 6, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
+	      // psLogMsg ("psphotFullForceSummary", 6, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
                 if (! readout->data_exists) { continue; }
 
Index: /trunk/psphot/src/psphotFullForceSummaryReadout.c
===================================================================
--- /trunk/psphot/src/psphotFullForceSummaryReadout.c	(revision 40808)
+++ /trunk/psphot/src/psphotFullForceSummaryReadout.c	(revision 40809)
@@ -1,4 +1,3 @@
 #include "psphotInternal.h"
-
 
 typedef struct {
@@ -637,5 +636,12 @@
                 // now save the model parameters in the model structure
 
-                pmPSF_AxesToModel (model->params->data.F32, axes, model->class->useReff);
+		// a crazy model can raise an assert in psEllise
+	        if ((fabs(axes.major) > 10000 * fabs(axes.minor)) || (fabs(axes.minor) > 10000 * fabs(axes.major))) {
+		  model->params->data.F32[PM_PAR_SXX] = NAN;
+		  model->params->data.F32[PM_PAR_SYY] = NAN;
+		  model->params->data.F32[PM_PAR_SXY] = NAN;
+		} else {
+		  pmPSF_AxesToModel (model->params->data.F32, axes, model->class->useReff);
+		}
 
                 model->mag = -2.5 * log10(flux0);