reflection_jd

 This function calculates the reflection on vessel

0001 function [P_ref,N_ref,E_ref] = reflection_jd(equil,P_end,N_end,E_end)
0002 %
0003 % This function calculates the reflection on vessel
0004 %
0005 if isfield(equil,'tokamak') && isfield(equil,'shotnum') && isfield(equil,'shotime'),
0006     tokamak = equil.tokamak;
0007     shotnum = equil.shotnum;
0008     shotime = equil.shotime;
0009 else
0010     istru = find(equil.id == '_');
0011     if length(istru) >= 2,
0012         tokamak = equil.id(1:istru(1)-1);
0013         shotnum = equil.id(istru(1)+1:istru(2)-1);
0014         shotime = equil.id(istru(2)+1:end);
0015     else           
0016         error('-----> Shot number information missing');
0017     end
0018 end
0019 %
0020 if ~ischar(shotime),
0021     shotime = num2str(shotime);
0022 end
0023 %
0024 if strcmp(tokamak,'TCV'),
0025     %
0026     addpath('/home/gnesin')
0027     %
0028     %save('rayend.mat','shotnum','shotime','P_end','K_end','E_end');
0029     %
0030     [P_ref,N_ref,E_ref] = reflect_tcv(shotnum,shotime,P_end,N_end,E_end);
0031     %
0032 end
0033