Hi all,
I am trying to change the num_levels for CarpetRegrid2 using the trigger, inspired from the gallery example at ET website. The trigger seems to kick in at the right place, as the simulation stopped while regridding for new num_levels.
The dedicated lines in parfile are as follows:
*CarpetRegrid2::freeze_unaligned_levels = "yes"CarpetRegrid2::movement_threshold_1 = 0.10CarpetRegrid2::movement_threshold_2 = 0.10CarpetRegrid2::num_centres = 2CarpetRegrid2::num_levels_1 = 7 #10CarpetRegrid2::num_levels_2 = 7CarpetRegrid2::position_x_1 = -13.543564CarpetRegrid2::position_x_2 = 13.543564CarpetRegrid2::radius_1[1] = 270.0CarpetRegrid2::radius_1[2] = 160.0CarpetRegrid2::radius_1[3] = 90.0CarpetRegrid2::radius_1[4] = 40.0CarpetRegrid2::radius_1[5] = 25.0CarpetRegrid2::radius_1[6] = 9.0CarpetRegrid2::radius_1[7] = 6.0CarpetRegrid2::radius_1[8] = 3.0CarpetRegrid2::radius_1[9] = 1.5CarpetRegrid2::radius_2[1] = 270.0CarpetRegrid2::radius_2[2] = 160.0CarpetRegrid2::radius_2[3] = 90.0CarpetRegrid2::radius_2[4] = 40.0CarpetRegrid2::radius_2[5] = 25.0CarpetRegrid2::radius_2[6] = 9.0CarpetRegrid2::regrid_every = 32CarpetRegrid2::symmetry_rotating180 = "no"CarpetRegrid2::verbose = "yes"*
*Trigger::Trigger_Number = 1Trigger::Trigger_Checked_Variable[0]="ADMBase::alp"Trigger::Trigger_Reduction [0]="minimum"Trigger::Trigger_Relation [0]="<"Trigger::Trigger_Checked_Value [0]=0.1Trigger::Trigger_Reaction [0]="steerscalar"Trigger::Trigger_Steered_Scalar [0] = "CarpetRegrid2::num_levels[0]" # == num_levels_1Trigger::Trigger_Steered_Scalar_Value[0] = "10"*
The outfile reports the following error:
*INFO (CarpetRegrid2): Centre 1 is at position [0.000191809,0.000291781,-0.00018496] with 10 levelsINFO (CarpetRegrid2): Centre 2 is at position [0.000191809,0.000291781,-0.00018496] with 7 levelsINFO (CarpetRegrid2): Regridding level 0 map 0 at iteration 4129 time 806.25INFO (CarpetRegrid2): Regridding levels 1 and up [1mERROR from host astro-Z390-AORUS-ELITE process 0 in thorn CarpetRegrid2, file /home/astro/Documents/ET/Cactus/configs/sim/build/CarpetRegrid2/regrid.cc:89: -> [0m Region 1 has 10 levels active, which is larger than the maximum number of refinement levels 7*
I could not figure out where it went wrong. The minimum lapse at iteration 4128 is 0.095, so the condition in trigger seems to be working fine. Some insight into it would be helpful. Thanks in advance.
Regards Shamim Haque Senior Research Fellow (SRF) Department of Physics IISER Bhopal