make_launch_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_launch = '';%launch local ID
0013 path_launch = '';%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 %************************************************************************************************************************************
0019 %
0020 % LAUNCH parameters
0021 %
0022 id = 'NTM_32';
0023 freq_GHz = 170;
0024 %
0025 R_L = 6.4848;
0026 Z_L = 4.110;
0027 P_L = 20;% [MW]
0028 %
0029 cone_div = 1.08;% equivalent exp(-2) power cone divergence
0030 cone_ori = -0.5;% equivalent cone origin
0031 %
0032 z_L = 0.35;% focal point
0033 %
0034 alpha_prater = -63.5;% angle between N_P and horizontal plane
0035 beta_prater = 22;% angle between N and poloidal plane
0036 %
0037 ns = 1000;
0038 method = 'cubic';
0039 %
0040 %************************************************************************************************************************************
0041 %
0042 [dkepath] = load_structures_yp('dkepath',id_dkepath,path_dkepath);
0043 %
0044 % launch structure
0045 %
0046 launch.id = [id,locid_launch];
0047 launch.type = 'EC';
0048 launch.omega_rf = [freq_GHz]*2*pi*1e9;%wave angular frequency (rad/s)
0049 %
0050 launch.yR_L = R_L;
0051 launch.yZ_L = Z_L;
0052 launch.yphi_L = 0;
0053 launch.yalpha_L = sign(beta_prater)*pi - atan(tand(beta_prater)/cosd(alpha_prater));
0054 launch.ybeta_L = acos(cosd(beta_prater)*sind(alpha_prater));
0055 launch.yP_L = P_L*1e6;
0056 launch.dNpar0 = NaN;
0057 launch.ns = ns;
0058 launch.method = method;
0059 %
0060 launch.w0 = 0.3/(freq_GHz*pi*cone_div*pi/180);
0061 launch.z_L = z_L;
0062 %
0063 %************************************************************************************************************************************
0064 %
0065 filename = ['LAUNCH_',launch.id,'.mat'];
0066 %
0067 save([path_launch,filename],'launch');
0068 %