Index: /branches/eam_branches/ipp-20110213/Ohana/src/opihi/cmd.data/fit1d.c
===================================================================
--- /branches/eam_branches/ipp-20110213/Ohana/src/opihi/cmd.data/fit1d.c	(revision 31143)
+++ /branches/eam_branches/ipp-20110213/Ohana/src/opihi/cmd.data/fit1d.c	(revision 31144)
@@ -42,5 +42,5 @@
 
   if (argc != 4) {
-    gprint (GP_ERR, "USAGE: fit x y order [-dy wt] [-quiet/-q] [-clip Nsigma Niter]\n");
+    gprint (GP_ERR, "USAGE: fit1d x y order [-dy wt] [-quiet/-q] [-clip Nsigma Niter]\n");
     return (FALSE);
   }
Index: /branches/eam_branches/ipp-20110213/Ohana/src/opihi/cmd.data/fit2d.c
===================================================================
--- /branches/eam_branches/ipp-20110213/Ohana/src/opihi/cmd.data/fit2d.c	(revision 31143)
+++ /branches/eam_branches/ipp-20110213/Ohana/src/opihi/cmd.data/fit2d.c	(revision 31144)
@@ -48,5 +48,5 @@
 
   if (argc != 5) {
-    gprint (GP_ERR, "USAGE: fit x y z order [-dz wt]\n");
+    gprint (GP_ERR, "USAGE: fit2d x y z order [-dz wt] [-quiet/-q] [-clip Nsigma Niter]\n");
     return (FALSE);
   }
Index: /branches/eam_branches/ipp-20110213/Ohana/src/opihi/dvo/skycoverage.c
===================================================================
--- /branches/eam_branches/ipp-20110213/Ohana/src/opihi/dvo/skycoverage.c	(revision 31143)
+++ /branches/eam_branches/ipp-20110213/Ohana/src/opihi/dvo/skycoverage.c	(revision 31144)
@@ -5,5 +5,5 @@
 int skycoverage (int argc, char **argv) {
 
-  int WITH_MOSAIC, SOLO_MOSAIC;
+  int WITH_MOSAIC, SOLO_MOSAIC, ShowDensity;
   off_t i, Nimage;
   int N, status, TimeSelect, ByName, xs, ys;
@@ -60,4 +60,10 @@
   }
 
+  ShowDensity = FALSE;
+  if ((N = get_argument (argc, argv, "-density"))) {
+    remove_argument (N, &argc, argv);
+    ShowDensity = TRUE;
+  }
+
   ByName = FALSE;
   if ((N = get_argument (argc, argv, "-name"))) {
@@ -141,5 +147,5 @@
   if (argc != 3) {
     gprint (GP_ERR, "USAGE: skycoverage (buffer) (Npts)\n");
-    gprint (GP_ERR, "  options: [-scale pixscale] [-center ra dec] [-size Nx Nx] [-proj projection] [-time start range] [-trange start stop] [-name name] [-photcode name] [+mosaic] [-mosaic]\n");
+    gprint (GP_ERR, "  options: [-scale pixscale] [-center ra dec] [-size Nx Nx] [-proj projection] [-time start range] [-trange start stop] [-name name] [-photcode name] [+mosaic] [-mosaic] [-density]\n");
     gprint (GP_ERR, "       (buffer) saves bitmapped image\n");
     gprint (GP_ERR, "       (Npts) gives the number of test points per image in each dimension\n");
@@ -148,4 +154,5 @@
     gprint (GP_ERR, "       -size (Nx) (Ny)    : specifies the size of the image [360/scale, 180/scale]\n");
     gprint (GP_ERR, "       -proj (projection) : specifies the projection choice [AIT]\n");
+    gprint (GP_ERR, "       -density           : create image with relative density (else binary on/off)\n");
     gprint (GP_ERR, "       note: we need 64800 / (pixscale)^2 pixels to represent the sky\n");
     return (FALSE);
@@ -186,5 +193,4 @@
 
   V = (float *)buf[0].matrix.buffer;
-  bzero (V, Nx*Ny*sizeof(float));
 
   for (ys = 0; ys < Ny; ys++) {
@@ -194,5 +200,7 @@
       status &= (r <= 360);
       if (status) {
-	V[ys*Nx + xs] = 2;
+	V[ys*Nx + xs] = ShowDensity ?  0 : 2;
+      } else {
+	V[ys*Nx + xs] = ShowDensity ? -1 : 0;
       }
     }
@@ -249,5 +257,9 @@
 	  xs = (int)Xs;
 	  ys = (int)Ys;
-	  V[ys*Nx + xs] = 1;
+	  if (ShowDensity) {
+	      V[ys*Nx + xs] += 1;
+	  } else {
+	      V[ys*Nx + xs] = 1;
+	  }
 	}
       }