diff --git a/fd/spatial.py b/fd/spatial.py
index 6b147fe966275b9f675c8b4a3d5b3994e4f7f184..43e88eaea0a200e71ed4e4dff0e649b4a5d869e3 100644
--- a/fd/spatial.py
+++ b/fd/spatial.py
@@ -86,6 +86,8 @@ def discretize_spatial(expr, dx, stencil=fd_stencils_standard):
             stencil = fd_stencils_standard
         elif stencil == 'isotropic':
             stencil = fd_stencils_isotropic
+        elif stencil == 'isotropic_hd':
+            stencil = fd_stencils_isotropic_high_density_code
         else:
             raise ValueError("Unknown stencil. Supported 'standard' and 'isotropic'")