Index: /branches/pap_branch_080617/ppSub/src/ppSubReadout.c
===================================================================
--- /branches/pap_branch_080617/ppSub/src/ppSubReadout.c	(revision 18287)
+++ /branches/pap_branch_080617/ppSub/src/ppSubReadout.c	(revision 18288)
@@ -197,4 +197,5 @@
 #ifdef TESTING
     {
+        pmReadoutMaskApply(minuend, maskBlank);
         psFits *fits = psFitsOpen("minuend.fits", "w");
         psFitsWriteImage(fits, NULL, minuend->image, 0, NULL);
@@ -202,4 +203,5 @@
     }
     {
+        pmReadoutMaskApply(subtrahend, maskBlank);
         psFits *fits = psFitsOpen("subtrahend.fits", "w");
         psFitsWriteImage(fits, NULL, subtrahend->image, 0, NULL);
Index: /branches/pap_branch_080617/ppSub/test/fake.c
===================================================================
--- /branches/pap_branch_080617/ppSub/test/fake.c	(revision 18287)
+++ /branches/pap_branch_080617/ppSub/test/fake.c	(revision 18288)
@@ -7,8 +7,8 @@
 
 // PSF for images
-#define X_AXIS_1 3.0                    // Length of x axis (FWHM) for image 1
-#define Y_AXIS_1 5.0                    // Length of y axis (FWHM) for image 1
-#define X_AXIS_2 5.0                    // Length of x axis (FWHM) for image 1
-#define Y_AXIS_2 3.0                    // Length of y axis (FWHM) for image 1
+#define X_AXIS_1 2.0                    // Length of x axis (FWHM) for image 1
+#define Y_AXIS_1 3.0                    // Length of y axis (FWHM) for image 1
+#define X_AXIS_2 3.0                    // Length of x axis (FWHM) for image 1
+#define Y_AXIS_2 2.0                    // Length of y axis (FWHM) for image 1
 
 // Image parameters
@@ -22,5 +22,5 @@
 #define SOURCES_DENSITY 30000.0         // Source density at nominated magnitude (per deg^2)
 #define SOURCES_ZP 25.0                 // Magnitude ZP
-#define SOURCES_GAUSSIAN 1              // Gaussian sources (boolean)?
+#define SOURCES_GAUSSIAN 0              // Gaussian sources (boolean)?
 
 // Noise properties