This function returns LH wave launching parameters by J. Decker and J.-F. Artaud, 02/02/2011
0001 function launchs = make_wavelaunchs_metis2luke_lh(option,lhparam) 0002 % 0003 % This function returns LH wave launching parameters 0004 % 0005 % by J. Decker and J.-F. Artaud, 02/02/2011 0006 % 0007 if lhparam.plh < lhparam.plhmin,% not enough power to account for wave 0008 launchs = cell(1,0); 0009 return 0010 end 0011 % 0012 if lhparam.mode == 0, 0013 % 0014 % simple LUKE case with v1, v2, D, to match prescription or HXR profile 0015 % 0016 error('lhparam.mode = 0 not yet implemented') 0017 % 0018 elseif lhparam.mode == 1, 0019 % 0020 % C3PO/LUKE case with one antenna and midplane vertical launching position 0021 % 0022 launch.id = 'LH'; 0023 launch.type = 'LH'; 0024 launch.omega_rf = lhparam.freqlh*2*pi*1e9;%wave angular frequency (rad/s) [1,1] 0025 % 0026 if lhparam.lob,%negative lobe accounted for 0027 launch.bNpar0 = option.signe*[lhparam.npar0,-6];%wave parallel index of refraction [1,n_lobe] 0028 launch.bPlhtot = lhparam.plh*[lhparam.directivity,1 - lhparam.directivity];%wave initial power [1,n_lobe] 0029 else 0030 launch.bNpar0 = [lhparam.npar0]; 0031 launch.bPlhtot = lhparam.plh*lhparam.directivity; 0032 end 0033 % 0034 launch.bNpar0 = -option.signe*launch.bNpar0;% for co-current 0035 % 0036 launch.bdNpar0 = lhparam.dnpar0*ones(size(launch.bNpar0));%ray spectral width for Fokker-Planck calculations [1,n_lobe] 0037 % 0038 launch.rZ0 = lhparam.z0;%vertical launching positions (m) [1,m_pol] 0039 % 0040 launch.m0 = lhparam.m0;%initial poloidal mode number [1,1] 0041 launch.phi0 = lhparam.phi0;%toroidal angle position [1,1] 0042 launch.i_ref = lhparam.i_ref;%Initial radial position index [1,1] 0043 launch.LFS = lhparam.LFS;% (-1) high field side launching, (1) low field side launching [1,1] 0044 launch.mmode = lhparam.mmode;% (-1) for slow wave, (1) for fast wave [1,1] 0045 launch.nparp = lhparam.nparp;% spectral broadening factor [1,1] 0046 launch.a_sdNpar = lhparam.a_sdNpar; % initial spectral upshift [1,1] 0047 % 0048 launchs{1} = launch; 0049 % 0050 elseif lhparam.mode == 2, 0051 % 0052 % C3PO/LUKE case with detailed antenna structure and/or ALOHA 0053 % 0054 error('lhparam.mode = 2 not yet implemented') 0055 % 0056 % gene_lh.signe.b0 = 1; 0057 % gene_lh.signe.ip = option.signe; 0058 % % 0059 % par_lh.freqlh = lhparam.freqlh; 0060 % par_lh.aloha_mode = lhparam.aloha_mode; 0061 % par_lh.lob = lhparam.lob; 0062 % par_lh.conf = lhparam.conf; 0063 % par_lh.a_sdNpar = lhparam.a_sdNpar; 0064 % par_lh.nparp = lhparam.nparp; 0065 % % 0066 % cons_lh = lhparam.plh*exp(1i*lhparam.npar0); 0067 % % 0068 % launchs = make_wavelaunchs_cronos_lh(gene_lh,par_lh,cons_lh);%for the LH waves 0069 % 0070 end 0071 % 0072