diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index 32ee5b0ec47589824fcdb6cf342b9fc7ea7602a6..729b45ef00e3d8430afb94354a8d88fbd821476c 100644
--- a/crpptbx/TCV/gdat_tcv.m
+++ b/crpptbx/TCV/gdat_tcv.m
@@ -1864,11 +1864,13 @@ for it=1:length(gdat_data.t)
   else
     ii=find(~isnan(gdat_data.grids_1d.rhopolnorm(:,it)));
   end
-  if length(ii)==length(gdat_data.grids_1d.rhopolnorm)
-    if (nbdim_x == 1)
+  if (nbdim_x == 1)
+    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);
       gdat_data.grids_1d.rhovolnorm(:,it)=interpos(rhovol.x,rhovol.data(:,it_vol_eff),gdat_data.grids_1d.rhopolnorm);
-    else
+    end
+  else
+    if length(ii)==size(gdat_data.grids_1d.rhopolnorm,1)
       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