c------------------------------------------------------
      double precision function HNose(ek,s,e,pi,Q,t_bath,dimen)
      implicit none
      double precision ek,s,e,pi,Q,t_bath,Rb
      integer dimen
      Rb=0.001986d0
      HNose=ek+e+pi**2/(2*Q)+dimen*Rb*t_bath*log(s)
c      print '(6f15.5,i5,a2,2f15.5)',ek,s,e,pi,Q,t_bath,dimen,"--",
c     &      pi**2/(2*Q),dimen*Rb*t_bath*log(s)
      return
      end