diff --git a/crpptbx/TCV/loadTCVdata.m b/crpptbx/TCV/loadTCVdata.m
index 4af7b51a6f41eacc12d4b0207a118e35bb58cbf9..c519945f7ed35850dd17e97678b517c7ac6af74e 100644
--- a/crpptbx/TCV/loadTCVdata.m
+++ b/crpptbx/TCV/loadTCVdata.m
@@ -1,4 +1,4 @@
- function [trace,error,varargout]=loadTCVdata(shot,data_type,varargin)
+ function [trace,error_status,varargout]=loadTCVdata(shot,data_type,varargin)
 %
 % Added option to load shot=-1 or >=100000
 % Added option shot=-9 to list keywords
@@ -89,12 +89,12 @@
 % function needed: mds functions-xtomo_geometry-get_xtomo_data (Furno's routines)
 %                  VsxrTCVradius
 % Example:
-%         [ip,error]=loadTCVdata(shot,'Ip',1);
-%         [sxr,error,R]=loadTCVdata(shot,'sxR',1);
+%         [ip,error_status]=loadTCVdata(shot,'Ip',1);
+%         [sxr,error_status,R]=loadTCVdata(shot,'sxR',1);
 %
 
 varargout{1}=cell(1,1);
-error=1;
+error_status=1;
 
 % To allow multiple ways of writing a specific keyword, use data_type_eff within this routine
 if exist('data_type') && ~isempty(data_type)
@@ -502,7 +502,7 @@ switch TCVkeywrdcase{index}
     else
       mdsopen(shot);
       % test if node exists
-      error=1;
+      error_status=1;
       ij=findstr(TCVsiglocation{index},'[');
       if isempty(ij); ij=length(TCVsiglocation{index})+1; end
       if eval(['~mdsdata(''node_exists("\' TCVsiglocation{index}(1:ij-1) '")'')'])
@@ -614,7 +614,7 @@ switch TCVkeywrdcase{index}
     end
     
     trace.name=[num2str(shot) ';' nodenameeff];
-    error=0;
+    error_status=0;
 
   %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
   case {'ne','te'}
@@ -980,7 +980,7 @@ switch TCVkeywrdcase{index}
   case {'profnerho','profterho'}
     % vol from psitbx
     mdsopen(shot);
-    error=1;
+    error_status=1;
     if strcmp(TCVkeywrdcase{index},'profnerho')
       nodenameeff=['\results::THOMSON.PROFILES.AUTO:ne'];
       avers=tdi('\results::THOMSON.PROFILES.AUTO:ne:version_num');
@@ -1001,9 +1001,9 @@ switch TCVkeywrdcase{index}
         if ~isempty(tracetdi.dim) && ~ischar(tracetdi.data)
           trace.x=tracetdi.dim{1};
           trace.t=tracetdi.dim{2};
-          error=0;
+          error_status=0;
         else
-         error=2;
+         error_status=2;
           trace.x=[];
           trace.t=[];
         end
@@ -1013,11 +1013,11 @@ switch TCVkeywrdcase{index}
           disp('assumes dim{2} for x in THOMSON.PROFILES.AUTO')
           trace.x=tracetdi.dim{2};
           trace.t=tracetdi.dim{1};
-          error=0;
+          error_status=0;
         else
           trace.x=[];
           trace.t=[];
-         error=2;
+         error_status=2;
         end
       end
     else
@@ -1211,7 +1211,7 @@ switch TCVkeywrdcase{index}
     trace.dim=[{trace.x} ; {trace.t}];
     trace.dimunits=[{'channel #'} ; {'time [s]'}];
     trace.name=[num2str(shot) ';' 'get_xtomo_data'];
-    error=0;
+    error_status=0;
 
   %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
   case 'ece'
@@ -1242,7 +1242,7 @@ switch TCVkeywrdcase{index}
     radius.data=trace.R;
     radius.t=trace.t;
     varargout{1}={radius};
-    error=0;
+    error_status=0;
     % isfield does not work since tracetdi is not a 'struct' but a tdi object, thus isfield using isa does not work
     if any(strcmp(fieldnames(tracetdi),'units'))
       trace.units=tracetdi.units;
@@ -1273,7 +1273,7 @@ switch TCVkeywrdcase{index}
     trace.name=[num2str(shot) ';' 'get_mds_mio(MPX)'];
     [xchord,ychord]=mpx_geometry;
     varargout{1}={VsxrTCVradius(zmag.data,xchord,ychord)};
-    error=0;
+    error_status=0;
 
   %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
   case 'IOH'
@@ -1292,7 +1292,7 @@ switch TCVkeywrdcase{index}
     end
     trace.name=[num2str(shot) ';' nodenameeff{1} ',' nodenameeff{2}];
     mdsclose;
-    error=0;
+    error_status=0;
     
   %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
   case 'vloop'
@@ -1311,7 +1311,7 @@ switch TCVkeywrdcase{index}
     end
     trace.name=[num2str(shot) ';' nodenameeff{1} ',' nodenameeff{2}];
     mdsclose;
-    error=0;
+    error_status=0;
     
   %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
   case 'pgyro'
@@ -1330,7 +1330,7 @@ switch TCVkeywrdcase{index}
     end
     trace.name=[num2str(shot) ';' nodenameeff];
     mdsclose;
-    error=0;
+    error_status=0;
     
   %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
   case 'jtor'
@@ -1349,7 +1349,7 @@ switch TCVkeywrdcase{index}
     end
     trace.name=[num2str(shot) ';' nodenameeff];
     mdsclose;
-    error=0;
+    error_status=0;
     
   %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
   case {'neft:trial','teft:trial','neftav:trial','teftav:trial'}
@@ -1371,7 +1371,7 @@ switch TCVkeywrdcase{index}
     end
     trace.name=[num2str(shot) ' ; ' nodenameeff{:} ' ; trialindx=' num2str(trialindx) ];
     mdsclose;
-    error=0;
+    error_status=0;
     
   %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
  case {'vi_tor', 'vi_torfit', 'vi_pol', 'vi_polfit', 'Ti', 'Tifit', 'ni', 'nifit', 'zeffcxrs', 'zeffcxrsfit'}
@@ -1414,7 +1414,7 @@ switch TCVkeywrdcase{index}
       trace.std_t = [];
     end
     mdsclose;
-    error=0;
+    error_status=0;
     
   %&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
   otherwise