def METHOD_NAME(draw, header=grid_heads(), zcorn=zcorns):
    grid_head = draw(header)
    dims = (grid_head.num_x, grid_head.num_y, grid_head.num_z)
    corner_size = (dims[0] + 1) * (dims[1] + 1) * 6
    coord = arrays(
        shape=corner_size,
        dtype=np.float32,
        elements=finites,
    )
    actnum = st.one_of(
        st.just(None),
        arrays(
            shape=dims[0] * dims[1] * dims[2],
            dtype=np.int32,
            elements=st.integers(min_value=0, max_value=3),
        ),
    )
    return xtge.GlobalGrid(
        coord=draw(coord),
        zcorn=draw(zcorn(dims)),
        actnum=draw(actnum),
        grid_head=grid_head,
        coord_sys=draw(map_axes),
        boxorig=draw(st.tuples(indices, indices, indices)),
        corsnum=draw(
            arrays(elements=indices, dtype="int32", shape=indices),
        ),
    )