mdsconnect('tokamak-profiledb.ccfe.ac.uk');
[a,b]=mdsopen('pr98_jet',19649); % shot in PR98 PDB (freely accessible)
%[a,b]=mdsopen('jet',53521);     % shot in working PDB (registration required for access)
%[a,b]=mdsopen('itb_tftr',94607);% shot in ITB PDB (registration required for access)

% comment example
[y,status]=mdsvalue('\top.comments:shot');

% 0D variable example
[y,status]=mdsvalue('\top.zerod:zeff');

% 1D variable example
[ip,status]=mdsvalue('\top.oned:ip');      % gets the 1D variable
[tip,status]=mdsvalue('dim_of(\top.oned:ip)');    % gets time of the 1D
variable
[ip_unit,status]=mdsvalue('units_of(\top.oned:ip)');    % gets unit of the
1D variable
[tip_units,status]=mdsvalue('units_of(dim_of(\top.oned:ip))');    % gets
unit of time of the 1D variable

% 2D variable example
[q,status]=mdsvalue('\top.twod:q');      % gets the 2D variable
[tq,status]=mdsvalue('dim_of(\top.twod:q,1)');    % gets time of the 2D
variable
[rq,status]=mdsvalue('dim_of(\top.twod:q,0)');    % gets radius of the 2D
variable
[q_unit,status]=mdsvalue('units_of(\top.twod:q)');      % gets unit of the
2D variable
[tq_unit,status]=mdsvalue('units_of(dim_of(\top.twod:q,1))');    % gets
unit of time of the 2D variable
[rq_unit,status]=mdsvalue('units_of(dim_of(\top.twod:q,0))');    % gets
unit of radius of the 2D variable


% other instructions ...
%[y,status]=mdsvalue('getnci(".comments:*","nid_number")'); % list of all
signals (node numbers) available in the "comments" branch
%[y,status]=mdsvalue('getnci(".zerod:*","nid_number")')% % list of all
signals (node number) available in the 0D branch

%[y,status]=mdsvalue('getnci(123,"node_name")'); % gets node name