0001 function output = iluke_output(output)
0002
0003 if ~isfield(output,'data_proc'),
0004
0005 output.data_proc = proc_luke_jd(output,struct);
0006
0007 end
0008
0009 flag_rf = output.data_proc.scalar.p0_rf_2piRp ~= 0;
0010
0011 output.info.txt = {...
0012 ['Total current : ',num2str(output.data_proc.scalar.I_tot),' MA (positive for co-current)'],...
0013 ['Current peak r/a : ',num2str(output.data_proc.scalar.jrhoT),' (square root of toroidal flux)'],...
0014 ['Current FWHM : ',num2str(output.data_proc.scalar.jdrhoT),' (square root of toroidal flux)']};
0015 if flag_rf,
0016 output.info.txt = [output.info.txt,...
0017 ['RF conv. param. : ',num2str(output.data_proc.simul.residue_rf(end))]];
0018 end
0019
0020 if isinf(output.data_proc.scalar.Rp),
0021 if flag_rf,
0022 output.info.txt = [output.info.txt,...
0023 ['Launched RF power : ',num2str(output.data_proc.scalar.p0_rf_2piRp),' MW/m (P/2piRp)'],...
0024 ['Absorbed RF power : ',num2str(output.data_proc.scalar.p_rf_2piRp),' MW/m ---> ',num2str(100*output.data_proc.scalar.p_rf_2piRp/output.data_proc.scalar.p0_rf_2piRp),' %. (P/2piRp)']];
0025 for in_rf = 1:length(output.data_proc.scalar.n_rf_list),
0026 output.info.txt = [output.info.txt,...
0027 ['RF power on harm. ',num2str(output.data_proc.scalar.n_rf_list(in_rf)),' : ',num2str(output.data_proc.scalar.np_rf_2piRp(in_rf)),' MW/m ---> ',num2str(100*output.data_proc.scalar.np_rf_2piRp(in_rf)/output.data_proc.scalar.p0_rf_2piRp),' %. (P/2piRp)']];
0028 end
0029 output.info.txt = [output.info.txt,...
0030 ['Linear RF power : ',num2str(output.data_proc.scalar.p_rf_2piRp_lin),' MW/m ---> ',num2str(100*output.data_proc.scalar.p_rf_2piRp_lin/output.data_proc.scalar.p0_rf_2piRp),' %. (P/2piRp)']];
0031 end
0032
0033 output.info.txt = [output.info.txt,...
0034 ['Ohmic power : ',num2str(output.data_proc.scalar.p_e_2piRp),' MW/m (P/2piRp)'],...
0035 ['Total power : ',num2str(output.data_proc.scalar.p_tot_2piRp),' MW/m (P/2piRp)'],...
0036 ['CD efficiency : ',num2str(output.data_proc.scalar.EtaRp),' A.m/W (I*Rp/P)']];
0037 else
0038 if flag_rf,
0039 output.info.txt = [output.info.txt,...
0040 ['Launched RF power : ',num2str(output.data_proc.scalar.p0_rf_2piRp*(2*pi*output.data_proc.scalar.Rp)),' MW'],...
0041 ['Absorbed RF power : ',num2str(output.data_proc.scalar.p_rf_2piRp*(2*pi*output.data_proc.scalar.Rp)),' MW ---> ',num2str(100*output.data_proc.scalar.p_rf_2piRp/output.data_proc.scalar.p0_rf_2piRp),' %.']];
0042 for in_rf = 1:length(output.data_proc.scalar.n_rf_list),
0043 output.info.txt = [output.info.txt,...
0044 ['RF power on harm. ',num2str(output.data_proc.scalar.n_rf_list(in_rf)),' : ',num2str(output.data_proc.scalar.np_rf_2piRp(in_rf)*(2*pi*output.data_proc.scalar.Rp)),' MW ---> ',num2str(100*output.data_proc.scalar.np_rf_2piRp(in_rf)/output.data_proc.scalar.p0_rf_2piRp),' %.']];
0045 end
0046 output.info.txt = [output.info.txt,...
0047 ['Linear RF power : ',num2str(output.data_proc.scalar.p_rf_2piRp_lin*(2*pi*output.data_proc.scalar.Rp)),' MW ---> ',num2str(100*output.data_proc.scalar.p_rf_2piRp_lin/output.data_proc.scalar.p0_rf_2piRp),' %.']];
0048 end
0049
0050 output.info.txt = [output.info.txt,...
0051 ['Ohmic power : ',num2str(output.data_proc.scalar.p_e_2piRp*(2*pi*output.data_proc.scalar.Rp)),' MW'],...
0052 ['Total power : ',num2str(output.data_proc.scalar.p_tot_2piRp*(2*pi*output.data_proc.scalar.Rp)),' MW'],...
0053 ['CD efficiency : ',num2str(output.data_proc.scalar.EtaRp/output.data_proc.scalar.Rp),' A/W (I/P)']];
0054 end
0055
0056 if flag_rf,
0057 output.info.txt = [output.info.txt,...
0058 ['RF efficiency : ',num2str(output.data_proc.scalar.EtaRFRp),' 10^19 A/(W.m^2) (n_19*I*Rp/PRF)']];
0059 end