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