Parameters for testing the effect of 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_JETliketest_purerippple 0002 % 0003 % Parameters for testing the effect of 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 % 0009 id_fluct = 'pureripple';%Plasma fluctuations identifier 0010 % 0011 % Path parameters 0012 % 0013 id_dkepath = '';%For all paths used by DKE solver 0014 path_dkepath = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0015 % 0016 % Equilibrium parameters 0017 % 0018 id_equil = 'JETliketest';%For plasma equilibrium 0019 path_equil = '../EQUIL/';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0020 % 0021 % Non axisymmetric perturbation of the toroidal MHD equilibrium for any rf waves : toroidal magnetic field ripple 0022 % 0023 id_ripple = 'rippletest';%Name of the machine 0024 % 0025 % ------------------------------------------------------------------------- 0026 % 0027 % Load structures 0028 % 0029 if exist('id_ripple'), 0030 if ~isempty(id_ripple), 0031 [equil,ripple] = load_structures_yp('equil',id_equil,path_equil,'ripple',id_ripple,''); 0032 fluctuations.naequilp.ripple = ripple; 0033 else 0034 [equil] = load_structures_yp('equil',id_equil,path_equil); 0035 end 0036 else 0037 [equil] = load_structures_yp('equil',id_equil,path_equil); 0038 end 0039 % 0040 [fluct] = fluctstruct_yp(equil,id_fluct,fluctuations); 0041 % 0042 % Data saving 0043 % 0044 if isempty(id_fluct), 0045 save_str = ['FLUCT_',id_equil,'.mat']; 0046 else 0047 save_str = ['FLUCT_',id_equil,'_',id_fluct,'.mat']; 0048 end 0049 eval(['save ',save_str,' fluct']); 0050 % 0051 info_dke_yp(2,'Fluctuations structure saved'); 0052 % 0053 % load('../EQUIL/EQUILV_JETliketest.mat'); 0054 % 0055 % [fluct] = fluctphase_yp(fluct); 0056 % [fluct_fit] = fitfluct_yp(fluct); 0057 % testfitfluct_yp(equil_fit,fluct,fluct_fit,0.8); 0058