If I understood that correctly you have 3 axis (x,y,z) and a 'datacube' full of values d(x,y,z).
If you wan to pot a surface the easiest way is
surf(X,X,Z,C)
Here Z values are plottet over the XY plane and C defines to color.
Now your 4D data set is actually not one surface but 93 flat 'layers' stacked on top of each other over the XY plane so it's probably the easiest to just plot the Z layers individually. Here is a minimum example plotting your data that way:
clear all; close all; clc
load profileNethreed;
x=[36:1:42];
y=[26:1:45];
z=[80:10:1000];
d=profileNethreed;
d = permute(d,[2 1 3]); %matchning dimensions of d to the meshgrid
[LAT,LON, HEIGHT] = meshgrid(x,y,z);
for n = 1:length(z)
surf(LAT(:,:,n),LON(:,:,n),HEIGHT(:,:,n),d(:,:,n))
hold on
end
An alterntive good way of viualizing you data is imshow3D
clear all; close all; clc
load profileNethreed;
d=profileNethreed;
LOW = min(d(:), [], 'all');
HIGH = max(d(:), [], 'all');
imshow3D (d, [LOW HIGH])
