#!/bin/sh
# extracts double occupancy from .mag file
# parameters: skip  (default: skip=3)
START=`expr ${1:-3} + 1`
TMP=`mktemp`

cat > ${TMP}


BETA=`cat ${TMP} | grep beta | printcol.sh 4`
U=`cat ${TMP} | grep U | grep "=" | printcol.sh 4`
DT0=`grep -v "#" ${TMP} | head -1 | printcol.sh 1`
DT=`echo ${BETA} ${DT0} | awk '{printf "%s  %1.0f", $1, $1/$2}' | awk '{printf "%8.6f\n", 1.0*$1/$2}'`

grep -v "#" ${TMP} | printcol.sh 5 | tail -n +$START | stats -s | awk '{print "'$BETA'  '$U'  '$DT'  ", $1,"  " $2 }' 
rm ${TMP}