diff --git a/matlab/IMAS/gdat_imas.m b/matlab/IMAS/gdat_imas.m
index 88b72b2f647f2f8ca0bb8b09601439ac62e11e9d..d2b06fefcbdf8d2fcc5fc3b3289cc933cdc28ddd 100644
--- a/matlab/IMAS/gdat_imas.m
+++ b/matlab/IMAS/gdat_imas.m
@@ -633,7 +633,7 @@ end
 function [gdat_data] = get_thomson_raw_data(shot,data_request_eff,gdat_data,doedge,nverbose);
 %
 try
-  time=mdsdata('\results::thomson:times');
+  time=mdsvalue('\results::thomson:times');
 catch
   gdat_data.error_bar = [];
   if strcmp(data_request_eff(1:2),'ne')
@@ -695,7 +695,7 @@ if strcmp(data_request_eff(1:2),'ne')
     disp('***********************************************************************')
   end
 end
-z=mdsdata(['\diagz::thomson_set_up' edge_str_dot ':vertical_pos']);
+z=mdsvalue(['\diagz::thomson_set_up' edge_str_dot ':vertical_pos']);
 gdat_data.dim=[{z};{time}];
 gdat_data.dimunits=[{'Z [m]'} ; {'time [s]'}];
 gdat_data.x=z;
@@ -716,7 +716,7 @@ function [psiscatvol,psi_max] = psi_scatvol_zshift(gdat_data, zshift, system_str
 
 % Use psitbx
 t_th = gdat_data.t;
-th_z = mdsdata(['dim_of(\thomson' system_str_dot ':te,1)']); % TODO: Isn't this gdat_data.dim{1}?
+th_z = mdsvalue(['dim_of(\thomson' system_str_dot ':te,1)']); % TODO: Isn't this gdat_data.dim{1}?
 th_r = 0.9*ones(size(th_z));
 ntt = numel(t_th);
 nch = numel(th_z); % number of chords