<div dir="ltr">Dear Roland, Gabriele,<div>thank you, that's a good point. Indeed, I let the system evolve for well more than 32 iterations after the refinement level to be enlarged pops out (at the time of checkpoint, before merger, that level is not yet present; it appears at the time of merger). But the level was not enlarged by simply steering the parameter with <span style="font-family:monospace">Trigger::Trigger_Reaction[0]="steerparam"</span><font face="arial, sans-serif">.</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">One thing that may be important (but I am not sure) is, as I just mentioned, that the refinement level I want to enlarge is not yet present at the time of checkpoint, it appears after ~2000 iterations.</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Sincerely,</font></div><div><font face="arial, sans-serif">Lorenzo</font></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno lun 14 giu 2021 alle ore 18:02 Gabriele Bozzola <<a href="mailto:bozzola.gabriele@gmail.com">bozzola.gabriele@gmail.com</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Lorenzo,<div><br></div><div>In my experience, it is not enough to change the value of the parameter:</div><div>you also need to regrid. If you don't regrid, the grid configuration will not</div><div>change. (At least, this is what I found using Trigger to change the number</div><div>of refinement levels.) In your par file you regrid every 32 iterations, have </div><div>you tried checking the output after regridding?</div><div><br></div><div>Gabriele</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jun 14, 2021 at 8:37 AM Lorenzo Ennoggi <<a href="mailto:lorenzo.ennoggi@gmail.com" target="_blank">lorenzo.ennoggi@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear Erik,<div>thank you very much. I am now trying the way Frank suggested before and, if again the finest refinement level does not get enlarged, I will probably try to edit the checkpoint files.</div><div><br></div><div>Sincerely,</div><div>Lorenzo</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno lun 14 giu 2021 alle ore 15:45 Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">On Mon, Jun 14, 2021 at 5:45 AM Lorenzo Ennoggi<br>
<<a href="mailto:lorenzo.ennoggi@gmail.com" target="_blank">lorenzo.ennoggi@gmail.com</a>> wrote:<br>
><br>
> Dear all,<br>
> I would like to restart a simulation from a checkpoint enlarging the finest refinement level in my grid. In particular, I would like to change CarpetRegrid2::radius_1[5] = 20.1 to CarpetRegrid2::radius_1[5] = 25. I tried to do that by making the parameter CarpetRegrid2::radius_1 always steerable and adding the following lines to my parameter file:<br>
><br>
><br>
> ActiveThorns = "Trigger"<br>
> Trigger::Trigger_Number = 1<br>
><br>
> Trigger::Trigger_Checked_Variable[0] = "hydrobase::rho"<br>
> Trigger::Trigger_Relation[0] = ">"<br>
> Trigger::Trigger_Checked_Value[0] = 1.e-4<br>
> Trigger::Trigger_Reduction[0] = "maximum"<br>
><br>
> Trigger::Trigger_Reaction[0] = "steerparam"<br>
> Trigger::Trigger_Steered_Parameter_Thorn[0] = "CarpetRegrid2"<br>
> Trigger::Trigger_Steered_Parameter_Name[0] = "radius_1[5]"<br>
> Trigger::Trigger_Steered_Parameter_Value[0] = "25." #New value of the above parameter<br>
> Trigger::Trigger_Once[0] = 1 #It means "yes"<br>
> Trigger::Trigger_Debug = 1 #It means "yes"<br>
><br>
><br>
> In this way, the trigger should act immediately, as the maximum of rho is ~1e-3 in my case. Indeed, it looks like the trigger is working fine, as the standard output says<br>
><br>
> [...]<br>
> INFO (Trigger): Testing triggers<br>
> INFO (Trigger): last_checked: -1<br>
> INFO (Trigger): reducing trigger 0 red_handle 2 varindex 428<br>
> INFO (Trigger): reducing was ok<br>
> INFO (Trigger): trigger nr. 0 fullfilled for rho (0.000804>0.000100)<br>
> INFO (Trigger): Steering parameter<br>
> INFO (Trigger): Parameter steered<br>
> [...]<br>
><br>
> but unfortunately, looking at the 2D output with VisIt, the finest refinement level does not get enlarged (i.e. radius_1[5] does not appear to have been steered at all). Do you have any suggestions on how to get this trigger to work? Or, are you aware of other ways of obtaining the same result without using the thorn Trigger?<br>
<br>
Changing these parameters has no effect, even if they were steerable.<br>
<br>
CarpetRegrid2 uses these parameters to set up the initial grid<br>
structure, by setting the grid scalars "positions", "radius", etc.<br>
During regridding, only these grid scalars are examined; the<br>
parameters are ignored. Other mechanisms to update the grid structure<br>
(e.g. by tracking horizons or punctures) update these grid scalars,<br>
which then determines the grid structure during regridding. These grid<br>
scalars are checkpointed and recovered, as all grid variables are.<br>
<br>
To modify the grid structure you need to modify these grid scalars.<br>
You can do this e.g. via Python by doing surgery on the checkpoint<br>
files, or by scheduling a function that runs right after recovery.<br>
Carpet performs a regridding step after recovering, before the first<br>
evolution step.<br>
<br>
-erik<br>
<br>
-- <br>
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/" rel="noreferrer" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a><br>
</blockquote></div>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote></div>
</blockquote></div>