Visualize fold change over multidimensional array

Hello,

I have several large 3D arrays (for example, of dimensions around 1560x21x48 (double)) containing flux data, where 1560 would be the reaction over 21 instances (according to a function) across 48 conditions. I would like to visualize the change (like a heatmap) in flux for each reaction according to how they change over the 48 conditions. So, the table would be 1560x21, and the colors would be the amount of change. I'm not dearly attached to any specific way of doing this, as I simply want a visual representation to get an overview of how much change is taking place.

So far I'm just using a bare bones fraction approach

z=fluxes.AT_spent(:,:,1)./fluxes.AT_spent(:,:,2:end)

and then taking the mean for each i,j.

This does indeed result in a chart where of course 1 is neutral, in a 1560x21 matrix where values are the changes around 1, but the heatmap is just a black image, even though the scale uses blue for 1. I must be missing something, clearly.

I'm sure there's a more elegant approach though, since my approach clearly lacks finesse.

Any advice is appreciated!