diff --git a/matlab/TCV_IMAS/tcv2ids.m b/matlab/TCV_IMAS/tcv2ids.m
index da618f38938e129e92e8816db3f204cc0a92b980..5e5f717dd1c139d4930a64c850147e9fd3d94a21 100644
--- a/matlab/TCV_IMAS/tcv2ids.m
+++ b/matlab/TCV_IMAS/tcv2ids.m
@@ -23,6 +23,7 @@ function [ids_from_tcv,varargout] = tcv2ids(shot,varargin);
 %           'time_out': if 2 values provided: get all time data within that time interval
 %                       otherwise get values at these times provided in time_out (with linear interpolation and cst extrapolation)
 %           'trialindx': trial_indx for relevant nodes, in particular CONF kinetic profiles nodes
+%           'liuqe': liuqe indx for equilibrium, only 1,2,3 are available
 %
 % Outputs
 % varargout{1}: return also the ids in array of structure with the names, to allow easy use of plotallids
@@ -42,6 +43,7 @@ p.addOptional('b0sign_out', 0, @(x) isempty(x) || (isscalar(x) && (x==0 || x==-1
 p.addOptional('nverbose', 1, @(x) isempty(x) || isnumeric(x) );
 p.addOptional('time_out', [], @(x) isempty(x) || isnumeric(x) );
 p.addOptional('trialindx', [], @(x) isempty(x) || isnumeric(x) );
+p.addOptional('liuqe', [], @(x) isempty(x) || isnumeric(x) );
 
 params_not_for_gdat_params = {'shot','ids_names'};
 
@@ -129,4 +131,4 @@ if nargout>=2
     varargout{1} = [];
     disp('problems to fill in varargout')
   end
-end
+end
\ No newline at end of file