anchor.diff:          call n_minus_el%append(n_minus)
           call link_angle_el%append(link_angle)
           if (modulo(i, colloid_sampling)==0) then
  -           call dimer_io%position%append(colloids%pos)
  -           call dimer_io%velocity%append(colloids%vel)
  -           call dimer_io%image%append(colloids%image)
  +           call enzyme_io%position%append(colloids%pos)
  +           call enzyme_io%velocity%append(colloids%vel)
  +           call enzyme_io%image%append(colloids%image)
           end if
   
        end if
  
positive.code:        link_angle, ior(H5MD_LINEAR, H5MD_STORE_TIME), step=N_MD_steps, &
         time=N_MD_steps*dt)
  
    call h5gcreate_f(dimer_io%group, 'box', box_group, error)
    call h5md_write_attribute(box_group, 'dimension', 3)
    call dummy_element%create_fixed(box_group, 'edges', solvent_cells%edges)
    call h5gclose_f(box_group, error)
  
negative.code:   call h5md_write_dataset(timers_group, main%name, main%total)
    call h5gclose_f(timers_group, error)
  
    call enzyme_io%close()
    call n_solvent_el%close()
    call n_plus_el%close()
    call n_minus_el%close()
  