def METHOD_NAME(data):
    etot = data["athena", "pressure"] / (data.ds.gamma - 1.0)
    etot += data["gas", "kinetic_energy_density"]
    if ("athena", "cell_centered_B_x") in self.field_list:
        etot += data["gas", "magnetic_energy_density"]
    return etot