anchor.diff:    call e_solvent% close()
     call e_solvent_v% close()
     call elem_tz% close()
  +  call elem_rhoz% close()
  +  call elem_T% close()
   
     call h5gclose_f(solvent_group, error)
   
  
positive.code: 
    call datafile% create('data_setup_simple_fluid.h5', 'RMPCDMD:setup_simple_fluid', '0.0 dev', 'Pierre de Buyl')
  
    call tz% init(&
         solvent_cells% origin(3), &
         solvent_cells% origin(3) + solvent_cells% L(3) * solvent_cells% a, &
         solvent_cells% L(3))
    call h5gcreate_f(datafile% id, 'observables', datafile% observables, error)
  
    call h5gcreate_f(datafile% particles, 'solvent', solvent_group, error)
  
negative.code:   type(particle_system_t) :: solvent
  
    type(profile_t) :: tz
    type(histogram_t) :: rhoz
  
    type(h5md_file_t) :: datafile
    type(h5md_element_t) :: elem
    type(h5md_element_t) :: e_solvent, e_solvent_v
    type(h5md_element_t) :: elem_tz, elem_tz_count
    type(h5md_element_t) :: elem_rhoz
    type(h5md_element_t) :: elem_T
    integer(HID_T) :: box_group, solvent_group
  
    integer, parameter :: N = 10800
  