phlow=5
phhigh=12
nsets=4
delph=.1
ndata=`echo "($phhigh - $phlow) / $delph * $nsets" | hoc`
datafile=`echo "$1" | sed -e "s/\.out\./.data./"`
outfile="$1"


cat $datafile |
sed	-e "1s/^/Theory: /" \
	-e "s/nvert=[ 	][ 	]*[0-9]*[0-9]/nvert= 1/" \
	-e "s/ndata=[ 	][ 	]*[0-9]*[0-9]/ndata= $ndata/" \
	-e "/^starting/q"

cat $outfile |
sed -n	-e "/^function/,/^mean/{
		/^[-+ 	]*[0-9]/p
	}"

echo |
nawk '
	{
	print
}
END	{
	print "data"
	print
	count = 0
	for (ident = 0; ident < nsets && count < ndata; ident++) {
		for (pH = phlow; pH < (phhigh - delph/2) && count < ndata; pH += delph) {
			count++
			printf "0\t0\t0\t%g\t%d\n", pH, ident
		}
	}
	if (count != ndata) print "ERROR", count, ndata
}' ndata="$ndata" phlow="$phlow" phhigh="$phhigh" nsets="$nsets" delph="$delph" -