def getPointRange(iFace, dims):
    """Return the correct point range for face iFace on a block with
    dimensions given in dims"""
    il = dims[0]
    jl = dims[1]
    kl = dims[2]
    if iFace == 0:
        return [[1, il], [1, jl], [1, 1]]
    elif iFace == 1:
        return [[1, il], [1, jl], [kl, kl]]
    elif iFace == 2:
        return [[1, 1], [1, jl], [1, kl]]
    elif iFace == 3:
        return [[il, il], [1, jl], [1, kl]]
    elif iFace == 4:
        return [[1, il], [1, 1], [1, kl]]
    elif iFace == 5:
        return [[1, il], [jl, jl], [1, kl]]