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