Support for 4d Data

[oelbert]

Description
This PR adds support for 4D data structures in NDSL, especially the translate tests, and also fixes a bug where TranslateFortranData2Py would cast everything to Float instead of respecting types.

How Has This Been Tested?
The PBL and surface layer translate tests have been run with these updates

Checklist:

• My code follows the style guidelines of this project
• I have performed a self-review of my own code
• I have commented my code, particularly in hard-to-understand areas
• I have made corresponding changes to the documentation
• My changes generate no new warnings
• Any dependent changes have been merged and published in downstream modules
• New check tests, if applicable, are included

[bensonr]

Remove the empty line.

I may have missed it, but there is no test to ensure ndims <=4. I'm guessing the gt4py.storage.from_array function will crash at that point since data is undefined when ndims >4.

[FlorianDeconinck]

Can we generalize this with something like

def _make_storage_data_Nd(
    data: Field,
    shape: Tuple[int, ...],
    start: Tuple[int, ...],
    *,
    backend: str,
) -> Field:
    buffer = zeros(shape, backend=backend)
    idx = tuple([slice(start, start + shape_dim) for shape_dim in data.shape])
    buffer[idx] = asarray(data, type(buffer))
    return buffer

and we just have to pass a default start if it doesn't exist (that can be build with something like tuple([0]*data.ndim)