diff --git a/matlab/TCV/gdat_tcv.m b/matlab/TCV/gdat_tcv.m
index 4885f998e7e92569d0111ca435cb213e4b6f2b7a..564eebba09b91327ede9d50930b34cf729cf8ae8 100644
--- a/matlab/TCV/gdat_tcv.m
+++ b/matlab/TCV/gdat_tcv.m
@@ -907,6 +907,24 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
        gdat_data.Btor.x = {L.rx,L.zx};
        gdat_data.Btor.data_fullpath = datapath;
        gdat_data.Btor.label = 'Toroidal magnetic field map in (R,Z)';
+       
+       % total magnetic field
+       gdat_data.Btot.data = sqrt(gdat_data.Btor.data.^2 + gdat_data.Bz.data.^2 + gdat_data.Br.data.^2);
+       gdat_data.Btot.units = 'T';
+       gdat_data.Btot.dim = {L.rx,L.zx,LY.t};
+       gdat_data.Btot.dimunits = {'m','m','s'};
+       gdat_data.Btot.t = LY.t;
+       gdat_data.Btot.x = {L.rx,L.zx};
+       gdat_data.Btot.data_fullpath = datapath;
+       gdat_data.Btot.label = 'Total magnetic field map in (R,Z)';
+       
+       gdat_data.data = gdat_data.Btot.data;
+       gdat_data.units = 'T';
+       gdat_data.dim = {L.rx,L.zx,LY.t};
+       gdat_data.dimunits = {'m','m','s'};
+       gdat_data.t = gdat_data.Btot.t;
+       gdat_data.x = {L.rx,L.zx};
+       gdat_data.label = 'Total magnetic field map in (R,Z)';
 
      else
 
@@ -914,6 +932,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
        gdat_data.Btor = [];
        gdat_data.Br   = [];
        gdat_data.Bz   = [];
+       gdat_data.Btot = [];
 
      end