[CIG-SHORT] An error when running Pylith .
faqidiao
faqidiao at asch.whigg.ac.cn
Wed Oct 26 18:39:35 PDT 2011
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
--
\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\
Dr. Carl W. Gable, Staff Scientist
Acting Group Leader, EES-16, Computational Earth Science
MS T003, Los Alamos National Laboratory, Los Alamos NM 87545
Email:gable at lanl.gov Voice:505-665-3533 Fax:505-665-8737
Schedule B
Homepage: http://www.ees.lanl.gov/staff/gable
LaGriT Mesh Generation Software: http://lagrit.lanl.gov
Mesh Generation, Geological Apps: http://meshing.lanl.gov
Focus:Geophysics, Hydrology, Mesh Generation, Computational Geometry
/\/\/\/\/\/\/\/\/\/\/Correspondence / TSPA\/\/\/\/\/\/\/\/\/\/\/\/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-short/attachments/20111027/fa565356/attachment-0001.htm
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wenc1.zip
Type: application/octet-stream
Size: 78337 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-short/attachments/20111027/fa565356/attachment-0001.obj
More information about the CIG-SHORT
mailing list