Commit 489cf838 authored by Lore, Jeremy's avatar Lore, Jeremy
Browse files

Update adas plotting routines

parent c416208c

matlab/bfield_library_jdl/ADAS_routines/get_adas_adf11_data.m

+1 −1
Original line number Diff line number Diff line
@@ -6,7 +6,7 @@ url_base = 'http://open.adas.ac.uk/download/adf11/'; 


types = {'acd','plt','scd','ccd','prb','prc'};

% species = {'c','h','he','li','be','ar','o','ne','n','w','fe','xe'};

species = {'he'};

species = {'ne'};

years = {'89','93','96','12'};



outname_base = 'C:\Work\ADAS\adf11_all';

matlab/bfield_library_jdl/ADAS_routines/plot_adas_data_recomb_ioniz_hydrogen.m

+2 −1
Original line number Diff line number Diff line 
clearvars;



suffix = '96'; element='h';

% suffix = '96'; element='h';

% suffix = '96'; element='c';

% suffix = '50'; element='w';

suffix = '96'; element='ne';



fprintf('Reading scd\n');

fname = ['C:\Work\ADAS\adf11_all\scd',suffix,'\','scd',suffix,'_',element,'.dat'];  % Effective ionization coefficients (cm^-3/s)

matlab/bfield_library_jdl/ADAS_routines/plot_iz_coeff_interp.m

+29 −5
Original line number Diff line number Diff line 
clearvars;



% suffix = '12'; element='h';

suffix = '96'; element='he';

% suffix = '96'; element='he';

suffix = '96'; element='ne';

fname = ['C:\Work\ADAS\adf11_all\scd',suffix,'\','scd',suffix,'_',element,'.dat'];  % Effective ionization coefficients (cm^-3/s)

scd = read_adas_adf11_file(fname); 



scd1 = squeeze(scd.coeff(:,:,1));  % Choose charge 



Te_test = linspace(.1,20,100);

ne_test = logspace(log10(1e11),log10(1e13),3);  % cm3

% Te_test = linspace(.1,20,100);

Te_test = logspace(0,3,1000);

% ne_test = logspace(log10(1e11),log10(1e13),3);  % cm3



ne_test = 1e14;



for i=1:length(ne_test)

    sv_iz_adas(:,i) = interp_adas_rate_coefficient(Te_test,ne_test(i),scd.te,scd.ne,scd1);
@@ -17,7 +21,8 @@ end 


s = styflipper(length(ne_test));



figure; hold on; box on

figure; hold on; box on; grid on



for i=1:length(ne_test)

    plot(Te_test,sv_iz_adas(:,i)./1e6,'b','linewidth',2,'linestyle',char(s{i}))

end
@@ -27,4 +32,23 @@ xlabel('T_e (eV)','fontsize',12) 
ylabel('<\sigmav> (m^3/s)','fontsize',12)

set(gca,'fontsize',12)

axis([.1,10,10^-20,10^-14])

legend('n_e = 1e17','n_e = 1e18','n_e = 1e19')
 
 No newline at end of file
 
% legend('n_e = 1e17','n_e = 1e18','n_e = 1e19')'





pc = phys_const;

sigma_v = sv_iz_adas./1e6;

% vp = sqrt(2.*0.026*pc.eV/(20*pc.mp));

% vp = sqrt(2.*0.1*pc.eV/(20*pc.mp));

vp = sqrt(2.*100*pc.eV/(20*pc.mp));

data_plot = vp./(ne_test*1e6.*sigma_v);



figure; hold on; box on; grid on

plot(Te_test,data_plot)

set(gca,'xscale','log')

set(gca,'yscale','log')

xlabel('T_e (eV)','fontsize',12)

ylabel('\lambda_{mfp} (m)','fontsize',12)

set(gca,'fontsize',12)



set(gcf,'color','w')