[CIG-LONG] Isostasy in GALE

[Walter Landry]

Robert Gray <graywacke at gmail.com> wrote:
> Hey Walter,
> Thanks for the help.
> I'm having trouble getting this input file going in Underworld.
>  Specifically, I don't know how to declare the equation to use for isostasy
> in the .xml file and not get a segmentation fault when I run the script.
>  I've attached the .xml file that I've set up for the simulation.  In this
> model, which includes temperature effects, I've set the coefficient of
> thermal expansion to zero for both materials (which have the same
> densities) to ensure an accurate isostatic compensation. Again, thanks for
> your help.

Underworld does not include an equation parser.  You have to use a
StandardConditionFunction instead.  I think that you can use
Velocity_SimpleShear to get the same effect.  Their StressBC may or
may not be correct.  I fixed a few things for the latest Gale release.

Cheers,
Walter Landry