script make_wave_TScirc_LHtest Parameters for test mode LHCD calculations This function has a benchmarking purpose only Approximations: - Circular equilibrium required - The polarization is purely electrostatic: in that case, only the parallel polarization term is non-zero. - Small FLR limit: J0(z) -> 1. by J. Decker (RLE/MIT) <jodecker@mit.edu> and Y. Peysson (DRFC/DSM/CEA) <yves.peysson@cea.fr>
0001 % script make_wave_TScirc_LHtest 0002 % 0003 % Parameters for test mode LHCD calculations 0004 % This function has a benchmarking purpose only 0005 % Approximations: 0006 % - Circular equilibrium required 0007 % - The polarization is purely electrostatic: in that case, only the parallel polarization term is non-zero. 0008 % - Small FLR limit: J0(z) -> 1. 0009 % 0010 % by J. Decker (RLE/MIT) <jodecker@mit.edu> and Y. Peysson (DRFC/DSM/CEA) <yves.peysson@cea.fr> 0011 % 0012 clear all 0013 % 0014 id_wave = 'LH_karney_offaxis_2'; 0015 % 0016 id_equil = 'TScirc';%For plasma equilibrium 0017 path_equil = '../EQUIL_files/';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0018 % 0019 [equil] = load_structures_yp('equil',id_equil,path_equil); 0020 % 0021 %************************************************************************************************************************************ 0022 % 0023 wavestruct.omega_lh = [4]*2*pi*1e9; %(GHz -> rad/s). Wave frequency [1,1] Indicative, no effect in small FLR limit opt_lh > 0 0024 %Option parameter for cross-comparison between old LH code: 0025 % - (1): 1/vpar dependence 0026 % - (2): no 1/vpar dependence and old grid technique for Dlh calculations (Karney, Shoucri) (see rfdiff_dke_jd) 0027 wavestruct.opt_lh = 2; % [1,1] 0028 % 0029 % Choose (vparmin_lh,vparmax_lh) or (Nparmin_lh,Nparmax_lh) for square n// LH wave power spectrum, 0030 % or (Npar_lh,dNpar_lh) for Gaussian shape 0031 % 0032 wavestruct.norm_ref = 1;%Normalization procedure for the LH quasilinear diffusion coefficient and spectrum boundaries 0033 % 0034 wavestruct.yvparmin_lh = [4];%LH wave square N// Spectrum: Lower limit of the plateau (vth_ref or vth) [1,n_scenario_lh] 0035 wavestruct.yvparmax_lh = [7];%LH wave square N// Spectrum: Upper limit of the plateau (vth_ref or vth) [1,n_scenario_lh] 0036 % 0037 wavestruct.yNparmin_lh = [NaN];%LH wave square N// Spectrum: Lower limit [1,n_scenario_lh] 0038 wavestruct.yNparmax_lh = [NaN];%LH wave square N// Spectrum: Upper limit [1,n_scenario_lh] 0039 wavestruct.yNpar_lh = [NaN];%LH wave Gaussian N// Spectrum: peak [1,n_scenario_lh] 0040 wavestruct.ydNpar_lh = [NaN];%LH wave Gaussian N// Spectrum: width [1,n_scenario_lh] 0041 % 0042 % Note: this diffusion coefficient is different from the general QL D0. It has a benchmarking purpose only 0043 wavestruct.yD0_in_c_lh = [2];%Central LH QL diffusion coefficient (nhuth_ref*pth_ref^2 or nhuth*pth^2) [1,n_scenario_lh] 0044 % 0045 wavestruct.yD0_in_lh_prof = [1];%Quasilinear diffusion coefficient radial profile: (0) uniform, (1) gaussian radial profile [1,n_scenario_lh] 0046 wavestruct.ypeak_lh = [0.4];%Radial peak position of the LH quasi-linear diffusion coefficient (r/a on midplane) [1,n_scenario_lh] 0047 wavestruct.ywidth_lh = [0.1];%Radial width of the LH quasi-linear diffusion coefficient (r/a on midplane) [1,n_scenario_lh] 0048 % 0049 wavestruct.ythetab_lh = [0]*pi/180;%(deg -> rad). Poloidal location of LH beam [0..2pi] [1,n_scenario_lh] 0050 % (0) from local values Te and ne, (1) from central values Te0 and ne0 0051 % 0052 %************************************************************************************************************************************ 0053 % 0054 wave = make_idealLHwave_jd(equil,wavestruct); 0055 % 0056 save_str = ['WAVE_',id_equil,'_',id_wave,'.mat']; 0057 eval(['save ',save_str,' wave']);