<div dir="ltr">Hi Roland,<div><br></div><div>Excellent, that makes sense and if I deactivate "rho" as an extra variable it works. I can create a ticket in a couple of days.</div><div><br></div><div>Thanks a lot.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El jue, 21 oct 2021 a las 14:04, Roland Haas (<<a href="mailto:rhaas@illinois.edu" target="_blank">rhaas@illinois.edu</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello Luciano,<br>
<br>
> Great points. I wasn't really sure how to use the backtrace, but now<br>
> following the addresses there I found that the problem was in the Outflow<br>
> thorn. In fact, paying closer attention to the parfile you shared a<br>
> while back, I see that you don't use the Outflow thorn. Is it possible to<br>
> use it within Llama maybe tweaking the "coord_system" parameter or the<br>
> interpolator? Or should I consider alternatives?<br>
I would have expected that Outflow should work with Llama since it only<br>
has fairly basics use of the interpolator.<br>
<br>
I can see that this may be triggered in your parfile since you set:<br>
<br>
Outflow::extra_variables = "<br>
HydroBase::rho<br>
HydroBase::Y_e<br>
HydroBase::entropy<br>
HydroBase::temperature<br>
GRHydro::poynting_scalar<br>
"<br>
<br>
and "rho" is already in the set of variables that Outflow interpolates<br>
all the time i.e. this matches the description ("same interpolation to<br>
be done into different output arrays.") of the trigger in the comment.<br>
<br>
So you could try what happens if you make sure that you are not<br>
requesting any of the variables that Outflow already interpolates then<br>
this may work. I think this is only "rho" right now. You can see the<br>
full list of default variables in Outflow's source code<br>
einsteinanalysis/Outflow/src/outflow.c:<br>
<br>
CCTK_STRING input_array_names[NUM_INPUT_ARRAYS]<br>
= { "ADMBase::gxx",<br>
"ADMBase::gxy",<br>
"ADMBase::gxz",<br>
"ADMBase::gyy",<br>
"ADMBase::gyz",<br>
"ADMBase::gzz",<br>
<br>
"HydroBase::vel[0]",<br>
"HydroBase::vel[1]",<br>
"HydroBase::vel[2]",<br>
"HydroBase::rho",<br>
<br>
"ADMBase::betax",<br>
"ADMBase::betay",<br>
"ADMBase::betaz",<br>
"ADMBase::alp",<br>
};<br>
<br>
If you could create a ticket with this parfile and linking to this<br>
email conversation (via the email archive<br>
<a href="http://lists.einsteintoolkit.org/pipermail/users/2021-October/date.html" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/pipermail/users/2021-October/date.html</a>)<br>
that would be great. This should be supported by Llama.<br>
<br>
Yours,<br>
Roland<br>
<br>
-- <br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://pgp.mit.edu" rel="noreferrer" target="_blank">http://pgp.mit.edu</a> .<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="trebuchet ms, sans-serif"><b>Luc Combi</b></font></div><div dir="ltr"><div style="font-family:sans-serif"><div style="color:rgb(117,117,117)">Instituto Argentino de Radioastronomía </div><div style="color:rgb(117,117,117)">CONICET, Buenos Aires, Argentina</div><div style="color:rgb(117,117,117)"><a href="http://sites.google.com/view/lucianocombi" target="_blank">web page.</a><br></div></div><div style="color:rgb(33,33,33);font-family:sans-serif"><br><br></div></div></div></div></div></div></div></div></div></div></div></div>