diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index c4e39ca04ea6257a63e3e31af7062098294c5a6b..5623d3e25c9250edadde1d2c5e04a34602d36739 100644
--- a/crpptbx/TCV/gdat_tcv.m
+++ b/crpptbx/TCV/gdat_tcv.m
@@ -1148,7 +1148,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
         itpsitdi=iround_os(psitdi.dim{3},gdat_data.t(it));
         psirz=psitdi.data(:,:,itpsitdi);
         %psiscatvol0=interp2(rmesh,zmesh,psirz',0.9*ones(size(zthom)),zthom-zeffshift(itpsitdi),'spline'); % faster with interpos
-        psiscatvol0=interpos2Dcartesian(rmesh,zmesh,psirz,0.9*ones(size(zthom)),zthom-zeffshift(itpsitdi),-0.1,-0.1);
+        %psiscatvol0=interpos2Dcartesian(rmesh,zmesh,psirz,0.9*ones(size(zthom)),zthom-zeffshift(itpsitdi),-0.1,-0.1);
         psiscatvol0=interp2(rmesh,zmesh,psirz',0.9*ones(size(zthom)),zthom-zeffshift(itpsitdi),'linear');
         psiscatvol.data(it,:)=psiscatvol0;
         % since take closest psi(R,Z) from psitdi, should also take closest for psi_max and not interpolating