test_bananatip_yp

 Global parameters for the vectorial magnetic equilibrium

0001 %
0002 % Global parameters for the vectorial magnetic equilibrium
0003 %
0004 clear all
0005 close all
0006 clear functions
0007 %
0008 %load('/Users/yves/Projets/LUKE/Simulations/Examples/EAST/EAST_040235/040235_004000_ne_av=2p5/EQUIL_EAST_040235_004000.mat');
0009 %load('/Users/yves/Projets/LUKE/Simulations/Examples/EAST/EAST_040235/040235_004000_ne_av=2p5/LUKE_DATA_EAST_040235_004000_C3PO_antenna_2p45Ghz_VLOOP.mat','output');
0010 %load('LUKE_RESULTS_WEST_149293_39.9683_ES_METIS_C3PO_C4_Npar_2.1.mat')
0011 %output = lukestructs.output;
0012 %equil = lukestructs.equil;
0013 %
0014 load('/Volumes/Triolet/WEST/WEST_52702/LUKE_DATA_WEST_52702_5.0442_C3PO_C3_C3PO_C4.mat')
0015 %
0016 XXf0 = output.dke_out.XXf0{1};
0017 psin = output.radialDKE.xpsin_f;
0018 dpsin = output.radialDKE.xdpsin_f_dke;
0019 mhu = output.momentumDKE.mhu;
0020 pn = output.momentumDKE.pn;
0021 mksa = output.mksa;
0022 %
0023 pn2dpn = output.momentumDKE.pn(:).*output.momentumDKE.pn(:).*output.momentumDKE.dpn(:);
0024 %
0025 ntheta = 71;
0026 %
0027 ripple = bananatip_yp(ntheta,psin,dpsin,pn,mhu,equil,XXf0,pn2dpn,mksa);
0028 %
0029 nR = 101;
0030 nZ = 51;
0031 nlevel = 11;
0032 Ec = [mksa.Te_ref,20,100,200];%keV
0033 %
0034 R_tip = linspace(min(min(ripple.R_tip)),max(max(ripple.R_tip)),nR);
0035 Z_tip = linspace(min(min(ripple.Z_tip)),max(max(ripple.Z_tip)),nZ);
0036 %
0037 XR_tip = R_tip(:)*ones(1,nZ);
0038 XZ_tip = ones(nR,1)*Z_tip(:)';
0039 %
0040 maxDNeDphi_tip = max(max(max(ripple.DNeDphi_tip)));
0041 %
0042 for iEc = 1:length(Ec),
0043     XDNeDphi_tip = griddata(ripple.R_tip,ripple.Z_tip,squeeze(ripple.DNeDphi_tip(sum(ripple.Ec_tip < Ec(iEc)),:,:)),XR_tip,XZ_tip)';
0044     subplot(2,2,iEc),ax = graph2D_jd(XR_tip,XZ_tip,real(XDNeDphi_tip),'R (m)','Z(m)',['Number of electron banana tips @ ',num2str(ripple.Ec_tip(sum(ripple.Ec_tip < Ec(iEc)))),' (keV)'],NaN,NaN,1,NaN,nlevel);colorbar
0045 end
0046 %
0047 disp(['Max number of electron banana tips (electron/rd) : ',num2str(maxDNeDphi_tip)])