anchor.diff:         n_solvent, H5MD_LINEAR, step=N_MD_steps, &
          time=N_MD_steps*dt)
   
  +  call q_el%create_time(hfile%observables, 'q', &
  +       rigid_janus%q, H5MD_LINEAR, step=1, &
  +       time=dt)
  +
  +  call janus_pos_el%create_time(hfile%observables, 'janus_pos', &
  +       rigid_janus%pos, H5MD_LINEAR, step=1, &
  +       time=dt)
  +
  +  call janus_vel_el%create_time(hfile%observables, 'janus_vel', &
  +       rigid_janus%vel, H5MD_LINEAR, step=1, &
  +       time=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)
  
positive.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)
          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)
          end if
  
          call time_flag%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)
  
  