anchor.diff:    write(*,*) '    e co so     |   e co co     |   kin co      |   kin so      |   total       |   temp        |'
     write(*,*) ''
   
  -  do i = 1, 100
  +  do i = 1, 350
        so_max = 0
        co_max = 0
        md: do j = 1, N_MD_steps
  
positive.code: program setup_simple_colloids
    use common
    use cell_system
    use particle_system
  
negative.code: 
    call h5close_f(error)
  
  end program setup_simple_colloid
  