#!/bin/sh

# Gather therodynamics results -> td.dat
gathertd_nosrt

# Copy the reference results w/o impurity for subtraction
cp ../td_0/td.dat td-ref.dat

# Average over z: td.dat -> td-avg.dat
tdavg -c

# Split according to columns
report td-avg.dat td

# Charge susceptibility
./susc td-Q2.dat td-Q.dat td-Qimp.dat
