Main Page
src/flux2d.f

c $Header: /home/proj/vtfng/vtf/amroc/rim/src/flux2d.f,v 1.2 2004/12/01 20:51:35 ralf Exp $
      subroutine fluxI(nx,ny,nbc,neqn,dx,dy,dt,fx)
      implicit double precision (a-h,o-z)

      include 'rim2d.i'

      common/Bldir/dn1,dn2,tandn1,tandn2

      dimension fx(neqn,1-nbc:nx+nbc,1-nbc:ny+nbc)
      
      ifc=1
      do j=1,ny
         do i=1,nx+1        
            dn1 = 1.d0
            dn2 = 0.d0
            tandn1 =-dn2
            tandn2 = dn1
      
            call traceback(i,j,ifc,pf,denf,uf,vf,zf,gammaf,dx,dy,dt)
      
            Tf = pf/denf 
            stbx(i,j) =-rcon*zf*dexp(-eact/Tf)*denf
            
            fx(1,i,j) = dy*denf*uf
            fx(2,i,j) = dy*pf + fx(1,i,j)*uf
            fx(3,i,j) = fx(1,i,j)*vf
            ett = Tf/(gammaf-1.d0) + 0.5d0*(uf*uf+vf*vf)+q0*zf
            fx(4,i,j) = fx(1,i,j)*(ett+Tf)
            fx(5,i,j) = fx(1,i,j)*zf
            fx(6,i,j) = gammaf
         enddo
      enddo
      
      return
      end
c     
c------------------------------------------------------------------
c
      subroutine fluxJ(nx,ny,nbc,neqn,dx,dy,dt,fy)
      implicit double precision (a-h,o-z)

      include 'rim2d.i'

      common/Bldir/dn1,dn2,tandn1,tandn2
      
      dimension fy(neqn,1-nbc:nx+nbc,1-nbc:ny+nbc)
      
      ifc=0
      do j=1,ny+1
         do i=1,nx        
            dn1 = 0.d0
            dn2 = 1.d0
            tandn1 =-1.d0
            tandn2 = 0.d0
      
            call traceback(i,j,ifc,pf,denf,uf,vf,zf,gammaf,dx,dy,dt)
      
            Tf = pf/denf 
            stby(i,j) =-rcon*zf*dexp(-eact/Tf)*denf 
            
            fy(1,i,j) = dx*denf*vf
            fy(2,i,j) = fy(1,i,j)*uf
            fy(3,i,j) = dx*pf + fy(1,i,j)*vf 
            ett = Tf/(gammaf-1.d0)+0.5d0*(uf*uf+vf*vf)+q0*zf
            fy(4,i,j) = fy(1,i,j)*(ett+Tf)
            fy(5,i,j) = fy(1,i,j)*zf
            fy(6,i,j) = gammaf
         enddo
      enddo
      
      return
      end