From 5809481873ec392f8f869056612f20d64ab75bca Mon Sep 17 00:00:00 2001
From: Olivier Sauter <Olivier.Sauter@epfl.ch>
Date: Wed, 18 Sep 2019 17:31:05 +0200
Subject: [PATCH] add ids ec_antennas
---
matlab/TCV_IMAS/tcv_get_ids_ec_antennas.m | 76 +++++++++++++++++++++++
1 file changed, 76 insertions(+)
create mode 100644 matlab/TCV_IMAS/tcv_get_ids_ec_antennas.m
diff --git a/matlab/TCV_IMAS/tcv_get_ids_ec_antennas.m b/matlab/TCV_IMAS/tcv_get_ids_ec_antennas.m
new file mode 100644
index 00000000..06d38aa4
--- /dev/null
+++ b/matlab/TCV_IMAS/tcv_get_ids_ec_antennas.m
@@ -0,0 +1,76 @@
+function [ids_ec_antennas,ids_ec_antennas_description,varargout] = tcv_get_ids_ec_antennas(shot,ids_ec_antennas_empty, gdat_params,varargin);
+%
+% [ids_ec_antennas,ids_ec_antennas_description,varargout] = tcv_get_ids_ec_antennas(shot,ids_ec_antennas_empty,varargin);
+%
+%
+% gdat_params: gdat_data.gdat_params to get all params passed from original call, in particular error_bar options
+%
+
+[ids_ec_antennas, params_ec_antennas] = tcv_ids_headpart(shot,ids_ec_antennas_empty,'ec_antennas','homogeneous_time',0,varargin{:});
+
+% As a general rule, for a new substructure under the main ids, construct a local structure like:
+% "global_quantities" with subfields being the relevant data to get and a local structure:
+% "global_quantities_desc" which contains the same subfields themselves containing the gdat string aftre shot used
+%
+
+pow=gdat_tcv(shot,'powers');
+
+if shot>=63000
+ warning('not sure if set-up correct for new launchers, beams, gyrotrons. Check with O. Sauter')
+ return
+end
+
+nb_antennas = 9; % at this stage keep same numbering as standard TCV launchers before mid-2019
+
+ids_ec_antennas.antenna(1:nb_antennas) = ids_ec_antennas.antenna(1); % copy empty structure for all units, then fill in
+% $$$ ids_ec_antennas.time = pow.ec.t; if homogeneous time
+% $$$ ids_ec_antennas_description.time = 'time from power';
+[r_phi0,z_phi0,r1,z1,rc,zc,thetphi_L,thetphi_tor,raysbylauncher,powers,fname,z_axis,launch_params] = ...
+ toray_raygeom_TCV(shot,pow.ec.t(1:10:end));
+
+for iant=1:nb_antennas
+ ids_ec_antennas.antenna{iant}.identifier = ['L' num2str(iant)];
+ switch iant
+ case {1,4}
+ ids_ec_antennas.antenna{iant}.name = 'LFS equatorial launcher';
+ case {2, 3, 5, 6}
+ ids_ec_antennas.antenna{iant}.name = 'LFS top launcher';
+ otherwise
+ ids_ec_antennas.antenna{iant}.name = 'TOP launcher';
+ end
+ if any(isfinite(pow.ec.data(:,iant)))
+ ids_ec_antennas.antenna{iant}.power_launched.data = pow.ec.data(:,iant);
+ ids_ec_antennas.antenna{iant}.power_launched.time = pow.ec.t;
+ ids_ec_antennas_description.antenna{iant}.power_launched = 'from gdat powers .ec';
+ ids_ec_antennas_description.antenna{iant}.launching_position = 'from launch_params obtained from toray_raygeom_TCV';
+ else
+ ids_ec_antennas_description.antenna{iant}.power_launched = 'no power for this launcher from gdat powers .ec';
+ end
+ if ~isempty(launch_params{iant})
+ % find 1st non-empty data
+ it_ok{iant} = [];
+ for it=1:length(launch_params{iant})
+ if ~isempty(launch_params{iant}{it})
+ it_ok{iant}(end+1) = it;
+ end
+ end
+ % non time-dependent quantities, take 1st ok values
+ ids_ec_antennas.antenna{iant}.frequency =launch_params{iant}{it_ok{iant}(1)}.freq;
+ ids_ec_antennas.antenna{iant}.mode.time = [pow.ec.t(1) pow.ec.t(end)];
+ ids_ec_antennas.antenna{iant}.mode.data = [-1 -1]; % at this stage assume X mode always, to change when available
+ for i=1:length(it_ok{iant})
+ r0 = sqrt(launch_params{iant}{it_ok{iant}(i)}.x0.^2 + launch_params{iant}{it_ok{iant}(i)}.y0.^2) / 100.; % in [m]
+ time_launch = launch_params{iant}{it_ok{iant}(i)}.time;
+ ids_ec_antennas.antenna{iant}.launching_position.r.data(i) = r0;
+ ids_ec_antennas.antenna{iant}.launching_position.z.data(i) = launch_params{iant}{it_ok{iant}(i)}.z0/100.;
+ ids_ec_antennas.antenna{iant}.launching_position.phi.data(i) = atan2(launch_params{iant}{it_ok{iant}(i)}.y0/100,r0);
+ ids_ec_antennas.antenna{iant}.launching_position.r.time(i) = time_launch;
+ ids_ec_antennas.antenna{iant}.launching_position.z.time(i) = time_launch;
+ ids_ec_antennas.antenna{iant}.launching_position.phi.time(i) = time_launch;
+ ids_ec_antennas.antenna{iant}.launching_angle_pol.data(i) = launch_params{iant}{it_ok{iant}(i)}.theta_toray * pi/180.;
+ ids_ec_antennas.antenna{iant}.launching_angle_pol.time(i) = time_launch;
+ ids_ec_antennas.antenna{iant}.launching_angle_tor.data(i) = launch_params{iant}{it_ok{iant}(i)}.phi_toray * pi/180.;
+ ids_ec_antennas.antenna{iant}.launching_angle_tor.time(i) = time_launch;
+ end
+ end
+end
--
GitLab