def METHOD_NAME(self, locs):
    """Compute bulk modulus field at locations.
    """
    (npts, dim) = locs.shape
    METHOD_NAME = (p_lambda + 2.0 / 3.0 * p_mu) * numpy.ones((1, npts, 1), dtype=numpy.float64)
    return METHOD_NAME