Ok, thank you so much Federico. I think I will try to enlarge the grid up to a radius of 25 following the suggestion Frank just gave me and, in case I still have problems, I will try enlarging it a bit more.
Thank you again, Lorenzo
Il giorno lun 14 giu 2021 alle ore 12:12 Federico Cipolletta < cipo87@gmail.com> ha scritto:
That could be one easy workaround, but it specifically depends on your needs (i.e. the region in which you need more accuracy, the resolution you are interested in and most of all the resources that are available to you).
Best, Federico Cipolletta.
Il lun 14 giu 2021, 12:06 Lorenzo Ennoggi lorenzo.ennoggi@gmail.com ha scritto:
Hi Federico, great to hear from you! And thank you for your quick answer. So you are suggesting to enlarge the finest refinement level a bit more, right?
Thank you, Lorenzo
Il giorno lun 14 giu 2021 alle ore 11:56 Federico Cipolletta < cipo87@gmail.com> ha scritto:
Hello Lorenzo,
I recall that I obtained something similar to what you are describing, when attempting a regrinding for very small radii of the finest refinement level (let's say a most refinement radius near 5.0).
I wasn't able to control precisely this setting and I concluded that it could be related to the number of buffer/ghost zones that are required (of course it also depends on the most refined resolution, though).
Adopting greater radii, instead, it seems that I could have better control.
I am sorry of not being of any more help and I wonder if someone could add more.
Best regards, Federico Cipolletta.
Il lun 14 giu 2021, 11:45 users-request@einsteintoolkit.org ha scritto:
Send Users mailing list submissions to users@einsteintoolkit.org
To subscribe or unsubscribe via the World Wide Web, visit http://lists.einsteintoolkit.org/mailman/listinfo/users or, via email, send a message with subject or body 'help' to users-request@einsteintoolkit.org
You can reach the person managing the list at users-owner@einsteintoolkit.org
When replying, please edit your Subject line so it is more specific than "Re: Contents of Users digest..."
Today's Topics:
- Issue with thorn Trigger (Lorenzo Ennoggi)
Message: 1 Date: Mon, 14 Jun 2021 11:45:37 +0200 From: Lorenzo Ennoggi lorenzo.ennoggi@gmail.com Subject: [Users] Issue with thorn Trigger To: users@einsteintoolkit.org Message-ID: <CAG9ZJAyU0akW7iOrViD5OSHENw5_soDVNS= T1N+TgFtoV9u_aw@mail.gmail.com> Content-Type: text/plain; charset="utf-8"
Dear all, 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:
ActiveThorns = "Trigger" Trigger::Trigger_Number = 1
Trigger::Trigger_Checked_Variable[0] = "hydrobase::rho" Trigger::Trigger_Relation[0] = ">" Trigger::Trigger_Checked_Value[0] = 1.e-4 Trigger::Trigger_Reduction[0] = "maximum"
Trigger::Trigger_Reaction[0] = "steerparam" Trigger::Trigger_Steered_Parameter_Thorn[0] = "CarpetRegrid2" Trigger::Trigger_Steered_Parameter_Name[0] = "radius_1[5]" Trigger::Trigger_Steered_Parameter_Value[0] = "25." #New value of the above parameter Trigger::Trigger_Once[0] = 1 #It means "yes" Trigger::Trigger_Debug = 1 #It means "yes"
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
[...] INFO (Trigger): Testing triggers INFO (Trigger): last_checked: -1 INFO (Trigger): reducing trigger 0 red_handle 2 varindex 428 INFO (Trigger): reducing was ok INFO (Trigger): trigger nr. 0 fullfilled for rho (0.000804>0.000100) INFO (Trigger): Steering parameter INFO (Trigger): Parameter steered [...]
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?
I attach the parameter file in case you need it.
Thank you very much in advance for your help.
Sincerely, Lorenzo Ennoggi