Harmonise Field interpolation order to consistently be (time, depth, lat, lon) throughout the code

[erikvansebille]

Since Parcels v2, the order for Field interpolation has been (time, depth, lat, lon). However, under the hood, there are still multiple cases (especially in JIT-mode) where that order is different. E.g.

https://github.com/OceanParcels/parcels/blob/b6dd5c852134ac78d9e19af575dd8793fbce0fc7/parcels/field.py#L1236-L1240

Or
https://github.com/OceanParcels/parcels/blob/b6dd5c852134ac78d9e19af575dd8793fbce0fc7/parcels/field.py#L1552

For consistency, it might be good to harmonise the order of variables throughout the code-base

[erikvansebille]

Another place where the order is wrong is in the FieldOutOfBoundSurfaceError message (and similar messages):

Parcels/parcels/tools/statuscodes.py

Line 57 in 5fff42c

message = f"{field.name if field else 'Field'} sampled out-of-bound, at ({self.x}, {self.y}, {self.z})"