Homec4science
Diffusion COSOlver (new_coulomb) Browse

cosolver/matlab/COSOlver_ToLorentz.m
41e7a23b16aenew_coulomb

Last Change
Raw File

COSOlver_ToLorentz.m
View Options

% sqrt(me/me) - scan
%% Clear current directory
clear all; close all;
%% Define COSOlver inputs
% BASIC Namelist
BASIC.eecolls = '.false.';
BASIC.iicolls = '.false.' ;
BASIC.eicolls = '.true.';
BASIC.iecolls = '.false.';
% Hermite-Laguerre Resolution
BASIC.Pmaxi = 5;
BASIC.Jmaxi = 5;
BASIC.Pmaxe = 10;
BASIC.Jmaxe = 4;
% Coulomb collision sum bound
BASIC.pamaxx = 6 ;
BASIC.jamaxx = 6 ;
BASIC.pbmaxx = 6;
BASIC.jbmaxx = 6 ;
% 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 = 1; % ... T_i/T_e
MODEL_PAR.kperp = 0;
% OPERATOR_MODEL Namelist, load Coulomb collision operator
OPERATOR_MODEL.ETEST = 1;
OPERATOR_MODEL.ITEST = 0 ;
OPERATOR_MODEL.EBACK = 0 ;
OPERATOR_MODEL.IBACK = 0 ;
OPERATOR_MODEL.ESELF = 0;
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;
%% Compute Lorentz matrix
OPERATOR_MODEL.ETEST =2;
INPUT.MODEL_PAR= MODEL_PAR;
resfile_SUFFIX = COSOlver_create_input(INPUT);
commandCOSOlver = './bin/CO';
[status,cmdout] = system(commandCOSOlver,'-echo');
% Load matrix from resfile *.h5
Lorentz= h5read(['./iCa/ei_',resfile_SUFFIX],'/Ceipj/CeipjT');
%% Define mu scan
mus = logspace(-8,0,5);
%% Perform scan using DK full-Coulomb
OPERATOR_MODEL.ETEST = 1;
normCeipjT = zeros(length(mus),1);
for ii=1:length(mus)
MODEL_PAR.mu = mus(ii);
INPUT.MODEL_PAR= MODEL_PAR;
resfile_SUFFIX = COSOlver_create_input(INPUT);
commandCOSOlver = './bin/CO';
[status,cmdout] = system(commandCOSOlver,'-echo');
% Load matrix from resfile *.h5
CeipjT(:,:,ii) = h5read(['./iCa/ei_',resfile_SUFFIX],'/Ceipj/CeipjT');
end
%% Plot norm of collisional matrices
%figure;
%loglog(mus,normCeipjT,'*r');
%xlabel('$\sqrt{m_e/m_i}$')
%ylabel('|| C_{ei}^T - C_{ei}^{L0} ||')
c4science · Help