diff --git a/matlab/prior_posterior_statistics.m b/matlab/prior_posterior_statistics.m
index 9c6cffcf0dcd1e59a0dc8e226b0f17fb4935c92c..a492b80f77beb35631b9208cfedfd9f091521665 100644
--- a/matlab/prior_posterior_statistics.m
+++ b/matlab/prior_posterior_statistics.m
@@ -228,14 +228,18 @@ else
         ifil(2,j+1) =ifil(2,j)+nfiles;
         nfiles = ceil(nBlockPerCPU(j)/MAX_nerro);
         ifil(3,j+1) =ifil(3,j)+nfiles;
-        nfiles = ceil(nBlockPerCPU(j)/MAX_naK);
-        ifil(4,j+1) =ifil(4,j)+nfiles;
+        if naK
+            nfiles = ceil(nBlockPerCPU(j)/MAX_naK);
+            ifil(4,j+1) =ifil(4,j)+nfiles;
+        end
         nfiles = ceil(nBlockPerCPU(j)/MAX_nruns);
         ifil(5,j+1) =ifil(5,j)+nfiles;
-        nfiles = ceil(nBlockPerCPU(j)/MAX_nforc1);
-        ifil(6,j+1) =ifil(6,j)+nfiles;
-        nfiles = ceil(nBlockPerCPU(j)/MAX_nforc2);
-        ifil(7,j+1) =ifil(7,j)+nfiles;
+        if horizon
+            nfiles = ceil(nBlockPerCPU(j)/MAX_nforc1);
+            ifil(6,j+1) =ifil(6,j)+nfiles;
+            nfiles = ceil(nBlockPerCPU(j)/MAX_nforc2);
+            ifil(7,j+1) =ifil(7,j)+nfiles;
+        end
         %       nfiles = ceil(nBlockPerCPU(j)/MAX_momentsno);
         %       ifil(8,j+1) =ifil(8,j)+nfiles;
     end