def METHOD_NAME(self, locs):
    """Compute initial traction field at locations.
    """
    (npts, dim) = locs.shape
    traction = numpy.zeros((1, npts, self.SPACE_DIM), dtype=numpy.float64)
    traction[:,:, 2] = T0
    return traction