Fabien -
I think the code below should get you going.
The changes I made were:
- `xVelocity` and `zVelocity` changed to rank-0 CellVariables. FiPy *always*
solves at cell centers.
- Created a rank-1 FaceVariable to hold the velocity. The cell components must
be interpolated and manually
Thanks a lot Jon.
I appreciate your help.
I will try to apply your remarks but I will certainly still need your help
again ;-) to develop the code.
Thanks again.
Fabien
> Le 15 août 2018 à 19:46, Guyer, Jonathan E. Dr. (Fed)
> a écrit :
>
> Fabien -
>
> I think the code below should get
