From 506176ad657727a3a4a0d61423a598b1e768a723 Mon Sep 17 00:00:00 2001
From: Luke Simons <luke.simons@epfl.ch>
Date: Wed, 23 Apr 2025 10:48:04 +0200
Subject: [PATCH 1/8] Fixes to tcv_get_ids_nbi for position, comments and COCOS
dependent on IMAS version
---
matlab/TCV_IMAS/tcv2ids.m | 2 +-
matlab/TCV_IMAS/tcv_get_ids_nbi.m | 13 +++++++------
2 files changed, 8 insertions(+), 7 deletions(-)
diff --git a/matlab/TCV_IMAS/tcv2ids.m b/matlab/TCV_IMAS/tcv2ids.m
index 312cb4bc..aa881ef4 100644
--- a/matlab/TCV_IMAS/tcv2ids.m
+++ b/matlab/TCV_IMAS/tcv2ids.m
@@ -16,7 +16,7 @@ function [ids_from_tcv,varargout] = tcv2ids(shot,varargin);
% empty or 'delta' (default): meaning difference in upper is set (standard error_bar
% 'added': errorbar is added: upper=data+delta and lower=data-delta
% 'delta_with_lower': as 'delta' but lower also set
-% 'cocos_out': (default 11) cocos to transform ids from TCV cocos_in=17 to cocos_out
+% 'cocos_out': (DD3 default 11, DD4 default 17) cocos to transform ids from TCV cocos_in=17 to cocos_out
% 'ipsign_out': if a specific sign fo Ip is desired in output within the cocos_out system (default 0=no specific sign)
% 'b0sign_out': if a specific sign fo B0 is desired in output within the cocos_out system (default 0=no specific sign)
% 'nverbose': (default 1), set it to 3 to have more messages, for example about not fully valid nodes when doing transformation (empty or Nans)
diff --git a/matlab/TCV_IMAS/tcv_get_ids_nbi.m b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
index 9d5cf6a5..0b394b77 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_nbi.m
+++ b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
@@ -5,6 +5,7 @@ function [ids_nbi,ids_nbi_description,varargout] = tcv_get_ids_nbi(shot,ids_nbi_
%
% gdat_params: gdat_data.gdat_params to get all params passed from original call, in particular error_bar options
%
+imas_version_number=getenv('IMAS_VERSION');
if exist('gdat_params')
[ids_nbi, params_nbi] = tcv_ids_headpart(shot,ids_nbi_empty,'nbi','homogeneous_time',0,'gdat_params',gdat_params,varargin{:});
@@ -22,8 +23,7 @@ ids_nbi_description='';
% "global_quantities_desc" which contains the same subfields themselves containing the gdat string aftre shot used
%
-
-nb_units = 3; % assume 2 units: 1st NBH and DNBI
+nb_units = 3; % assume 3 units: 1st NB1. 2nd NB2 and DNBI
ids_nbi.unit(1:nb_units) = ids_nbi.unit(1); % copy empty structure for all units, then fill in
% create lists of what is different for each units so that can scan through units
@@ -76,9 +76,9 @@ beamlets_group.position(1).z = 0.;
beamlets_group.position(2).phi = 58.8255*pi/180.;
beamlets_group.position(2).r = 4.5889;
beamlets_group.position(2).z = 0.;
-beamlets_group.position(2).phi = 295.2416*pi/180.;
-beamlets_group.position(2).r = 4.9274;
-beamlets_group.position(2).z = 0.;
+beamlets_group.position(3).phi = 295.2416*pi/180.;
+beamlets_group.position(3).r = 4.9274;
+beamlets_group.position(3).z = 0.;
params_eff = params_eff_ref;
for iunit=1:nb_units
@@ -172,7 +172,8 @@ end
% cocos automatic transform
-if exist('ids_generic_cocos_nodes_transformation_symbolic') == 2
+if exist('ids_generic_cocos_nodes_transformation_symbolic') == 2 && ...
+ strcmp(imas_version_number(1),'3')
[ids_nbi,cocoscoeff]=ids_generic_cocos_nodes_transformation_symbolic(ids_nbi,'nbi',gdat_params.cocos_in, ...
gdat_params.cocos_out,gdat_params.ipsign_out,gdat_params.b0sign_out,gdat_params.ipsign_in,gdat_params.b0sign_in, ...
gdat_params.error_bar,gdat_params.nverbose);
--
GitLab
From efe1e6263508a684281055b0ac1487a88058c3fa Mon Sep 17 00:00:00 2001
From: Luke Simons <luke.simons@epfl.ch>
Date: Wed, 23 Apr 2025 11:57:54 +0200
Subject: [PATCH 2/8] Added comment about data source
---
matlab/TCV_IMAS/tcv_get_ids_nbi.m | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/matlab/TCV_IMAS/tcv_get_ids_nbi.m b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
index 0b394b77..e9ba38c1 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_nbi.m
+++ b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
@@ -5,6 +5,10 @@ function [ids_nbi,ids_nbi_description,varargout] = tcv_get_ids_nbi(shot,ids_nbi_
%
% gdat_params: gdat_data.gdat_params to get all params passed from original call, in particular error_bar options
%
+%
+% Most information taken from NBI and NB_Model wiki pages:
+% https://spcwiki.epfl.ch/wiki/NBI
+% https://spcwiki.epfl.ch/wiki/NB_Model
imas_version_number=getenv('IMAS_VERSION');
if exist('gdat_params')
--
GitLab
From 85b59b7926ee51c8bc050199b8e3bc956d416fef Mon Sep 17 00:00:00 2001
From: Luke Simons <luke.simons@epfl.ch>
Date: Wed, 23 Apr 2025 13:37:58 +0200
Subject: [PATCH 3/8] Changed handling of COCOS transformation in
tcv_get_ids_nbi to use gdat cocos_out convention
---
matlab/TCV_IMAS/tcv2ids.m | 8 +++++++-
1 file changed, 7 insertions(+), 1 deletion(-)
diff --git a/matlab/TCV_IMAS/tcv2ids.m b/matlab/TCV_IMAS/tcv2ids.m
index aa881ef4..a78490fb 100644
--- a/matlab/TCV_IMAS/tcv2ids.m
+++ b/matlab/TCV_IMAS/tcv2ids.m
@@ -29,6 +29,12 @@ function [ids_from_tcv,varargout] = tcv2ids(shot,varargin);
% varargout{1}: return also the ids in array of structure with the names, to allow easy use of plotallids
%
+imas_version_number=str2num(getenv('IMAS_VERSION'));
+cocos_out_default=11;
+if imas_version_number(1) >= 4
+ cocos_out_default=17;
+end
+
% initialize input parser
p = inputParser;
% no required inputs here so one can call with empty args and get defaults parameters
@@ -37,7 +43,7 @@ ids_names = tcv_available_ids;
p.addOptional('shot', [], @(x) (isnumeric(x) && isscalar(x) && (x == round(x)))); % integer
p.addOptional('ids_names', ids_names, @(x) isempty(x) || iscell(x) ); % char or cell array
p.addOptional('error_bar', 'delta', @(x) isempty(x) || ischar(x) ); % char or cell array
-p.addOptional('cocos_out', 11, @(x) isempty(x) || isnumeric(x) ); % char
+p.addOptional('cocos_out', cocos_out_default, @(x) isempty(x) || isnumeric(x) ); % char
p.addOptional('ipsign_out', 0, @(x) isempty(x) || (isscalar(x) && (x==0 || x==-1 || x==+1)) );
p.addOptional('b0sign_out', 0, @(x) isempty(x) || (isscalar(x) && (x==0 || x==-1 || x==+1)) );
p.addOptional('nverbose', 1, @(x) isempty(x) || isnumeric(x) );
--
GitLab
From e711fb1f9c537b5c7ae6e2e19c43545b724e5937 Mon Sep 17 00:00:00 2001
From: Luke Simons <luke.simons@epfl.ch>
Date: Wed, 23 Apr 2025 13:54:11 +0200
Subject: [PATCH 4/8] Changed handling of COCOS transformation in
tcv_get_ids_nbi to use gdat cocos_out convention
---
matlab/TCV_IMAS/tcv2ids.m | 4 ++--
matlab/TCV_IMAS/tcv_get_ids_nbi.m | 4 +---
2 files changed, 3 insertions(+), 5 deletions(-)
diff --git a/matlab/TCV_IMAS/tcv2ids.m b/matlab/TCV_IMAS/tcv2ids.m
index a78490fb..9e034a8a 100644
--- a/matlab/TCV_IMAS/tcv2ids.m
+++ b/matlab/TCV_IMAS/tcv2ids.m
@@ -29,9 +29,9 @@ function [ids_from_tcv,varargout] = tcv2ids(shot,varargin);
% varargout{1}: return also the ids in array of structure with the names, to allow easy use of plotallids
%
-imas_version_number=str2num(getenv('IMAS_VERSION'));
+imas_version_number=getenv('IMAS_VERSION');
cocos_out_default=11;
-if imas_version_number(1) >= 4
+if str2num(imas_version_number(1)) >= 4
cocos_out_default=17;
end
diff --git a/matlab/TCV_IMAS/tcv_get_ids_nbi.m b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
index e9ba38c1..e4d15d27 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_nbi.m
+++ b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
@@ -9,7 +9,6 @@ function [ids_nbi,ids_nbi_description,varargout] = tcv_get_ids_nbi(shot,ids_nbi_
% Most information taken from NBI and NB_Model wiki pages:
% https://spcwiki.epfl.ch/wiki/NBI
% https://spcwiki.epfl.ch/wiki/NB_Model
-imas_version_number=getenv('IMAS_VERSION');
if exist('gdat_params')
[ids_nbi, params_nbi] = tcv_ids_headpart(shot,ids_nbi_empty,'nbi','homogeneous_time',0,'gdat_params',gdat_params,varargin{:});
@@ -176,8 +175,7 @@ end
% cocos automatic transform
-if exist('ids_generic_cocos_nodes_transformation_symbolic') == 2 && ...
- strcmp(imas_version_number(1),'3')
+if exist('ids_generic_cocos_nodes_transformation_symbolic') == 2
[ids_nbi,cocoscoeff]=ids_generic_cocos_nodes_transformation_symbolic(ids_nbi,'nbi',gdat_params.cocos_in, ...
gdat_params.cocos_out,gdat_params.ipsign_out,gdat_params.b0sign_out,gdat_params.ipsign_in,gdat_params.b0sign_in, ...
gdat_params.error_bar,gdat_params.nverbose);
--
GitLab
From eb7e363b1a1b26f2bca0f52e8e3a44f0a9d9eb2a Mon Sep 17 00:00:00 2001
From: Luke Simons <luke.simons@epfl.ch>
Date: Wed, 23 Apr 2025 13:59:56 +0200
Subject: [PATCH 5/8] Added try catch block incase IMAS_VERSION is undefined
---
matlab/TCV_IMAS/tcv2ids.m | 10 +++++++---
1 file changed, 7 insertions(+), 3 deletions(-)
diff --git a/matlab/TCV_IMAS/tcv2ids.m b/matlab/TCV_IMAS/tcv2ids.m
index 9e034a8a..1cc251db 100644
--- a/matlab/TCV_IMAS/tcv2ids.m
+++ b/matlab/TCV_IMAS/tcv2ids.m
@@ -30,9 +30,13 @@ function [ids_from_tcv,varargout] = tcv2ids(shot,varargin);
%
imas_version_number=getenv('IMAS_VERSION');
-cocos_out_default=11;
-if str2num(imas_version_number(1)) >= 4
- cocos_out_default=17;
+cocos_out_default=17;
+try
+ if str2num(imas_version_number(1)) < 4
+ cocos_out_default=11;
+ end
+catch
+ warning('IMAS_VERSION undefined, assuming cocos_out=17');
end
% initialize input parser
--
GitLab
From a13d4a989b3be50fcad5a47cbd2b6759ba4d1dae Mon Sep 17 00:00:00 2001
From: Luke Simons <luke.simons@epfl.ch>
Date: Wed, 23 Apr 2025 14:07:30 +0200
Subject: [PATCH 6/8] Changed warning to disp
---
matlab/TCV_IMAS/tcv2ids.m | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/matlab/TCV_IMAS/tcv2ids.m b/matlab/TCV_IMAS/tcv2ids.m
index 1cc251db..7c0a293f 100644
--- a/matlab/TCV_IMAS/tcv2ids.m
+++ b/matlab/TCV_IMAS/tcv2ids.m
@@ -36,7 +36,7 @@ try
cocos_out_default=11;
end
catch
- warning('IMAS_VERSION undefined, assuming cocos_out=17');
+ disp('IMAS_VERSION undefined, assuming cocos_out=17');
end
% initialize input parser
--
GitLab
From 0d5ad4a6874d0e087dc10c8b86ea1c0ff8265e1d Mon Sep 17 00:00:00 2001
From: Luke Simons <luke.simons@epfl.ch>
Date: Wed, 23 Apr 2025 14:56:48 +0200
Subject: [PATCH 7/8] Changes to have dynamically allocated NBI systems in DD
---
matlab/TCV_IMAS/tcv_get_ids_nbi.m | 86 ++++++++++++++-----------------
1 file changed, 40 insertions(+), 46 deletions(-)
diff --git a/matlab/TCV_IMAS/tcv_get_ids_nbi.m b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
index e4d15d27..6eae1073 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_nbi.m
+++ b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
@@ -9,6 +9,7 @@ function [ids_nbi,ids_nbi_description,varargout] = tcv_get_ids_nbi(shot,ids_nbi_
% Most information taken from NBI and NB_Model wiki pages:
% https://spcwiki.epfl.ch/wiki/NBI
% https://spcwiki.epfl.ch/wiki/NB_Model
+i_nbi1=1; i_nbi2=2; i_dnbi=3;
if exist('gdat_params')
[ids_nbi, params_nbi] = tcv_ids_headpart(shot,ids_nbi_empty,'nbi','homogeneous_time',0,'gdat_params',gdat_params,varargin{:});
@@ -27,15 +28,20 @@ ids_nbi_description='';
%
nb_units = 3; % assume 3 units: 1st NB1. 2nd NB2 and DNBI
-ids_nbi.unit(1:nb_units) = ids_nbi.unit(1); % copy empty structure for all units, then fill in
-
% create lists of what is different for each units so that can scan through units
unit_identifier = {'NB1', 'NB2', 'DNBI'};
unit_name = {'25keV 1st NBH source', '50keV 2nd NBH source', 'diagnostic NBI'};
results_subname = {'nb1', 'nb2', 'dnbi'};
+
if shot<70811
- results_subname = {'nbh', 'nb2', 'dnbi'};
+ unit_identifier = {'NB1', 'DNBI'};
+ unit_name = {'25keV 1st NBH source', 'diagnostic NBI'};
+ results_subname = {'nbh','dnbi'};
+ nb_units = 2; % assume 2 units: 1st NBI and DNBI
end
+nbi_indices=1:nb_units;
+ids_nbi.unit(nbi_indices) = ids_nbi.unit(1); % copy empty structure for all units, then fill in
+
species.a = [2., 2., 1.];
species.z_n = [1., 1., 1.];
species.label = {'D', 'D','H'};
@@ -46,45 +52,45 @@ beamlets_group.width_horizontal = [250, 250, 87.2]*1e-3;
beamlets_group.width_vertical = [250, 250, 87.2]*1e-3;
beamlets_group.focus(1:nb_units)=struct('focal_length_horizontal',[],'focal_length_vertical',[],'width_min_horizontal',[],'width_min_vertical',[]);
-beamlets_group.focus(1).focal_length_horizontal = 3.76;
-beamlets_group.focus(1).focal_length_vertical = 3.98;
-beamlets_group.focus(1).width_min_horizontal = 21.6*1e-2;
-beamlets_group.focus(1).width_min_vertical = 9.4*1e-2;
+beamlets_group.focus(i_nbi1).focal_length_horizontal = 3.76;
+beamlets_group.focus(i_nbi1).focal_length_vertical = 3.98;
+beamlets_group.focus(i_nbi1).width_min_horizontal = 21.6*1e-2;
+beamlets_group.focus(i_nbi1).width_min_vertical = 9.4*1e-2;
% So far NB2 parameters are merely a copy of NB1 parameters
-beamlets_group.focus(2).focal_length_horizontal = 3.76;
-beamlets_group.focus(2).focal_length_vertical = 3.98;
-beamlets_group.focus(2).width_min_horizontal = 21.6*1e-2;
-beamlets_group.focus(2).width_min_vertical = 9.4*1e-2;
-beamlets_group.focus(3).focal_length_horizontal = 1.8;
-beamlets_group.focus(3).focal_length_vertical = 1.8;
-beamlets_group.focus(3).width_min_horizontal = 12.1*1e-2;
-beamlets_group.focus(3).width_min_vertical = 12.1*1e-2;
+beamlets_group.focus(i_nbi2).focal_length_horizontal = 3.76;
+beamlets_group.focus(i_nbi2).focal_length_vertical = 3.98;
+beamlets_group.focus(i_nbi2).width_min_horizontal = 21.6*1e-2;
+beamlets_group.focus(i_nbi2).width_min_vertical = 9.4*1e-2;
+beamlets_group.focus(i_dnbi).focal_length_horizontal = 1.8;
+beamlets_group.focus(i_dnbi).focal_length_vertical = 1.8;
+beamlets_group.focus(i_dnbi).width_min_horizontal = 12.1*1e-2;
+beamlets_group.focus(i_dnbi).width_min_vertical = 12.1*1e-2;
beamlets_group.divergence(1:nb_units) = struct('particle_fraction',[],'vertical',[],'horizontal',[]);
-beamlets_group.divergence(1).particle_fraction = 1.;
-beamlets_group.divergence(1).vertical = 0.59 *pi/180.;
-beamlets_group.divergence(1).horizontal = 1.4 *pi/180.;
-beamlets_group.divergence(2).particle_fraction = 1.;
-beamlets_group.divergence(2).vertical = 0.59 *pi/180.;
-beamlets_group.divergence(2).horizontal = 1.4 *pi/180.;
-beamlets_group.divergence(3).particle_fraction = 1.;
-beamlets_group.divergence(3).vertical = 0.53 *pi/180.;
-beamlets_group.divergence(3).horizontal = 0.53 *pi/180.;
+beamlets_group.divergence(i_nbi1).particle_fraction = 1.;
+beamlets_group.divergence(i_nbi1).vertical = 0.59 *pi/180.;
+beamlets_group.divergence(i_nbi1).horizontal = 1.4 *pi/180.;
+beamlets_group.divergence(i_nbi2).particle_fraction = 1.;
+beamlets_group.divergence(i_nbi2).vertical = 0.59 *pi/180.;
+beamlets_group.divergence(i_nbi2).horizontal = 1.4 *pi/180.;
+beamlets_group.divergence(i_dnbi).particle_fraction = 1.;
+beamlets_group.divergence(i_dnbi).vertical = 0.53 *pi/180.;
+beamlets_group.divergence(i_dnbi).horizontal = 0.53 *pi/180.;
%dcd_NBH = psitbxdcd(4.5889, 0.0, 211.9535*pi/180, 0.0, -9.2308*pi/180);
beamlets_group.position(1:nb_units) = struct('phi',[],'r',[],'z',[]);
-beamlets_group.position(1).phi = 211.9535*pi/180.;
-beamlets_group.position(1).r = 4.5889;
-beamlets_group.position(1).z = 0.;
-beamlets_group.position(2).phi = 58.8255*pi/180.;
-beamlets_group.position(2).r = 4.5889;
-beamlets_group.position(2).z = 0.;
-beamlets_group.position(3).phi = 295.2416*pi/180.;
-beamlets_group.position(3).r = 4.9274;
-beamlets_group.position(3).z = 0.;
+beamlets_group.position(i_nbi1).phi = 211.9535*pi/180.;
+beamlets_group.position(i_nbi1).r = 4.5889;
+beamlets_group.position(i_nbi1).z = 0.;
+beamlets_group.position(i_nbi2).phi = 58.8255*pi/180.;
+beamlets_group.position(i_nbi2).r = 4.5889;
+beamlets_group.position(i_nbi2).z = 0.;
+beamlets_group.position(i_dnbi).phi = 295.2416*pi/180.;
+beamlets_group.position(i_dnbi).r = 4.9274;
+beamlets_group.position(i_dnbi).z = 0.;
params_eff = params_eff_ref;
-for iunit=1:nb_units
+for iunit=nbi_indices
ids_nbi.unit{iunit}.identifier = unit_identifier{iunit};
ids_nbi.unit{iunit}.name = unit_name{iunit};
%% power
@@ -131,21 +137,9 @@ for iunit=1:nb_units
% https://spcwiki.epfl.ch/wiki/NB_Model
ids_nbi.unit{iunit}.beamlets_group{1}.direction = beamlets_group.direction(iunit); %clockwise
ids_nbi.unit{iunit}.beamlets_group{1}.tangency_radius = beamlets_group.tangency_radius(iunit);
- % ids_nbi.unit{iunit}.beamlets_group{1}.tangency_radius_error_index: -999999999
- % ids_nbi.unit{iunit}.beamlets_group{1}.tangency_radius_error_lower: -9.0000e+40
- % ids_nbi.unit{iunit}.beamlets_group{1}.tangency_radius_error_upper: -9.0000e+40
ids_nbi.unit{iunit}.beamlets_group{1}.angle = beamlets_group.angle(iunit); %injection parallel to midplane
- % ids_nbi.unit{iunit}.beamlets_group{1}.angle_error_index: -999999999
- % ids_nbi.unit{iunit}.beamlets_group{1}.angle_error_lower: -9.0000e+40
- % ids_nbi.unit{iunit}.beamlets_group{1}.angle_error_upper: -9.0000e+40
ids_nbi.unit{iunit}.beamlets_group{1}.width_horizontal = beamlets_group.width_horizontal(iunit);
- % ids_nbi.unit{iunit}.beamlets_group{1}.width_horizontal_error_index: -999999999
- % ids_nbi.unit{iunit}.beamlets_group{1}.width_horizontal_error_lower: -9.0000e+40
- % ids_nbi.unit{iunit}.beamlets_group{1}.width_horizontal_error_upper: -9.0000e+40
ids_nbi.unit{iunit}.beamlets_group{1}.width_vertical = beamlets_group.width_vertical(iunit);
- % ids_nbi.unit{iunit}.beamlets_group{1}.width_vertical_error_index: -999999999
- % ids_nbi.unit{iunit}.beamlets_group{1}.width_vertical_error_lower: -9.0000e+40
- % ids_nbi.unit{iunit}.beamlets_group{1}.width_vertical_error_upper: -9.0000e+40
% Should always copy "leaves" only to avoid deleting non-filled in values like error_bars
fields_to_copy = fieldnames(beamlets_group.focus(iunit));
--
GitLab
From a39c66b1e744fa426e0c4dbb0b8e33e9c8137e2f Mon Sep 17 00:00:00 2001
From: Luke Simons <luke.simons@epfl.ch>
Date: Wed, 23 Apr 2025 13:13:03 +0000
Subject: [PATCH 8/8] Small suggested changes to comments
---
matlab/TCV_IMAS/tcv2ids.m | 4 ++--
matlab/TCV_IMAS/tcv_get_ids_nbi.m | 1 -
2 files changed, 2 insertions(+), 3 deletions(-)
diff --git a/matlab/TCV_IMAS/tcv2ids.m b/matlab/TCV_IMAS/tcv2ids.m
index 7c0a293f..83662207 100644
--- a/matlab/TCV_IMAS/tcv2ids.m
+++ b/matlab/TCV_IMAS/tcv2ids.m
@@ -16,7 +16,7 @@ function [ids_from_tcv,varargout] = tcv2ids(shot,varargin);
% empty or 'delta' (default): meaning difference in upper is set (standard error_bar
% 'added': errorbar is added: upper=data+delta and lower=data-delta
% 'delta_with_lower': as 'delta' but lower also set
-% 'cocos_out': (DD3 default 11, DD4 default 17) cocos to transform ids from TCV cocos_in=17 to cocos_out
+% 'cocos_out': (DD3 and before default 11, DD4 default 17) cocos to transform ids from TCV cocos_in=17 to cocos_out
% 'ipsign_out': if a specific sign fo Ip is desired in output within the cocos_out system (default 0=no specific sign)
% 'b0sign_out': if a specific sign fo B0 is desired in output within the cocos_out system (default 0=no specific sign)
% 'nverbose': (default 1), set it to 3 to have more messages, for example about not fully valid nodes when doing transformation (empty or Nans)
@@ -36,7 +36,7 @@ try
cocos_out_default=11;
end
catch
- disp('IMAS_VERSION undefined, assuming cocos_out=17');
+ disp('IMAS_VERSION undefined, assuming cocos_out_default=17');
end
% initialize input parser
diff --git a/matlab/TCV_IMAS/tcv_get_ids_nbi.m b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
index 6eae1073..a9baab2b 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_nbi.m
+++ b/matlab/TCV_IMAS/tcv_get_ids_nbi.m
@@ -5,7 +5,6 @@ function [ids_nbi,ids_nbi_description,varargout] = tcv_get_ids_nbi(shot,ids_nbi_
%
% gdat_params: gdat_data.gdat_params to get all params passed from original call, in particular error_bar options
%
-%
% Most information taken from NBI and NB_Model wiki pages:
% https://spcwiki.epfl.ch/wiki/NBI
% https://spcwiki.epfl.ch/wiki/NB_Model
--
GitLab