<div dir="ltr">Miguel<div><br></div><div>Yes, this is possible. The trick is to set up a domain that is still 3D, but which is only 1 grid point wide in the z direction.</div><div><br></div><div>Clearly, this is not possible with cell-centred refinement, as refinement always splits cells into two. However, as vertex-centred refinement is the default, this might not be relevant here.</div><div><br></div><div>I also assume that you have a true 2+1 code that doesn't require any ghost zones in the z direction. You would then explicitly set the number of ghost zones and boundary zones to be 0 in the z direction.</div><div><br></div><div>The "usual" problem I encounter when setting up such a domain is that it's slightly inconvenient to specify the coordinates and grid spacings in the z direction. I prefer to define the domain extent via CoordBase, so I assume you're doing the same.</div><div><br></div><div>Usually, you would specify zmin, zmax, and ncells_z to define the domain. However, if you specify zmin=0, zmax=0, and then set ncells_z=0 (0 cells correspond to 1 point), then Carpet calculates dz = (zmax - zmin) / ncells_z, which is undefined since you divide zero by zero. So instead of specifying zmin, zmax, ncells_z, you have to specify zmin, zmax, and dz. Carpet then calculates ncells_k = (zmax - zmin) / dz, which is fine. In short, setting CoordBase::spacing = "numcells" will fail, and you have to use CoordBase::spacing = "gridspacing" instead.</div><div><br></div><div>That's about it. Of course, it's easy to make off-by-one errors in the parameter file, so looking at 1D ASCII output to see whether the coordinates and domain sizes you specified are correct is important.</div><div><br></div><div>Unfortunately I do not have a sample parameter file.</div><div><br></div><div>-erik</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 12, 2018 at 8:14 PM, Miguel Zilhão <span dir="ltr"><<a href="mailto:miguel.zilhao.nogueira@tecnico.ulisboa.pt" target="_blank">miguel.zilhao.nogueira@tecnico.ulisboa.pt</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">hi all,<br>
<br>
i have a 2+1 simulation that i would like to run with the ET, and i would like to use Carpet for the<br>
mesh refinement in the spatial (x and y) directions. is this possible? is there an example parameter<br>
file that i could have a look at?<br>
<br>
many thanks,<br>
Miguel<br>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.<wbr>org/mailman/listinfo/users</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>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><br></div></div></div>
</div>