vtf-logo

src/equations/cles_source2d.f


c     =====================================================
      SUBROUTINE cles_source(ux,vx,work,ncomps,ixhi,ixlo,
     $     iyhi,iylo,rhs,nvars,nx,ny,dx,dy)
c     =====================================================
      IMPLICIT none
  
      Integer ncomps, nvars
      INTEGER ixhi,ixlo,iyhi,iylo
      INTEGER nx,ny
      DOUBLE PRECISION dx,dy
      DOUBLE PRECISION ux(ncomps,ixlo:ixhi,iylo:iyhi)
      DOUBLE PRECISION vx(nvars,ixlo:ixhi,iylo:iyhi)
      DOUBLE PRECISION work(ixlo:ixhi,iylo:iyhi)
      DOUBLE PRECISION rhs(nvars,nx,ny)
      
      INTEGER i,j

c     ---- default to no source term
c     do j=1, ny
c        do i=1, nx
c           rhs(2,i,j) = rhs(2,i,j)-ux(1,i,j) 
c           rhs(5,i,j) = rhs(5,i,j)-ux(2,i,j) 
c        enddo
c     enddo
      
      return
      END

<