diff --git a/sources/io.F90 b/sources/io.F90
index f3002cb..c6a2c02 100644
--- a/sources/io.F90
+++ b/sources/io.F90
@@ -3904,7 +3904,7 @@ module io
 #if NDIMS == 3
         l         = rank(faces)
         dims(1:l) = shape(faces)
-        call store_dataset_h5(grp_id, 'faces', &
+        call read_dataset_h5(grp_id, 'faces', &
                             H5T_NATIVE_INTEGER, dims(1:l), faces, status)
 #endif /* NDIMS == 3 */
         l         = rank(edges)