Merge commit '339420dd5dd874c41f6bab5808291fb4036dd022' as 'funtools'

1 files changed, 171 insertions, 0 deletions

diff --git a/funtools/funcnts_plot.tmpl b/funtools/funcnts_plot.tmpl
new file mode 100755
index 0000000..dace624
--- /dev/null
+++ b/funtools/funcnts_plot.tmpl
@@ -0,0 +1,171 @@
+#!/bin/sh
+# set -x
+
+LC_NUMERIC=en_US.ISO8859-1
+export LC_NUMERIC
+
+files="None"
+while [ x"$1" != x ]; do
+    case $1 in
+	-file)
+            files="$2"
+	    shift
+	    shift
+	    continue;;
+	*)  
+            break;;
+    esac
+done
+
+if [ x"$1" = xgnuplot ]; then
+
+  if [ x`which gnuplot 2>/dev/null` = x ]; then
+    echo "ERROR: gnuplot not available"
+    exit 1
+  fi
+  AWK '
+  BEGIN{
+    mode=1
+    xlabel="unknown"
+    ylabel="unknown"
+    prev="noprev"
+    ind=0
+    plots=0
+    counts = 0
+    initparam=0
+  }
+  mode==1{
+    if( initparam == 0 ){
+      if( files == "None" ){
+        files=""
+        dofile="true"
+      }
+      else{
+        dofile="false"
+      }
+      initparam=1
+    }
+    if( dofile == "true" && $1 == "#" && $2 == "data" && $3 == "file:" ){
+      if( files != "" ) files = files ","
+      files = files $4
+    }
+    else if( $1 == "#" && $2 == "radii:" ){
+      xlabel = $3
+    }
+    else if( $1 == "#" && $2 == "surf_bri:" ){
+      ylabel = $3
+    }
+    else if( $1 == "----" ){
+      head = sprintf("set title \"funcnts(%s)\"\n", files)
+      head = head sprintf("set xlabel \" radius(%s)\"\n", xlabel)
+      head = head sprintf("set ylabel \"surf_bri(%s)\"\n", ylabel)
+      head = head "plot "
+      mode = 2
+      next
+    }
+  }
+  mode==2{
+    if( NF == 12 ){
+      if( $11 == prev ){
+	ind += 1
+      }
+      else{
+	ind = 1
+        plots += 1
+        i = (plots - 1) * 4 + 2
+        if( plots > 1 )
+          head = head ", "
+        head = head sprintf("\"-\" using 1:%d:%d:%d:%d title \"profile #%d\" with boxerrorbars", i, i+1, i+2, i+3, plots)
+      }
+      if( data[ind] == "" )
+        data[ind] = sprintf("%s", ($9+$10)/2)
+      data[ind] = data[ind] sprintf(" %s %s %s %s", $7, $7-$8, $7+$8, $10-$9)
+      counts += $7
+      prev = $11
+    }
+    else{
+      mode = 3
+    }
+  }
+  END{
+    if( counts == 0 ){
+      exit 0
+    }
+    print head
+    for(i=1; i<=plots; i++){
+      for(j=1; j<=ind; j++){
+        print data[j]
+      }
+      print "e"
+    }    
+  }
+  ' files="$files" | gnuplot -persist - 1>/dev/null 2>&1
+  exit 0
+
+elif [ x"$1" = xds9 ]; then
+
+  AWK '
+  BEGIN{
+    mode=1
+    xlabel="unknown"
+    ylabel="unknown"
+    initparam=0
+  }
+  mode==1{
+    if( initparam == 0 ){
+      if( files == "None" ){
+        files=""
+        dofile="true"
+      }
+      else{
+        dofile="false"
+      }
+      initparam=1
+    }
+    if( dofile == "true" && $1 == "#" && $2 == "data" && $3 == "file:" ){
+      if( files != "" ) files = files ","
+      files = files $4
+    }
+    else if( $1 == "#" && $2 == "radii:" ){
+      xlabel = $3
+    }
+    else if( $1 == "#" && $2 == "surf_bri:" ){
+      ylabel = $3
+    }
+    else if( $1 == "----" ){
+      ds9header="'"${DS9_PLOT_HEADER}"'"
+      if( ds9header != "" ) {
+	  printf "%s 3\n", ds9header
+      }
+      else{
+          printf "funcnts(%s) avg_radius(%s) surf_bri(%s) 3\n", files, xlabel, ylabel
+      }
+      mode = 2
+      next
+    }
+  }
+  mode==2{
+    if( NF == 12 ){
+      if( $9 != "NA" ){
+        print ($9+$10)/2, $7, $8
+      }
+      else{
+         print "ERROR: annulus, circle, panda regions only"
+	 exit 1
+#        print $7, "non-annulus-region", "specified!"
+#        mode = 3
+      }
+    }
+    else{
+      mode = 3
+    }
+  }
+  ' files="$files"
+  exit 0
+else
+
+  echo "ERROR: unknown argument: " $1
+  echo "funcnts -r ... | funcnts.plot [ds9|gnuplot]"
+  exit 1
+
+fi