fix q_rho again

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@5871 d63d8f72-b253-0410-a779-e742ad2e26cf

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);