c-----------------------------------------------------------------
      double precision function HNose_nh(eki,e)
      implicit real*8 (a-h,o-z)
      include 'DIMENSIONS'
      include 'COMMON.MD'
      HNose_nh=eki+e+dimen*Rb*t_bath*xlogs(1)+qmass(1)*vlogs(1)**2/2
      do i=2,nnos
        HNose_nh=HNose_nh+qmass(i)*vlogs(i)**2/2+Rb*t_bath*xlogs(i)
      enddo
c      write(4,'(5e15.5)') 
c     &       vlogs(1),xlogs(1),HNose,eki,e
      return
      end