def METHOD_NAME(self):
    pn = op.network.Cubic(shape=[3, 3, 3])
    pn.add_model(propname='pore.z',
                 model=op.models.network.coordination_number)
    assert np.all(pn['pore.z'][:5] == np.array([3, 4, 3, 4, 5]))