<div dir="ltr">Luca<div><br></div><div>This check ensures that there is one active grid point in the y direction, and that the number of symmetry (i.e. boundary) points is the same as the number of ghost points.</div><div><br></div><div>This needs to hold independent of Carpet. If CarpetRegrid2 chooses a different domain size, then this is wrong, and you will need to correct the parameter file.</div><div><br></div><div>I have just added a test file to Cartoon2D. Look in particular at the CoordBase parameters there.</div><div><br></div><div>-erik<br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 28, 2015 at 3:26 AM, Luca Baiotti <span dir="ltr"><<a href="mailto:baiotti@ipc.phys.sci.osaka-u.ac.jp" target="_blank">baiotti@ipc.phys.sci.osaka-u.ac.jp</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
I have a question about CactusNumerical/Cartoon2D.<br>
<br>
In the source file ApplyCartoon.c, there is a condition for the grid setup:<br>
<br>
if(cctk_gsh[1] != 2*cctk_nghostzones[1]+1)<br>
{<br>
CCTK_WARN(0, "Grid size in y direction inappropriate.");<br>
}<br>
<br>
As far as I could tell, this is generally violated in multilevel<br>
Carpet (CarpetRegrid2) runs, with standard settings that determine the<br>
minimum extension of the refined grids, like<br>
Carpet::use_buffer_zones = yes<br>
mol::ode_method = rk4<br>
<br>
That condition for Cartoon2D can be achieved on the coarsest level by<br>
setting CoordBase::ymax equal to the spacing dy and CoordBase::ymin<br>
equal to -dy. But it cannot be set for refined levels, because<br>
CarpetRegrid2 adds as many points as required for the buffer zones.<br>
Indeed, Carpet adds enough buffer points for its own purposes that<br>
Cartoon2D seems to work fine.<br>
<br>
Should the above condition be changed to<br>
<br>
if(cctk_gsh[1] >= 2*cctk_nghostzones[1]+1)<br>
<br>
?<br>
<br>
Yours,<br>
Luca<br>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>><br><a href="http://www.perimeterinstitute.ca/personal/eschnetter/" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a></div>
</div>