def METHOD_NAME(self, locs):
    """Compute shear modulus field at locations.
    """
    (npts, dim) = locs.shape
    METHOD_NAME = p_mu * numpy.ones((1, npts, 1), dtype=numpy.float64)
    return METHOD_NAME