!     Heating of the material in a less violent manner then
!     just by velocities distribution (as in Vel.f [KFLAG=2])
!     MSE 6270, Leonid Zhigilei

      SUBROUTINE HEATING(TEMPTR,HDONE)
      INCLUDE 'common.h'
      Parameter (SKF=0.001d0)
      IF(TEMPTR.GE.QTEM) THEN
        HDONE=1
        RETURN
      ELSE
        SC=SQRT(1.0d0+SKF*(QTEM-TEMPTR)/(QTEM))
        t_loop: DO I=1,NAN
!         Sometimes you may want to heat up only a part of your system
!         IF(XD(3,I).GT.(ZCENTR+ZL/4.0d0)) CYCLE t_loop
!         IF(KHIST(J).NE.2) CYCLE t_loop
          Q1D(1,I)=Q1D(1,I)*SC
          Q1D(2,I)=Q1D(2,I)*SC
          IF(NDIM.EQ.3) Q1D(3,I)=Q1D(3,I)*SC
        End do t_loop
      ENDIF

      RETURN
      END