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