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 % Note : tn is adjusted for a direct comparison with the fluctuation case 0014 % 0015 tn = [94.9783];% Final simulation time (normalized to the central thermal collision time) 0016 dtn = [100i];% internal LUKE time step. As an imaginary number, dkeparam.dtn prescribes the number of internal time steps 0017 % 0018 % *********************** Specify LUKE structures ***************************** 0019 % 0020 locid_simul = '';%Simulation ID 0021 path_simul = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0022 % 0023 psin_S = [];%Normalized poloidal flux grid where calculations are performed (0 < psin_S < 1) (If one value: local calculation only, not used if empty) 0024 rho_S = [0:0.1:0.6,0.61:0.01:0.79,0.80:0.05:1]; 0025 % 0026 id_dkepath = '';%For all paths used by DKE solver 0027 path_dkepath = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0028 % 0029 id_equil = 'ITER_Scen2_200103121816_430_129';%For plasma equilibrium 0030 path_equil = 'EQUIL/';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0031 % 0032 id_fluct = '';%For fast electron radial transport 0033 path_fluct = 'FLUCT/';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0034 % 0035 id_dkeparam = 'EC_RT';%For DKE code parameters 0036 path_dkeparam = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0037 % 0038 id_dkedisplay = 'NO_DISPLAY';%For output code display 0039 path_dkedisplay = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0040 % 0041 id_ohm = '';%For Ohmic electric contribution 0042 path_ohm = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0043 % 0044 ids_wavestruct = {'C3PO_NTM_32_24rays'};%For RF waves contribution (put all the type of waves needed) 0045 paths_wavestruct = {''};%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0046 % 0047 id_transpfaste = '';%For fast electron radial transport 0048 path_transpfaste = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0049 % 0050 id_ripple = '';%For fast electron magnetic ripple losses 0051 path_ripple = '';%if nothing is specified, the working directory is first used and then MatLab is looking in all the path 0052 % 0053 %************************************************************************************************************************************ 0054 % 0055 [dkepath,equil,dkeparam,dkedisplay,ohm,wavestructs,transpfaste,ripple,equil.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); 0056 % 0057 %************************************************************************************************************************************ 0058 % 0059 simul.id = make_luke_simulid_jd(equil,ohm,transpfaste,ripple,'',wavestructs,locid_simul); 0060 simul.path = path_simul; 0061 % 0062 dkeparam.timevol = 1;%to calculate time evolution of moments 0063 dkeparam.tn = tn;% as an imaginary number, dkeparam.tn prescribes the number of fluctuation iterations 0064 dkeparam.dtn = dtn;% internal LUKE time step. As an imaginary number, dkeparam.dtn prescribes the number of internal time steps 0065 dkeparam.psin_S = psin_S; 0066 dkeparam.rho_S = rho_S; 0067 % 0068 lukestructs.simul = simul; 0069 lukestructs.equil = equil; 0070 lukestructs.dkeparam = dkeparam; 0071 lukestructs.dkedisplay = dkedisplay; 0072 lukestructs.ohm = ohm; 0073 lukestructs.wavestructs = wavestructs; 0074 lukestructs.transpfaste = transpfaste; 0075 lukestructs.ripple = ripple; 0076 % 0077 lukestructs = run_lukert(lukestructs); 0078 % 0079 % ---------------------------------------------- Save results ---------------------------------------------- 0080 % 0081 filename = [path_simul,'LUKE_RESULTS_',simul.id,'.mat']; 0082 save(filename); 0083 info_dke_yp(2,['LUKE Results saved in ',filename]); 0084 % 0085 opt.ir_display = 0; 0086 opt.rho = 'g'; 0087 opt.peakmode = 0; 0088 opt.peak = 'n'; 0089 opt.wave = 'n'; 0090 opt.spec = 0; 0091 opt.diaryname = ['res_',lukestructs.simul.id,'.log']; 0092 % 0093 proc_luke_jd(lukestructs.output,opt); 0094 % 0095 0096 0097