[aspect-devel] free surf in 3D

Timo Heister timo.heister at gmail.com
Sat Mar 28 04:33:38 PDT 2015

Can you please provide a simple setup that we can run to reproduce it?

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

>  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 )
> 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.
> _______________________________________________
> Aspect-devel mailing list
> 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/20150328/635830ac/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2015-03-27 at 18.27.32.png
Type: image/png
Size: 19361 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20150328/635830ac/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: v_end.png
Type: image/png
Size: 1384185 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20150328/635830ac/attachment-0003.png>

More information about the Aspect-devel mailing list