Index: branches/eam_branches/ipp-20130207/psModules/src/imcombine/pmStack.c
===================================================================
--- branches/eam_branches/ipp-20130207/psModules/src/imcombine/pmStack.c	(revision 35125)
+++ branches/eam_branches/ipp-20130207/psModules/src/imcombine/pmStack.c	(revision 35236)
@@ -1418,5 +1418,7 @@
       }
       combined->image->data.F32[y][x] = stats->robustMedian;
-      combined->mask->data.PS_TYPE_IMAGE_MASK_DATA[y][x] = 0;
+      if (combined->mask) {
+	combined->mask->data.PS_TYPE_IMAGE_MASK_DATA[y][x] = 0;
+      }
 #if (0)
       if ((x == 59)&&(y > 40)&&(y < 50)) {
Index: branches/eam_branches/ipp-20130207/psModules/src/objects/models/pmModel_TRAIL.c
===================================================================
--- branches/eam_branches/ipp-20130207/psModules/src/objects/models/pmModel_TRAIL.c	(revision 35125)
+++ branches/eam_branches/ipp-20130207/psModules/src/objects/models/pmModel_TRAIL.c	(revision 35236)
@@ -256,5 +256,7 @@
     if (!isfinite(axes.theta)) return false;
 
-    float size = (axes.major > source->moments->Mrf) ? axes.major : source->moments->Mrf;
+    float size = (axes.major > sqrt(source->moments->Mrf)) ? axes.major : sqrt(source->moments->Mrf);
+    //    if (size > psfAxes.major) { size -= psfAxes.major; }
+    //else { size = psfAxes.major; }
 
     // axes.major is a sigma in the major direction; scale to