File: id-5.f

package info (click to toggle)
gcc-arm-none-eabi 15%3A12.2.rel1-1
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 959,712 kB
  • sloc: cpp: 3,275,382; ansic: 2,061,766; ada: 840,956; f90: 208,513; makefile: 76,132; asm: 73,433; xml: 50,448; exp: 34,146; sh: 32,436; objc: 15,637; fortran: 14,012; python: 11,991; pascal: 6,787; awk: 4,779; perl: 3,054; yacc: 338; ml: 285; lex: 201; haskell: 122
file content (19 lines) | stat: -rw-r--r-- 669 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
      subroutine shell(Re,Pr,nx,ny,nz,
     $nuim,nuex2,nuex4,cfl,scheme,conf,ni,maxit) 
      real*8 q(5,nx,ny,nz),dq(5,nx,ny,nz),rhs(5,nx,ny,nz),e(5,nx,ny,nz),
     1     f(5,nx,ny,nz),g(5,nx,ny,nz),ev(5,nx,ny,nz),fv(5,nx,ny,nz),
     2     gv(5,nx,ny,nz),diss(5,nx,ny,nz)
      do k=1,nz
         do j=1,ny
            do i=1,nx
               do l=1,5
                  t1= -0.5d0*dt*(
     3            (g(l,i,j,kp1)-g(l,i,j,km1))/dz) +
     4            dt/Re*((ev(l,i,j,k)-ev(l,im1,j,k))/dx +
     6                  (gv(l,i,j,k)-gv(l,i,j,km1))/dz)
                  rhs(l,i,j,k)=t1+t2
               enddo
            enddo
         enddo
      enddo
      end