diff --git a/wk/parameters/lin_AUG.m b/wk/parameters/lin_AUG.m
index 2d3e19bc32aaabeb1981c5568b7de5cb9d66851f..bfa0eddbe45e7d4fb8e8c666937c00d632e8abbe 100644
--- a/wk/parameters/lin_AUG.m
+++ b/wk/parameters/lin_AUG.m
@@ -46,12 +46,11 @@ PB_PHASE= 0;
 %% TIME PARAMETERS
 TMAX     = 5;  % Maximal time unit
 DT       = 2e-3;   % Time step
-DTSAVE0D = 1;      % Sampling per time unit for 0D arrays
-DTSAVE2D = -1;     % Sampling per time unit for 2D arrays
-DTSAVE3D = 3;      % Sampling per time unit for 3D arrays
-DTSAVE5D = 100;     % Sampling per time unit for 5D arrays
+DTSAVE0D = 0.1;      % Sampling time for 0D arrays
+DTSAVE2D = -1;     % Sampling time for 2D arrays
+DTSAVE3D = 0.1;      % Sampling time for 3D arrays
+DTSAVE5D = 1;     % Sampling time for 5D arrays
 JOB2LOAD = -1;     % Start a new simulation serie
-
 %% OPTIONS
 LINEARITY = 'linear';   % activate non-linearity (is cancelled if KXEQ0 = 1)
 CO        = 'DG';       % Collision operator (LB:L.Bernstein, DG:Dougherty, SG:Sugama, LR: Lorentz, LD: Landau)
diff --git a/wk/parameters/lin_DTT_AB_rho85_PT.m b/wk/parameters/lin_DTT_AB_rho85_PT.m
index 4fa60133f5ffb1811c686a3f74ed2fc6900cee96..262c352546fca43c1995c27eadd08318c2505c42 100644
--- a/wk/parameters/lin_DTT_AB_rho85_PT.m
+++ b/wk/parameters/lin_DTT_AB_rho85_PT.m
@@ -45,10 +45,10 @@ PB_PHASE = 0;
 %% TIME PARAMETERS
 TMAX     = 15;  % Maximal time unit
 DT       = 1e-3;   % Time step
-DTSAVE0D = 1;      % Sampling per time unit for 0D arrays
-DTSAVE2D = -1;     % Sampling per time unit for 2D arrays
-DTSAVE3D = 2;      % Sampling per time unit for 3D arrays
-DTSAVE5D = 100;     % Sampling per time unit for 5D arrays
+DTSAVE0D = 0.5;      % Sampling time for 0D arrays
+DTSAVE2D = -1;     % Sampling time for 2D arrays
+DTSAVE3D = 0.5;      % Sampling time for 3D arrays
+DTSAVE5D = 100;     % Sampling time for 5D arrays
 JOB2LOAD = -1;     % Start a new simulation serie
 
 %% OPTIONS
diff --git a/wk/parameters/lin_ETG_adiab_i.m b/wk/parameters/lin_ETG_adiab_i.m
index b3b3d7ed66d0eca85839704488aea817953df80b..219bcce7af0d850ef57aa7987bf89f25a8efd115 100644
--- a/wk/parameters/lin_ETG_adiab_i.m
+++ b/wk/parameters/lin_ETG_adiab_i.m
@@ -47,10 +47,10 @@ PB_PHASE= 0;
 %% TIME PARAMETERS
 TMAX     = 50;  % Maximal time unit
 DT       = 1e-4;   % Time step
-DTSAVE0D = 1;      % Sampling per time unit for 0D arrays
-DTSAVE2D = -1;     % Sampling per time unit for 2D arrays
-DTSAVE3D = 1;      % Sampling per time unit for 3D arrays
-DTSAVE5D = 100;     % Sampling per time unit for 5D arrays
+DTSAVE0D = 1.0;    % Sampling time for 0D arrays
+DTSAVE2D = -1;     % Sampling time for 2D arrays
+DTSAVE3D = 1.0;    % Sampling time for 3D arrays
+DTSAVE5D = 100;    % Sampling time for 5D arrays
 JOB2LOAD = -1;     % Start a new simulation serie
 
 %% OPTIONS
diff --git a/wk/parameters/lin_Entropy.m b/wk/parameters/lin_Entropy.m
index 18a3c2ad714091e40d46d2afcf811db471ab4952..ea0b655db7e4e66f06ff53dea77bfe562c9ac8ce 100644
--- a/wk/parameters/lin_Entropy.m
+++ b/wk/parameters/lin_Entropy.m
@@ -47,10 +47,10 @@ PB_PHASE= 0;
 %% TIME PARAMETERS
 TMAX     = 100;  % Maximal time unit
 DT       = 1e-3;   % Time step
