anchor.diff:         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 h5gcreate_f(enzyme_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)
  
positive.code:         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)
          end if
  
       end if
  
negative.code:   double precision :: temperature, kin_e
    double precision :: v_com(3)
    double precision :: tmp_vec(3), normal_vec(3), tmp_q(4)
    type(particle_system_io_t) :: enzyme_io
    double precision :: bulk_rate(2)
    logical :: bulk_rmpcd
  
  