anchor.diff:    seed = clock + 37 * [ (i - 1, i = 1, seed_size) ]
     call random_seed(put = seed)
     deallocate(seed)
  +
  +  call h5open_f(error)
   
     L = [8, 5, 5]
   
  
positive.code: 
    call solvent_cells%count_particles(solvent% pos)
  
    print *, sum(solvent_cells%cell_count)
    print *, solvent_cells%cell_count
    print *, solvent_cells%cell_start(1), solvent_cells%cell_start(solvent_cells%N)
  
    call solvent% sort(solvent_cells)
  
  
negative.code:   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, error
    integer, allocatable :: seed(:)
  
    call random_seed(size = seed_size)
  