# Do the measurement of the specific heat for the HiPCiP project -
# Two-Dimensional XY Model
# Coded by Chris
# 03/06/1999

echo  >h4heat.awk "/Temperature/ {T=\$2;}"
echo  >>h4heat.awk "/Specific heat/,/Ergodicity/ {print T,\$3; T=T+0.02;}"

echo >h4heat.out "Temperature 0.20"
echo -e >h4heat.dat "3\nf\n20\n5000\n0.20\n0.1\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 0.30"
echo -e >h4heat.dat "3\nf\n20\n5000\n0.30\n0.1\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 0.40"
echo -e >h4heat.dat "3\nf\n20\n5000\n0.40\n0.2\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 0.50"
echo -e >h4heat.dat "3\nf\n20\n5000\n0.50\n0.2\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 0.60"
echo -e >h4heat.dat "3\nf\n20\n5000\n0.60\n0.2\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 0.70"
echo -e >h4heat.dat "3\nf\n20\n5000\n0.70\n0.2\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out

echo >>h4heat.out "Temperature 0.80"
echo -e >h4heat.dat "3\nf\n20\n10000\n0.80\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n0.82\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n0.84\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n0.86\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n0.88\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out

echo -e >h4heat.dat "3\nf\n20\n10000\n0.90\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n0.92\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n0.94\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n0.96\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n0.98\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out

echo -e >h4heat.dat "3\nf\n20\n10000\n1.00\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.02\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.04\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.06\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.08\n0.3\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out

echo -e >h4heat.dat "3\nf\n20\n10000\n1.10\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.12\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.14\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.16\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.18\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out

echo -e >h4heat.dat "3\nf\n20\n10000\n1.20\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.22\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.24\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.26\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo -e >h4heat.dat "3\nf\n20\n10000\n1.28\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out

echo -e >h4heat.dat "3\nf\n20\n5000\n1.30\n0.4\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 1.40"
echo -e >h4heat.dat "3\nf\n20\n5000\n1.40\n0.5\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 1.50"
echo -e >h4heat.dat "3\nf\n20\n5000\n1.50\n0.5\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 1.60"
echo -e >h4heat.dat "3\nf\n20\n5000\n1.60\n0.6\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 1.70"
echo -e >h4heat.dat "3\nf\n20\n5000\n1.70\n0.6\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 1.80"
echo -e >h4heat.dat "3\nf\n20\n5000\n1.80\n0.7\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 1.90"
echo -e >h4heat.dat "3\nf\n20\n5000\n1.90\n0.7\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out
echo >>h4heat.out "Temperature 2.00"
echo -e >h4heat.dat "3\nf\n20\n5000\n2.00\n0.8\nf\n1000"
./HiPCiP4 <h4heat.dat >>h4heat.out

awk -f h4heat.awk h4heat.out >heat.dat
rm h4heat.dat h4heat.out h4heat.awk


