<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Dear all<div class=""><br class=""></div><div class="">I am trying to write a simple first order upwind FVM to solve linear advection in 2d on triangular meshes using dmplex.</div><div class=""><br class=""></div><div class="">Based on some suggestions in recent discussions, I cooked up this simple example</div><div class=""><br class=""></div><div class=""><a href="https://github.com/cpraveen/cfdlab/tree/master/petsc/dmplex_convect2d" class="">https://github.com/cpraveen/cfdlab/tree/master/petsc/dmplex_convect2d</a></div><div class=""><br class=""></div><div class="">A gaussian profile must rotate around origin but in my code, there is some errors developing in just one or two cells, see picture.</div><div class=""><br class=""></div><div class="">I am looking for any suggestions/hints/tips on how to identify the problem. Since I have never written a dmplex code before, I am not able to ask a specific question.</div><div class=""><br class=""></div><div class="">Thanks</div><div class="">praveen</div><div class=""><br class=""></div><div class=""><img apple-inline="yes" id="92317BB5-11AB-4781-AE04-4670C0A6360D" src="cid:F64D5419-2817-418F-A3C4-605781AA1CE3" class=""></div></body></html>