-DTSAVE0D = 1;      % Sampling per time unit for 0D arrays
-DTSAVE2D = -1;     % Sampling per time unit for 2D arrays
-DTSAVE3D = 1;      % Sampling per time unit for 3D arrays
-DTSAVE5D = 100;     % Sampling per time unit for 5D arrays
+DTSAVE0D = 1.0;    % Sampling time for 0D arrays
+DTSAVE2D = -1;     % Sampling time for 2D arrays
+DTSAVE3D = 2.0;    % Sampling time for 3D arrays
+DTSAVE5D = 100;    % Sampling time for 5D arrays
 JOB2LOAD = -1;     % Start a new simulation serie
 
 %% OPTIONS
diff --git a/wk/parameters/lin_ITG.m b/wk/parameters/lin_ITG.m
index 6519a58f31a5d685e45ac6330efffde09de31414..7e2167a3f8646b08f56338a158002d9e295063ae 100644
--- a/wk/parameters/lin_ITG.m
+++ b/wk/parameters/lin_ITG.m
@@ -42,10 +42,10 @@ NPOL   = 1;       % Number of poloidal turns
 %% TIME PARAMETERS
 TMAX     = 50;  % Maximal time unit
 DT       = 1e-2;   % Time step
-DTSAVE0D = 1;      % Sampling per time unit for 0D arrays
-DTSAVE2D = -1;     % Sampling per time unit for 2D arrays
-DTSAVE3D = 2;      % Sampling per time unit for 3D arrays
-DTSAVE5D = 100;     % Sampling per time unit for 5D arrays
+DTSAVE0D = 1;      % Sampling time for 0D arrays
+DTSAVE2D = -1;     % Sampling time for 2D arrays
+DTSAVE3D = 1;      % Sampling time for 3D arrays
+DTSAVE5D = 100;     % Sampling time for 5D arrays
 JOB2LOAD = -1;     % Start a new simulation serie
 
 %% OPTIONS
diff --git a/wk/parameters/lin_Ivanov.m b/wk/parameters/lin_Ivanov.m
index b1b7e7eb78be86b697d965abddffb1946594e2be..a638d755eca8e64b372504c7a5f5d6abb581518c 100644
--- a/wk/parameters/lin_Ivanov.m
+++ b/wk/parameters/lin_Ivanov.m
@@ -43,10 +43,10 @@ NPOL   = 1;       % Number of poloidal turns
 %% TIME PARAMETERS
 TMAX     = 50;  % Maximal time unit
 DT       = 1e-2;   % Time step
-DTSAVE0D = 1;      % Sampling per time unit for 0D arrays
-DTSAVE2D = -1;     % Sampling per time unit for 2D arrays
-DTSAVE3D = 1;      % Sampling per time unit for 3D arrays
-DTSAVE5D = 20;     % Sampling per time unit for 5D arrays
+DTSAVE0D = 1;      % Sampling time for 0D arrays
+DTSAVE2D = -1;     % Sampling time for 2D arrays
+DTSAVE3D = 1;      % Sampling time for 3D arrays
+DTSAVE5D = 25;     % Sampling time for 5D arrays
 JOB2LOAD = -1;     % Start a new simulation serie
 
 %% OPTIONS
diff --git a/wk/parameters/lin_KBM.m b/wk/parameters/lin_KBM.m
index dab2fe08308a13edcecfe55bd7658366c7aec5c2..fc2000de7a634cb90ad735bfe510dc2240bec4a0 100644
--- a/wk/parameters/lin_KBM.m
+++ b/wk/parameters/lin_KBM.m
@@ -42,10 +42,10 @@ NPOL   = 1;       % Number of poloidal turns
 %% TIME PARAMETERS
 TMAX     = 15;  % Maximal time unit
 DT       = 5e-3;   % Time step
-DTSAVE0D = 1;      % Sampling per time unit for 0D arrays
-DTSAVE2D = -1;     % Sampling per time unit for 2D arrays
-DTSAVE3D = 2;      % Sampling per time unit for 3D arrays
-DTSAVE5D = 100;     % Sampling per time unit for 5D arrays
+DTSAVE0D = 1.0;    % Sampling time for 0D arrays
+DTSAVE2D = -1;     % Sampling time for 2D arrays
+DTSAVE3D = 0.5;    % Sampling time for 3D arrays
+DTSAVE5D = 100;    % Sampling time for 5D arrays
 JOB2LOAD = -1;     % Start a new simulation serie
 
 %% OPTIONS