[aspect-devel] VectorTools::integrate_difference fails to return an error that decreases under grid refinement (Elbridge Gerry Puckett) (Harsha Lokavarapu)

[Ying He]

On 2016-10-25 19:55, aspect-devel-request at geodynamics.org wrote:

> Today's Topics:
> 
>    1. Re: VectorTools::integrate_difference fails to return an
>       error that decreases under grid refinement (Elbridge Gerry
>       Puckett) (Harsha Lokavarapu)
> 
> 

> <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161025/0443df83/attachment.cc>
> -------------- next part --------------
> set Additional shared libraries            = ./libegp_hvl_box.so
> 
> ############### Global parameters
> 
> set Use years in output instead of seconds     = false
> 

Hi Gerry and Harsha,

There global parameter 'set Use years in output instead of seconds' was 
set up differently between the box and the annulus prm files. In your 
annulus model, the parameter wasn't chosen to assign a value, so by 
default, it uses years in stead of seconds. Therefore, for the final 
output solutions, there is a re-scaling back to the unit m/year. however 
your exact solution is still assumed with the unit m/s. That's why when 
you tried to compute the difference, the error is fixed no matter how 
the mesh refined.

Hope this help.

---
Best,
Ying