Let's suppose "plate" is the temperature distribution.
[M,N] = size(plate);
%Pay attention:
% M --> raws N --> coloumns
%BUT for meshgrid you must introduce (x,y) coordinates, so:
% M --> y N--> x
[x,y] = meshgrid(1:N,1:M);
surf(x,y,plate);
If plate would exceed [0-100] bounds, you could use a support variable:
plate_2 = plate;
plate_2(plate_2>100) = 100;
plate_2(plate_2<0) = 0;
surf(x,y,plate_2);
