Invert residuals into data-model

src/o_chires.c

@@ -250,6 +250,13 @@ int o_chi_lhood0(double *chi_ext, double *lhood0_ext, double *np_b0)
             if(M.output)
             {
                 fprintf(stderr, "Write output residuals with same header as input - residuals.fits\n");
+
+		//invert residuals to produce data - model instead of model - data (within contours)
+                for( i = 0; i < imFrame.ny; i++ )
+			for( j = 0; j < imFrame.nx; j++)
+				ero[i][j]=-ero[i][j];
+
+
                 wrf_from_header(imFrame.pixfile,"residuals.fits", imFrame.nx, imFrame.ny, ero);
             }