anchor.diff:    call h5gclose_f(timers_group, error)
   
     call dimer_io%close()
  -  call solvent_io%close()
     call n_solvent_el%close()
     call link_angle_el%close()
     call hfile%close()
  
positive.code:         end do
          call n_solvent_el%append(n_solvent)
          call link_angle_el%append(link_angle)
       end if
  
    end do
  
negative.code:         end do
          call varia%tac()
  
       end do md_loop
  
       call solvent_cells%random_shift(state(1))
  