Added "integration_example" Nek setup files

Author: Ronald Rahaman

examples/integration_example/3D_sideset.exo

@@ -0,0 +1,118 @@
+C     Dimension file to be included
+C
+C     HCUBE array dimensions
+C
+      parameter (ldim=3)
+      parameter (lx1=5,ly1=lx1,lz1=lx1,lelt=960,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-2)
+      parameter (lx3=lx1)
+      parameter (ly3=ly1)
+      parameter (lz3=lz1)
+
+      parameter (lp = 32)
+      parameter (lelg = 1000)
+c
+c     parameter (lpelv=lelv,lpelt=lelt,lpert=3)  ! perturbation
+c     parameter (lpx1=lx1,lpy1=ly1,lpz1=lz1)     ! array sizes
+c     parameter (lpx2=lx2,lpy2=ly2,lpz2=lz2)
+c
+      parameter (lpelv=1,lpelt=1,lpert=1)        ! perturbation
+      parameter (lpx1=1,lpy1=1,lpz1=1)           ! array sizes
+      parameter (lpx2=1,lpy2=1,lpz2=1)
+c
+c     parameter (lbelv=lelv,lbelt=lelt)          ! MHD
+c     parameter (lbx1=lx1,lby1=ly1,lbz1=lz1)     ! array sizes
+c     parameter (lbx2=lx2,lby2=ly2,lbz2=lz2)
+c
+      parameter (lbelv=1,lbelt=1)                ! MHD
+      parameter (lbx1=1,lby1=1,lbz1=1)           ! array sizes
+      parameter (lbx2=1,lby2=1,lbz2=1)
+ 
+C     LX1M=LX1 when there are moving meshes; =1 otherwise
+      parameter (lx1m=lx1,ly1m=ly1,lz1m=lz1)
+      parameter (ldimt= 4)                       ! 3 passive scalars + T
+      parameter (ldimt1=ldimt+1)
+      parameter (ldimt3=ldimt+3)
+c
+c     Note:  In the new code, LELGEC should be about sqrt(LELG)
+c
+      PARAMETER (LELGEC = 1)
+      PARAMETER (LXYZ2  = 1)
+      PARAMETER (LXZ21  = 1)
+ 
+      PARAMETER (LMAXV=LX1*LY1*LZ1*LELV)
+      PARAMETER (LMAXT=LX1*LY1*LZ1*LELT)
+      PARAMETER (LMAXP=LX2*LY2*LZ2*LELV)
+      PARAMETER (LXZ=LX1*LZ1)
+      PARAMETER (LORDER=3)
+      PARAMETER (MAXOBJ=4,MAXMBR=LELT*6)
+      PARAMETER (lhis=100)         ! # of pts a proc reads from hpts.in
+                                   ! Note: lhis*np > npoints in hpts.in
+C
+C     Common Block Dimensions
+C
+      PARAMETER (LCTMP0 =2*LX1*LY1*LZ1*LELT)
+      PARAMETER (LCTMP1 =4*LX1*LY1*LZ1*LELT)
+C
+C     The parameter LVEC controls whether an additional 42 field arrays
+C     are required for Steady State Solutions.  If you are not using
+C     Steady State, it is recommended that LVEC=1.
+C
+      PARAMETER (LVEC=1)
+C
+C     Uzawa projection array dimensions
+C
+      parameter (mxprev = 40)
+      parameter (lgmres = 40)
+C
+C     Split projection array dimensions
+C
+      parameter(lmvec = 1)
+      parameter(lsvec = 1)
+      parameter(lstore=lmvec*lsvec)
+c
+c     NONCONFORMING STUFF
+c
+      parameter (maxmor = lelt)
+C
+C     Array dimensions
+C
+      COMMON/DIMN/NELV,NELT,NX1,NY1,NZ1,NX2,NY2,NZ2
+     $,NX3,NY3,NZ3,NDIM,NFIELD,NPERT,NID
+     $,NXD,NYD,NZD
+
+c automatically added by makenek
+      parameter(lxo   = lx1) ! max output grid size (lxo>=lx1)
+
+c automatically added by makenek
+      parameter(lpart = 62000) ! max number of particles
+
+c automatically added by makenek
+      integer ax1,ay1,az1,ax2,ay2,az2
+      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
+
+c automatically added by makenek
+      parameter (lfdm=0)  ! == 1 for fast diagonalization method
+
+c automatically added by makenek
+      common/IOFLAG/nio  ! for logfile verbosity control
+
+c automatically added by makenek
+      integer toteq
+      parameter(toteq = 1  ) ! Number of conserved variables 
+c IFCMT=FALSE  then toteq=1
+
+c automatically added by makenek
+      integer lelcmt
+      parameter(lelcmt = 1  ) ! # of cmt elements 
+c IFCMT=FALSE  then lelcmt=1