From c1b7920fce4888ed54ad5e899fd5ea8b19dc46df Mon Sep 17 00:00:00 2001
From: Olivier Sauter <Olivier.Sauter@epfl.ch>
Date: Sat, 2 May 2020 14:34:23 +0200
Subject: [PATCH] rm blanks then same as CXRS analysis
---
matlab/TCV/CXRS_fit_profies.m | 16 ++++++++--------
1 file changed, 8 insertions(+), 8 deletions(-)
diff --git a/matlab/TCV/CXRS_fit_profies.m b/matlab/TCV/CXRS_fit_profies.m
index 26c1c32e..6244844a 100644
--- a/matlab/TCV/CXRS_fit_profies.m
+++ b/matlab/TCV/CXRS_fit_profies.m
@@ -6,15 +6,15 @@ function profiles=CXRS_fit_profies(V,Times,param_prof,param_plot);
sel=find((V.t>=Times(1))&(V.t<=Times(2)));
x = V.x(sel); dx = V.dx(sel); y = V.y(sel); dy = V.dy(sel); t = V.t(sel); dt = V.dt(sel); sys = V.sys(sel);
-% shift of duplicated points in rho
+% shift of duplicated points in rho
Xb=x;
for ip=1:length(Xb)
sss=1:length(Xb); sss=find(sss~=ip);
iff=(find(Xb(ip)==Xb(sss)));
if length(iff)>0, Xb(iff(1))=Xb(iff(1))+1.e-8; end
-end
+end
% sort points in rho
-[x,iii] = sort(Xb);
+[x,iii] = sort(Xb);
if length(iii)==0 % No data
profiles.t_lim=Times;
profiles.x_prof=[];
@@ -32,12 +32,12 @@ if length(iii)==0 % No data
end
return;
end
-dx=dx(iii); y=y(iii); dy=dy(iii); t=t(iii); dt=dt(iii);sys=sys(iii);
-if isempty(param_prof.xout), param_prof.xout=linspace(min(x)*0.98,max(x)*1.02,101); end
-if ((min(x)-min(param_prof.xout))>0.45),
+dx=dx(iii); y=y(iii); dy=dy(iii); t=t(iii); dt=dt(iii);sys=sys(iii);
+if isempty(param_prof.xout), param_prof.xout=linspace(min(x)*0.98,max(x)*1.02,101); end
+if ((min(x)-min(param_prof.xout))>0.45),
param_prof.xout=linspace(min(x)*0.85,max(param_prof.xout),length(param_prof.xout));
- %disp(sprintf('Redefine rho for %4.4f-%4.4fsec',Times));
-end
+ %disp(sprintf('Redefine rho for %4.4f-%4.4fsec',Times));
+end
x_prof = param_prof.xout;
XXmax=max([max(x)+0.01 1.3]);
if length(iii)>=10,
--
GitLab