diff --git a/docs/documentation/source/ats_demos b/docs/documentation/source/ats_demos
index 8297c37a..bde4885e 160000
--- a/docs/documentation/source/ats_demos
+++ b/docs/documentation/source/ats_demos
@@ -1 +1 @@
-Subproject commit 8297c37a3e39885ea2931c3791b86b3aa68fe06f
+Subproject commit bde4885e4a0c4defd7a83027a7b03941ac4dedf8
diff --git a/tools/utils/ats_xdmf.py b/tools/utils/ats_xdmf.py
index 3c7ec4dd..afd50841 100644
--- a/tools/utils/ats_xdmf.py
+++ b/tools/utils/ats_xdmf.py
@@ -340,7 +340,7 @@ def meshElemPolygons(etype, coords, conn):
     if etype != 'HEX':
         raise RuntimeError("Only works for Hexs")
 
-    y_mean = np.array([c[1] for c in coords.values()]).mean()
+    y_mean = np.array([c[1] for c in coords]).mean()
     
     coords2 = np.array([[coords[i][0::2] for i in c[1:] if coords[i][1] > y_mean] for c in conn])
     try: