def defineBlock(ccMesh,p0,p1,vals=[0,1]):
    sigma = np.zeros(ccMesh.shape[0]) + vals[1]
    ind   = getIndicesBlock(p0,p1,ccMesh)

    sigma[ind] = vals[0]

    return mkvc(sigma)
Build a block with the conductivity specified by condVal.  Returns an array.
dataset/ETHPy150Open simpeg/simpeg/SimPEG/Utils/ModelBuilder.py/defineBlock