[aspect-devel] spherical models with S40RTS

[Timo Heister]

Hey Scott,

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.
>
> Scott
>
>
>
>
> 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
>>> setting:
>>>
>>> "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
>> http://www.math.clemson.edu/~heister/
>
>



-- 
Timo Heister
http://www.math.clemson.edu/~heister/