Index: branches/eam_branches/relphot.20210521/src/GridOps.c
===================================================================
--- branches/eam_branches/relphot.20210521/src/GridOps.c	(revision 41632)
+++ branches/eam_branches/relphot.20210521/src/GridOps.c	(revision 41637)
@@ -60,4 +60,6 @@
     GridCorr[i] = NULL;
   }
+
+  fprintf (stderr, "Generating grid corrections for %d photcodes\n", NGridCorr);
   return;
 }
@@ -102,4 +104,5 @@
     result = GridCorr[i];
     *Nlast = i;
+    break;
   }
   return result;
@@ -127,6 +130,8 @@
   Image *images = getimages (&Nimages, NULL);
 
+  int NGridReal = 0;
+
   for (int i = 0; i < Nimages; i++) {
-    int code = images[0].photcode;
+    int code = images[i].photcode;
     myAssert (code >= 0, "oops");
     myAssert (code < NGridCorr, "oops");
@@ -137,4 +142,6 @@
 
     if (GridCorr[code]) continue; // already created this one
+
+    NGridReal ++;
 
     ALLOCATE(GridCorr[code], GridCorrectionType, 1); 
@@ -157,4 +164,5 @@
   }
   resetMgrid(); // start with values of 0
+  fprintf (stderr, "Init grid corrections for %d photcodes\n", NGridReal);
 }
 
@@ -274,6 +282,5 @@
 	GridCorr[code]-> Mgrid[ix][iy] = Mgrid;
 	GridCorr[code]->dMgrid[ix][iy] = sqrt(r*(Mgrid2 - Mgrid*Mgrid)); // sample stdev
-	fprintf (stderr, "grid code %d, %d x %d : %f +/- %f : %d\n", code, ix, iy,
-		 GridCorr[code]-> Mgrid[ix][iy], GridCorr[code]->dMgrid[ix][iy], GridCorr[code]->nMgrid[ix][iy]);
+	// fprintf (stderr, "grid code %d, %d x %d : %f +/- %f : %d\n", code, ix, iy, GridCorr[code]-> Mgrid[ix][iy], GridCorr[code]->dMgrid[ix][iy], GridCorr[code]->nMgrid[ix][iy]);
       }
     }