[CIG-SHORT] No output

[Birendra jha]

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
> > 
> 
>