Index: branches/eam_branches/20100225/Ohana/src/relastro/src/relastro.c
===================================================================
--- branches/eam_branches/20100225/Ohana/src/relastro/src/relastro.c	(revision 27421)
+++ branches/eam_branches/20100225/Ohana/src/relastro/src/relastro.c	(revision 27509)
@@ -1,3 +1,9 @@
 # include "relastro.h"
+
+# define MARKTIME(MSG,...) { \
+  float dtime; \
+  gettimeofday (&stop, (void *) NULL); \
+  dtime = DTIME (stop, start); \
+  fprintf (stderr, MSG, __VA_ARGS__); }
 
 int main (int argc, char **argv) {
@@ -6,6 +12,9 @@
   Catalog *catalog;
   FITS_DB db;
+  struct timeval start, stop;
 
   SkyList *skylist = NULL;
+
+  gettimeofday (&start, (void *) NULL);
 
   /* get configuration info, args */
@@ -26,14 +35,20 @@
   if (db.dbstate == LCK_EMPTY) Shutdown ("ERROR: No images in catalog %s (1)", db.filename);
   if (!dvo_image_load (&db, VERBOSE, FALSE)) Shutdown ("can't read image catalog %s", db.filename);
+  MARKTIME("load image data: %f sec\n", dtime);
 
   /* load regions and images based on specified sky patch (default depth) */
   skylist = load_images (&db, &UserPatch);
+  MARKTIME("load images: %f sec\n", dtime);
 
   /* load catalog data from region files : subselect high-quality measurements */
   catalog = load_catalogs (skylist, &Ncatalog, TRUE);
+  MARKTIME("load catalog data: %f sec\n", dtime);
 
   /* match measurements with images */
   initImageBins (catalog, Ncatalog);
+  MARKTIME("make image bins: %f sec\n", dtime);
+
   findImages (catalog, Ncatalog);
+  MARKTIME("set up image indexes: %f sec\n", dtime);
 
   if (PLOTSTUFF) {
Index: branches/eam_branches/20100225/Ohana/src/relphot/src/relphot.c
===================================================================
--- branches/eam_branches/20100225/Ohana/src/relphot/src/relphot.c	(revision 27421)
+++ branches/eam_branches/20100225/Ohana/src/relphot/src/relphot.c	(revision 27509)
@@ -1,4 +1,10 @@
 # include "relphot.h"
 # define USE_DIRECT 0
+
+# define MARKTIME(MSG,...) { \
+  float dtime; \
+  gettimeofday (&stop, (void *) NULL); \
+  dtime = DTIME (stop, start); \
+  fprintf (stderr, MSG, __VA_ARGS__); }
 
 int main (int argc, char **argv) {
@@ -7,6 +13,8 @@
   Catalog *catalog;
   FITS_DB db;
+  struct timeval start, stop;
+  SkyList *skylist = NULL;
 
-  SkyList *skylist = NULL;
+  gettimeofday (&start, (void *) NULL);
 
   /* get configuration info, args */
@@ -39,4 +47,5 @@
     if (!dvo_image_load (&db, VERBOSE, FALSE)) Shutdown ("can't read image catalog %s", db.filename);
   }
+  MARKTIME("load image data: %f sec\n", dtime);
 
   /* load regions and images based on specified sky patch */
@@ -45,4 +54,5 @@
   // XXX this is fairly lame: argv[1] is photcode if UserPatchSelect is true
   skylist = load_images (&db, argv[1], &UserPatch, UserPatchSelect);
+  MARKTIME("load images: %f sec\n", dtime);
 
   /* unlock, if we can (else, unlocked below) */
@@ -51,4 +61,5 @@
   /* load catalog data from region files */
   catalog = load_catalogs (skylist, &Ncatalog);
+  MARKTIME("load catalog data: %f sec\n", dtime);
   
   /* add in a loop over the catalogs calling dvo_catalog_chipcoords */
@@ -56,4 +67,6 @@
   /* match measurements with images, mosaics */
   initImageBins  (catalog, Ncatalog);
+  MARKTIME("make image bins: %f sec\n", dtime);
+
   initMosaicBins (catalog, Ncatalog);
   initGridBins   (catalog, Ncatalog);
@@ -61,5 +74,9 @@
 
   findImages (catalog, Ncatalog);
+  MARKTIME("set up image indexes: %f sec\n", dtime);
+
   findMosaics (catalog, Ncatalog);  /* also sets Grid values */
+  MARKTIME("set up mosaic indexes: %f sec\n", dtime);
+
   SAVEPLOT = FALSE;