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