anchor.diff:    type(h5md_file_t) :: hfile
     type(h5md_element_t) :: dummy_element
     integer(HID_T) :: fields_group, params_group
  +  integer(HID_T) :: kinetics_group
     integer(HID_T) :: box_group
     integer(HID_T) :: dummy_id
     type(thermo_t) :: thermo_data
  
positive.code:     colloids%mass(enz_2) = colloids%mass(enz_2) + 1
  
      bound_molecule_id(enzyme_idx) = solvent%id(idx)
      solvent%species(idx) = 0
  
      ! compute conformational energy difference when changing the angle
  
negative.code:   type(histogram_t), allocatable :: radial_hist(:)
    double precision, allocatable :: dummy_hist(:,:,:)
  
    !! Kinetics data
    type(enzyme_kinetics_t), allocatable :: kinetics_data(:)
  
    integer :: equilibration_loops
    integer :: colloid_sampling
    logical :: sampling
  