diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index 6360949d4ff46c344a450dc23d7e6d66b8378142..ecebc95f05cc64893563dd44e0f33f3e04e66ba5 100644
--- a/crpptbx/TCV/gdat_tcv.m
+++ b/crpptbx/TCV/gdat_tcv.m
@@ -1856,10 +1856,19 @@ for it=1:length(gdat_data.t)
   % do an interpolation on closest point to avoid 2D interp
   it_rt_eff = it_rt(it);
   it_vol_eff = it_vol(it);
-  ii=find(~isnan(gdat_data.grids_1d.rhopolnorm(:,it)));
+  if (nbdim_x == 1)
+    ii=find(~isnan(gdat_data.grids_1d.rhopolnorm));
+  else
+    ii=find(~isnan(gdat_data.grids_1d.rhopolnorm(:,it)));
+  end
   if length(ii)==length(gdat_data.grids_1d.rhopolnorm)
-    gdat_data.grids_1d.rhotornorm(:,it)=interpos(rhotor.x,rhotor.data(:,it_rt_eff),gdat_data.grids_1d.rhopolnorm(:,it));
-    gdat_data.grids_1d.rhovolnorm(:,it)=interpos(rhovol.x,rhovol.data(:,it_vol_eff),gdat_data.grids_1d.rhopolnorm(:,it));
+    if (nbdim_x == 1)
+      gdat_data.grids_1d.rhotornorm(:,it)=interpos(rhotor.x,rhotor.data(:,it_rt_eff),gdat_data.grids_1d.rhopolnorm);
+      gdat_data.grids_1d.rhovolnorm(:,it)=interpos(rhovol.x,rhovol.data(:,it_vol_eff),gdat_data.grids_1d.rhopolnorm);
+    else
+      gdat_data.grids_1d.rhotornorm(:,it)=interpos(rhotor.x,rhotor.data(:,it_rt_eff),gdat_data.grids_1d.rhopolnorm(:,it));
+      gdat_data.grids_1d.rhovolnorm(:,it)=interpos(rhovol.x,rhovol.data(:,it_vol_eff),gdat_data.grids_1d.rhopolnorm(:,it));
+    end
   end
 end
 gdat_data.grids_1d.rhotor_edge=interpos(rhotor.t',rhotor.rhotor_edge,gdat_data.t',-0.01);