c-----------------------------------------------------------------------
c  nek5000 user-file template
c
c  user specified routines:
c     - uservp  : variable properties
c     - userf   : local acceleration term for fluid
c     - userq   : local source term for scalars
c     - userbc  : boundary conditions
c     - useric  : initial conditions
c     - userchk : general purpose routine for checking errors etc.
c     - userqtl : thermal divergence for lowMach number flows 
c     - usrdat  : modify element vertices 
c     - usrdat2 : modify mesh coordinates
c     - usrdat3 : general purpose routine for initialization
c     
c-----------------------------------------------------------------------
      subroutine uservp (ix,iy,iz,eg)

      integer ix,iy,iz,eg
     
      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'

      integer e
c     e = gllel(eg)

      udiff  = 0.0
      utrans = 0.0

      return
      end
c-----------------------------------------------------------------------
      subroutine userf(ix,iy,iz,eg)

      integer ix,iy,iz,eg

      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'

      integer e
c     e = gllel(eg)
      
      ffx = 0.0
      ffy = 0.0
      ffz = temp*param(35)

      return
      end
c-----------------------------------------------------------------------
      subroutine userq  (ix,iy,iz,eg)

      integer ix,iy,iz,eg

      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'

      integer e
c     e = gllel(eg)

      qvol   = 0.0

      return
      end
c-----------------------------------------------------------------------
      subroutine userbc(ix,iy,iz,iside,eg)

      integer ix,iy,iz,iside,eg

      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'

      ux   = 0.0
      uy   = 0.0
      uz   = 0.0
      temp = 0.0
      
      if(x.eq.0)	temp =  0.5
      if(x.eq.0.25)	temp = -0.5
      
      return
      end
c-----------------------------------------------------------------------
      subroutine useric(ix,iy,iz,eg)

      integer ix,iy,iz,eg

      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'

      ux	= 0.0
      uy	= 0.0
      uz	= 0.0
      temp	= 0.0

      return
      end
c-----------------------------------------------------------------------
      subroutine userchk()
      
      include 'SIZE'
      include 'TOTAL'
      
      return
      end
c-----------------------------------------------------------------------
      subroutine userqtl ! Set thermal divergence
      call userqtl_scig 

      return
      end
c-----------------------------------------------------------------------
      subroutine usrdat()

      include 'SIZE'
      include 'TOTAL'

      return
      end
c-----------------------------------------------------------------------
      subroutine usrdat2()

      include 'SIZE'
      include 'TOTAL'
      
      param(2) = param(35)/sqrt(param(34))
      param(8) = 1./sqrt(param(34))
      
      return
      end
c-----------------------------------------------------------------------
      subroutine usrdat3()

      include 'SIZE'
      include 'TOTAL'

      end