anchor.diff:    type(lj_params_t) :: solvent_colloid_lj
     type(lj_params_t) :: colloid_lj
   
  +  type(h5md_file_t) :: datafile
  +  type(h5md_element_t) :: elem
  +  integer(HID_T) :: colloids_group, box_group
  +
     integer, parameter :: N = 1000
     integer, parameter :: N_colloids = 3
   
     double precision :: epsilon, sigma, sigma_cut
   
  -  integer :: i, L(3), seed_size, clock
  +  integer :: i, L(3), seed_size, clock, error
     integer, allocatable :: seed(:)
   
     call random_seed(size = seed_size)
  
positive.code:   end do
    close(12)
  
  end program try_all
  
negative.code:   seed = clock + 37 * [ (i - 1, i = 1, seed_size) ]
    call random_seed(put = seed)
    deallocate(seed)
  
    call h5open_f(error)
  
    L = [8, 5, 5]
  
  