def METHOD_NAME(self):
    net = op.network.Cubic(shape=[5, 5, 5])
    # Create first isolated cluster
    Ps = net['pore.coords'][:, 2] == 2.5
    Ps = Ps*(net['pore.coords'][:, 1] > 2.5)
    Ts = net.find_neighbor_throats(pores=Ps, mode='exclusive_or')
    trim(network=net, throats=Ts)
    # Create a second isolated cluster
    Ps = net['pore.coords'][:, 2] == 3.5
    Ps = Ps*(net['pore.coords'][:, 1] < 2.5)
    Ts = net.find_neighbor_throats(pores=Ps, mode='exclusive_or')
    trim(network=net, throats=Ts)
    a = op.utils.check_network_health(net)
    # assert len(a['disconnected_clusters']) == 3
    assert len(a['disconnected_pores']) == 20