adding rod bundle case

examples/coupled_heat_conduction/SIZE

@@ -2,22 +2,22 @@ C     Dimension file to be included
 C
 C     HCUBE array dimensions
 C
-      parameter (ldim=2)
-      parameter (lx1=8,ly1=lx1,lz1=1,lelt=300,lelv=lelt)
-      parameter (lxd=12,lyd=lxd,lzd=1)
-      parameter (lelx=20,lely=20,lelz=1)
+      parameter (ldim=3)
+      parameter (lx1=5,ly1=lx1,lz1=lx1,lelt=240,lelv=lelt)
+      parameter (lxd=8,lyd=lxd,lzd=lxd)
+      parameter (lelx=1,lely=1,lelz=1)
  
       parameter (lzl=3 + 2*(ldim-3))
  
       parameter (lx2=lx1-2)
       parameter (ly2=ly1-2)
-      parameter (lz2=lz1  )
+      parameter (lz2=lz1-2)
       parameter (lx3=lx1)
       parameter (ly3=ly1)
       parameter (lz3=lz1)
 
-      parameter (lp = 512)
-      parameter (lelg = 4100)
+      parameter (lp = 32)
+      parameter (lelg = 500)
 c
 c     parameter (lpelv=lelv,lpelt=lelt,lpert=3)  ! perturbation
 c     parameter (lpx1=lx1,lpy1=ly1,lpz1=lz1)     ! array sizes
@@ -36,8 +36,8 @@ c
       parameter (lbx2=1,lby2=1,lbz2=1)
  
 C     LX1M=LX1 when there are moving meshes; =1 otherwise
-      parameter (lx1m=1,ly1m=1,lz1m=1)
-      parameter (ldimt= 2)                       ! 2 passive scalars + T
+      parameter (lx1m=lx1,ly1m=ly1,lz1m=lz1)
+      parameter (ldimt= 4)                       ! 3 passive scalars + T
       parameter (ldimt1=ldimt+1)
       parameter (ldimt3=ldimt+3)
 c
@@ -69,8 +69,8 @@ C
 C
 C     Uzawa projection array dimensions
 C
-      parameter (mxprev = 20)
-      parameter (lgmres = 30)
+      parameter (mxprev = 40)
+      parameter (lgmres = 40)
 C
 C     Split projection array dimensions
 C
@@ -92,11 +92,11 @@ c automatically added by makenek
       parameter(lxo   = lx1) ! max output grid size (lxo>=lx1)
 
 c automatically added by makenek
-      parameter(lpart = 1  ) ! max number of particles
+      parameter(lpart = 62000) ! max number of particles
 
 c automatically added by makenek
       integer ax1,ay1,az1,ax2,ay2,az2
-      parameter (ax1=lx1,ay1=ly1,az1=lz1,ax2=lx2,ay2=ly2,az2=lz2) ! running averages
+      parameter (ax1=1,ay1=1,az1=1,ax2=1,ay2=1,az2=1) ! running averages
 
 c automatically added by makenek
       parameter (lxs=1,lys=lxs,lzs=(lxs-1)*(ldim-2)+1) !New Pressure Preconditioner

examples/coupled_heat_conduction/rod12.usr

+c-----------------------------------------------------------------------
+C
+C  USER SPECIFIED ROUTINES:
+C
+C     - boundary conditions
+C     - initial conditions
+C     - variable properties
+C     - local acceleration for fluid (a)
+C     - forcing function for passive scalar (q)
+C     - general purpose routine for checking errors etc.
+C
+c-----------------------------------------------------------------------
+      subroutine uservp (ix,iy,iz,eg)
+      include 'SIZE'
+      include 'TOTAL'
+      include 'NEKUSE'
+
+      integer e,f,eg
+c     e = gllel(eg)
+
+      udiff =0.
+      utrans=0.
+      return
+      end
+c-----------------------------------------------------------------------
+      subroutine userf  (ix,iy,iz,eg)
+      include 'SIZE'
+      include 'TOTAL'
+      include 'NEKUSE'
+
+      integer e,f,eg
+c     e = gllel(eg)
+
+
+c     Note: this is an acceleration term, NOT a force!
+c     Thus, ffx will subsequently be multiplied by rho(x,t).
+      ffx=0.0
+      ffy=0.0
+      ffz=0.0     
+
+      return
+      end
+c-----------------------------------------------------------------------
+      subroutine userq  (ix,iy,iz,eg)
+      include 'SIZE'
+      include 'TOTAL'
+      include 'NEKUSE'
+
+      integer e,f,eg
+c     e = gllel(eg)
+
+      qvol   = 0.0
+
+      return
+      end
+c-----------------------------------------------------------------------
+      subroutine userchk
+      include 'SIZE'
+      include 'TOTAL'
+
+      n1=nelt*lx1*ly1*lz1
+      n2=nelt*lx2*ly2*lz2
+
+      pmin=glmin(pr,n2)
+      pmax=glmax(pr,n2)
+
+      wmin=glmin(vz,n1)
+      wmax=glmax(vz,n1)
+
+      tmin=glmin(t,n1)
+      tmax=glmax(t,n1)
+
+      ifflow=.false. 
+      if (nid.eq.0) then
+         write(6,*)"*** Pressure: ",pmin," - ",pmax        
+         write(6,*)"*** Velocity: ",wmin," - ",wmax
+         write(6,*)"*** Temperature: ",tmin," - ",tmax
+      endif
+
+      return
+      end
+c-----------------------------------------------------------------------
+      subroutine userbc (ix,iy,iz,iside,ieg)
+      include 'SIZE'
+      include 'TOTAL'
+      include 'NEKUSE'
+
+      integer e,ieg
+      real ucx, ucy, ucz, ucy_e, yy 
+
+      e=gllel(ieg)
+
+      zz  = zm1(ix,iy,iz,e)
+      xx  = xm1(ix,iy,iz,e)
+      yy  = ym1(ix,iy,iz,e)
+      rr  = sqrt(xx**2 + yy**2)
+
+      ux = 0.0
+      uy = 0.0
+      uz = 1.0
+      temp = 0.0
+      flux = 1.0
+
+      return
+      end
+c-----------------------------------------------------------------------
+      subroutine useric (ix,iy,iz,ieg)
+      include 'SIZE'
+      include 'TOTAL'
+      include 'NEKUSE'
+
+      integer idum, e
+      save    idum
+      data    idum / 0 /
+
+      if (idum.eq.0) idum = 99 + nid
+      eps = .35
+
+      uz=0.0
+      uy=0.0
+      ux=0.0
+      temp=0
+
+      return
+      end
+c-----------------------------------------------------------------------
+      subroutine usrdat
+      include 'SIZE'
+      include 'TOTAL'
+
+      return
+      end
+c-----------------------------------------------------------------------
+      subroutine usrdat2
+      include 'SIZE'
+      include 'TOTAL'
+
+      param(66) = 4.   ! These give the std nek binary i/o and are 
+      param(67) = 4.   ! good default values
+      ifuservp=.false.
+
+      return
+      end
+c-----------------------------------------------------------------------
+      subroutine usrdat3
+      include 'SIZE'
+      include 'TOTAL'
+c
+      return
+      end
+