cluster_tools/mergeDosePerEnergyFile.sh

   1 #!/bin/bash
   2 set -u
   3
   4 function usage {
   5         echo "$0 -i <file1> -j <file2> -o <result>"
   6         exit 1
   7 }
   8
   9 if [ $# != 6 ]
  10 then
  11         usage
  12 fi
  13
  14 IN1=$2
  15 IN2=$4
  16 RESULT=$6
  17
  18 test -f ${IN1} && test -f ${IN2} || usage
  19
  20 TMP="$(mktemp)"
  21 echo "merging dose file"
  22 for PARAM in `awk '$1 == "#" {print $3}' ${IN1}`
  23 do
  24         echo "merging ${PARAM}"
  25         V1=`awk -v P=${PARAM} '$3 == P {print $4}' ${IN1} `
  26         V2=`awk -v P=${PARAM} '$3 == P {print $4}' ${IN2} `
  27         V1=`echo ${V1} | sed -e 's/[eE]+*/\*10\^/'`
  28         V2=`echo ${V2} | sed -e 's/[eE]+*/\*10\^/'`
  29         R=`echo "scale=30; ${V1} + ${V2}" | bc -l`
  30         test -z "${R}" && continue
  31         echo "# energydose ${PARAM} ${R}" >> ${TMP}
  32 done
  33 mv -f ${TMP} ${RESULT}