[CIG-SHORT] An error when running Pylith .

Brad Aagaard baagaard at usgs.gov
Thu Oct 27 08:52:24 PDT 2011 

Faqi,

It looks like you are trying to share vertices/nodes across multiple 
fault surfaces. Nodes/vertices cannot be shared by multiple fault 
surfaces. We have an explicit check for this when using nodesets/psets 
to specify faults, but not when using a fault mesh from LaGriT. Using a 
fault mesh to specify faults is an old experimental feature that is 
fragile and untested.

All fault surfaces that are simply connected should be merged into a 
single fault surface.

Regards,
Brad


On 10/26/2011 06:39 PM, faqidiao wrote:
> Carl&  pylith developers,
>
>    Thank your all for your kindly help!
>
>    This problem perhaps induced by the fault mesh, and you may can help.
> Let me give some simple introduction about my mesh.
>
>    I attached my script,  from which I want to:
>
> 1) read in four faults and merge them into the background mesh.
> 2) extract the fault information from the combined mesh and make it as the
>      input file of Pylith.
>
>     From the output files of LaGrit, it seems that I have extracted the fault information
> correctly. I said it is correct because the extracted fault information is consistent
> with the original fault information (tmp_extractc_1.inp&  fault_info_slip_1 .table, etc).
>
>    I also checked the tmp_extractc_1.gmv, and didn't find any' hole' or 'silver' element
> on the fault mesh.
>
>    It should be mentioned that all faults in my mesh is curved, which means that the dip
> decreased with depth. And the strike of the faults are not uniform.
>
>   Please have a check of my mesh and give some suggestions to improve it.
>
>    Thank you for your help in advance!
>    Best wishes,
>
>    Faqi Diao
>
>
>
>
>
>
>
>
>
> From: Carl W. Gable
> Date: 2011-10-27 04:22
> To: Faqi Diao
> CC: Charles Williams
> Subject: Re: [CIG-SHORT] An error when running Pylith .
> Faqi,
> Let me know if you want any help analyzing the mesh to insure what you
> think you have is actually have. For example, I can tell you how to
> query the mesh to check for 3 faces (4 nodes) on fault, etc.
> Carl
>
>
> On 10/26/2011 01:49 PM, Charles Williams wrote:
>> It can also happen if there is too much fault curvature. For this to
>> happen, there must be 3 faces on the fault, which means pretty extreme
>> curvature or your cell size is not fine enough for the fault geometry
>> you are representing.
>>
>> Charles
>>
>>
>> On 27/10/2011, at 4:29 AM, Matthew Knepley wrote:
>>
>>> On Wed, Oct 26, 2011 at 10:13 AM, faqidiao<faqidiao at asch.whigg.ac.cn
>>> <mailto:faqidiao at asch.whigg.ac.cn>>  wrote:
>>>
>>>      __
>>>      Brad,
>>>      Thanks for your rapid response. I try my problem with new version
>>>      (pylith-1.6.1). However, error still occurred.
>>>      Please find the attached error information. I agree with you, this
>>>      perhaps induced by the fault topology. BUT,
>>>      I have checked my Lagrit running and it's result , all seems well.
>>>      So, I still puzzled about this.
>>>      I also attached my problem, and I don't know whether you have time
>>>      to run my model and check where the
>>>      problem is?
>>>
>>>
>>> Here is the problem. Somewhere in your mesh, there is a terahedron
>>> with all 4 vertices on a fault. This can happen
>>> because
>>>
>>> a) there is a problem in your fault description
>>>
>>> b) You have two faults which get too close to each other
>>>
>>> Thanks,
>>>
>>> Matt
>>>
>>>      Thank you again!
>>>      Faqi Diao
>>>      ------------------------------------------------------------------------
>>>      *From:* Brad Aagaard<mailto:baagaard at usgs.gov>
>>>      *Date:* 2011-10-26 12:45
>>>      *To:* faqidiao<mailto:faqidiao at asch.whigg.ac.cn>
>>>      *CC:* cig-short<mailto:cig-short at geodynamics.org>
>>>      *Subject:* Re: [CIG-SHORT] An error when running Pylith .
>>>      Faqi,
>>>      The error message suggests that the nodeset/pset you are using to
>>>      create the fault has a problem. This can happen if the fault has
>>>      holes in it or there is a problem with the topology of the mesh.
>>>      Another possibility is that there is a bug. We have made a lot of
>>>      improvements since v1.4.2, so I recommend that you use the current
>>>      release (v.1.6.1) and see if you still get the same error message.
>>>      Regards,
>>>      Brad
>>>      -----faqidiao<faqidiao at asch.whigg.ac.cn
>>>      <mailto:faqidiao at asch.whigg.ac.cn>>  wrote: -----
>>>      =======================
>>>      To: cig-short<cig-short at geodynamics.org
>>>      <mailto:cig-short at geodynamics.org>>
>>>      From: faqidiao<faqidiao at asch.whigg.ac.cn
>>>      <mailto:faqidiao at asch.whigg.ac.cn>>
>>>      Date: 10/25/2011 07:07PM
>>>      Subject: [CIG-SHORT] An error when running Pylith .
>>>      =======================
>>>      Hi all,
>>>      Recently, I met a problem when running Pylith1.4.2, and I don't
>>>      know how to fix
>>>      it, though I have paid much time.
>>>      Attached is the error message, can you have a look and give some
>>>      suggestions
>>>      on how to solve the problem?
>>>      Thank you in advance!
>>>      Best regards,
>>>      Faqi Diao_______________________________________________
>>>      CIG-SHORT mailing list
>>>      CIG-SHORT at geodynamics.org<mailto:CIG-SHORT at geodynamics.org>
>>>      http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>>>      [attachment(s) "error message.png" removed by Brad
>>>      Aagaard/GD/USGS/DOI]=
>>>
>>>      _______________________________________________
>>>      CIG-SHORT mailing list
>>>      CIG-SHORT at geodynamics.org<mailto:CIG-SHORT at geodynamics.org>
>>>      http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>>>
>>>
>>>
>>>
>>> --
>>> What most experimenters take for granted before they begin their
>>> experiments is infinitely more interesting than any results to which
>>> their experiments lead.
>>> -- Norbert Wiener
>>> _______________________________________________
>>> CIG-SHORT mailing list
>>> CIG-SHORT at geodynamics.org<mailto:CIG-SHORT at geodynamics.org>
>>> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>>
>> Charles A. Williams
>> Scientist
>> GNS Science
>> 1 Fairway Drive, Avalon
>> PO Box 30368
>> Lower Hutt 5040
>> New Zealand
>> ph (office): 0064-4570-4566
>> fax (office): 0064-4570-4600
>> C.Williams at gns.cri.nz<mailto:C.Williams at gns.cri.nz>
>>
>>
>>
>> _______________________________________________
>> CIG-SHORT mailing list
>> CIG-SHORT at geodynamics.org
>> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>
>

More information about the CIG-SHORT mailing list