Thanks for your quick reply ...
On May 31, 2011, at 10:08 AM, Marc Saudreau wrote:
1. In GMSH you can group some mesh elements and name them (usually
with a code number).
Is it possible to get this code number within FiPy to be able to work
on these elements specifically.
For instance, a method like mesh.getFaceColor() ?
FiPy does not presently read that information.
Is there a possibility to get it in a future version of FiPy ?
Or is there a way for me to get this information by modifying some
*.py files ... with your help ?
This does not seem to difficult to achieve (but may be I'm wrong) since
this information is in the gmsh file,and you read this file. The
next step should be to store this information and to get it via a
getFaceColor() method.
If not, I guess that the unique solution is to used face coordinates
as in many examples you gave?
Correct.
This is not a good new since my geometry can not be obtain from
coordinates.
My question is about the variable outCoeffTop which is, in this
particular case, related to the y direction. But
how to define the same variable when the flux direction is imposed by
the normal direction of faces.
You can use mesh._getFaceNormals() to get that.
I thought to use the mesh._getFaceNormals() method to get that but how
to use it. Is this the best way to proceed ??
outCoeffTop = (maskmyfaces* (h*mesh._getFaceNormals)).getDivergence()
Thanks again for your help
Marc
