geodynamics / aspect

A parallel, extensible finite element code to simulate convection in both 2D and 3D models.
https://aspect.geodynamics.org/
Other
217 stars 232 forks source link

Improve continent cookbook #5888

Closed jdannberg closed 3 weeks ago

jdannberg commented 3 weeks ago

Because of #5882 I had a look through the cookbook to try and see what the problem was. These are unrelated improvements (unfortunately they do not fix the SIGFPE).

The only change in the setup is that the model had particles, it tracked the initial composition with these particles, but it did not interpolated them back onto the compositional field and instead solved the equation with the field method. I did not think this was intended, so I changed the particles to be active instead of passive, which also meant I had to change the generator to avoid empty cells. Also see #5882.

tjhei commented 3 weeks ago

This should address the first point of #5882

jdannberg commented 3 weeks ago

I have changed

set Number of particles per cell per direction = 4

and that avoids the SIGFPE.