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:   call h5md_write_dataset(timers_group, change_timer%name, buffer_timer%total)
    call h5md_write_dataset(timers_group, buffer_timer%name, buffer_timer%total)
    call h5md_write_dataset(timers_group, neigh%time_update%name, neigh%time_update%total)
    call h5md_write_dataset(timers_group, varia%name, varia%total)
    call h5md_write_dataset(timers_group, neigh%time_force%name, neigh%time_force%total)
  
    call h5md_write_dataset(timers_group, '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 + &
         varia%total + neigh%time_force%total)
  
    call h5gclose_f(timers_group, error)
  
  