anchor.diff:    call flag_timer%init('flag')
     call change_timer%init('change')
     call buffer_timer%init('buffer')
  +  call randomize_timer%init('randomize')
     call varia%init('varia')
   
     n_threads = omp_get_max_threads()
  
positive.code: 
    call neigh% update_list(colloids, solvent, max_cut+skin, solvent_cells)
  
    e1 = compute_force(colloids, solvent, neigh, solvent_cells% edges, solvent_colloid_lj)
    e2 = compute_force_n2(colloids, solvent_cells% edges, colloid_lj)
    e_wall = lj93_zwall(colloids, solvent_cells% edges, walls_colloid_lj, 3)
  
negative.code:   integer :: i, L(3),  n_threads
    integer :: j, k, m
    
    type(timer_t) :: flag_timer, change_timer, buffer_timer, varia, randomize_timer
    integer(HID_T) :: timers_group
  
    integer, allocatable :: rho_xy(:,:,:)
  