anchor.diff:    colloid_sampling = PTread_i(config, 'colloid_sampling', loc=params_group)
     dt = tau / N_MD_steps
     N_loop = PTread_i(config, 'N_loop', loc=params_group)
  +  N_thermo_loop = PTread_i(config, 'N_thermo_loop', loc=params_group)
   
     N_colloids = PTread_i(config, 'N_colloids', loc=params_group)
     epsilon = PTread_d(config, 'epsilon', loc=params_group)
  
positive.code: !! \param N_MD             number MD steps occuring in tau
  !! \param colloid_sampling interval (in MD steps) of sampling the colloid position and velocity
  !! \param N_loop           number of MPCD timesteps
  !! \param N_colloids       number of colloids
  !! \param epsilon          solvent-colloid epsilon
  !! \param sigma            radius of the colloids
  
negative.code:   double precision :: e1, e2
    double precision :: temperature, local_mass, total_mass, kin_e, v_com(3)
    double precision :: tau, dt, T, T_init, T_final
    integer :: N_MD_steps, N_loop, N_thermo_loop
    integer :: colloid_sampling
    integer :: n_extra_sorting, loop_i_last_sort
    integer :: n_threads
  