Homec4science
Diffusion COSOlver (llechot) Browse

cosolver/matlab_scripts/COSOlver_PamaxxScan.m
b3f94e22e87cllechot

Last Change
Raw File

COSOlver_PamaxxScan.m
View Options

% Pamaxx-scan
%% Clear current directory
clear all; close all;
%% Define COSOlver inputs
% BASIC Namelist
BASIC.eecolls = '.true.';
BASIC.iicolls = '.false.' ;
BASIC.eicolls = '.false.';
BASIC.iecolls = '.false.';
% Hermite-Laguerre Resolution
BASIC.Pmaxi = 5;
BASIC.Jmaxi = 5;
BASIC.Pmaxe = 5;
BASIC.Jmaxe = 5;
% Coulomb collision sum bound
BASIC.pamaxx = 5 ;
BASIC.jamaxx = 5 ;
BASIC.pbmaxx = 5;
BASIC.jbmaxx = 5 ;
% FLR sum bound
BASIC.nemaxxFLR = 0;
BASIC.nimaxxFLR = 0;
BASIC.OVERWRITE = '.true.'; % OVERWRITE existing resfile if true
% MODEL_PAR Namelist
MODEL_PAR.nu = 1; % ... normalized collision freq.
MODEL_PAR.mu = 0.023; % ... sqrt(m_e/m_i)
MODEL_PAR.tau = 0.01; % ... T_i/T_e
MODEL_PAR.kperp = 0;
% OPERATOR_MODEL Namelist, load Coulomb collision operator
OPERATOR_MODEL.ETEST = 0;
OPERATOR_MODEL.ITEST = 0 ;
OPERATOR_MODEL.EBACK = 0 ;
OPERATOR_MODEL.IBACK = 0 ;
OPERATOR_MODEL.ESELF = 1;
OPERATOR_MODEL.ISELF = 0;
OPERATOR_MODEL.GKE = 0;
OPERATOR_MODEL.GKI = 0;
% Create input structure
INPUT.BASIC = BASIC;
INPUT.MODEL_PAR = MODEL_PAR;
INPUT.OPERATOR_MODEL = OPERATOR_MODEL;
COSOlver_input = INPUT;
%% Define pamaxx scan
pamaxxs = 4:1:20;
%% Perform scan
normCee = zeros(length(pamaxxs),1);
for ip=1:length(pamaxxs)
BASIC.pamaxx = pamaxxs(ip) ;
BASIC.jamaxx = pamaxxs(ip) ;
BASIC.pbmaxx = pamaxxs(ip);
BASIC.jbmaxx = pamaxxs(ip) ;
INPUT.BASIC = BASIC;
resfile_SUFFIX = COSOlver_create_input(INPUT);
commandCOSOlver = './bin/CO';
[status,cmdout] = system(commandCOSOlver,'-echo');
% Load matrix from resfile *.h5
Cabpj = h5read(['./iCa/self_',resfile_SUFFIX],'/Caapj/Ceepj');
normCee(ip) = norm(Cabpj);
end
%% Plot convergence plot
figure;
plot(pamaxxs,normCee,'*b');
xlabel('J_{amaxx}');
ylabel('||C_{ee}||');
c4science · Help