make_fluct_ne_sigma0p1_delta0p02

PURPOSE ^

Parameters for testing the effect of density or magnetic field fluctuations on ray trajectories

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 Parameters for testing the effect of density or magnetic field fluctuations on ray trajectories

 by Y. Peysson (DRFC/DSM/CEA) <yves.peysson@cea.fr> and J. Decker (DRFC/DSM/CEA) <joan.decker@cea.fr>

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 %
0002 % Parameters for testing the effect of density or magnetic field fluctuations on ray trajectories
0003 %
0004 % by Y. Peysson (DRFC/DSM/CEA) <yves.peysson@cea.fr> and J. Decker (DRFC/DSM/CEA) <joan.decker@cea.fr>
0005 %
0006 clear all
0007 close all
0008 clc
0009 %
0010 id_fluct = 'ne_sigma0p1_delta0p02';%Plasma fluctuations identifier
0011 %
0012 % Path parameters
0013 %
0014 id_dkepath = '';%For all paths used by DKE solver
0015 path_dkepath = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0016 %
0017 % Equilibrium parameters
0018 %
0019 id_equil = 'ITER_Scen2_200103121816_430_129';%For plasma equilibrium
0020 path_equil = '../EQUIL/';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0021 %
0022 % Core density fluctuations
0023 %
0024 naequilp.core.fluct.type = {'ne'};%Type of fluctuations or modulations ('ne': electron density, 'B': magnetic field) {1,nfluct_types}
0025 naequilp.core.fluct.model = [5];%Statistical ne fluctuation model (>= 1) : [1,nfluct_types]
0026                         %    - (1) -> 3-D Gaussian model (rho,theta,phi), relative epsi values (benchmark of C3PO)
0027                         %    - (2) -> 3-D Gaussian model (rho,theta,phi), absolute epsi values (m)
0028                         %    - (3) -> 2-D Gaussian drift-wave like model global (rho,curtheta), absolute epsi values (m)
0029                         %    - (4) -> 1-D Gaussian drift-wave like model local (curtheta), absolute epsi values (m)
0030                         %    - (5) -> 1-D Gaussian drift-wave like model local (cn,cm), absolute epsi values (m)
0031                         %    - (6) -> 2-D Gaussian drift-wave like model local (rho,cn,cm), absolute epsi values (m)
0032                         %    - (7) -> *** NOT IMPLEMENTED *** 2-D Gaussian drift-wave like model local (rho,lperp), absolute epsi values (m)
0033 naequilp.core.fluct.sigmar_max = [0.1];%Maximum value of the relative fluctuations variance at the poloidal angle theta = 0 [1,nfluct_types]
0034 naequilp.core.fluct.sigmar_hwhm = [0.02];%Radial half width at half maximum of the relative fluctuations variance at the poloidal angle theta = 0 [1,nfluct_types]
0035 naequilp.core.fluct.sigmar_rho = [1];%Radial location where the relative fluctuations variance peaks at the poloidal angle theta = 0 [1,nfluct_types]
0036 naequilp.core.fluct.polmode = [0.1];%HFS/LFS relative fluctuations variance ratio. No poloidal dependency corresponds to 1 [1,nfluct_types]
0037 %
0038 naequilp.core.fluct.epsi_rho = [1];%
0039 naequilp.core.fluct.epsi_theta = [0.01];%theta direction is perp direction for models 3 & 4
0040 naequilp.core.fluct.epsi_phi = [0.01];%useless for models 3 & 4
0041 %
0042 naequilp.core.fluct.lmin = [1];%
0043 naequilp.core.fluct.mmin = [1];%for perp direction for models 3 & 4
0044 naequilp.core.fluct.nmin = [1];%useless for models 3 & 4
0045 %
0046 naequilp.core.fluct.lmax = [2];%
0047 naequilp.core.fluct.mmax = [1000];%for perp direction for models 3 & 4
0048 naequilp.core.fluct.nmax = [200];%useless for models 3 & 4
0049 %
0050 % Edge plasma fluctuations outside the separatrix (for LH wave only, run with run_lukert or a specific script)
0051 %
0052 naequilp.edge.fluct.dtn = 1;%Gaussian time correlation model, edge_dtn = 2.14*t_correlation (if edge_dtn is Inf, no fluctuations)
0053 naequilp.edge.fluct.hwhm_npar = 0.1;%Gaussian distribution (if hwhm is NaN or empty, no fluctuations)
0054 %
0055 % -------------------------------------------------------------------------
0056 %
0057 % Load structures
0058 %
0059 [equil] = load_structures_yp('equil',id_equil,path_equil);
0060 % =========================================================================
0061 %
0062 [fluct] = fluctstruct_yp(equil,id_fluct,naequilp);
0063 fluct.naequilp = naequilp;
0064 %
0065 % Data saving
0066 %
0067 save_str = ['FLUCT_',id_equil,'_',id_fluct,'.mat'];
0068 eval(['save ',save_str,' fluct']);
0069 %
0070 info_dke_yp(2,'Plasma fluctuations saved');

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