def METHOD_NAME(self):
    fd = op.algorithms.FickianDiffusion(network=self.pn, phase=self.air)
    # Remove given type from given pores
    fd['pore.bc.rate'][1] = 1.0
    fd['pore.bc.value'][0] = 1.0
    fd.set_value_BC(pores=[0, 1], mode='remove')
    mask = np.isfinite(fd['pore.bc.value'])
    assert mask.sum() == 0
    assert np.isfinite(fd['pore.bc.rate']).sum() == 1
    # Remove given type from all pores
    fd['pore.bc.rate'][1] = 1.0
    fd['pore.bc.value'][0] = 1.0
    fd.set_BC(pores=None, bctype='rate', mode='remove')
    assert np.isfinite(fd['pore.bc.value']).sum() == 1
    assert np.isfinite(fd['pore.bc.rate']).sum() == 0
    # Remove all types from all pores
    fd['pore.bc.rate'][1] = 1.0
    fd['pore.bc.value'][0] = 1.0
    fd.set_BC(bctype=fd['pore.bc'].keys(), mode='remove')
    assert np.isfinite(fd['pore.bc.value']).sum() == 0
    assert np.isfinite(fd['pore.bc.rate']).sum() == 0
    # Remove all types from all pores given bctype is []
    fd['pore.bc.rate'][1] = 1.0
    fd['pore.bc.value'][0] = 1.0
    fd.set_BC(mode='remove')
    assert np.isfinite(fd['pore.bc.value']).sum() == 0
    assert np.isfinite(fd['pore.bc.rate']).sum() == 0
    # Remove all types from given pores
    fd['pore.bc.rate'][1] = 1.0
    fd['pore.bc.value'][0] = 1.0
    fd.set_BC(pores=[0, 1], mode='remove')
    assert np.isfinite(fd['pore.bc.value']).sum() == 0
    assert np.isfinite(fd['pore.bc.rate']).sum() == 0