anchor.diff:    !e_wall = colloid_wall_interaction(colloids, walls_colloid_lj,solvent_cells% edges)
     solvent% force_old = solvent% force
     colloids% force_old = colloids% force
  -  write(*,*) colloids% force
  -  write(*,*) ''
  -  write(*,*) '    i           |    e co so     |   e co co     |   kin co      |   kin so      |   total       |   temp        |'
  -  write(*,*) ''
   
     call vx% init(0.d0, solvent_cells% edges(3), L(3))
   
     i = 0
  -  kin_co = (mass(1)*sum(colloids% vel(:,1)**2)+mass(2)*sum(colloids% vel(:,2)**2))/2
  -  call thermo_write
  +
  +  write(*,*) 'Running for', N_loop, 'loops'
     !start RMPCDMD
     setup: do i = 1, N_loop
  +     if (modulo(i,20) == 0) write(*,'(i09)',advance='no') i
        md: do j = 1, N_MD_steps
           call mpcd_stream_zwall_light(solvent, solvent_cells, dt,g)
   
  
positive.code:   double precision :: skin, co_max, so_max
    integer :: N_MD_steps, N_loop
    integer :: n_extra_sorting
    double precision :: kin_co
  
    double precision :: conc_z(400)
    double precision :: colloid_pos(3,2)
  
    type(PTo) :: config
    integer(c_int64_t) :: seed
  
negative.code:   write(*,'(a16,f8.3)') buffer_timer%name, buffer_timer%total
    write(*,'(a16,f8.3)') neigh%time_update%name, neigh%time_update%total
    write(*,'(a16,f8.3)') neigh%time_force%name, neigh%time_force%total
  
    write(*,'(a16,f8.3)') 'total                          ', &
         solvent%time_stream%total + solvent%time_step%total + solvent%time_count%total +&
         solvent%time_sort%total + solvent%time_ct%total + solvent%time_md_vel%total + solvent%time_max_disp%total + &
         flag_timer%total + change_timer%total + buffer_timer%total + neigh%time_update%total + neigh%time_force%total
    
  contains
  
    subroutine flag_particles
    double precision :: dist_to_C_sq
  