[aspect-devel] free surf in 3D

[Thieulot, C. (Cedric)]

Hi Ian,

I have been unfortunately disctracted (derailed ?) by administrative tasks recently.
I remember creating such a minimal case a while back. I will revive it and will make sure to have it ready by the hackathon.
Thanks for checking,

C.



On 12 May 2015, at 01:52, Ian Rose <ian.rose at berkeley.edu<mailto:ian.rose at berkeley.edu>> wrote:

Hi Cedric,

I realize this came up a while ago, but haven't been able to look into it in detail (sorry).  Any further progress on this?  I have not been able to reproduce the issue, so some kind of minimal example would be useful.  Also, are you certain about the line that is hanging?  I am aware of an issue that I need to address in project_normal_velocity_onto_boundary(), so if that is the offending line things may be simpler...

Best,
Ian


On Sat, Mar 28, 2015 at 4:33 AM, Timo Heister <timo.heister at gmail.com<mailto:timo.heister at gmail.com>> wrote:
Can you please provide a simple setup that we can run to reproduce it? Thanks.

On Fri, Mar 27, 2015 at 1:30 PM, Thieulot, C. (Cedric) <c.thieulot at uu.nl<mailto:c.thieulot at uu.nl>> wrote:
Hi all,

I have been trying to do some simple crustal 3D runs with my visco-plastic rheology and I am now using the free surface.
Since my rheology is nonlinear I use Iterated Stokes, so I had to ‘force’ Ian’s routines to accept this solver.
When I use no AMR, things go rather well, I get enormous mountains, very large (surface) deformation, very happy (see pic).
When I turn AMR on, the model does some tens of time steps and then stalls: the code does not crash and
my job still ‘runs’ on the cluster and uses 10’s of cores but their activity is zero and nothing happens for hours.
(I have checked the mesh, the fields, etc … everything looks as expected, no visually suspicious clue).
I have been putting flags all over the place to figure out where precisely it seems to stop and I have narrowed it down to
the highlighted line on the attached pic (taken from Simulator<dim>::FreeSurfaceHandler::make_constraints() in  source/simulator/freesurface.cc<http://freesurface.cc/> )
The code does not seem to pass this line and stalls there.
Would anybody have an idea of what could go wrong there ? and how to fix it ? :)

Cheers,

C.



<Screen Shot 2015-03-27 at 18.27.32.png>

<v_end.png>

_______________________________________________
Aspect-devel mailing list
Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel


_______________________________________________
Aspect-devel mailing list
Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

_______________________________________________
Aspect-devel mailing list
Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20150512/087fc7fc/attachment.html>