c =====================================================
SUBROUTINE cles_source(ux,vx,work,ncomps,ixhi,ixlo,
$ iyhi,iylo,izhi,izlo,rhs,nvars,nx,ny,nz,dx,dy,dz)
c =====================================================
IMPLICIT none
Integer ncomps, nvars
INTEGER ixhi,ixlo,iyhi,iylo,izhi,izlo
INTEGER nx,ny,nz
DOUBLE PRECISION dx,dy,dz
DOUBLE PRECISION ux(ncomps,ixlo:ixhi,iylo:iyhi,izlo:izhi)
DOUBLE PRECISION vx(nvars,ixlo:ixhi,iylo:iyhi,izlo:izhi)
DOUBLE PRECISION work(ixlo:ixhi,iylo:iyhi,izlo:izhi)
DOUBLE PRECISION rhs(nvars,nx,ny,nz)
INTEGER i,j,k
c ---- default to no source term
c do k=1, nz
c do j=1, ny
c do i=1, nx
c rhs(2,i,j,k) = rhs(2,i,j,k)-ux(1,i,j,k)
c rhs(5,i,j,k) = rhs(5,i,j,k)-ux(2,i,j,k)
c enddo
c enddo
c enddo
return
END