diff --git a/src/cmfd_utils.F90 b/src/cmfd_utils.F90
index 90dc1253e0..7ec2c1a264 100644
--- a/src/cmfd_utils.F90
+++ b/src/cmfd_utils.F90
@@ -519,9 +519,23 @@ contains
 
      character(LEN=8), parameter :: filename = "filef.h5" ! File name
      character(LEN=4), parameter :: grpname = "cmfd"      ! Group name
+     character(LEN=9), parameter :: dsetname = "cmfd/dset" ! dataset name
      integer(HID_T) :: file_id                            ! File identifier
-     integer(HID_T) :: group_id                           ! Group identifieer
+     integer(HID_T) :: group_id                           ! Group identifier
+     integer(HID_T) :: dataspace_id                       ! Data space identifier
+     integer(HID_T) :: dataset_id                         ! Dataset identifier
      integer        :: error                              ! Error flag
+     integer        :: i,j
+     integer, dimension(3,3) :: dset_data
+     integer(HSIZE_T), dimension(2) :: dims = (/3,3/)
+     integer :: rank = 2
+
+     ! Create example dataset
+     do i = 1,3
+       do j = 1,3
+         dset_data(i,j) = j
+       end do
+     end do
 
      ! Initialize FORTRAN interface.
      call h5open_f(error)
@@ -532,6 +546,21 @@ contains
      ! Create the CMFD group
      call h5gcreate_f(file_id, grpname, group_id, error)
 
+     ! Create the data space for the dataset
+     call h5screate_simple_f(rank, dims, dataspace_id, error)
+
+     ! Create the data set in CMFD group
+     call h5dcreate_f(file_id, dsetname, H5T_NATIVE_INTEGER,dataspace_id,dataset_id,error)
+
+     ! Write the dataset
+     call h5dwrite_f(dataset_id, H5T_NATIVE_INTEGER, dset_data, dims, error)
+
+     ! Close the data space
+     call h5sclose_f(dataspace_id, error)
+
+     ! Close the dataset
+     call h5dclose_f(dataset_id, error)
+ 
      ! Close the CMFD group
      call h5gclose_f(group_id, error)