<div dir="ltr">Hello,<br>
I am running a simulation with Daniel Siegel where we use a base level<br>
plus 6 refinement levels, like the following:<br>
CarpetRegrid2::radius_1[1] = 336.0 # dx/64<br>
CarpetRegrid2::radius_1[2] = 168.0 # dx/64<br>
CarpetRegrid2::radius_1[3] = 84.0 # dx/64<br>
CarpetRegrid2::radius_1[4] = 42.0 # dx/64<br>
CarpetRegrid2::radius_1[5] = 21.0 # dx/64<br>
CarpetRegrid2::radius_1[6] = 10.5 # dx/64<br>
<br>
We ran simulations of a neutrino-driven wind from a Neutron star, out<br>
to ~150 ms, and now we want to turn on a high magnetic field. We have<br>
tried B-fields of varying strengths, but if we try a B-field strength<br>
that is too high, we find errors generated on the innermost refinement<br>
level boundary. We believe if we enlarge the innermost box size ie.<br>
CarpetRegrid2::radius_1[6] = 10.5 --> ~19-20, these errors would<br>
disappear. But we are having trouble doing this. In our attempt, we<br>
activate the trigger thorn with the following settings:<br>
<br>
ActiveThorns = "Trigger"<br>
Trigger::Trigger_Debug = 1<br>
Trigger::Trigger_Once[0]=1<br>
Trigger::Trigger_Number =1<br>
Trigger::Trigger_Checked_Variable[0]="GRHydro::GRHydro_rho_min"<br>
Trigger::Trigger_Reduction [0]=""<br>
Trigger::Trigger_Relation [0]="<"<br>
Trigger::Trigger_Checked_Value [0]= 1.1<br>
Trigger::Trigger_Reaction [0]="steerparam"<br>
Trigger::Trigger_Steered_Parameter_Value [0] = "19.32"<br>
Trigger::Trigger_Steered_Parameter_Name [0] = "radius_1[6]"<br>
Trigger::Trigger_Steered_Parameter_Thorn [0] = "CarpetRegrid2"<br>
<br>
(We use a trivially satisfied trigger condition so that the parameter<br>
value is reset immediately upon restart.)<br>
Since the CarpetRegrid2 variable "radius_1[6]" is not steerable, we<br>
had to modify this to make it such. After turning debug on, we see the<br>
following output:<br>
<br>
INFO (Trigger): Testing triggers<br>
INFO (Trigger): last_checked: -1<br>
INFO (Trigger): trigger nr. 0 fulfilled for GRHydro_rho_min (0.000000<1.100000)<br>
INFO (Trigger): Steering parameter<br>
INFO (Trigger): Parameter steered<br>
<br>
Which supposedly successfully resets the value correctly. We believe<br>
that even though this may be the case, the actual grid is not<br>
changing. We plot the first 2 boundary levels and show them in the<br>
attached file.<br>
<br>
There must be interpolation of the second innermost grid onto the<br>
innermost grid, and we believe this is not actually handled<br>
mid-simulation and thus doesn't not result in the change we expect.<br>
<br>
Is there any simple way of changing the inner grid easily?<br><div>
Any help is much appreciated!</div><div><table class="gmail-ajC" cellpadding="0"><tbody><tr class="gmail-ajv"><td colspan="2" tabindex="0" class="gmail-gG"><span class="gmail-gI"></span><br><img src="cid:ii_l2kycqzs0" alt="image(1).png" width="357" height="558"><br><br></td><td colspan="2" tabindex="0" class="gmail-gL"><br></td></tr></tbody></table></div>Thanks,<br>
Dhruv Desai</div>