Index: trunk/Ohana/src/opihi/cmd.data/gridify.c
===================================================================
--- trunk/Ohana/src/opihi/cmd.data/gridify.c	(revision 2876)
+++ trunk/Ohana/src/opihi/cmd.data/gridify.c	(revision 2938)
@@ -3,5 +3,5 @@
 int gridify (int argc, char **argv) {
 
-  int i, Nx, Ny, Xb, Yb;
+  int i, Nx, Ny, Xb, Yb, Normalize, N;
   float Xmin, Xmax, dX, Ymin, Ymax, dY;
   float *buf, *val, *x, *y, *z;
@@ -9,4 +9,10 @@
   Buffer *bf;
   Vector *vx, *vy, *vz;
+
+  Normalize = TRUE;
+  if ((N = get_argument (argc, argv, "-raw"))) {
+    remove_argument (N, &argc, argv);
+    Normalize = FALSE;
+  }
 
   if (argc != 11) {
@@ -66,5 +72,5 @@
     if (Xb >= Nx) continue;
     if (Yb >= Ny) continue;
-    val[Xb + Yb*Nx] = *z;
+    val[Xb + Yb*Nx] += *z;
     Nval[Xb + Yb*Nx]++;
   }
@@ -72,9 +78,13 @@
   buf = (float *) bf[0].matrix.buffer;
   for (i = 0; i < Nx*Ny; i++) {
-    if (Nval[i] == 0) {
-      buf[i] = 0;
-      continue;
+    if (Normalize) {
+      if (Nval[i] == 0) {
+	buf[i] = 0;
+	continue;
+      }
+      buf[i] = val[i] / Nval[i];
+    } else {
+      buf[i] = val[i];
     }
-    buf[i] = val[i] / Nval[i];
   }
 
Index: trunk/Ohana/src/opihi/cmd.data/imsmooth.c
===================================================================
--- trunk/Ohana/src/opihi/cmd.data/imsmooth.c	(revision 2876)
+++ trunk/Ohana/src/opihi/cmd.data/imsmooth.c	(revision 2938)
@@ -1,12 +1,18 @@
 # include "data.h"
-# define NSIGMA 3
 
 int imsmooth (int argc, char **argv) {
   
-  int i, j, n, Nx, Ny, Ns, Ngauss;
+  int i, j, n, N, Nx, Ny, Ns, Ngauss;
   float *vi, *vo, *gauss, *gaussnorm;
-  float g, s, sigma;
+  float g, s, sigma, Nsigma;
   Buffer *in;
   float *temp;
+
+  Nsigma = 3;
+  if ((N = get_argument (argc, argv, "-Nsigma"))) {
+    remove_argument (N, &argc, argv);
+    Nsigma = atof (argv[N]);
+    remove_argument (N, &argc, argv);
+  }
 
   if (argc != 3) {
@@ -23,5 +29,5 @@
 
   /* build a 1D gaussian */
-  Ns = (int) (NSIGMA*sigma + 0.5);
+  Ns = (int) (Nsigma*sigma + 0.5);
   Ngauss = 2*Ns + 1;
   ALLOCATE (gaussnorm, float, Ngauss);