c ===================================================== SUBROUTINE cles_source(ux,vx,work,ncomps,ixhi,ixlo, $ rhs,nvars,nx,dx) c ===================================================== IMPLICIT none Integer ncomps, nvars INTEGER ixhi,ixlo INTEGER nx DOUBLE PRECISION dx DOUBLE PRECISION ux(ncomps,ixlo:ixhi) DOUBLE PRECISION vx(nvars,ixlo:ixhi) DOUBLE PRECISION work(ixlo:ixhi) DOUBLE PRECISION rhs(nvars,nx) INTEGER i c ---- default to no source term c do i=1, nx c rhs(2,i) = rhs(2,i)-ux(1,i) c rhs(5,i) = rhs(5,i)-ux(2,i) c enddo return END