<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Sorry about that, Birendra. I didn't actually look at the test before sending it. It seems to me this one may not be what you need. The examples/2d/subduction folder may be more what you need.<div><br></div><div>How curved is the fault in your problem?</div><div><br></div><div>Cheers,</div><div>Charles</div><div><br></div><div><br><div><div>On 6/08/2013, at 11:09 AM, Birendra jha wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Thank you Charles. Two comments/questions:<br>1. There is nothing about solver settings in the pylithapp.cfg, no petsc settings either.<br>2. Density is specified as 2500.0*m/2.<br><br>regards<br>Birendra<br><br>--------------------------------------------<br>On Tue, 8/6/13, Charles Williams <<a href="mailto:willic3@gmail.com">willic3@gmail.com</a>> wrote:<br><br> Subject: Re: [CIG-SHORT] Fault opening error<br> To: <a href="mailto:cig-short@geodynamics.org">cig-short@geodynamics.org</a><br> Date: Tuesday, August 6, 2013, 2:53 AM<br><br> Hi Birendra,<br> It's possible the tests don't get<br> distributed with the binaries. I'm attaching a<br> tarball with the directory Brad mentioned.<br> Cheers,Charles<br><br> On 6/08/2013, at 5:48 AM, Birendra jha<br> wrote:<br> Hi Brad<br><br> I did not find the tests/2d/nonplanar folder in the pylith<br> directory. Did you mean examples/2d/subduction?<br><br> Also, can you please send me the pylithapp.cfg with modified<br> 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>><br> 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<br> 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<br> and<br> didn't get an <br> error. I am out of the office for another week and<br> 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> Hi Brad<br><br> I did that. It did not<br> solve the fault opening<br> problem.<br><br> Do you have some time<br> to run this case? It should<br> 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<br> these<br> nodes is assigned some non zero value in the normal<br> direction during numerical iteration. Since fault is<br> under<br> compression everywhere, it gives an error from<br> FaultCohesiveDyn::integrateResidual. Is it possible<br> that<br> after the iterations have converged, there is no<br> slip<br> component in the normal direction? If so, then this<br> error<br> check in inetgrateResidual is probably too strict?<br><br> Do you have a simple<br> case where spontaneous fault<br> slip<br> is modeled on a curved fault surface?<br><br> Thanks and best<br> regards<br> Birendra<br><br><br> --------------------------------------------<br> On Sun, 8/4/13, Brad<br> Aagaard <<a href="mailto:baagaard@usgs.gov">baagaard@usgs.gov</a>><br> wrote:<br><br> Subject: Re: [CIG-SHORT] Fault<br> opening<br> error<br> To: <a href="mailto:cig-short@geodynamics.org">cig-short@geodynamics.org</a><br> Date:<br> Sunday, August 4, 2013, 8:23 PM<br><br> Birendra,<br><br> Your<br> solver tolerances need<br> adjustment. The FaultCohesiveDyn<br> object in a<br> quasi-static simulation<br> requires use<br> of absolute tolerances<br> in order to<br> properly detect stick/slip<br> with the<br> iterative solver (see<br> the friction<br> session in the CDM2013<br> tutorial for<br> details). You should<br> lower the<br> ksp_rtol and snes_rtol<br> tolerances to<br> 1.0e-20 to force use of<br> the<br> absolute tolerances.<br><br> Brad<br><br> On<br> 8/2/13 2:35 PM, Birendra jha<br> wrote:<br> > Hi everyone,<br> ><br> ><br> I am trying to run a model with a<br> fault<br> (FaultCohesiveDyn)<br> in it. The fault<br> surface is not planar<br> but is curved. I<br> am getting fault<br> opening error at some<br> node, which stops<br> the simulation. Why<br> does this happen? What<br> should I do to<br> avoid getting this<br> error?<br> ><br> ><br> Please see the attached files<br> which should run.<br> ><br> ><br> Here is what I have tried so<br> far:<br> > 1. different<br> fault dips<br> ><br> 2. different fault traces (plane<br> to curved)<br> > 3. different<br> fault friction<br> coefficients<br> ><br> ><br> My objective is to get some<br> reverse slip on this curved<br> fault surface. The<br> slip does not have<br> to be on the curved<br> part, it can be on<br> the part that<br> strikes parallel to the<br> y_pos boundary.<br> ><br> ><br> Thank you very much for any<br> help.<br> ><br> ><br> Birendra<br> ><br> ><br> ><br> ><br><br> _______________________________________________<br> > CIG-SHORT<br> mailing list<br> ><br> <a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br> ><br> <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br> ><br><br> _______________________________________________<br> CIG-SHORT mailing list<br> <a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br> <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br><br><br><br><br><br> _______________________________________________<br> CIG-SHORT mailing list<br> <a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br> <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><br><br><br> Charles A. WilliamsScientist<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): <a href="mailto:0064-4570-4600C.Williams@gns.cri.nz">0064-4570-4600C.Williams@gns.cri.nz</a><br><br><br><br> -----Inline Attachment Follows-----<br><br> _______________________________________________<br> CIG-SHORT mailing list<br> <a href="mailto:CIG-SHORT@geodynamics.org">CIG-SHORT@geodynamics.org</a><br> <a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short</a><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="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><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></span>
</div>
<br></div><div>
<div style="border-top: solid 1px black; border-bottom: solid 1px black;
padding: 10px 0; margin: 20px 0; font-size: 8pt;
font-family: Verdana, Arial, Helvetica, sans-serif;">Notice: This email and any attachments are confidential. If received in error please destroy and immediately notify us. Do not copy or disclose the contents.<br></div>
</div>
</body></html>