gcc/testsuite/gfortran.dg/graphite/id-5.f

   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