[CIG-SHORT] Modifying volcano_chamberdike example

Brad Aagaard baagaard at usgs.gov
Tue Sep 6 12:14:23 PDT 2016


On 09/05/2016 08:17 PM, Nuria Lumbreras wrote:
> Dear Sirs,
>
> My name is Nuria Lumbreras and I am working with PyLith for my master's
> thesis.
>
> With volcano_chamberdike example for Ruapehu volcano as a starting
> point, I've changed the topography and geometry trying to model the
> Tagoro volcano (which erupted in 2011 for first time) in El Hierro
> (Canary Islands). I am using a spherical magma chamber as well, and a
> 12km-length horizontal dike instead of a vertical one.
>
> However, I am finding some troubles when trying to set a time dependence
> for this last one. What I want to get is a temporally advancing dike,
> which is, no dike at time zero, and the whole dike (12km) after three
> months, evolving in a constant rate (6km at time 1.5 months). Would you
> please send me more information about this?

Please provide a more detailed explanation of what your difficulties 
are. Are you able to run the corresponding static problem?

> I also would appreciate if you could explain to me how to run the model
> with an spheroid. When I scale the sphere (in the z axis in my case), I
> can visualize the spheroid in CUBIT, but when I run the Pylith script I
> get the following message error:
>
> RuntimeError: Could not find group of points 'magma_chamber' in PETSc DM
> mesh.
> application called MPI_Abort(MPI_COMM_WORLD, -1) - process 0
> /home/nuria/pylith/bin/nemesis: mpirun: exit 255
> /home/nuria/pylith/bin/pylith: /home/nuria/pylith/bin/nemesis: exit 1

You are getting this because PyLith could not find a nodeset called 
"magma_chamber" in your mesh file. Make sure you create nodesets in your 
mesh for all groups of vertices where you want boundary conditions, etc.

Regards,
Brad






More information about the CIG-SHORT mailing list