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:   use hilbert
    use neighbor_list
    use interaction
    implicit none
  
    type(cell_system_t) :: solvent_cells
  
negative.code:   end do
    close(12)
  
    call h5gclose_f(colloids_group, error)
  
    call h5close_f(error)
  
  end program try_all
  