From 92ede64e801b015f565fe6ca4c6eead2779b27c2 Mon Sep 17 00:00:00 2001
From: Marc Vila-Oliva <vilaoliva@linux15.dg.creatis.insa-lyon.fr>
Date: Thu, 4 Jul 2013 19:06:58 +0200
Subject: [PATCH] new script for merging dose per energy files

---
 cluster_tools/mergeDosePerEnergyFile.sh | 33 +++++++++++++++++++++++++
 1 file changed, 33 insertions(+)
 create mode 100755 cluster_tools/mergeDosePerEnergyFile.sh

diff --git a/cluster_tools/mergeDosePerEnergyFile.sh b/cluster_tools/mergeDosePerEnergyFile.sh
new file mode 100755
index 0000000..f849016
--- /dev/null
+++ b/cluster_tools/mergeDosePerEnergyFile.sh
@@ -0,0 +1,33 @@
+#!/bin/bash
+set -u
+
+function usage {
+	echo "$0 -i <file1> -j <file2> -o <result>"
+	exit 1
+}
+
+if [ $# != 6 ]
+then
+	usage
+fi
+
+IN1=$2
+IN2=$4
+RESULT=$6
+
+test -f ${IN1} && test -f ${IN2} || usage
+
+TMP="$(mktemp)"
+echo "merging dose file"
+for PARAM in `awk '$1 == "#" {print $3}' ${IN1}`
+do
+        echo "merging ${PARAM}"
+	V1=`awk -v P=${PARAM} '$3 == P {print $4}' ${IN1} `
+	V2=`awk -v P=${PARAM} '$3 == P {print $4}' ${IN2} `
+        V1=`echo ${V1} | sed -e 's/[eE]+*/\*10\^/'`
+        V2=`echo ${V2} | sed -e 's/[eE]+*/\*10\^/'`
+	R=`echo "${V1} + ${V2}" | bc -l`
+        test -z "${R}" && continue
+        echo "# energydose ${PARAM} ${R}" >> ${TMP}
+done
+mv -f ${TMP} ${RESULT}
-- 
2.45.2