diff --git a/matlab/TCV/CXRS_fit_profies.m b/matlab/TCV/CXRS_fit_profies.m index 26c1c32ec0b885ac839a70a191d9342b85a39274..6244844aba40cde3ad5b2349925ed06011bf8dfa 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,