Index: /trunk/Ohana/src/addstar/include/loadgalphot.h
===================================================================
--- /trunk/Ohana/src/addstar/include/loadgalphot.h	(revision 39023)
+++ /trunk/Ohana/src/addstar/include/loadgalphot.h	(revision 39024)
@@ -1,2 +1,4 @@
+
+int IMAGE_ID;
 
 typedef struct {
Index: /trunk/Ohana/src/addstar/src/args_loadgalphot.c
===================================================================
--- /trunk/Ohana/src/addstar/src/args_loadgalphot.c	(revision 39023)
+++ /trunk/Ohana/src/addstar/src/args_loadgalphot.c	(revision 39024)
@@ -83,4 +83,11 @@
   }
 
+  IMAGE_ID = 0;
+  if ((N = get_argument (*argc, argv, "-image-id"))) {
+    remove_argument (N, argc, argv);
+    IMAGE_ID = atoi (argv[N]);
+    remove_argument (N, argc, argv);
+  }
+
   /* extra error messages */
   VERBOSE = FALSE;
Index: /trunk/Ohana/src/addstar/src/loadgalphot_readstars.c
===================================================================
--- /trunk/Ohana/src/addstar/src/loadgalphot_readstars.c	(revision 39023)
+++ /trunk/Ohana/src/addstar/src/loadgalphot_readstars.c	(revision 39024)
@@ -232,4 +232,5 @@
     stars[i].galphot.detID     = ID_gal[i];
     stars[i].galphot.photcode  = options->photcode;
+    stars[i].galphot.imageID   = IMAGE_ID;
 
     // I have a grid of measurements with (Flux, dFlux, Chisq) at each point
Index: /trunk/Ohana/src/addstar/src/mkcmf.c
===================================================================
--- /trunk/Ohana/src/addstar/src/mkcmf.c	(revision 39023)
+++ /trunk/Ohana/src/addstar/src/mkcmf.c	(revision 39024)
@@ -267,5 +267,5 @@
   coords.crpix2 = CRPIX2;
   
-  if (isStack || isForcedWarp || isDiff) {
+  if (isStack || isForcedWarp || isForcedGalaxy || isDiff) {
     coords.cdelt1 = 0.25/3600.0;
     coords.cdelt2 = 0.25/3600.0;
@@ -350,5 +350,7 @@
 
   gfits_modify (&header, "NSTARS",   "%d", 1, Nstars);
-  gfits_modify (&header, "PHOTCODE", "%s", 1, photcode);
+  if (!isForcedGalaxy) {
+    gfits_modify (&header, "PHOTCODE", "%s", 1, photcode);
+  }
   if (isfinite(mjd)) {
     gfits_modify (&header, "MJD-OBS", "%lf", 1, mjd);