Skip to content
Snippets Groups Projects
Commit d2439390 authored by Olivier Sauter's avatar Olivier Sauter
Browse files

add freq etc from ECS param set, add list of keywords in help loadAUGdata

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@4384 d63d8f72-b253-0410-a779-e742ad2e26cf
parent 1e238903
No related branches found
No related tags found
No related merge requests found
......@@ -52,6 +52,38 @@ function [trace,error,varargout]=loadAUGdata(shot,data_type,varargin)
% [ip,error]=loadAUGdata(shot,'MAG/Ipi');
% [n2,error]=loadAUGdata(shot,'MOD/EvenNAmp');
%
% List of keywords (can be used in adition to 'DIAG/tracenam'), with comments when not obvious:
%
% Ip :
% b0 :
% zmag :
% rmag :
% rgeo :
% zgeo :
% q0 :
% q95 :
% kappa :
% delta :
% deltatop :
% deltabot :
% neint :
% neterho :
% cxrs : CXRS structure from CEZ with vrot, Ti, ...
% cxrs_rho : same as cxrs but project on rho as well (using 'equil' information)
% equil : equilibrium structure, rhopol, rhotor, rhovol, etc from EQI
% equil_fpp : as equil but from FPP
% equil_eqm : as equil but from EQM
% equil_eqr : as equil but from EQR
% sxr : from SXR/A or B (from old stuff, not sure still OK)
% sxR : from SXR/A or B adding R of chords (from old stuff, not sure still OK)
% sxb : 'SXB/J' chords
% ece :
% eced :
% Halpha :
% pgyro : for each gyrotrons, power, freq, etc (ask for more)
% powers : power traces for each sources
%
varargout={cell(1,1)};
error=1;
......@@ -901,12 +933,21 @@ switch AUGkeywrdcase{index}
trace_all.ecrh1(i) = a;
trace_all.pgyro(:,i) = reshape(a.data,nb_timepoints,1);
end
tracename_eff = ['f-G' num2str(i)];
[a,e]=rdaAUG_eff(shot,ppftype,tracename_eff,shotfile_exp_eff);
if isempty(a) || e~=0
a = sf2par('ECS',shot,'gyr_freq',['P_sy1_g' num2str(i)]);
if isempty(a)
else
trace_all.freq_ecrh1(i) = a;
trace_all.freq_ech(i) = a.data(1);
trace_all.freq_ech(i) = a.value;
end
a = sf2par('ECS',shot,'GPolPos',['P_sy1_g' num2str(i)]);
if isempty(a)
else
trace_all.polpos_ecs(i) = a.value;
end
a = sf2par('ECS',shot,'GTorPos',['P_sy1_g' num2str(i)]);
if isempty(a)
else
trace_all.torpos_ecs(i) = a.value;
end
% "new" ECRH2 gyrotrons: gyro 5 to 8 in pgyro
tracename_eff = ['PG' num2str(i) 'N'];
......@@ -916,12 +957,21 @@ switch AUGkeywrdcase{index}
trace_all.ecrh2(i) = a;
trace_all.pgyro(:,i+4) = reshape(a.data,nb_timepoints,1);
end
tracename_eff = ['f-G' num2str(i) 'N'];
[a,e]=rdaAUG_eff(shot,ppftype,tracename_eff,shotfile_exp_eff);
if isempty(a) || e~=0
a = sf2par('ECS',shot,'gyr_freq',['P_sy2_g' num2str(i)]);
if isempty(a)
else
trace_all.freq_ecrh2(i) = a;
trace_all.freq_ech(i+4) = a.data(1);
trace_all.freq_ech(i+4) = a.value;
end
a = sf2par('ECS',shot,'GPolPos',['P_sy2_g' num2str(i)]);
if isempty(a)
else
trace_all.polpos_ecs(i+4) = a.value;
end
a = sf2par('ECS',shot,'GTorPos',['P_sy2_g' num2str(i)]);
if isempty(a)
else
trace_all.torpos_ecs(i+4) = a.value;
end
end
if ~isempty(trace_all)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment