proc_waves

0001 close all
0002 clear all
0003 %
0004 %Wave display procedure. Data are in /local cephee disk
0005 %
0006 if strcmp(computer,'GLNXA64'),%
0007     root_dir = '/local/peysson/Simulations/Yves/DKE/FLUCT/';%if on cephee
0008 elseif strcmp(computer,'MACI64'),
0009     root_dir = '/IRFM-LOCAL/Simulations/Yves/DKE/FLUCT/';%Macbook pro Yves
0010     root_dir = '';
0011 else
0012     root_dir = '';
0013 end
0014 %
0015 scenario = input_dke_yp('Scenario name ?','');
0016 itnmax = input_dke_yp('Number of iterations considered ?',1);
0017 %
0018 all_waves = [];
0019 for itn = 1:itnmax
0020     load([root_dir,'wave_results/waves_',scenario,'_',int2str(itn),'.mat']); 
0021     %
0022     for iw = 1:length(waves)
0023         for jr = 1:length(waves{iw}.rays),
0024             waves{iw}.rays{jr}.P0_2piRp = waves{iw}.rays{jr}.P0_2piRp/itnmax;
0025             Npar0(itn,iw,jr) = real(waves{iw}.rays{jr}.sNpar(1));
0026         end
0027     end
0028     %
0029     all_waves = [all_waves,waves];
0030 end
0031 %
0032 for jr = 1:length(waves{1}.rays),
0033     figure,hist(Npar0(:,iw,jr))
0034 end
0035 %
0036 p_opt = input_dke_yp('Option for output figures : -1 (nothing), 0 (print), 1 (print and save), 2 (save)',-1,-1:2,'',[1,1]);
0037 g_opt.abs_max = input_dke_yp('Option for ray trajectories : 0 (full trajectory), 1 (cut when damped)',0,0:1,'',[1,1]);
0038 %
0039 g_opt.mode = input_dke_yp('Option for showing ray polarization (n/y)','n',{'n','y'});
0040 %xvaropt = input_dke_yp('Option for rho definition : 0 (geometric), 1 (poloidal flux), 2 (toroidal flux)',0,0:2,'',[1,1]);
0041 xvar = input_dke_yp('Option for rho definition : ''g'' (geometric), ''p'' (poloidal flux), ''t'' (toroidal flux), ''v'' (volumic rho)','g',{'g','p','t','v'},'',[1,1]);
0042 %
0043 g_opt.print = p_opt;
0044 g_opt.TeNpar = 1;
0045 nr_dep = 40;
0046 rho_disp = 0:0.1:1;
0047 ymask = NaN;
0048 %
0049 graph_wave_jd(wavestructs{1}.equil,all_waves,nr_dep,g_opt,rho_disp,ymask,xvar);
0050 %