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

@@ -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