anchor.diff:    integer, dimension(N_species) :: n_solvent
     type(h5md_element_t) :: n_solvent_el
   
  +  type(h5md_element_t) :: q_el, janus_pos_el, janus_vel_el
  +
     type(histogram_t) :: radial_hist
     type(h5md_element_t) :: radial_hist_el
     double precision :: polar_r_max
  
positive.code:   end if
    solvent% force_old = solvent% force
    colloids% force_old = colloids% force
    if (do_quaternion) call rigid_janus%compute_force_torque(colloids, solvent_cells%edges)
  
    write(*,*) 'Running for', equilibration_loops, '+', N_loop, 'loops'
    call main%tic()
  
negative.code:   call solvent_io%close()
    call radial_hist_el%close()
    call n_solvent_el%close()
    call q_el%close()
    call janus_pos_el%close()
    call janus_vel_el%close()
    call hfile%close()
    call h5close_f(error)
  
  