anchor.diff:    !! Histogram variable
     type(histogram_t), allocatable :: radial_hist(:)
     double precision, allocatable :: dummy_hist(:,:,:)
  -
  -  integer, parameter :: block_length = 8
  -  type(axial_correlator_t) :: axial_cf
  -  integer(HID_T) :: correlator_group
   
     integer :: equilibration_loops
     integer :: colloid_sampling
  
positive.code:   call h5gclose_f(params_group, error)
    call PTkill(config)
  
    call axial_cf%init(block_length, N_loop, N_loop*N_MD_steps)
  
    do enzyme_i = 1, N_enzymes
       colloids% species(3*(enzyme_i-1)+1) = 2
       colloids% species(3*(enzyme_i-1)+2) = 1
  
negative.code:   call dummy_element%create_fixed(box_group, 'edges', solvent_cells%edges)
    call h5gclose_f(box_group, error)
  
    call solvent% random_placement(solvent_cells% edges, colloids, solvent_colloid_lj, state(1))
  
    call solvent% sort(solvent_cells)
  