anchor.diff:    call time_change%init('change')
     call bulk_reac_timer%init('bulk_reac')
   
  -  call timer_list%init(20)
  +  call timer_list%init(22)
     call timer_list%append(varia)
     call timer_list%append(time_flag)
     call timer_list%append(time_refuel)
  
positive.code:   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(colloids%time_rattle_pos)
    call timer_list%append(colloids%time_rattle_vel)
    call timer_list%append(colloids%time_elastic)
    call timer_list%append(neigh%time_update)
    call timer_list%append(neigh%time_force)
    call timer_list%append(colloids%time_max_disp)
  
negative.code: 
          call thermo_data%append(hfile, temperature, e1+e2+e3, kin_e, e1+e2+e3+kin_e, v_com)
          call axial_cf%add(i-equilibration_loops, com_pos, unit_r)
          call varia%tac()
       end if
  
       call solvent_cells%random_shift(state(1))
  
       call apply_pbc(solvent, solvent_cells% edges)
       call apply_pbc(colloids, solvent_cells% edges)
  