%session 30-11-04

cd E:\temp
dir
load p1lg1THERMOall_sansentete.txt
yd=p1lg1THERMOall_sansentete;
plot(yd(:,1))
plot(diff(yd(:,1)),'.')
plot(diff(yd(:,1))*1440,'.')
diff(yd(4154:4155,1))*1440*10000
diff(yd(4154:4155,1))*1440*100
diff(yd(4153:4154,1))*1440*100
load p2lg1THERMOall_sansentete.txt
yd2=p2lg1THERMOall_sansentete;
load p3lg1THERMOall_sansentete.txt
yd3=p3lg1THERMOall_sansentete;
figure
plot(diff(yd2(:,1))*1440,'.')
figure
plot(diff(yd3(:,1))*1440,'.')
1/1440
load P1L1_carbinorgdissous_jourjulien_170703.txt
zd=P1L1_carbinorgdissous_jourjulien_170703;
figure
plot(zd(:,21),'.
plot(zd(:,21),'.')
title('cinorg')
figure
plot(zd(:,8),'.')
help julien
julien(02,14,00,01,05)
julien(02,14,24,01,05)
indh=find(zd(:,8)==24);
zdt=zd;
zdt(indh,21)=zdt(indh,21)-1;
figure
plot(zdt(:,21),'.')
plot(diff(zdt(:,21))'.')
inddi=find(diff(zdt(:,21))<1/1440);
min(diff(zdt(:,21)))
min(diff(zdt(:,21)))*1440
for k=1:2651
end
nli=ones(2651,1)*NaN;
for k=1:2651
nli(k)=find(abs(yd(:,1)-zdt(k,21))<0.5/1440);
end
k
nli(k)
find(abs(yd(:,1)-zdt(k,21))<0.5/1440);
find(abs(yd(:,1)-zdt(k,21))<0.5/1440)
disp(yd(1715,1),yd(1716,1))
disp([yd(1715,1),yd(1716,1)])
ztd(k,21)
zdt(k,21)
zdt(k,21)-yd(1715,1)
zdt(k,21)-yd(1716,1)
(zdt(k,21)-yd(1716,1))*1440
(zdt(k,21)-yd(1716,1))*1440
(zdt(k,21)-yd(1715,1))*1440
[minm,indm]=find(abs(yd(:,1)-zdt(k,21))<0.5/1440)
indm=find(abs(yd(:,1)-zdt(k,21))<0.5/1440)
[mi,indi]=find(abs(yd(indm,1)-zdt(k,21))
[mi,indi]=find(abs(yd(indm,1)-zdt(k,21)))
[mi,indi]=min(abs(yd(indm,1)-zdt(k,21)))
nli=ones(2651,1)*NaN;
for k=1:2651
indm=find(abs(yd(:,1)-zdt(k,21))<0.5/1440);
if(length(indm)>1;[mi,indi]=min(abs(yd(indm,1)-zdt(k,21)));nli(k)=indm(indi);end
nli=ones(2651,1)*NaN;
for k=1:2651
indm=find(abs(yd(:,1)-zdt(k,21))<0.5/1440);
if length(indm)>1;[mi,indi]=min(abs(yd(indm,1)-zdt(k,21)));nli(k)=indm(indi);end
end
nomb=find(isnan(nli))
nomb=find(isnan(nli));
nli=ones(2651,1)*NaN;
for k=1:2651
indm=find(abs(yd(:,1)-zdt(k,21))<0.5/1440);
if length(indm)>1;[mi,indi]=min(abs(yd(indm,1)-zdt(k,21)));nli(k)=indm(indi);end
if length(indm)==1;nli(k)=indm;end
end
nomb=find(isnan(nli));
nomb
zdt(70,21)
nli(70)
nli(69))
nli(69)
nli(71)
ydt=round(yd*1440);
figuree
figure
plot(ydt;'.')
plot(ydt,'.')
ydt=round(yd(:,1)*1440);
plot(ydt,'.')
plot(diff(ydt),'.')
nli=ones(2651,1)*NaN;
for k=1:2651
nli(k)=find(abs(ydt(:,1)-zdt(k,21))<0.5/1440);
end
k
nli=ones(2651,1)*NaN;
for k=1:2651
nli(k)=find(abs(ydt-zdt(k,21))<0.5/1440);
end
find(abs(ydt-zdt(1,21))<0.5/1440)
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21))<0.5/1440);
if length(indo==1);nli(k)=indo;end
end
k
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21))<0.5/1440);
if length(indo>0);nli(k)=indo;end
end
nomb=find(isnan(nli));
help indo
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<0.5/1440);
if length(indo>0);nli(k)=indo;end
end
nomb
nomb=find(isnan(nli));
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<0.5);
if length(indo>0);nli(k)=indo;end
end
nomb=find(isnan(nli));
nomb
zdt(305,21)
zdt(305,21)*1440
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<0.5);
if length(indo>0);nli(k)=indo;end
p
end
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<0.5);
if length(indo>1);nli(k)=indo(1);end
if length(indo>0);nli(k)=indo;end
end
nomb=find(isnan(nli));
indo=find(abs(ydt-zdt(305,21)*1440)<0.5);
indo
indo=find(abs(ydt-zdt(306,21)*1440)<0.5);
indo
indo=find(abs(ydt-zdt(304,21)*1440)<0.5);
indo
37.2281-37.22778
(37.2281-37.22778)*1440
find(yd(:,1)==37.22778)
ydt(4798,1)
ydt(4798,1)-37.2281*1440
ydt(4799,1)-37.2281*1440
ydt(4797,1)-37.2281*1440
ydt(4798,1)-zdt(305,21)*1440
ydt(4798,1)-zdt(305,21)*1440+0.5
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<=0.5);
if length(indo>1);nli(k)=indo(1);end
if length(indo>0);nli(k)=indo;end
end
k
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<=0.5);
if length(indo)>1;nli(k)=indo(1);end
if length(indo)>0;nli(k)=indo;end
end
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<0.5);
if length(indo)>1;nli(k)=indo(1);end
if length(indo)>0;nli(k)=indo;end
end
nomb=find(isnan(nli));
length(indo)
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<=0.5);
if length(indo)>1;nli(k)=indo(1);end
if length(indo)>0;nli(k)=indo;end
end
length(indo)
nli(k)=indo(1);
if length(indo)>1;nli(k)=indo(1);end
if length(indo)>0;nli(k)=indo;end
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<=0.5);
if length(indo)>1;nli(k)=indo(1);end
if length(indo)==1;nli(k)=indo;end
end
nomb=find(isnan(nli));
nli=ones(2651,1)*NaN;
for k=1:2651
indo=find(abs(ydt-zdt(k,21)*1440)<=0.5);
%if length(indo)>1;nli(k)=indo(1);end
if length(indo)>0;nli(k)=indo(1);end
end
nomb=find(isnan(nli));