Commit c416208c authored by Lore, Jeremy's avatar Lore, Jeremy
Browse files

Update adas routines

parent bd532885

matlab/bfield_library_jdl/ADAS_routines/get_adas_adf11_data.m

+2 −1
Original line number Diff line number Diff line
@@ -5,7 +5,8 @@ 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 = {'c','h','he','li','be','ar','o','ne','n','w','fe','xe'};

species = {'he'};

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



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

matlab/bfield_library_jdl/ADAS_routines/interp_adas_rate_coefficient.m

+8 −1
Original line number Diff line number Diff line 
function sigma_v = interp_adas_rate_coefficient(Te,ne,te_adas,ne_adas,coeff_adas)

% Te in eV, ne in cm^-3

% coefficient cm^3/s  

% coefficient cm^3/s   (Wcm^3 for plt)



nTeAdas = length(te_adas);

nNeAdas = length(ne_adas);



if ~all(size(coeff_adas) == [nNeAdas,nTeAdas])

    error('sizes do not match!')

end



log_interp = 1;

% my_method = 'spline';

matlab/bfield_library_jdl/ADAS_routines/plot_Lz_adas.m

0 → 100644
+53 −0
Original line number Diff line number Diff line 
clearvars;





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

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

fname = ['C:\Work\ADAS\adf11_all\plt',suffix,'\','plt',suffix,'_',element,'.dat'];  % Radiated power (W cm^3)

plt = read_adas_adf11_file(fname);

fname = ['C:\Work\ADAS\adf11_all\prb',suffix,'\','prb',suffix,'_',element,'.dat'];  % Radiated power (W cm^3)

prb = read_adas_adf11_file(fname);







% Te_test = linspace(0.1,5,100);

Te_test = logspace(-1,4,100);  % eV

ne_test = 1e13;   % cm^-3





nCharge = plt.charge;

for iCharge = 1:nCharge

    for i=1:length(ne_test)

        plt_eval(:,iCharge) = interp_adas_rate_coefficient(Te_test,ne_test,plt.te,plt.ne,plt.coeff(:,:,iCharge));

        prb_eval(:,iCharge) = interp_adas_rate_coefficient(Te_test,ne_test,prb.te,prb.ne,prb.coeff(:,:,iCharge));

    end

end





figure; hold on; box on; grid on;

plot(Te_test,plt_eval/1e6)

% plot(Te_test,sum(plt_eval,2)/1e6,'k')

plot(Te_test,prb_eval/1e6)

% plot(Te_test,(prb_eval+prb_eval)/1e6)

% plot(Te_test,sum(prb_eval,2)/1e6,'k')

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

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

ylabel('<\sigmav> (Wm^3)','fontsize',12)

legend('0','1','2','3','4','5','tot')

axis([0.1,1e4,1e-35,1e-30])



% asdfasd

% s = styflipper(length(ne_test));

% for i=1:length(ne_test)

%     plot(Te_test,sv_iz(:,i),'r','linewidth',2,'linestyle',char(s{i}))

%     plot(Te_test,sv_rc(:,i),'b','linewidth',2,'linestyle',char(s{i}))

% end

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

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

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

% ylabel('<\sigmav> (cm^3/s)','fontsize',12)

% set(gca,'fontsize',12)

% % axis([1,2e4,10^-12,10^-7])

% % axis([0,5,10^-16,10^-8])

% axis([1,2e4,10^-16,10^-7])

matlab/bfield_library_jdl/ADAS_routines/plot_iz_coeff_interp.m

+5 −3
Original line number Diff line number Diff line
@@ -3,7 +3,9 @@ clearvars; 
% suffix = '12'; element='h';

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

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));

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
@@ -15,11 +17,11 @@ end 


s = styflipper(length(ne_test));



% figure; hold on; box on

figure; hold on; box on

for i=1:length(ne_test)

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

end

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

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

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

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

ylabel('<\sigmav> (m^3/s)','fontsize',12)

matlab/bfield_library_jdl/ADAS_routines/test_adas_interp.m

+9 −9
Original line number Diff line number Diff line 
% test_adas_interp

clearvars;



% This is not general for non-Hydrogen cases, must then choose/squeeze coefficients!



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

[te_scd,ne_scd,scd] = read_adas_adf11_file(fname);

scd = read_adas_adf11_file(fname);

fname = ['C:\Work\ADAS\adf11_all\acd96\','acd96_h.dat'];  % Effective recombination coefficients (cm^-3/s)

[te_acd,ne_acd,acd] = read_adas_adf11_file(fname);

acd = read_adas_adf11_file(fname);

fname = ['C:\Work\ADAS\adf11_all\ccd96\','ccd96_h.dat'];  % Effective cx coefficients (cm^-3/s)

[te_ccd,ne_ccd,ccd] = read_adas_adf11_file(fname);

ccd = read_adas_adf11_file(fname);





% Te_test = linspace(0.1,5,100);

Te_test = logspace(-1,4,100);

ne_test = logspace(13,15,3);



scd1 = squeeze(scd(1,:,:));

acd1 = squeeze(acd(1,:,:));

ccd1 = squeeze(ccd(1,:,:));



for i=1:length(ne_test)

    sv_iz(:,i) = interp_adas_rate_coefficient(Te_test,ne_test(i),te_scd,ne_scd,scd1);

    sv_rc(:,i) = interp_adas_rate_coefficient(Te_test,ne_test(i),te_acd,ne_acd,acd1);

    sv_cx(:,i) = interp_adas_rate_coefficient(Te_test,ne_test(i),te_ccd,ne_ccd,ccd1);

    sv_iz(:,i) = interp_adas_rate_coefficient(Te_test,ne_test(i),scd.te,scd.ne,scd.coeff);

    sv_rc(:,i) = interp_adas_rate_coefficient(Te_test,ne_test(i),acd.te,acd.ne,acd.coeff);

    sv_cx(:,i) = interp_adas_rate_coefficient(Te_test,ne_test(i),ccd.te,ccd.ne,ccd.coeff);

end



s = styflipper(length(ne_test));