diff --git a/gdat.m b/gdat.m
index 202bb029384b2827fac2f94756936a079c92e8c4..17c3a4eacda5999dc235495f06a8591c66d67c2d 100644
--- a/gdat.m
+++ b/gdat.m
@@ -11,8 +11,26 @@ function [trace,error,varargout] = gdat(shot,data_type,varargin)
 % 'sxR'  =  soft x-ray emission with varargout{1} option (requires varargin{5}!)
 %
 % specific to TCV (see help loadTCVdata for more information)
-% 'MPX'  =  wire chambers with varargout{1} option (requires varargin{5}!)
 % 'xx_2 or xx_3' for Liuqe2 or 3: same as above for xx related to equilibrium
+% 'rcont' =  R of plama flux surfaces
+% 'zcont' =  Z of plama flux surfaces
+% 'vol' =  volume of flux surfaces
+% 'qrho' =  q profile on rho mesh
+% 'neint' =  line-integrated electron density [m/m^3]
+% 'ne'= ne raw profile on (z,t). ADD error bars in .std
+% 'te'= Te raw profile on (z,t). ADD error bars in .std
+% 'nerho'= ne profile on (rho=sqrt(psi),time) mesh.Note rho is a 2D array as depends on time. ADD error bars in .std
+% 'terho'= Te profile on (rho=sqrt(psi),time) mesh.Note rho is a 2D array as depends on time. ADD error bars in .std
+% 'profnerho' =  ne smoothed or fitted , vs (rho,t) (from Thomson fit)
+% 'profterho' =  te smoothed or fitted , vs (rho,t) (from Thomson fit)
+% 'neft' =  ne fitted from data on rho mesh (from proffit.local_time:neft)
+% 'teft' =  te fitted from data on rho mesh (from proffit.local_time:teft)
+% 'neftav' =  ne fitted from averaged over time data on rho mesh (from proffit.avg_time:neft)
+% 'teftav' =  te fitted from averaged over time data on rho mesh (from proffit.avg_time:teft)
+% 'ece'  =  electron cyclotron emission
+% 'sxr'  =  soft x-ray emission
+% 'sxR'  =  soft x-ray emission with varargout{1} option (requires varargin{4}!)
+% 'MPX'  =  soft x-ray from wire chambers
 %
 % JET
 % Special case compatible with old gdat.m allows (JET related): gdat(51994,'ppf','efit/xip',1)