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:

   1. 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20210614/96210e68/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: BNS_IF_1.5_1.5_Ab2.2_lorenz_damp_fac_0.89_dx0125_MissingLink.par
Type: application/octet-stream
Size: 30523 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20210614/96210e68/attachment.obj

------------------------------

_______________________________________________
Users mailing list
Users@einsteintoolkit.org
http://lists.einsteintoolkit.org/mailman/listinfo/users


End of Users Digest, Vol 135, Issue 16
**************************************
_______________________________________________
Users mailing list
Users@einsteintoolkit.org
http://lists.einsteintoolkit.org/mailman/listinfo/users