From 2a9d08a61a5f9a86f2de31b080cdad397a771f95 Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Mon, 13 May 2019 07:53:32 +0000
Subject: [PATCH] add min_power_on in input and change default to 2e3
git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@11873 d63d8f72-b253-0410-a779-e742ad2e26cf
---
crpptbx/AUG/read_ech_py.m | 7 ++++++-
1 file changed, 6 insertions(+), 1 deletion(-)
diff --git a/crpptbx/AUG/read_ech_py.m b/crpptbx/AUG/read_ech_py.m
index 5945b5db..56018c79 100644
--- a/crpptbx/AUG/read_ech_py.m
+++ b/crpptbx/AUG/read_ech_py.m
@@ -6,6 +6,7 @@ function [torbeam] = read_ech_py(shot,varargin);
% transforms ECS/ECN into theta,phi of Torbeam
%
% varargin{1}: time array for Torbeam relevant data output
+% varargin{2}: min_power_on to decide if the beam was on at a given time (default P>=2e3)
%
torbeam.rfmod = - ones(1,8); % 1 for O-mode, -1 for X-mode
torbeam.tb_par.('xxb') = [238., 238., 231.1, 231.1, 236.4, 236.4, 236.4 ,236.4];
@@ -154,7 +155,11 @@ for igy=1:ecsystems(isys)
end
end
-min_power_on = 10;
+if nargin>=3 && ~isempty(varargin{2})
+ min_power_on = varargin{2};
+else
+ min_power_on = 2e3;
+end
torbeam.beam_on = false.*ones(size(torbeam.power_t));
torbeam.beam_on(torbeam.power_t >= min_power_on) = true;
torbeam.t_some_beam_on = torbeam.t(any(torbeam.beam_on,1));
--
GitLab