<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Walter,<div><br></div><div>This seems to have resolved the TimeIntegrand error issue – thanks!</div><div><br></div><div>However, I have another question that relates to one of the modifications you and George made to the model. The way the model ends are set up now, there is effectively no inflow of material, because the two ends of the model move at a rate equal to the velocity of the inflowing material. As a result, the models end up shortening more and more the longer they run, and with some of the smaller/shorter models, we practically run out of material in the late stages of the run. (Obviously you can't try to simulate 70 km of displacement on a 60 km model if the model is going to shorten by the full amount of the displacement.)</div><div><br></div><div>George and I experimented a bit with the Front/BackEquations, but we were unable to modify them to allow the model to actually retain its dimensions or change them more slowly; we always got the "wall is moving in the wrong direction" error. Is there any workaround that will allow us to keep the models at a constant size, but have material flow in and out of the two ends? Or should we try something else, like just defining the models to be 70+ km longer than we need them to be at the start, and allow that material to be gradually lost as the simulation progresses?</div><div><br></div><div>Thanks,</div><div>Maria</div><div><br></div><div><br></div><div><br></div><div><div><div>On May 9, 2012, at 9:29 AM, María Helga Guðmundsdóttir wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Great, thanks for catching that. I'll take those lines out and run a test model to see what happens.<div><br></div><div>-M</div><div><br><div apple-content-edited="true">
<span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; -webkit-text-decorations-in-effect: none; text-indent: 0px; -webkit-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>María Helga Guðmundsdóttir</div><div>MS Student, Stanford University</div><div>Geological and Environmental Sciences</div></div></span></div></span><span class="Apple-style-span" style="font-size: 12px; "><a href="mailto:mariahg@stanford.edu">mariahg@stanford.edu</a></span><br class="Apple-interchange-newline"></span>
</div>
<br><div><div>On May 8, 2012, at 6:10 PM, Walter Landry wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Walter Landry <<a href="mailto:walter@geodynamics.org">walter@geodynamics.org</a>> wrote:<br><blockquote type="cite">María Helga Guðmundsdóttir <<a href="mailto:mariahg@stanford.edu">mariahg@stanford.edu</a>> wrote:<br></blockquote><blockquote type="cite"><blockquote type="cite">Hi Walter,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">I've been running the restraining bend models that George Hilley has<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">been troubleshooting with you lately, with mixed results. Two of<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">them have run fine upto the amount of displacement we want, but<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">another four models have crashed in an identical fashion at<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">completely different points during the run (timesteps 8, 11, 90, and<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">289 respectively, with total displacements of 4, 10, 17, and<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">37km).<br></blockquote></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The problem is the mesh distortion. I am attaching a picture taken<br></blockquote><blockquote type="cite">from the inside, and you can see that the height of the mesh almost<br></blockquote><blockquote type="cite">goes to zero. I am not sure what to do about this. What do you change<br></blockquote><blockquote type="cite">between the runs that work and the runs that do not?<br></blockquote><br>One thing I noticed about your input file is that in EulerDeform, you<br>have the lines<br><br> "staticTopFront": "True",<br> "staticTopBack": "True",<br><br>Those should be False or just not there. That will induce grid<br>distortion at outflow, because while the bulk is thinning, the<br>boundary is kept high.<br><br>Cheers,<br>Walter Landry<br></div></blockquote></div><br></div></div>_______________________________________________<br>CIG-LONG mailing list<br><a href="mailto:CIG-LONG@geodynamics.org">CIG-LONG@geodynamics.org</a><br>http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long<br></blockquote></div><br></div></body></html>