No this should not be. Can  you select you single epoch after computing the noise csd?

Also instead `np.zeros(np.shape(stcs[0].data))` use `np.zeros_like(stcs[0].data)`