From 5d081669343ce47d83bb5ba4b19e91ec269d9155 Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Wed, 8 Jun 2016 10:11:13 +0000
Subject: [PATCH] fix further grids_1d for 1d cases
git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@5901 d63d8f72-b253-0410-a779-e742ad2e26cf
---
crpptbx/TCV/gdat_tcv.m | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index 32ee5b0e..729b45ef 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
--
GitLab