[aspect-devel] Initial conditions for free-surface

[Ian Rose]

Hi Wolfgang,

It is possible to provide initial conditions, and without modifying the
FreeSurfaceHandler. What you would need to do is create a new GeometryModel
that has the initial shape that you require, and then flag the appropriate
boundary indicator as a free surface. Generally speaking, you can start
from a geometry model that has close to the right shape, then perturb it as
necessary.

For an example, you can look at what I have done with the ReboundBox
geometry model that is in
`benchmark/crameri_et_al/case_1/crameri_benchmark_1.cc` which sets up a
Cartesian box with sinusoidal topography for a relaxation benchmark.

Hope this helps!
Best,
Ian

On Tue, Jan 12, 2016 at 12:51 AM, Wolfgang Szwillus <
szwillus at geophysik.uni-kiel.de> wrote:

> Hi!
> Is it possible to provide initial conditions for the free surface? I would
> like to give an ascii file describing real topography. I guess what I have
> to do is to write a new class that works essentially like
> FreeSurfaceHandler, but has a differential initialization. Is it enough to
> make changes in FreeSurfaceHandler::execute() to achieve this effect?
>
> Cheers
> Wolfgang
> _______________________________________________
> 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/20160113/025f8bcd/attachment.html>