+++ /dev/null
-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