make_wave_C3POparam_EC

PURPOSE ^

C3PO RAY TRACING DEFAULT PARAMETERS

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 C3PO RAY TRACING DEFAULT PARAMETERS


 Equilibrium parameters for the ray-tracing

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 %
0002 % C3PO RAY TRACING DEFAULT PARAMETERS
0003 %
0004 %
0005 % Equilibrium parameters for the ray-tracing
0006 %
0007 fitparam.mode_equil = 1;%Magnetic equilibrium grid type: (1): (psi-theta), (2): (x-y)
0008 fitparam.method = 'spline';%nearest,spline,pchip
0009 fitparam.nharm = NaN;%Number of harmonics in the magnetic equilibrium interpolation (less than ntheta_equil/2)
0010 fitparam.ngridresample = 1001;%Number of grid points for resampling the radial profile of magnetic equilibrium parameters
0011 fitparam.opt_load = 0;%Reload existing vectorial magnetic equilibrium (1) or overwrite it (0).
0012 %
0013 % Global parameters for the ray-tracing
0014 %
0015 rayparam.testmode = 0;
0016 rayparam.tensortype = 0;%(0:cold,1:warm,2:hot;3:weak realtivistic,4:full relativistic)
0017 rayparam.t0 = 0;
0018 rayparam.tfinal = 50000;
0019 rayparam.dt0 = 1.e-4;
0020 rayparam.dS = 1.e-3;
0021 rayparam.tol = 1e-8;%when tolerance is increased (less accurate calculation of D=0), tfinal must be increased accordingly
0022 rayparam.kmax = 50000;
0023 rayparam.ncyclharm = 3;%number of cyclotron harmonics (just for hot and relativistic dielectric tensors)
0024 rayparam.reflection = 0;%1:Enforce wave reflection at plasma boundary, 0: the code calculates itself if the ray must leave of not the plasma
0025 rayparam.rel_opt = 1;%option for (1) relativistic or (0) non-relativistic calculations
0026 rayparam.nperp = 1000;%number of points in pperp integration for damping calculations
0027 rayparam.pperpmax = 10;%maximum value of pperp in damping calculations
0028 rayparam.tau_lim = Inf;%value of optical depth beyond which the wave is considered absorbed
0029 rayparam.kextra = 1000;%extra points if QL absorption is delayed
0030 rayparam.colldamp = 0;%non-resonant collisional damping -> (1) actived (default), (0) not activated
0031 %
0032 % C3PO computing parameters
0033 %
0034 C3POparam.clustermode.main_C3PO_jd.scheduler.mode = 0;%MatLab distributed computing environment disabled (0), enabled with the dedicated toolbox (1), enabled with a private method (2) for the function coll_dke_jd.m (MDC toolbox must be installed for option 1)
0035 C3POparam.clustermode.main_C3PO_jd.scheduler.memory = 500;%required memory (in mb)
0036 %
0037 % C3PO Display parameters
0038 %
0039 C3POdisplay.ray = 0;
0040 C3POdisplay.equilibrium = 0;
0041 C3POdisplay.p_opt = -1;%Printing or saving option of the figures
0042 %
0043 % Wave parameters
0044 %
0045 waveparam.mmode = 0.5;%cold plasma mode [1] : (-1) m (1) p, p is the slow mode when kperp > 0 (slow/fast LH wave: -1/+1, O/X mode: -1/+1)
0046 waveparam.kmode = 0;%(0:cold,1:warm,2:hot;3:weak realtivistic,4:full relativistic)
0047 %
0048 %Option parameter for FLR effects and cross-comparison between old FP code:
0049 %    - (0): all FLR effects
0050 %    - (1): small FLR effects and 1/vpar dependence
0051 %    - (2): small FLR effects and no 1/vpar dependence and old grid technique for DQL calculations (Karney, Shoucri) (see rfdiff_dke_jd)
0052 %
0053 waveparam.opt_rf = 1;
0054 %
0055 waveparam.dsmin = 0;%minimum size for ray fragments
0056 %
0057 waveparam.nd = 4;%Number of transverse distance positions for beamlets
0058 waveparam.nchi = 6;%Number of angular positions for beamlets
0059 %
0060 waveparam.n_rf_list = 1:3;
0061 waveparam.ns = 1;%ray smoothing parameter
0062 %
0063 waveparam.rhomin = 0.01;%minimum rho allowed in C3PO (to avoid spurious kicks)
0064 %
0065 save('C3POPARAM_EC.mat','fitparam','rayparam','C3POparam','C3POdisplay','waveparam');
0066

Community support and wiki are available on Redmine. Last update: 18-Apr-2019.