Dear List, i am trying to implement a source term based on the gradient of an auxiliary field in the APESolver. My problem is quite general, so i in terms of the simpler ADR solver, a minimal example would look like this: du / dt + v * grad(u) = c * grad(v) with u being a transported variable, v the advection velocity and c a constant. The problem is the gradient on the right hand side. I tried computing it by using ExpList::PhysDeriv(), but since the projection is not C1-continuous, this fails at the element edges. I attached an example plot to illustrate the problem. p0/rho0 uses the left axis, the gradient the right hand side one. Do you have any suggestions on how to handle the gradient term? Is there some other method i should use? Cheers, Kilian -- Kilian Lackhove, M.Sc. Fachgebiet für Energie- und Kraftwerkstechnik L1|08 114 Technische Universität Darmstadt Jovanka-Bontschits-Straße 2 D-64287 Darmstadt Germany Tel.: +49 6151 16 - 28915 Fax: +49 6151 16 - 6555 e-mail: lackhove@ekt.tu-darmstadt.de