Skip to content
Snippets Groups Projects

Check cxrs

Merged Olivier Sauter requested to merge check_cxrs into master
1 file
+ 8
8
Compare changes
  • Side-by-side
  • Inline
@@ -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,
Loading