COSOlver_get_started.m
b3f94e22e87cmaster

COSOlver_get_started.m
View Options

			%% COSOlver : Get started
			% -----------------------
			%
			% Loads the basic input parameters for default case.
			%
			%
			%
			fprintf('COSOlver Get started ... \n')

			% Input parameters

			% BASIC Namelist
			BASIC.eecolls = '.true.';
			BASIC.iicolls = '.true. ' ;
			BASIC.eicolls = '.true.';
			BASIC.iecolls = '.true.';

			% Hermite-Laguerre Resolution

			BASIC.Pmaxi = 5;
			BASIC.Jmaxi = 5;
			BASIC.Pmaxe = 5;
			BASIC.Jmaxe = 5;

			% Coulomb collision sum bound
			BASIC.pamaxx = 5 ; % ! upper bound collision operator first sum first species
			BASIC.jamaxx = 5 ; %! upper bound collision operator second sum first species
			BASIC.pbmaxx = 5; % ! upper bound collision operator first sum second species
			BASIC.jbmaxx = 5 ; % ! upper bound collision operator second sum second species

			% 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 = 1;
			OPERATOR_MODEL.ITEST= 1 ;
			OPERATOR_MODEL.EBACK= 1 ;
			OPERATOR_MODEL.IBACK= 1 ;
			OPERATOR_MODEL.ESELF = 1;
			OPERATOR_MODEL.ISELF = 1;
			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;
			clear('BASIC','MODEL_PAR','OPERATOR_MODEL','INPUT');
			fprintf('COSOlver default case loaded ... \n');