Integrated energy density is not equal to the total energy

**Describe the bug**
I might be doing something wrong, but I can't get the integrated energy density to be equal to the total energy. Specifically, I am comparing the VMC total energy from the `scalar.dat` file against the integrated energy density from the `stat.h5` file. The integration is carried out as `np.sum(self.total) * domain_volumes / np.sum(self.weights)` in the attached `get_edens.py` script.

VMC total: -84.891(5) Ha
Integrated Total Energy: -355.2907763843961

**To Reproduce**
An example is attached.
[vmc_J123.zip](https://github.com/user-attachments/files/17722718/vmc_J123.zip)

**Expected behavior**
VMC total energy and the integrated edens should be equal.

**System:**
NERSC