diff --git a/las_manipulation.py b/las_manipulation.py
index 2739730..2cb8d7b 100644
--- a/las_manipulation.py
+++ b/las_manipulation.py
@@ -159,7 +159,11 @@ def find_beach_slope(chain_list, z_ave_list, ref_index, check_value):
     while i>0 and z_ave_list[i]<check_value+1:
         i=i-1
 
+    # Hide numpy divide-by-zero warnings
+    np.seterr(all='ignore')
     slope=(z_ave_list[i]-z_ave_list[ref_index])/(chain_list[i]-chain_list[ref_index])
+    np.seterr(all=None)
+
 
     beach_ave=interpolate.interp1d([min(chain_list),max(chain_list)], [(min(chain_list)-chain_list[ref_index])*slope+z_ave_list[ref_index], (z_ave_list[ref_index]-(chain_list[ref_index]-max(chain_list))*slope)])