[aspect-devel] spherical models with S40RTS
heister at clemson.edu
Fri Apr 10 13:08:26 PDT 2015
I can reproduce the problem here and it seems to be related to hanging
nodes in combination with null space removal (and only in 3d). I will
keep investigating. Thanks for the testcase.
On Wed, Apr 8, 2015 at 5:34 PM, Scott King <sdk at vt.edu> wrote:
> Hi All;
> Here is the file we have been using. We still get a solid body rotation that completely swamps everything.
> On Apr 8, 2015, at 4:25 PM, Timo Heister <heister at clemson.edu> wrote:
>> (cc'ed to the mailing list)
>>> I thought maybe you might know who has been doing the modeling with S40RTS.
>> This was contributed by Jacky:
>> https://github.com/geodynamics/aspect/commit/31b758 and there is no
>> example but at least a test case: tests/initial_condition_S20RTS.prm
>> I would be great to have an example in cookbooks, though!
>>> We can not seem to get the rotation out of the calculation even when
>>> "Remove nullspace net rotation" as on page 102. The result had smaller
>>> velocities than without this, but it still looks like a solid body rotation
>>> more than anything.
>> Are you willing to share the prm so we can take a look?
>>> Shangxin sent e-mail to the
>>> aspect-revel mailing list but it keeps bouncing. Errrr.
>> That would be the place to ask. What is the error message? Is the
>> email address you are sending from subscribed?
>>> Also curious, if we use the vtk output, how do we plot the depth_average.vtu
>>> file? I don't get anything useful with paraview but I've never tried to
>>> plot 2D files. The solution-xxxxx.vtu files work fine.
>> The way to do it with gnuplot is described in the manual, but I think
>> we could write a section for this. I can just open the
>> depth_average.vtk in paraview and can look at the depth averages,
>> though. What happens if you do that?
>> Timo Heister
More information about the Aspect-devel