# prints out t, gamma(t), Sigma(t)
# the arrays below must be tuned to your specific computation

lat=1632
#beta=$2
#mass=$3

beta=4.0

for mass in 0.0025 0.02 0.035 0.04
do

rm sigma1632_${beta}_${mass}_1
for t in 1.5 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0
do

cat Out/fWflow_low_${lat}_${beta}_-0.25_${mass}_${t}1.[5-9][0-9] Out/fWflow_low_${lat}_${beta}_-0.25_${mass}_${t}1.[1-9][0-9][0-9] >tt
awk -v lat=$lat -v beta=$beta -v mass=$mass -v t=$t '/pbp:/{x+=($4+$5)/2;x2+=($4+$5)**2/4;i++}END{x=x/i;x2=x2/i-x**2;print t"1",x,sqrt(x2/(i-1)),beta,mass,lat}' tt >>sigma1632_${beta}_${mass}_1
done

done