From a67434b804e923ee14b4ffa455f3eb9bf85b560a Mon Sep 17 00:00:00 2001 From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch> Date: Fri, 10 Jul 2020 09:18:12 +0200 Subject: [PATCH] added kr parameter for MOLI --- matlab/MOLI_kperp_scan.m | 4 +++- matlab/MOLI_time_solver.m | 4 +++- 2 files changed, 6 insertions(+), 2 deletions(-) diff --git a/matlab/MOLI_kperp_scan.m b/matlab/MOLI_kperp_scan.m index 8697f988..c3c8d4ef 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 0bd7d118..b2f25e39 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 -- GitLab