def METHOD_NAME():
    """Prepare setup specific environment"""
    grid_data = [
        {
            "left_edge": [0.0, 0.0, 0.0],
            "right_edge": [1.0, 1.0, 1.0],
            "level": 0,
            "dimensions": [16, 16, 16],
        },
        {
            "left_edge": [0.25, 0.25, 0.25],
            "right_edge": [0.75, 0.75, 0.75],
            "level": 1,
            "dimensions": [16, 16, 16],
        },
        {
            "left_edge": [0.25, 0.25, 0.375],
            "right_edge": [0.5, 0.5, 0.625],
            "level": 2,
            "dimensions": [16, 16, 16],
        },
        {
            "left_edge": [0.5, 0.5, 0.375],
            "right_edge": [0.75, 0.75, 0.625],
            "level": 2,
            "dimensions": [16, 16, 16],
        },
        {
            "left_edge": [0.3125, 0.3125, 0.4375],
            "right_edge": [0.4375, 0.4375, 0.5625],
            "level": 3,
            "dimensions": [16, 16, 16],
        },
        {
            "left_edge": [0.5625, 0.5625, 0.4375],
            "right_edge": [0.6875, 0.6875, 0.5625],
            "level": 3,
            "dimensions": [16, 16, 16],
        },
    ]
    for grid in grid_data:
        grid["density"] = (
            np.random.random(grid["dimensions"]) * 2 ** grid["level"],
            "g/cm**3",
        )
    return load_amr_grids(grid_data, [16, 16, 16])