Test case:

Run any existing PMF case in PeleC, remove existing refinement criteria, and add, for example, the following refinement criteria:
tagging.refinement_indicators = gtemp
tagging.gtemp.adjacent_difference_greater = 10
tagging.gtemp.field_name = Temp
tagging.gtemp.max_level = 1

Result:

Running with AMR max level set to 0, the simulation runs fine. Setting the max level to 1 breaks the simulation with a segfault. Turning on debug mode shows that a 4 dimensional array (state?) is being referenced out of bounds (-1:16, -1:16, -1:16, 0:0).