Description

This is a regression in 5.0x. When evaluating a model with one output and specifying with_bounding_box=True
the output is truncated

Expected behavior

from astropy.modeling.models import *

x = [1,2]
y = [3,4]

model = Mapping((1,))
model(x, y, with_bounding_box = True)
[3, 4]

Actual behavior

model(x, y, with_bounding_box = True)
3