#!/bin/bash

# Gather therodynamics results -> td.dat
${NRGPREFIX} gathertd_nosrt

# Average over z: td.dat -> td-avg.dat, and subtract reference data at the
# same time.
cp -v ../ref/td.dat td-ref.dat
${NRGPREFIX} tdavg -c

# Split according to columns
${NRGPREFIX} report td-avg.dat td