[CIG-SHORT] No output

[Brad Aagaard]

Birendra,

I get vtk files in my output directory with the solution and state 
variables. The displacements are nonzero (~200m of vertical motion). You 
can check your solution and initial stress state using the solution. If 
the stresses match, then the displacements should be zero.

Regards,
Brad


On 03/13/2012 04:01 PM, Birendra jha wrote:
> Dear Brad, Charles,
>
> I clearly forgot to include the output directory in the zip file. But I have that on my computer so that's not the issue.
> I will try to reduce the length of the spatialdb files. Yes, it doesn't make sense to have 1 m resolution in spatialdb files when the grid cels are 100s of m. I may not be able to change to a tetrahedral grid but I'll try.
>
> Granted all that, I am sure it runs in few minutes. Are you able to see the vtk files in output folder? Are the displacements zero?
>
> I have run the same problem with constant density (just 1 line in density spatialdb and 2 lines in initial stress spatialdb) and get the desired result of zero displacement. Now I just want to see it with depth varying density. As you can probably imagine, I have density file for each of the 33 regions. I have their respective initial stress files too.
>
> regards
> Birendra
>
> --- On Wed, 3/14/12, Brad Aagaard<baagaard at usgs.gov>  wrote:
>
>> From: Brad Aagaard<baagaard at usgs.gov>
>> Subject: Re: No output
>> To: "Birendra jha"<bjha7333 at yahoo.com>
>> Cc: cig-short at geodynamics.org
>> Date: Wednesday, March 14, 2012, 4:17 AM
>> Birenda,
>>
>> I was able to run your test.cfg file after creating the
>> output directory (as Charles noted), but initializing the
>> integrators is taking a rather long time. Your spatial
>> databases that describe the variation in physical properties
>> are quite large. For every quadrature point in every cell
>> PyLith has to find which points to use for the
>> interpolation. Using a cloud of points with a much finer
>> spacing than the size of the cell will cause PyLith to run
>> much slower than necessary.
>>
>> The solution did converge, but when I looked at your mesh I
>> see some cells that have very poor aspect ratios so the
>> solution may not be very accurate. See the PyLith Wiki (http://www.geodynamics.org/cig/community/workinggroups/short/workarea/pylith-wiki/index_html)
>> for hints on checking the aspect ratio. Using tetehedral
>> cells will give you much greater flexibility in varying the
>> discretization size while keeping the aspect ratio of cells
>> close to 1. See examples/2d/subduction for an example of
>> using the bias scheme.
>>
>> Regards,
>> Brad
>>
>>
>>
>> On 03/13/2012 03:08 PM, Birendra jha wrote:
>>> Dear developers,
>>>
>>> I am trying to run a simple case like Step16 in the
>> tutorial except
>>> that density now varies with depth. I set the initial
>> geostatic
>>> stress from this density profile. I expect to see a
>> zero
>>> displacement. I have attached the files as zipped.
>>>
>>> Why I am not getting any output (It must be some small
>> mistake on my
>>> part. Just can't figure it out)?
>>>
>>> Thanks and regards Birendra
>>>
>>
>>
>