diff --git a/scripts/2-rom/5-hpod/hpod_run.py b/scripts/2-rom/5-hpod/hpod_run.py
index 705af56..b296aa6 100644
--- a/scripts/2-rom/5-hpod/hpod_run.py
+++ b/scripts/2-rom/5-hpod/hpod_run.py
@@ -178,6 +178,7 @@ def hilbert_axis0_mirror_right(a, axis=0, return_analytic=False):
 
     # Initialize coef to get the mass matrix
     bm = meshgrid["mass"]
+    bm[bm < 1e-8] = 1e-8 # Avoid zero mass entries
 
     # Instance io helper that will serve as buffer for the snapshots
     ioh = IoHelp(comm, number_of_fields = number_of_pod_fields, batch_size = pod_batch_size, field_size = bm.size, field_data_type=pod_dtype, mass_matrix_data_type=pod_dtype)
diff --git a/scripts/2-rom/5-hpod/inputs.json b/scripts/2-rom/5-hpod/inputs.json
index be1fecc..5b6ac75 100644
--- a/scripts/2-rom/5-hpod/inputs.json
+++ b/scripts/2-rom/5-hpod/inputs.json
@@ -4,7 +4,7 @@
 "batch_size"         : 20,
 "keep_modes"         : 20,
 "write_modes"        : 20,
-"fields"             : ["u","v"],
+"fields"             : ["u","v", "w"],
 "dtype"              : "double",
 "backend"            : "numpy",
 "hpod"               : true,