Script for running LUKE This script is automatically generated by irunluke_jd. by J. Decker (joan.decker@cea.fr) and Y. Peysson (yves.peysson@cea.fr)
0001 % Script for running LUKE 0002 % 0003 % This script is automatically generated by irunluke_jd. 0004 % by J. Decker (joan.decker@cea.fr) and Y. Peysson (yves.peysson@cea.fr) 0005 % 0006 clear all 0007 clear mex 0008 clear functions 0009 close all 0010 warning('off','all') 0011 pause(1) 0012 % 0013 sigmar_max = [0.7];% Maximum value of the relative fluctuations variance at the poloidal angle theta = 0 [1,nfluct_types] 0014 sigmar_hwhm = [0.02];% Radial half width at half maximum of the relative fluctuations variance at the poloidal angle theta = 0 [1,nfluct_types] 0015 Nfluct = [1];% Ratio of the LUKE evolution time (phase refresh) to the fluctuation time - also equal to the number of duplicated rays 0016 % 0017 tn = [100i];% Final simulation time (normalized to the central thermal collision time) 0018 dtn = [1i];% internal LUKE time step. As an imaginary number, dkeparam.dtn prescribes the number of internal time steps 0019 % 0020 p_opt = -1; 0021 opt_save = -2;%In case of crash, to restart simulation from an intermediate step stored in a *.mat file (<0 for launching simulation from backup, 0 no backup, >0 for performing backup) 0022 % 0023 % *********************** Specify LUKE structures ***************************** 0024 % 0025 locid_simul = '';%Simulation ID 0026 path_simul = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0027 % 0028 psin_S = [];%Normalized poloidal flux grid where calculations are performed (0 < psin_S < 1) (If one value: local calculation only, not used if empty) 0029 rho_S = [0:0.1:0.6,0.61:0.01:0.79,0.80:0.05:1]; 0030 % 0031 id_dkepath = '';%For all paths used by DKE solver 0032 path_dkepath = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0033 % 0034 id_equil = 'ITER_Scen2_200103121816_430_129';%For plasma equilibrium 0035 path_equil = 'EQUIL/';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0036 % 0037 id_fluct = 'fluctn';%For fast electron radial transport 0038 path_fluct = 'FLUCT/';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0039 % 0040 id_dkeparam = 'EC_RT';%For DKE code parameters 0041 path_dkeparam = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0042 % 0043 id_dkedisplay = 'NO_DISPLAY';%For output code display 0044 path_dkedisplay = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0045 % 0046 id_ohm = '';%For Ohmic electric contribution 0047 path_ohm = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0048 % 0049 ids_wavestruct = {'C3PO_NTM_32_24rays'};%For RF waves contribution (put all the type of waves needed) 0050 paths_wavestruct = {''};%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0051 % 0052 id_transpfaste = '';%For fast electron radial transport 0053 path_transpfaste = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0054 % 0055 id_ripple = '';%For fast electron magnetic ripple losses 0056 path_ripple = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0057 % 0058 %************************************************************************************************************************************ 0059 % 0060 [dkepath,equil,dkeparam,dkedisplay,ohm,wavestructs,transpfaste,ripple,fluct] = load_structures_yp('dkepath',id_dkepath,path_dkepath,'equil',id_equil,path_equil,'dkeparam',id_dkeparam,path_dkeparam,'dkedisplay',id_dkedisplay,path_dkedisplay,'ohm',id_ohm,path_ohm,'wavestructs',ids_wavestruct,paths_wavestruct,'transpfaste',id_transpfaste,path_transpfaste,'ripple',id_ripple,path_ripple,'fluct',id_fluct,path_fluct); 0061 % 0062 id_simul = [id_equil,'_',wavestructs{1}.id,'_',id_fluct,locid_simul]; 0063 % 0064 fluct_fitparam.method = 'pchip';%nearest,spline,pchip,cubic 0065 fluct_fitparam.nharm = 32;%Number of harmonics in the plasma fluctuations interpolation (less than ntheta_equil/2) 0066 fluct_fitparam.ngridresample = 201;%Number of grid points for resampling the radial profile of plasma fluctuations parameters (very slow if too big !!) 0067 % 0068 %************************************************************************************************************************************ 0069 % 0070 dkeparam.rt_mode = 1; 0071 dkeparam.rho_S = rho_S; 0072 dkeparam.psin_S = psin_S; 0073 dkeparam.timevol = 1;%to calculate time evolution of moments 0074 % 0075 dkedisplay.display_mode = 0; 0076 dkedisplay.display_time_mode = 0; 0077 % 0078 equil.fluct = fluct; 0079 equil.fluct.fitparam = fluct_fitparam; 0080 % 0081 %************************************************************************************************************************************ 0082 % 0083 luke_input.id_simul = id_simul; 0084 luke_input.dkeparam = dkeparam; 0085 luke_input.dkedisplay = dkedisplay; 0086 luke_input.equil = equil; 0087 luke_input.waves = ''; 0088 luke_input.wavestructs = wavestructs; 0089 luke_input.ohm = ohm; 0090 luke_input.transpfaste = transpfaste; 0091 luke_input.ripple = ripple; 0092 luke_input.Zf0_interp = ''; 0093 % 0094 frundke_fluctn(luke_input,dkepath,opt_save,p_opt,path_simul,dtn,tn,Nfluct,sigmar_max,sigmar_hwhm); 0095 %