Morning,
Wow, answers straight from the gurus. That you very much. Now I know I
can dig deeper in the literature and in the code and be confident that
it will lead somewhere. I really appreciate it. It will probably take me
a while to process everything though!
Best regards
Jonathan Guyer wrote:
On Dec 17, 2007, at 12:43 PM, I wrote:
You can also get the gradient at the face between cells, which
happens to sit exactly at r = R, with:
>>> y.getFaceGrad()[rF == 4][0]
array([ 8.])
Oops. Make that
>>> rF = mesh.getFaceCenters()[...,0]
>>> y.getFaceGrad()[rF == R][0]
array([ 8.])