[CIG-SHORT] Bug: Gravity with faults

[Brad Aagaard]

Chris-

Yes, this is a bug. This is related to a shortcoming of Pyrex/Pyrexembed 
that made it easy for this bug to occur. The bug prevents using gravity 
with faults. The Pyrex module for faults is missing the function to set 
the gravity field.

Since this bug disappears with the switch to SWIG and we will have a new 
release in the next month, I don't think it is worth our time to create 
a patch to fix this bug.

Brad


Charles Williams wrote:
> Hi Brad,
> 
> This does seem like a bug.  For our current fault implementation, I 
> think we should pass GravityField, but not do anything with it.  For 
> frictional faults, we may need it (I'm not sure).
> 
> Charles
> 
> 
> Begin forwarded message:
> 
>> *From: *Christopher DiCaprio <dicaprio at gps.caltech.edu 
>> <mailto:dicaprio at gps.caltech.edu>>
>> *Date: *14 May 2009 5:31:10 PM
>> *To: *cig-short at geodynamics.org <mailto:cig-short at geodynamics.org>
>> *Subject: **[CIG-SHORT] gravity problem*
>>
>> Thanks for the quick reply on the bug.
>>
>> I'm encountering a problem using gravity in a model containing a 
>> fault. It works fine for the same setup without a fault. I've attached 
>> the input files to this email.
>>
>> Thanks,
>> Chris
>>
>>
>> This is the error message:
>>
>> -- Pre-initializing problem.
>> >> 
>> /home/dicaprio/local/pylith-build/lib64/python2.5/site-packages/pylith/problems/Problem.py:141:verifyConfiguration
>> -- timedependent(info)
>> -- Verifying compatibility of problem configuration.
>> >> 
>> /home/dicaprio/local/pylith-build/lib64/python2.5/site-packages/pylith/problems/TimeDependent.py:113:initialize
>> -- timedependent(info)
>> -- Initializing problem.
>> Traceback (most recent call last):
>>  File "<string>", line 1, in <module>
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/mpi/__init__.py", 
>> line 42, in mpistart
>>    kwds = kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/pyre/applications/__init__.py", 
>> line 41, in start
>>    shell.run(**kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/pyre/applications/Shell.py", 
>> line 142, in run
>>    method(*args, **kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/pyre/applications/SuperScript.py", 
>> line 49, in execute
>>    self.main(*args, **kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/pyre/applications/SuperScript.py", 
>> line 58, in main
>>    self.runSubscript(*args, **kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/pyre/applications/SuperScript.py", 
>> line 62, in runSubscript
>>    self.subscript.run(*args, **kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/pyre/applications/Application.py", 
>> line 42, in run
>>    shell.run(*args, **kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/pyre/applications/Shell.py", 
>> line 142, in run
>>    method(*args, **kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/mpi/Application.py", 
>> line 118, in _onComputeNodes
>>    self.onComputeNodes(*args, **kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib/python2.5/site-packages/pythia-0.8.1.7-py2.5.egg/mpi/Application.py", 
>> line 124, in onComputeNodes
>>    self.main(*args, **kwds)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib64/python2.5/site-packages/pylith/PyLithApp.py", 
>> line 102, in main
>>    self.problem.initialize()
>>  File 
>> "/home/dicaprio/local/pylith-build/lib64/python2.5/site-packages/pylith/problems/TimeDependent.py", 
>> line 114, in initialize
>>    self.formulation.initialize(self.dimension)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib64/python2.5/site-packages/pylith/problems/Implicit.py", 
>> line 94, in initialize
>>    Formulation.initialize(self, dimension)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib64/python2.5/site-packages/pylith/problems/Formulation.py", 
>> line 186, in initialize
>>    integrator.initialize(totalTime, numTimeSteps)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib64/python2.5/site-packages/pylith/faults/FaultCohesiveKin.py", 
>> line 147, in initialize
>>    Integrator.initialize(self, totalTime, numTimeSteps)
>>  File 
>> "/home/dicaprio/local/pylith-build/lib64/python2.5/site-packages/pylith/feassemble/Integrator.py", 
>> line 98, in initialize
>>    self.cppHandle.gravityField = self.gravityField.cppHandle
>> AttributeError: 'faults.FaultCohesiveKin' object has no attribute 
>> 'gravityField'
>> --pyre-start: mpiexec: exit 1
>> /home/dicaprio/local/pylith-build/bin/pylith: 
>> /home/dicaprio/local/pylith-build/bin/nemesis: exit 1
>>
>>
>>
> 
> ------------------------------------------------------------------------
> 
>>
>>
> 
> ------------------------------------------------------------------------
> 
>>
>>
>>
>> _______________________________________________
>> 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>
> NOTE NEW E-MAIL ADDRESS
>