diff --git a/Sound Fields/MATLAB/Functions/ptSrcFieldTD.m b/Sound Fields/MATLAB/Functions/ptSrcFieldTD.m
index 422cb6ead77fc2e3aae39345dc04881c21a7b164..0c17c25af6c913f75c6016eeeecd0b0f667f50ee 100644
--- a/Sound Fields/MATLAB/Functions/ptSrcFieldTD.m	
+++ b/Sound Fields/MATLAB/Functions/ptSrcFieldTD.m	
@@ -170,7 +170,7 @@ function [rSig, rSigMean, rSigMtx, rSigMtxMean, Q] = ptSrcFieldTD(sPos, rPos, fs
     end
 
     % Sum source signals at each receiver position
-    rSig = squeeze(sum(rSigMtx, 3));
+    rSig = reshape(sum(rSigMtx, 3), size(rSigMtx, 1), size(rSigMtx, 2), size(rSigMtx, 4));
 
     % ====================================================
     % Calculate output arguments