anchor.diff:             n_solvent(j) = n_solvent(j) + 1
           end do
           call n_solvent_el%append(n_solvent)
  +        call n_plus_el%append(n_plus)
  +        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)
  
positive.code: 
    call dimer_io%close()
    call n_solvent_el%close()
    call link_angle_el%close()
    call hfile%close()
    call h5close_f(error)
  
negative.code:        n_solvent, ior(H5MD_LINEAR, H5MD_STORE_TIME), step=N_MD_steps, &
         time=N_MD_steps*dt)
  
    call n_plus_el%create_time(hfile%observables, 'n_plus', &
         n_plus, ior(H5MD_LINEAR, H5MD_STORE_TIME), step=N_MD_steps, &
         time=N_MD_steps*dt)
  
    call n_minus_el%create_time(hfile%observables, 'n_minus', &
         n_minus, ior(H5MD_LINEAR, H5MD_STORE_TIME), step=N_MD_steps, &
         time=N_MD_steps*dt)
  
    call link_angle_el%create_time(hfile%observables, 'link_angle', &
         link_angle, ior(H5MD_LINEAR, H5MD_STORE_TIME), step=N_MD_steps, &
         time=N_MD_steps*dt)
  