Thank you again Walter,<br><br>For the thermal problem, like you said in the previous mail, it is somehow related to my configuration.<br><br>Cheers<br><br><br><div class="gmail_quote">On Tue, Mar 6, 2012 at 1:29 AM, Walter Landry <span dir="ltr"><<a href="mailto:walter@geodynamics.org">walter@geodynamics.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Cc'ing the list.<br>
<div class="im"><br>
nicolas riel <<a href="mailto:a.nicolas.riel@gmail.com">a.nicolas.riel@gmail.com</a>> wrote:<br>
> Hello Walter,<br>
><br>
> I'm still working with the code and encountered some problems when running<br>
> a model with simple rheologies and without thermal properties. This model<br>
> (attached file) uses an inflow boundary on the left side with a wrapTop<br>
> condition (with or without a stressBC simulating isostasy on the bottom).<br>
><br>
> The step 0 is going well but the step 1 crashes and gives this error :<br>
><br>
> "side is moving in the wrong direction".<br>
><br>
> Walter, do you think that an air layer with a staticTop instead of a<br>
> wrapTop would fix this problem ?<br>
><br>
> I would like to try it, however I don't know how to set the "maxAirY" with<br>
> the "json" syntax. I think that could be a good idea to add it in the<br>
> manual. This way one can fix the Top boundary with vy=0 and stabilizes the<br>
> model.<br>
<br>
</div>I would not recommend using an air layer in general. An air layer is<br>
only useful if the height gets very small. For example, if you are<br>
simulating a wedge. If you do want to use an air layer, you do not<br>
have to do anything special. "maxAirY" is no longer needed because<br>
HydrostaticTerm is no longer needed.<br>
<br>
As for your current input file, you just need to fix the top left corner<br>
<br>
"staticLeftTop": "True"<br>
<div class="im"><br>
> Moreover when I try to add the thermal properties for example in the<br>
> Geomod2004 benshmark, the thermal field at step 0 is always rough with<br>
> irregular boundaries.<br>
<br>
</div>What do you mean by "rough"? The bottom is set to 873 while the bulk<br>
is set to 573. That is what I see. What do you see?<br>
<br>
Cheers,<br>
Walter Landry<br>
</blockquote></div><br><br clear="all"><br>-- <br>Nicolas RIEL<br>Doctorant / Phd student<br><br>Institut des Sciences de la Terre<br>(ISTerre-UMR CNRS/UJF 5275) Maison des Géosciences<br>1381, Rue de la Piscine<br>38400 - Saint Martin d'Hères<br>
<br><a href="http://nicolas.riel.free.fr" target="_blank">http://nicolas.riel.free.fr</a><br>
Tél. : 33.(0)4.76.51.40.57<br>Fax. : 33.(0)4.76.51.40.58<br>Mob. : 33.(0)6.65.58.97.99<br>