Commit 0a7fb2d6 authored by Olivier Sauter's avatar Olivier Sauter
Browse files

temp

parent 2e85bf0c
Pipeline #94213 passed with stages
in 3 minutes
function [gdat_data,gdat_params,error_status,varargout] = gdat_tcv_OS(shot,data_request,varargin) function [gdat_data,gdat_params,error_status,varargout] = gdat_tcv(shot,data_request,varargin)
% %
% function [gdat_data,gdat_params,error_status,varargout] = gdat(shot,data_request,varargin) % function [gdat_data,gdat_params,error_status,varargout] = gdat(shot,data_request,varargin)
% %
...@@ -2574,8 +2574,8 @@ elseif strcmp(mapping_for_tcv.method,'switchcase') ...@@ -2574,8 +2574,8 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
end end
if any(strmatch('dnbi',gdat_data.gdat_params.source)) if any(strmatch('dnbi',gdat_data.gdat_params.source))
% NB2 % NB2
nodenameeff = '\RESULTS::DNBI:POWR_TCV'; nodenameeff = '\RESULTS::DNBI:POWR_TCV';
nb2_data_tdi = tdi(nodenameeff); nb2_data_tdi = tdi(nodenameeff);
if ~isempty(nb2_data_tdi.data) && ~ischar(nb2_data_tdi.data) && ~isempty(nb2_data_tdi.dim) if ~isempty(nb2_data_tdi.data) && ~ischar(nb2_data_tdi.data) && ~isempty(nb2_data_tdi.dim)
nbi_neutral_power_tot = nb2_data_tdi.data.*1e6; % in W nbi_neutral_power_tot = nb2_data_tdi.data.*1e6; % in W
...@@ -2706,7 +2706,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase') ...@@ -2706,7 +2706,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
end end
case {'phi_tor', 'phitor', 'toroidal_flux'} case {'phi_tor', 'phitor', 'toroidal_flux'}
% Phi(LCFS) = int(Bphi dSphi), can use Eq.(11) of "Tokamak coordinate conventions: COCOS" paper: % Phi(LCFS) = int(Bphi dSphi), can use Eq.(11) of "Tokamak coordinate conventions: COCOS" paper:
% O. Sauter, S.Yu. Medvedev, Comput. Phys. Commun. 184 (2013) 293302 % O. Sauter, S.Yu. Medvedev, Comput. Phys. Commun. 184 (2013) 293???302
% since cocos=17 for LIUQE we get: % since cocos=17 for LIUQE we get:
% q = -dPhi/dpsi => Phi = - int(q*dpsi) which should always have the sign of B0 % q = -dPhi/dpsi => Phi = - int(q*dpsi) which should always have the sign of B0
% need to get q_rho but to avoid loop for rhotor in grids_1d, get q_rho explicitely here % need to get q_rho but to avoid loop for rhotor in grids_1d, get q_rho explicitely here
...@@ -2928,7 +2928,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase') ...@@ -2928,7 +2928,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
case {'rhotor_edge', 'rhotor', 'rhotor_norm'} case {'rhotor_edge', 'rhotor', 'rhotor_norm'}
% Phi(LCFS) = int(Bphi dSphi), can use Eq.(11) of "Tokamak coordinate conventions: COCOS" paper: % Phi(LCFS) = int(Bphi dSphi), can use Eq.(11) of "Tokamak coordinate conventions: COCOS" paper:
% O. Sauter, S.Yu. Medvedev, Comput. Phys. Commun. 184 (2013) 293302 % O. Sauter, S.Yu. Medvedev, Comput. Phys. Commun. 184 (2013) 293???302
% since cocos=17 for LIUQE we get: % since cocos=17 for LIUQE we get:
% q = -dPhi/dpsi => Phi = - int(q*dpsi) which should always have the sign of B0 % q = -dPhi/dpsi => Phi = - int(q*dpsi) which should always have the sign of B0
% need to get q_rho but to avoid loop for rhotor in grids_1d, get q_rho explicitely here % need to get q_rho but to avoid loop for rhotor in grids_1d, get q_rho explicitely here
...@@ -3315,7 +3315,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase') ...@@ -3315,7 +3315,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
gdat_data.x = [1:size(sig,1)]; gdat_data.x = [1:size(sig,1)];
gdat_data.dimunits = {'20 chords per camera'; 's'}; gdat_data.dimunits = {'20 chords per camera'; 's'};
else else
% keyboard keyboard
% extract only given channels % extract only given channels
gdat_data.data = sig(channel_xtomo,:); gdat_data.data = sig(channel_xtomo,:);
gdat_data.x = channel_xtomo; gdat_data.x = channel_xtomo;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment