<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body text="#000000" bgcolor="#FFFFFF">
    <div class="moz-cite-prefix">On 10/04/2012 10:41 AM, Matthew Knepley
      wrote:<br>
    </div>
    <blockquote
cite="mid:CAMYG4GkE2zUS2WdTN=pHciBE=Fe6gM=QuhDwovTDVbqVtxh9gg@mail.gmail.com"
      type="cite">On Thu, Oct 4, 2012 at 10:39 AM, Hongfeng Yang <span
        dir="ltr">&lt;<a moz-do-not-send="true"
          href="mailto:hyang@whoi.edu" target="_blank">hyang@whoi.edu</a>&gt;</span>
      wrote:<br>
      <div class="gmail_quote">
        <blockquote class="gmail_quote" style="margin:0 0 0
          .8ex;border-left:1px #ccc solid;padding-left:1ex">
          Hi Developers,<br>
          <br>
          I am tracking time history of variables on particular vertices
          on the<br>
          fault. However, I found the vertex ID in Pylith (v_fault) is
          larger than<br>
          the total number of vertices on the fault. How can I convert
          the value<br>
          of v_fault back to the vertex location on the fault surface?<br>
          <br>
          For example, I got the following fault opening error message<br>
          WARNING! Fault opening with nonzero traction., v_fault: 4154,
          opening:<br>
          1.03716e-10, normal traction: -0.0042952<br>
          mpinemesis: faults/FaultCohesiveDyn.cc:382: virtual void<br>
          pylith::faults::FaultCohesiveDyn::integrateResidual(const<br>
          pylith::topology::Field&lt;pylith::topology::Mesh&gt;&amp;,
          PylithScalar,<br>
          pylith::topology::SolutionFields*): Assertion
          `fabs(tractionNormal) &lt;<br>
          _zeroTolerance' failed.<br>
          [0]0:Return code = 0, signaled with Aborted<br>
          <br>
          v_fault points to vertex ID 4154. The total number of vertices
          on the<br>
          fault is 1472, however.<br>
        </blockquote>
        <div><br>
        </div>
        <div>The fault vertices are just vertices in the mesh, so this
          refers to the full</div>
        <div>mesh vertex numbering.</div>
        <div><br>
        </div>
      </div>
    </blockquote>
    Does Pylith output the vertex and topology in a re-numbering way? I
    read the output value of each cell on the fault and could not find
    the value of v_fault in any cell.<br>
    <br>
    Thanks,<br>
    <br>
    Hongfeng <br>
    <br>
  </body>
</html>