geodynamics/aspect — pull request #140 by gassmoeller (2014-06-20 19:38:36+00:00)
Title: Minimum refinement function for cell vertices.

When prescribing a minimum refinement level function, one often wants to refine boundary layers. If the current refinement however is to coarse to resolve this depth (e.g. a cell size of 300 km at a boundary layer width of 100 km) then even the midpoint of the cell at which the minimum refinement function is evaluated is not resolving this boundary layer. Therefore, I would argue to change the logic to evaluate the function at every cell vertex and as soon as one vertex requires action (stop coarsening or set refine flag) do this. It is still possible to define function that are not resolved at all, but at least every function of the kind: if depth < x is caught, because the outermost vertices have a depth of 0. I will add a line to changes.h in the next pull request, since I do not want to create conflicts.