<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div></div><div><br><div><div>On 6/08/2013, at 5:48 AM, Birendra jha wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Hi Brad<br><br>I did not find the tests/2d/nonplanar folder in the pylith directory. Did you mean examples/2d/subduction?<br><br>Also, can you please send me the pylithapp.cfg with modified tolerances?<br><br>Thanks a lot<br>Birendra<br><br>--------------------------------------------<br>On Mon, 8/5/13, Birendra jha <<a href="mailto:bjha7333@yahoo.com">bjha7333@yahoo.com</a>> wrote:<br><br> Subject: Re: [CIG-SHORT] Fault opening error<br> To: "Brad Aagaard" <<a href="mailto:baagaard@usgs.gov">baagaard@usgs.gov</a>><br> Date: Monday, August 5, 2013, 8:16 PM<br><br> Brad<br><br> Is it possible for you to send me the pylithapp.cfg, if<br> that's the one you modified?<br><br> I changed my ksp_rtol and snes_rtol to 1e-20 and I am still<br> getting the error.<br><br> Thank you<br> Birendra<br><br> --------------------------------------------<br> On Mon, 8/5/13, Brad Aagaard <<a href="mailto:baagaard@usgs.gov">baagaard@usgs.gov</a>><br> wrote:<br><br> Subject: Re: [CIG-SHORT] Fault opening error<br> To: "Birendra jha" <<a href="mailto:bjha7333@yahoo.com">bjha7333@yahoo.com</a>><br> Cc: <a href="mailto:cig-short@geodynamics.org">cig-short@geodynamics.org</a><br> Date: Monday, August 5, 2013, 7:44 PM<br><br> Birendra,<br><br> I ran a few time steps with the tighter tolerances and<br> didn't get an <br> error. I am out of the office for another week and don't<br> have time to <br> run anything until I return.<br><br> There is a curved fault dynamic simulation test in<br> tests/2d/nonplanar.<br><br> Regards,<br> Brad<br><br><br> On 8/5/13 7:11 AM, Birendra jha wrote:<br><blockquote type="cite">Hi Brad<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I did that. It did not solve the fault opening<br></blockquote> problem.<br><blockquote type="cite"><br></blockquote><blockquote type="cite">Do you have some time to run this case? It should<br></blockquote> take<br> 1-2 minutes. The main problem is that when slip is<br> non-zero<br> on the curved fault surface, slip vector on one of these<br> nodes is assigned some non zero value in the normal<br> direction during numerical iteration. Since fault is under<br> compression everywhere, it gives an error from<br> FaultCohesiveDyn::integrateResidual. Is it possible that<br> after the iterations have converged, there is no slip<br> component in the normal direction? If so, then this error<br> check in inetgrateResidual is probably too strict?<br><blockquote type="cite"><br></blockquote><blockquote type="cite">Do you have a simple case where spontaneous fault<br></blockquote> slip<br> is modeled on a curved fault surface?<br><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks and best regards<br></blockquote><blockquote type="cite">Birendra<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">--------------------------------------------<br></blockquote><blockquote type="cite">On Sun, 8/4/13, Brad Aagaard <<a href="mailto:baagaard@usgs.gov">baagaard@usgs.gov</a>><br></blockquote> wrote:<br><blockquote type="cite"><br></blockquote><blockquote type="cite"> Subject: Re: [CIG-SHORT] Fault opening<br></blockquote> error<br><blockquote type="cite"> To: <a href="mailto:cig-short@geodynamics.org">cig-short@geodynamics.org</a><br></blockquote><blockquote type="cite"> Date: Sunday, August 4, 2013, 8:23 PM<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> Birendra,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> Your solver tolerances need<br></blockquote> adjustment. The FaultCohesiveDyn<br><blockquote type="cite"> object in a<br></blockquote><blockquote type="cite"> quasi-static simulation requires use<br></blockquote> of absolute tolerances<br><blockquote type="cite"> in order to<br></blockquote><blockquote type="cite"> properly detect stick/slip with the<br></blockquote> iterative solver (see<br><blockquote type="cite"> the friction<br></blockquote><blockquote type="cite"> session in the CDM2013 tutorial for<br></blockquote> details). You should<br><blockquote type="cite"> lower the<br></blockquote><blockquote type="cite"> ksp_rtol and snes_rtol tolerances to<br></blockquote> 1.0e-20 to force use of<br><blockquote type="cite"> the<br></blockquote><blockquote type="cite"> absolute tolerances.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> Brad<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> On 8/2/13 2:35 PM, Birendra jha<br></blockquote> wrote:<br><blockquote type="cite"> > Hi everyone,<br></blockquote><blockquote type="cite"> ><br></blockquote><blockquote type="cite"> > I am trying to run a model with a<br></blockquote> fault<br><blockquote type="cite"> (FaultCohesiveDyn) in it. The fault<br></blockquote> surface is not planar<br><blockquote type="cite"> but is curved. I am getting fault<br></blockquote> opening error at some<br><blockquote type="cite"> node, which stops the simulation. Why<br></blockquote> does this happen? What<br><blockquote type="cite"> should I do to avoid getting this<br></blockquote> error?<br><blockquote type="cite"> ><br></blockquote><blockquote type="cite"> > Please see the attached files<br></blockquote> which should run.<br><blockquote type="cite"> ><br></blockquote><blockquote type="cite"> > Here is what I have tried so<br></blockquote> far:<br><blockquote type="cite"> > 1. different fault dips<br></blockquote><blockquote type="cite"> > 2. different fault traces (plane<br></blockquote> to curved)<br><blockquote type="cite"> > 3. different fault friction<br></blockquote> coefficients<br><blockquote type="cite"> ><br></blockquote><blockquote type="cite"> > My objective is to get some<br></blockquote> reverse slip on this curved<br><blockquote type="cite"> fault surface. The slip does not have<br></blockquote> to be on the curved<br><blockquote type="cite"> part, it can be on the part that<br></blockquote> strikes parallel to the<br><blockquote type="cite"> y_pos boundary.<br></blockquote><blockquote type="cite"> ><br></blockquote><blockquote type="cite"> > Thank you very much for any<br></blockquote> help.<br><blockquote type="cite"> ><br></blockquote><blockquote type="cite"> > Birendra<br></blockquote><blockquote type="cite"> ><br></blockquote><blockquote type="cite"> ><br></blockquote><blockquote type="cite"> ><br></blockquote><blockquote type="cite"> ><br></blockquote> _______________________________________________<br><blockquote type="cite"> > CIG-SHORT mailing list<br></blockquote><blockquote type="cite"> > <a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br></blockquote><blockquote type="cite"> > <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br></blockquote><blockquote type="cite"> ><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> _______________________________________________<br></blockquote><blockquote type="cite"> CIG-SHORT mailing list<br></blockquote><blockquote type="cite"> <a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br></blockquote><blockquote type="cite"> <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><br><br><br>_______________________________________________<br>CIG-SHORT mailing list<br><a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br>http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short<br></div></blockquote></div><br><div>
<span class="Apple-style-span" style="font-size: 12px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font face="Helvetica" size="3" style="font: normal normal normal 12px/normal Helvetica; ">Charles A. Williams</font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Scientist<br>GNS Science<br>1 Fairway Drive, Avalon<br>PO Box 30368<br>Lower Hutt 5040<br>New Zealand<br>ph (office): 0064-4570-4566<br>fax (office): 0064-4570-4600</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><a href="mailto:C.Williams@gns.cri.nz">C.Williams@gns.cri.nz</a></div></div></div></span>
</div>
<br></div></body></html>