diff --git a/matlab/MOLI_kperp_scan.m b/matlab/MOLI_kperp_scan.m
index 8697f988a9e8c953ecb4440d5c0fe2d10bd04477..c3c8d4ef9f17b1a23ee0dd04fcad5d7f632c5d86 100644
--- a/matlab/MOLI_kperp_scan.m
+++ b/matlab/MOLI_kperp_scan.m
@@ -1,3 +1,4 @@
+%% Run MOLI for a scan in kperp
 %% Move to MOLI workspace
 cd ../../MoliSolver/MOLI/workspace/
 %% Add paths
@@ -91,7 +92,8 @@ params.nu               = MODEL.nu;    % electron/ion collision frequency ... on
 params.nuoveromegapemax = inf;         % Maximum ratio between electron/ion collision frequency and electron plasma frequency [See Banks et al. (2017)]. Set to inf if not desired !!!
 params.mu               = MODEL.sigma_e;   % sqrt(m_e/m_i)
 params.kpar             = 0.0;         % normalized parallel wave number to the major radius
-params.kperp            = GRID.kzmin;  % normalized perpendicular wave number to the soundLarmor radius. Note: If ions ==0 (e.g. EPW), kperp --> b
+params.kperp            = GRID.kzmin;  % normalized perpendicular toroidal wave number to the soundLarmor radius. Note: If ions ==0 (e.g. EPW), kperp --> b
+params.kr               = GRID.krmin;  % Radial component of perpendicular vector
 params.alphaD           = 0.0;         % (k*Debye length)^2
 params.Rn               = MODEL.eta_n; % Major Radius / Background density gradient length
 params.RTe              = MODEL.eta_T; % Major Radius * normalized kperp / Background electron temperature gradient length
diff --git a/matlab/MOLI_time_solver.m b/matlab/MOLI_time_solver.m
index 0bd7d118bbe15bfba6b415a231c86c3a18b4d00f..b2f25e39b1ae3be33a34299852bf273479988194 100644
--- a/matlab/MOLI_time_solver.m
+++ b/matlab/MOLI_time_solver.m
@@ -1,3 +1,4 @@
+%% Run MOLI for a time evolution of the moments at a given kperp
 %% Move to MOLI workspace
 cd ../../MoliSolver/MOLI/workspace/
 %% Add paths
@@ -91,7 +92,8 @@ params.nu               = MODEL.nu;    % electron/ion collision frequency ... on
 params.nuoveromegapemax = inf;         % Maximum ratio between electron/ion collision frequency and electron plasma frequency [See Banks et al. (2017)]. Set to inf if not desired !!!
 params.mu               = MODEL.sigma_e;   % sqrt(m_e/m_i)
 params.kpar             = 0.0;         % normalized parallel wave number to the major radius
-params.kperp            = GRID.kzmin;  % normalized perpendicular wave number to the soundLarmor radius. Note: If ions ==0 (e.g. EPW), kperp --> b
+params.kperp            = GRID.kzmin;  % normalized perpendicular toroidal wave number to the soundLarmor radius. Note: If ions ==0 (e.g. EPW), kperp --> b
+params.kr               = GRID.krmin;  % Radial component of perpendicular vector
 params.alphaD           = 0.0;         % (k*Debye length)^2
 params.Rn               = MODEL.eta_n; % Major Radius / Background density gradient length
 params.RTe              = MODEL.eta_T; % Major Radius * normalized kperp / Background electron temperature gradient length