From fd4533df6d94a55024c38a280be5dd63282f6f57 Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Wed, 2 Dec 2015 14:47:54 +0000
Subject: [PATCH] add some n1 n2 signals from new magnetics
git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@5309 d63d8f72-b253-0410-a779-e742ad2e26cf
---
crpptbx/TCV/gdat_tcv.m | 16 +++++++++++++---
1 file changed, 13 insertions(+), 3 deletions(-)
diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index d0cd8ca8..8998ecd9 100644
--- a/crpptbx/TCV/gdat_tcv.m
+++ b/crpptbx/TCV/gdat_tcv.m
@@ -791,9 +791,19 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
case {'mhd'}
% load n=1, 2 and 3 Bdot from magnetic measurements
- n1=tdi('abs(mhdmode("LFS",1,1))');
- n2=tdi('abs(mhdmode("LFS",2,1))');
- n3=tdi('abs(mhdmode("LFS",3,1))');
+ if shot< 50926
+ n1=tdi('abs(mhdmode("LFS",1,1))');
+ n2=tdi('abs(mhdmode("LFS",2,1))');
+ n3=tdi('abs(mhdmode("LFS",3,1))');
+ else
+ aaLFSz23_sect3=tdi('\atlas::DT196_MHD_001:channel_067');
+ aaHFSz23_sect3=tdi('\atlas::DT196_MHD_002:channel_018');
+ n1 = aaLFSz23_sect3;
+ n1.data = aaLFSz23_sect3.data - aaHFSz23_sect3.data;
+ n2 = aaLFSz23_sect3;
+ n2.data = aaLFSz23_sect3.data + aaHFSz23_sect3.data;
+ n3=n1;
+ end
if ~isempty(n1.data)
gdat_data.data(:,1) = reshape(n1.data,length(n1.data),1);
if length(n2.data)==length(n1.data); gdat_data.data(:,2) = reshape(n2.data,length(n2.data),1); end
--
GitLab