def METHOD_NAME(self, data, w=1.0, axis=-1):
    """
    Evaluate a numerical volume integral of the given data
    on this grid.
    data: Data to integrate.
    w:    Optional weighting function.
    axis: Axis to integrate over.
    """
    return (self.VpVol*self.dr * w * data).sum(axis)