diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index d0cd8ca826357c1fc26e7ff063bd93e9ff1fd0c5..8998ecd94f21a7467f937695ed1e541724054aa6 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