anchor.diff:      dims(1)     = 0
       max_dims(1) = H5S_UNLIMITED_F
       call h5screate_simple_f(rank, dims, file_s, h5_error, max_dims)
  -    chunk_dims(1) = 1
  +    chunk_dims(1) = 128
       call h5pcreate_f(H5P_DATASET_CREATE_F, plist, h5_error)
       call h5pset_chunk_f(plist, rank, chunk_dims, h5_error)
       call h5dcreate_f(g_id, 'samples', H5T_NATIVE_DOUBLE, file_s, ID% d_id, h5_error, plist)
  
positive.code:     dims(1)     = 0
      max_dims(1) = H5S_UNLIMITED_F
      call h5screate_simple_f(rank, dims, file_s, h5_error, max_dims)
      chunk_dims(1) = 1
      call h5pcreate_f(H5P_DATASET_CREATE_F, plist, h5_error)
      call h5pset_chunk_f(plist, rank, chunk_dims, h5_error)
      call h5dcreate_f(g_id, 'samples', H5T_NATIVE_INTEGER, file_s, ID% d_id, h5_error, plist)
  
negative.code:     dims(3)     = 0
      max_dims(3) = H5S_UNLIMITED_F
      call h5screate_simple_f(rank, dims, file_s, h5_error, max_dims)
      chunk_dims(3) = 128
      call h5pcreate_f(H5P_DATASET_CREATE_F, plist, h5_error)
      call h5pset_chunk_f(plist, rank, chunk_dims, h5_error)
      call h5dcreate_f(g_id, 'samples', H5T_NATIVE_DOUBLE, file_s, ID% d_id, h5_error, plist)
  