# 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 beta in 4.0 4.5 5.0 
do
for mass in 0.0025 0.02 0.03 0.035 0.04 0.045 0.05 0.055 0.06 0.07 0.08 0.09
do

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

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

done

done