Andrew, >From this email, I am assuming things are copacetic.
Cheers On Jan 25, 2008 5:51 PM, asreeve <[EMAIL PROTECTED]> wrote: > > Daniel, > > The problem with the script I previously sent is in the boundary condition > term. fluxIn should be: > > fluxIn.setValue(101.*vel.getDivergence(),where=(Z==ZValue)) > > not > > fluxIn.setValue(101.*discharge.getDivergence(),where=(Z==ZValue)) > > I'm still using a negative DiffusionTerm in the tranport equation, which > doesn't make sense to me, but the output matches teh analytic solution. > > Thanks, > > Andy > > ---------------------- > Andrew Reeve > Dept. of Earth Science > University of Maine > Orono, ME 04469 > > > -- Daniel Wheeler
