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>
0001 function [] = make_fluct_JETh77 0002 % 0003 % Parameters for testing the effect of density or magnetic field fluctuations on ray trajectories 0004 % 0005 % by Y. Peysson (DRFC/DSM/CEA) <yves.peysson@cea.fr> and J. Decker (DRFC/DSM/CEA) <joan.decker@cea.fr> 0006 % 0007 close all 0008 clear fluct 0009 % 0010 id_fluct = 'test';%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 = 'JETh77';%For plasma equilibrium 0020 path_equil = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0021 % 0022 % Density fluctuations (kpar = 0 fluctuation spectrum model) 0023 % 0024 struct_ne.type = 'ne';%Type of fluctuations or modulations [1,1] 0025 struct_ne.model = [1,1];%Statistical ne fluctuation model (>= 1) : (1) -> Gaussian, (2) -> ... [1,nfluct_types] 0026 struct_ne.sigmar_max = [0.3,0.1];%Maximum value of the relative ne fluctuations variance at the poloidal angle theta = 0 [1,nfluct_types] 0027 struct_ne.sigmar_hwhm = [0.2,0.3];%Radial half width at half maximum of the relative ne fluctuations variance at the poloidal angle theta = 0 [1,nfluct_types] 0028 struct_ne.sigmar_rho = [1,1];%Radial location where the relative ne fluctuations variance peaks at the poloidal angle theta = 0 [1,nfluct_types] 0029 struct_ne.polmode = [0.2,1];%HFS/LFS relative ne fluctuations variance ratio. No poloidal dependency corresponds to 1 [1,nfluct_types] 0030 struct_ne.kperpc = [300,100];%Inverse of the correlation length of ne fluctuations perpendicular to the magnetic field (ap^(-1) unit) [1,nfluct_types] 0031 struct_ne.kperpref = [0,0];%ne fluctuations reference perpendicular wave vector with respect to the local magnetic field direction (ap^(-1) unit) [1,nfluct_types] 0032 struct_ne.raperp = [0.01,0.01];%Minimum relative amplitude considered in the Fourier series expansion [1,nfluct_types] 0033 struct_ne.nkperp = [200,600];%Number of points for describing the kperp power spectrum (odd number) [1,nfluct_types] 0034 % 0035 % Magnetic field fluctuations (kpar = 0 fluctuation spectrum model) 0036 % 0037 struct_B.type = 'Bperp';%Type of fluctuations or modulations [1,1] 0038 struct_B.model = [1,1];%Statistical Bperp fluctuation model (>= 1) : (1) -> Gaussian, (2) -> ... [1,nfluct_types] 0039 struct_B.sigmar_max = [5e-3,1e-4];%Maximum value of the relative Bperp fluctuations variance at the poloidal angle theta = 0 [1,nfluct_types] 0040 struct_B.sigmar_hwhm = [0.3,0.2];%Radial half width at half maximum of the relative Bperp fluctuations variance at the poloidal angle theta = 0 [1,nfluct_types] 0041 struct_B.sigmar_rho = [0,1];%Radial location where the relative Bperp fluctuations variance peaks at the poloidal angle theta = 0 [1,nfluct_types] 0042 struct_B.polmode = [1,0.3];%HFS/LFS relative Bperp fluctuations variance ratio. No poloidal dependency corresponds to 1 [1,nfluct_types] 0043 struct_B.kperpc = [300,200];%Inverse of the correlation length of Bperp fluctuations perpendicular to the magnetic field (ap^(-1) unit) [1,nfluct_types] 0044 struct_B.kperpref = [0,0];%Bperp fluctuations reference perpendicular wave vector with respect to the local magnetic field direction (ap^(-1) unit) [1,nfluct_types] 0045 struct_B.raperp = [0.01,0.01];%Minimum relative amplitude considered in the Fourier series expansion [1,nfluct_types] 0046 struct_B.nkperp = [300,200];%Number of points for describing the kperp power spectrum (odd number) [1,nfluct_types] 0047 % 0048 % ------------------------------------------------------------------------- 0049 % 0050 % Load structures 0051 % 0052 [equil] = load_structures_yp('equil',id_equil,path_equil); 0053 % 0054 % ========================================================================= 0055 % 0056 [fluct] = fluctstruct_yp(equil,id_fluct,struct_ne,struct_B); 0057 % 0058 % Data saving 0059 % 0060 save_str = ['FLUCT_',id_equil,'_',id_fluct,'.mat']; 0061 eval(['save ',save_str,' fluct']); 0062 % 0063 info_dke_yp(2,'Plasma fluctuations saved');