anchor.diff:    type(h5md_element_t) :: dummy_element
     integer(HID_T) :: fields_group, params_group
     integer(HID_T) :: box_group
  +  integer(HID_T) :: dummy_id
     type(thermo_t) :: thermo_data
     double precision :: temperature, kin_e
     double precision :: v_com(3)
  
positive.code:   call solvent_io%velocity%append(solvent%vel)
    call solvent_io%image%append(solvent%image)
    call solvent_io%species%append(solvent%species)
  
    ! Store timing data
    call timer_list%append(solvent%time_stream)
  
negative.code: 
    call neigh% update_list(colloids, solvent, max_cut+skin, solvent_cells, solvent_colloid_lj)
  
    e1 = compute_force(colloids, solvent, neigh, solvent_cells% edges, solvent_colloid_lj)
    e2 = compute_force_n2(colloids, solvent_cells% edges, colloid_lj)
    e3 = compute_bead_force()
  