Patchiness in initial composition from GWB in spherical shell

ASPECT
1

I am trying to use GWB to create initial conditions in a spherical shell geometry. When visualizing in GWB it is ok, but some patchy behavior is observed at the surface for the compositions of ASPECT output.

This is from visualization

Screenshot from 2023-01-03 10-03-11
Screenshot from 2023-01-03 10-03-111920×1080 319 KB

This is from the output

Screenshot from 2023-01-03 10-08-05
Screenshot from 2023-01-03 10-08-051920×1080 158 KB

I ran this model in debug mode. For mesh generation this was the script

subsection Mesh refinement
  set Initial global refinement                = 3
  set Initial adaptive refinement              = 3
  set Time steps between mesh refinement       = 5
  set Strategy                                 = temperature, composition, topography
  set Refinement criteria merge operation      = plus
end

Increasing or decreasing the refinement level only affects the resolution of the patches but they remain regardless. What can be done for this?

Thanks,
Mainak

2

Hello Mainak,

Do you set the parameter “min depth” for this region? In spherical geometry it might be necessary to set the min depth to a negative value so the region that is “inside” the region does not just miss the nodes at the surface. Try setting “min depth”:-10e3, like in this example:

“composition models”:[
{“model”:“uniform”, “compositions”:[0], "min depth”:-10e3, “max depth”:10e3}]},

Magali

3

Hi Magali,

Many thanks for the suggestion, it has worked. When declaring the “min depth” it had to be done both in general description as well as in “composition models” then it worked. Like this

   {"model":"continental plate", "name":"continental plate", "min depth":-15e3,
     "composition models":[{"model":"uniform","compositions":[0], "min depth":-15e3}]}

Thank you!

Mainak

4

So happy to hear that its fixed!
Cheers,
Magali