1      subroutine shell(Re,Pr,nx,ny,nz,
2     $nuim,nuex2,nuex4,cfl,scheme,conf,ni,maxit)
3      real*8 q(5,nx,ny,nz),dq(5,nx,ny,nz),rhs(5,nx,ny,nz),e(5,nx,ny,nz),
4     1     f(5,nx,ny,nz),g(5,nx,ny,nz),ev(5,nx,ny,nz),fv(5,nx,ny,nz),
5     2     gv(5,nx,ny,nz),diss(5,nx,ny,nz)
6      do k=1,nz
7         do j=1,ny
8            do i=1,nx
9               do l=1,5
10                  t1= -0.5d0*dt*(
11     3            (g(l,i,j,kp1)-g(l,i,j,km1))/dz) +
12     4            dt/Re*((ev(l,i,j,k)-ev(l,im1,j,k))/dx +
13     6                  (gv(l,i,j,k)-gv(l,i,j,km1))/dz)
14                  rhs(l,i,j,k)=t1+t2
15               enddo
16            enddo
17         enddo
18      enddo
19      end
20