C3PO RAY TRACING DEFAULT PARAMETERS Equilibrium parameters for the ray-tracing
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