anchor.diff:    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)
   
  
positive.code:        n_solvent, H5MD_LINEAR, step=N_MD_steps, &
         time=N_MD_steps*dt)
  
    call h5gcreate_f(janus_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)
  
negative.code:            unit_r = rel_pos(head_pos, com_pos, solvent_cells%edges)
             unit_r = unit_r / norm2(unit_r)
             call axial_cf%add_fast((i-equilibration_loops)*N_MD_steps+j-1, v_com, unit_r)
             call janus_pos_el%append(rigid_janus%pos)
             call janus_vel_el%append(rigid_janus%vel)
          end if
  
          if ((sampling) .and. (modulo(j, colloid_sampling)==0)) then
             call janus_io%position%append(colloids%pos)
             call janus_io%velocity%append(colloids%vel)
             call janus_io%image%append(colloids%image)
             call q_el%append(rigid_janus%q)
          end if
  
          call time_flag%tic()
  