<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 --&gt; ~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 = &quot;Trigger&quot;<br>
Trigger::Trigger_Debug = 1<br>
Trigger::Trigger_Once[0]=1<br>
Trigger::Trigger_Number =1<br>
Trigger::Trigger_Checked_Variable[0]=&quot;GRHydro::GRHydro_rho_min&quot;<br>
Trigger::Trigger_Reduction       [0]=&quot;&quot;<br>
Trigger::Trigger_Relation        [0]=&quot;&lt;&quot;<br>
Trigger::Trigger_Checked_Value   [0]= 1.1<br>
Trigger::Trigger_Reaction        [0]=&quot;steerparam&quot;<br>
Trigger::Trigger_Steered_Parameter_Value [0] = &quot;19.32&quot;<br>
Trigger::Trigger_Steered_Parameter_Name  [0] = &quot;radius_1[6]&quot;<br>
Trigger::Trigger_Steered_Parameter_Thorn [0] = &quot;CarpetRegrid2&quot;<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 &quot;radius_1[6]&quot; 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&lt;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&#39;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>