runhxr

0001 %
0002 permission = test_permissions_yp;
0003 %
0004 if ~permission 
0005     disp('Please move the script to a local folder where you have write permission before to run it')
0006     return;
0007 end
0008 %
0009 id_simul = 'LH_karney_ripple';%Simulation ID
0010 path_simul = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0011 %
0012 id_equil = 'TScirc';%For plasma equilibrium
0013 path_equil = [dkepath.luke_root,'Project_DKE/Database/EQUIL_files/'];%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0014 %
0015 id_hxr = 'TS99';%For HXR diagnostic configuration
0016 path_hxr = [dkepath.luke_root,'Project_DKE/Database/HXR_files/'];%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0017 %
0018 %************************************************************************************************************************************
0019 %************************************************************************************************************************************
0020 %************************************************************************************************************************************
0021 %
0022 if exist([path_hxr,'HXR_',id_hxr,'.mat']) == 0,%HXR
0023     info_dke_yp(2,[path_hxr,'HXR_',id_hxr,'.mat does not exist']);
0024     break
0025 else
0026     eval(['load ',path_hxr,'HXR_',id_hxr,'.mat']);
0027     info_dke_yp(2,'Structures ''hxr and hxrparam'' successfully loaded');
0028 end
0029 %
0030 if exist([path_simul,'DKE_RESULTS_',id_equil,'_',id_simul,'.mat']) == 0,%DKE results
0031     info_dke_yp(2,[path_simul,'DKE_RESULTS_',id_equil,'_',id_simul,'.mat does not exist']);
0032     break
0033 else
0034     eval(['load ',path_simul,'DKE_RESULTS_',id_equil,'_',id_simul,'.mat']);
0035     info_dke_yp(2,['DKE data results in file ',path_simul,'DKE_RESULTS_',id_equil,'_',id_simul,'.mat successfully loaded']);
0036 end
0037 %
0038 [Zbremchord,equilHXR] = bremchord_dke_yp(dkeparam,dkedisplay,equil,radialDKE,Zcurr,hxr,hxrparam);                        
0039 [Zbremplasma] = bremsstrahlung_dke_yp(dkeparam,display,equilHXR,radialDKE,mksa,momentumDKE,dke_out,Zbremchord,hxr,hxrparam);
0040 [Zbremdiag] = bremdiag_dke_yp(Zbremplasma,hxr,hxrparam);
0041 %
0042 %For display HXR
0043 %
0044 bremchorddisplay_dke_yp(Zbremchord,1);
0045 bremlocaldisplay_dke_yp(dkeparam,display,equilHXR,momentumDKE,radialDKE,dke_out,mksa,Zbremchord,Zbremplasma,0.3,180);
0046 bremdiagdisplay_dke_yp(hxr,Zbremchord,Zbremplasma,Zbremdiag,42,50);
0047 %
0048 eval(['save ',path_simul,'HXR_RESULTS_',id_equil,'_',id_simul,'_',id_hxr,'.mat Zbremchord Zbremplasma Zbremdiag equilHXR']);
0049 info_dke_yp(2,['Data saved in ',path_simul,'DKE_RESULTS_',id_equil,'_',id_simul,'_',id_hxr,'.mat']);