[CIG-LONG] [Gale] EulerDeform: coordinate out of bounds when increasing model resolution

[Walter Landry]

<Guillaume.Duclaux at csiro.au> wrote:
> G'day all,
> 
> I'm running compression models with similar geometries and
> mechanical properties at different resolution and a fixed timestep
> (dt=100 ka).
> An error with EulerDeform seems to occur when increasing the
> resolution. The higher the resolution, the earlier the error
> happens.

This is expected.  When a new mesh is created, Gale needs to
interpolate values onto the new mesh.  In your case, the mesh has
moved so far that it needs data on other processors which it does not
have.

So your time step is too big, and you need to reduce it.  With the
current timestep, your results are inaccurate anyway.  Note that the
timestep needs to get smaller for higher resolution grids.  The
simplest way is to let Gale choose the timestep for you.

Cheers,
Walter Landry
walter at geodynamics.org