diff --git a/src/main.rs b/src/main.rs
index 2696b9ed..77aab15d 100644
--- a/src/main.rs
+++ b/src/main.rs
@@ -88,7 +88,7 @@ enum Commands {
         #[arg(long, short = 'z', value_name = "NEL")]
         nelz: Option<usize>,
 
-        /// Voxel IDs to remove from the mesh [default: 0]
+        /// Voxel IDs to remove from the mesh
         #[arg(long, short, value_name = "ID")]
         remove: Option<Vec<u8>>,
 
diff --git a/src/voxel/mod.rs b/src/voxel/mod.rs
index 58688475..b2640ec1 100644
--- a/src/voxel/mod.rs
+++ b/src/voxel/mod.rs
@@ -86,7 +86,7 @@ impl Voxels {
 fn filter_voxel_data(data: &VoxelData, remove: Option<Vec<u8>>) -> (VoxelDataSized<NSD>, Blocks) {
     #[cfg(feature = "profile")]
     let time = Instant::now();
-    let mut removed_data = remove.unwrap_or(vec![0]);
+    let mut removed_data = remove.unwrap_or_default();
     removed_data.sort();
     removed_data.dedup();
     let filtered_voxel_data_combo: VoxelDataSized<4> = data
diff --git a/src/voxel/py.rs b/src/voxel/py.rs
index b54e45ae..3375a1a8 100644
--- a/src/voxel/py.rs
+++ b/src/voxel/py.rs
@@ -19,7 +19,7 @@ pub struct Voxels {
 #[pymethods]
 impl Voxels {
     /// Converts the voxels type into a finite elements type.
-    #[pyo3(signature = (remove=[0].to_vec(), scale=[1.0, 1.0, 1.0], translate=[0.0, 0.0, 0.0]))]
+    #[pyo3(signature = (remove=[].to_vec(), scale=[1.0, 1.0, 1.0], translate=[0.0, 0.0, 0.0]))]
     pub fn as_finite_elements(
         &self,
         remove: Option<Vec<u8>>,
diff --git a/src/voxel/test.rs b/src/voxel/test.rs
index ba77b5fa..bb825b3f 100644
--- a/src/voxel/test.rs
+++ b/src/voxel/test.rs
@@ -48,7 +48,7 @@ const VOXELS_GOLD: [[usize; 3]; NUM_ELEMENTS] = [
 #[test]
 fn filter() {
     let spn = Voxels::from_npy("tests/input/letter_f_3d.npy").unwrap();
-    let (filtered_voxel_data, element_blocks) = filter_voxel_data(spn.get_data(), None);
+    let (filtered_voxel_data, element_blocks) = filter_voxel_data(spn.get_data(), Some(vec![0]));
     assert_eq!(element_blocks.len(), NUM_ELEMENTS);
     BLOCKS_GOLD
         .iter()
diff --git a/tests/fem.rs b/tests/fem.rs
index c19623b4..342c1cbf 100644
--- a/tests/fem.rs
+++ b/tests/fem.rs
@@ -39,7 +39,7 @@ fn compare_files(
 ) {
     let voxels = Voxels::from_spn(spn_path, nel).unwrap();
     let fem = voxels
-        .into_finite_elements(None, &scale, &translate)
+        .into_finite_elements(Some(vec![0]), &scale, &translate)
         .unwrap();
     fem.write_inp(file_path).unwrap();
     let mut gold = String::new();
@@ -73,7 +73,7 @@ mod read_inp {
     fn letter_f_3d() {
         let voxels = Voxels::from_spn("tests/input/letter_f_3d.spn", [4, 5, 3]).unwrap();
         let fem = voxels
-            .into_finite_elements(None, &[1.0, 1.0, 1.0], &[0.0, 0.0, 0.0])
+            .into_finite_elements(Some(vec![0]), &[1.0, 1.0, 1.0], &[0.0, 0.0, 0.0])
             .unwrap();
         fem.write_inp("target/letter_f_3d.inp").unwrap();
         let read = FiniteElements::from_inp("target/letter_f_3d.inp").unwrap();
@@ -119,7 +119,7 @@ mod write_mesh {
     fn letter_f_3d() {
         let voxels = Voxels::from_spn("tests/input/letter_f_3d.spn", [4, 5, 3]).unwrap();
         let fem = voxels
-            .into_finite_elements(None, &[1.0, 1.0, 1.0], &[0.0, 0.0, 0.0])
+            .into_finite_elements(Some(vec![0]), &[1.0, 1.0, 1.0], &[0.0, 0.0, 0.0])
             .unwrap();
         fem.write_mesh("target/letter_f_3d.mesh").unwrap();
     }
@@ -131,7 +131,7 @@ mod write_vtk {
     fn letter_f_3d() {
         let voxels = Voxels::from_spn("tests/input/letter_f_3d.spn", [4, 5, 3]).unwrap();
         let fem = voxels
-            .into_finite_elements(None, &[1.0, 1.0, 1.0], &[0.0, 0.0, 0.0])
+            .into_finite_elements(Some(vec![0]), &[1.0, 1.0, 1.0], &[0.0, 0.0, 0.0])
             .unwrap();
         fem.write_vtk("target/letter_f_3d.vtk").unwrap();
     }