Index: trunk/psphot/src/psphotFullForceSummaryReadout.c
===================================================================
--- trunk/psphot/src/psphotFullForceSummaryReadout.c	(revision 40048)
+++ 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);