diff --git a/matlab/TCV/gdat_tcv.m b/matlab/TCV/gdat_tcv.m
index af9247e6d4ec348f3123ec1e14aef1e070e88a58..2101495af9bb379adf316dfc803a32e5cc7adc8a 100644
--- a/matlab/TCV/gdat_tcv.m
+++ b/matlab/TCV/gdat_tcv.m
@@ -2012,8 +2012,9 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
       r_edge = gdat(gdat_data.shot,params_eff);
       params_eff.data_request='z_contour_edge';
       z_edge = gdat(gdat_data.shot,params_eff);
-      R_thomson = gdat_tcv(66540,'\tcv_shot::top.diagz.measurements.thomson.set_up:radial_pos');
       Z_thomson = psiscatvol.dim{1};
+      %OS: not ok for old shots: R_thomson = gdat_tcv(gdat_data.shot,'\tcv_shot::top.diagz.measurements.thomson.set_up:radial_pos');
+      R_thomson.data = 0.90*ones(size(Z_thomson));
       psinorm_max = max(max(psi_norm));
       for it=1:numel(psiscatvol.dim{2})
         iteq=iround_os(r_edge.t,psiscatvol.dim{2}(it));