anchor.diff:    call flag_timer%init('flag')
     call change_timer%init('change')
     call buffer_timer%init('buffer')
  +  call varia%init('varia')
   
     n_threads = omp_get_max_threads()
     allocate(state(n_threads))
  
positive.code:         call vx% reset()
       end if
  
       call compute_rho_xy
       call rho_xy_el%append(rho_xy)
  
       call solvent% sort(solvent_cells)
  
negative.code:   integer :: i, L(3),  n_threads
    integer :: j, k, m
    
    type(timer_t) :: flag_timer, change_timer, buffer_timer, varia
    integer(HID_T) :: timers_group
  
    integer, allocatable :: rho_xy(:,:,:)
  