make_wavestruct_ITER_Scen2_C3PO_EC

 script make_wavestruct_Scen2_C3PO_EC

0001 % script make_wavestruct_Scen2_C3PO_EC
0002 %
0003 clear all
0004 clear mex
0005 clear functions
0006 close all
0007 warning('off','all')
0008 pause(1)
0009 %
0010 % *********************** Specify LUKE structures *****************************
0011 %
0012 locid_wavestruct = '_24rays';%wavestruct local ID
0013 path_wavestruct = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0014 %
0015 id_dkepath = '';%For all paths used by DKE solver
0016 path_dkepath = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0017 %
0018 id_equil = 'ITER_Scen2_200103121816_430_129';%For plasma equilibrium
0019 path_equil = 'EQUIL/';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0020 %
0021 id_rtparam = 'EC';%For C3PO code parameters
0022 path_rtparam = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0023 %
0024 ids_launch = {'NTM_32'};%For plasma equilibrium
0025 paths_launch = {''};%if nothing is specified, the working directory is first used and then MatLab is looking in all the path
0026 %
0027 %************************************************************************************************************************************
0028 %
0029 % C3PO computing parameters
0030 %
0031 C3POparam = struct;
0032 %
0033 % Display parameters for C3PO
0034 %
0035 C3POdisplay = struct;
0036 %
0037 % Wave parameters for C3PO
0038 %
0039 waveparam.mmode = -1;% cold plasma mode (1) m (-1) p, (slow/fast LH wave: -1/+1, O/X mode: -1/+1)
0040 waveparam.nd = 4;%Number of transverse distance positions for beamlets
0041 waveparam.nchi = 6;%Number of angular positions for beamlets
0042 %
0043 % Equilibrium parameters C3PO
0044 %
0045 fitparam.nharm = 32;%default value has since been changed to 64
0046 %
0047 % Global parameters for C3PO
0048 %
0049 rayparam = struct;
0050 %
0051 %************************************************************************************************************************************
0052 %
0053 [dkepath,equil,rtparam,launchs] = load_structures_yp('dkepath',id_dkepath,path_dkepath,'equil',id_equil,path_equil,'rtparam',id_rtparam,path_rtparam,'launchs',ids_launch,paths_launch);
0054 %
0055 fitparam = conc_struct_jd(rtparam.fitparam,fitparam);
0056 rayparam = conc_struct_jd(rtparam.rayparam,rayparam);
0057 C3POparam = conc_struct_jd(rtparam.C3POparam,C3POparam);
0058 C3POdisplay = conc_struct_jd(rtparam.C3POdisplay,C3POdisplay);
0059 waveparam = conc_struct_jd(rtparam.waveparam,waveparam);
0060 %
0061 % wavestruct
0062 %
0063 wavestruct.wavesolver = 'C3PO';
0064 wavestruct.launch = launchs{1};
0065 wavestruct.id = [wavestruct.wavesolver,'_',launchs{1}.id,locid_wavestruct];
0066 wavestruct.waveparam = waveparam;
0067 wavestruct.rayparam = rayparam;
0068 wavestruct.fitequilparam = fitparam;
0069 wavestruct.raydisplay = C3POdisplay;
0070 wavestruct.C3POparam = C3POparam;
0071 %
0072 %************************************************************************************************************************************
0073 %
0074 filename = ['WAVESTRUCT_',equil.id,'_',wavestruct.id,'.mat'];
0075 %
0076 save([path_wavestruct,filename],'wavestruct');
0077 %