anchor.diff:    !! dimension
     !! @param link_from is the name of another trajectory from which the time can be copied
     !! @param compress Switch to toggle SZIP compression.
  -  subroutine h5md_add_trajectory_data(file_id, trajectory_name, N, D, ID, group_name, species_react, link_from, compress)
  +  !! @param force_kind Force the supplied kind, "integer" or "double", for the dataset.
  +  subroutine h5md_add_trajectory_data(file_id, trajectory_name, N, D, ID, group_name, species_react, &
  +       link_from, compress, force_kind)
       integer(HID_T), intent(in) :: file_id
       character(len=*), intent(in) :: trajectory_name
       integer, intent(in) :: N, D
  
positive.code:        call h5pset_szip_f(plist, H5_SZIP_NN_OM_F, 8, h5_error)
      end if
      call h5pset_chunk_f(plist, rank, chunk_dims, h5_error)
      if (trajectory_name .ne. 'species' .and. trajectory_name .ne. 'jumps') then
         call h5dcreate_f(g_id, 'coordinates', H5T_NATIVE_DOUBLE, s_id, ID% d_id, h5_error, plist)
      else
         call h5dcreate_f(g_id, 'coordinates', H5T_NATIVE_INTEGER, s_id, ID% d_id, h5_error, plist)
      end if
      call h5pclose_f(plist, h5_error)
      call h5sclose_f(s_id, h5_error)
  
negative.code:     logical, intent(in), optional :: species_react
      character(len=*), intent(in), optional :: link_from
      logical, intent(in), optional :: compress
      character(len=*), intent(in), optional :: force_kind
      
      character(len=128) :: path
      integer :: rank
  