[CIG-LONG] Crash =?iso-2022-jp-3?B?GyQoUSN8GyhC?= Error in TimeIntegrand

[Walter Landry]

Walter Landry <walter at geodynamics.org> wrote:
> María Helga Guðmundsdóttir <mariahg at stanford.edu> wrote:
>> Hi Walter,
>> 
>> I've been running the restraining bend models that George Hilley has
>> been troubleshooting with you lately, with mixed results. Two of
>> them have run fine upto the amount of displacement we want, but
>> another four models have crashed in an identical fashion at
>> completely different points during the run (timesteps 8, 11, 90, and
>> 289 respectively, with total displacements of 4, 10, 17, and
>> 37km).
> 
> The problem is the mesh distortion.  I am attaching a picture taken
> from the inside, and you can see that the height of the mesh almost
> goes to zero.  I am not sure what to do about this. What do you change
> between the runs that work and the runs that do not?

One thing I noticed about your input file is that in EulerDeform, you
have the lines

  "staticTopFront": "True",
  "staticTopBack": "True",

Those should be False or just not there.  That will induce grid
distortion at outflow, because while the bulk is thinning, the
boundary is kept high.

Cheers,
Walter Landry