This is often required in analyzing climate data from netcdf files. Climate variables are stored in 3-d matrix (lat, lon, time). I want to find correlation coefficient between two 3-d matrices at each grid cell (i.e. correlation along time or 3rd dimension). Is there an easier way to do this, without reshaping, for geometric intuition of the grid map? Also how to get the p-value at each grid cell after calculating the correlation?