[Bug] Data dimensions: silent unroll

[FlorianDeconinck]

This code fails because the data dims get unrolled

import numpy as np
from gt4py.cartesian import gtscript
from gt4py.cartesian.gtscript import (
    PARALLEL,
    computation,
    interval,
)
from gt4py import storage as gt_storage

backend = "numpy"
full_shape = (1, 1, 1)
aligned_index = (0, 0, 0)

datadim_len = 4

FLOAT64_NDDIM = (np.float64, (datadim_len,))


@gtscript.stencil(backend=backend)
def stencil(
    out_vec_field: gtscript.Field[FLOAT64_NDDIM],
    vec_field: gtscript.Field[FLOAT64_NDDIM],
):
    with computation(PARALLEL), interval(...):
        # The indexation doesn't work, it will silently unrolls for all dims
        # instead of failing because out_vect_field is not indexed properly
        out_vec_field[0, 0, 0] = vec_field[0, 0, 0][2]


out_vec_field = gt_storage.ones(
    full_shape, backend=backend, aligned_index=aligned_index, dtype=FLOAT64_NDDIM
)

vec_field = gt_storage.ones(
    full_shape, backend=backend, aligned_index=aligned_index, dtype=FLOAT64_NDDIM
)
vec_field[:] = 2.0

stencil(
    out_vec_field,
    vec_field,
    origin=(0, 0, 0),
    domain=full_shape,
)

Parent: GEOS-ESM/SMT-Nebulae#36

---

• Disallow this silent unroll and feedback like [Bug] Data dimensions: wrong rhs has no feedback #44 that the indexation is not compatible