[aspect-devel] The input for limiting the mesh refinement area

[Nan Zhang]

Hi all,

I added the "maximum refinement function" into shell_simple_3d.prm like:
**************************************************************
...
subsection Mesh refinement
  set Initial global refinement          = 2
  set Initial adaptive refinement        = 4
  set Time steps between mesh refinement = 15
  set Strategy                           = temperature, maximum refinement
function

  subsection Maximum refinement function
    set Coordinate system = spherical
    set Function constants =
    set Variable names = r,phi,theta
    set Function expression = theta<1.0472      #60 degree
  end
end
...
**************************************************************
I tried to limit mesh refinement into a region of colatitude smaller than
60 degree. The mesh I got is attached. It seems that the north 60o region
has been refined, or other region is coarsen? What I don't understand is
the (1) how to connect the initial global refinement 2 with this figure?
(2) as well as the initial adaptive refinement = 4? (3) how to see the
"maximum refinement level" in the north 60o region?

Bests,
Nan

On Thu, Jan 28, 2016 at 11:08 AM, Juliane Dannberg <dannberg at gfz-potsdam.de>
wrote:

>
> Hi Nan,
>
> I think there are plugins called "Minimum refinement function" and
> "Maximum refinement function" that allow you to set a minimum and
> maximum level of refinement as function of, for example, the radius.
> Maybe this is what you are looking for? The syntax is:
>
> subsection Mesh refinement
>   set Strategy                                 = temperature, maximum
> refinement function
>
>   subsection Maximum refinement function
>     set Coordinate system      = spherical
>     set Function constants     =
>     set Function expression    = 6 +
> 2*sin((r+250000)*2*pi/250000)*sin(12*phi)
>     set Variable names         = r,phi
>   end
> end
>
> Best,
> Juliane
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20160129/5ee5f230/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MeshArea.PNG
Type: image/png
Size: 221039 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20160129/5ee5f230/attachment-0001.png>