# Do the measurement of the average energy per site for the HiPCiP project -
# Two-Dimensional XY Model
# Coded by Chris
# 26/05/1999

echo  >h4epsilon.awk "/Average energy over/,/sweeps/{T=T+0.2;print T,\$6,\$8;}"

echo -e >h4epsilon.dat "1\nf\n20\n5000\n0.2\n0.1\nf\n1000"
./HiPCiP4 <h4epsilon.dat >h4epsilon.out
echo -e >h4epsilon.dat "1\nf\n20\n5000\n0.4\n0.2\nf\n1000"
./HiPCiP4 <h4epsilon.dat >>h4epsilon.out
echo -e >h4epsilon.dat "1\nf\n20\n5000\n0.6\n0.3\nf\n1000"
./HiPCiP4 <h4epsilon.dat >>h4epsilon.out
echo -e >h4epsilon.dat "1\nf\n20\n5000\n0.8\n0.3\nf\n1000"
./HiPCiP4 <h4epsilon.dat >>h4epsilon.out
echo -e >h4epsilon.dat "1\nf\n20\n5000\n1.0\n0.4\nf\n1000"
./HiPCiP4 <h4epsilon.dat >>h4epsilon.out
echo -e >h4epsilon.dat "1\nf\n20\n5000\n1.2\n0.4\nf\n1000"
./HiPCiP4 <h4epsilon.dat >>h4epsilon.out
echo -e >h4epsilon.dat "1\nf\n20\n5000\n1.4\n0.5\nf\n1000"
./HiPCiP4 <h4epsilon.dat >>h4epsilon.out
echo -e >h4epsilon.dat "1\nf\n20\n5000\n1.6\n0.5\nf\n1000"
./HiPCiP4 <h4epsilon.dat >>h4epsilon.out
echo -e >h4epsilon.dat "1\nf\n20\n5000\n1.8\n0.6\nf\n1000"
./HiPCiP4 <h4epsilon.dat >>h4epsilon.out
echo -e >h4epsilon.dat "1\nf\n20\n5000\n2.0\n0.6\nf\n1000"
./HiPCiP4 <h4epsilon.dat >>h4epsilon.out

awk -f h4epsilon.awk h4epsilon.out >epsilon.dat
rm h4epsilon.dat h4epsilon.out h4epsilon.awk


