anchor.diff:    call time_refuel%init('refuel')
     call time_change%init('change')
   
  +  call timer_list%init(16)
  +  call timer_list%append(varia)
  +  call timer_list%append(time_flag)
  +  call timer_list%append(time_refuel)
  +  call timer_list%append(time_change)
  +
     call h5open_f(error)
   
     prob = PTread_d(config,'probability')
  
positive.code: 
    integer :: i, L(3)
    integer :: j, k
    type(timer_t) :: varia, main
    type(timer_t) :: time_flag, time_refuel, time_change
    type(threefry_rng_t), allocatable :: state(:)
    integer :: n_threads
    type(h5md_file_t) :: hfile
  
negative.code:   call solvent_io%image%append(solvent%image)
    call solvent_io%species%append(solvent%species)
  
    ! Store timing data
    call timer_list%append(solvent%time_stream)
    call timer_list%append(solvent%time_step)
    call timer_list%append(solvent%time_count)
    call timer_list%append(solvent%time_sort)
    call timer_list%append(solvent%time_ct)
    call timer_list%append(solvent%time_md_pos)
    call timer_list%append(solvent%time_md_vel)
    call timer_list%append(solvent%time_max_disp)
    call timer_list%append(colloids%time_self_force)
    call timer_list%append(neigh%time_update)
    call timer_list%append(neigh%time_force)
    call timer_list%append(colloids%time_max_disp)
  
    call h5gcreate_f(hfile%id, 'timers', timers_group, error)
    call timer_list%write(timers_group, total_time)
    call h5md_write_dataset(timers_group, 'total', total_time)
    call h5md_write_dataset(timers_group, main%name, main%total)
    call h5gclose_f(timers_group, error)
  
    call janus_io%close()
    call solvent_io%close()
    call radial_hist_el%close()
    call n_solvent_el%close()
    call hfile%close()
    call h5close_f(error)
  
  contains
  
  