4-dimensional interpolation of .inp files

Hello, I want to make an 4 dimensional Interpolation (x,y,z,t) of data extracted from .inp files. So my data looks like this: I have a n*3 double which represents the x,y,z coordinates of each Node

-2          0          -99.5000
-1.7321     -1.0000    -99.5000
0           -2         -99.5000
...
-1.6194     -0.6500    -30.9843
33.9000     20.6922    0.9872
...

Than I have a Struct with my Values per time and node looking like:

data.t1_1234:          data.t2_3456:         data.4_5678:
123456                 134562                187635
123455                 134533                187601
120321                 99882                 35403
...
2465990                1000000               3409330
100000                 100000                123456
...

Where the each entry of each data.t[...] is the value at the equvalent position in node in that given time.

e.g.: the Node with the coordiantes (x,y,z)

-2    0    -99.5000

has the value 123456 at time t=1.1234 and the value 134562 at time t=2.3456.

Now I want to give that function the coordinates and a time and want an interpolated value for it as the result.

So as Input I want to have (x,y,z,t) and interpolate my data in 4 dimensions as the coordinates and time most likely won't be exact.

I thought about using a 4D Interpolation with interpn, but I honestly didn't understood how to use it in 4 dimensions and I don't even know if I can use my data like that or if I need to format it.

An other thought was to do two 3D Interpolation. One after my given time and one before. Than just interpolate the two values. I hope that I can save some time with that method. Would this one be smart to save some time?

Any suggestions how I can interpolate my data in 4 Dimensions?

Best regards,