pro analogico_N2_scatter
device,decomposed=0
loadct,38
dove='d:\compatti\'

z=fltarr(2000)
for iz=1,2000 do z(iz-1)=iz
arra1=fltarr(2013,30)
arra2=fltarr(2013,30)
arra3=fltarr(816,30)
!p.multi=[0,3,2]

window,0,xsize=1000,ysize=700;,title=cchan(ich)
;plot,[0,0],[0,0],yrange=[1,20000],xrange=[1,20000],/xlog,/ylog,ystyle=1,color=0,background=255,xstyle=1,charsize=1.8,ytitle='D06',xtitle='D04'

;plot,arra1(32:62),arra2(32:62),color=0,psym=3,background=255 ;3
;nome='10\100824nocho'
nome='10\100914_satu'
openr,1,dove+nome+'D05.out',error=err1
openr,2,dove+nome+'D04.out',error=err2
openr,3,dove+nome+'A04.out',error=err2
readf,1,arra1,format='(i4,5i3,i4,f6.1,1x,f7.3,f4.1,3i5,2000i6)'
readf,2,arra2,format='(i4,5i3,i4,f6.1,1x,f7.3,f4.1,3i5,2000i6)'
readf,3,arra3

close,1
close,2
close,3

plot,[0,0],[0,0],color=0,psym=3,background=255,xrange=[1,10000],yrange=[1,10000],/xlog,/ylog,xtitle='D04',ytitle='D06'
oplot,arra1(13:31,*),arra2(13:31,*),color=0,psym=1
oplot,arra1(32:72,*),arra2(32:72,*),color=20,psym=1
oplot,arra1(72:112,*),arra2(72:112,*),color=40,psym=1
oplot,arra1(113:2012,*),arra2(113:2012,*),color=80,psym=1
oplot,arra1(20:27,*),arra2(20:27,*),color=120,psym=1

;plot,[0,0],[0,0],color=0,psym=3,background=255,xrange=[100,9000],yrange=[50,600],ystyle=1,xstyle=1,/xlog,/ylog,xtitle='D06',ytitle='A06',charsize=2
plot,[0,0],[0,0],color=0,psym=3,background=255,xrange=[100,10000],yrange=[50,600],ystyle=1,xstyle=1,xtitle='D06',ytitle='A06',charsize=2
oplot,arra2(13:23,*),arra3(16:26,*),color=0,psym=1
oplot,arra2(24:34,*),arra3(27:37,*),color=20,psym=1
oplot,arra2(35:45,*),arra3(38:48,*),color=40,psym=1
oplot,arra2(46:56,*),arra3(49:59,*),color=60,psym=3
oplot,arra2(57:67,*),arra3(60:70,*),color=80,psym=3
oplot,arra2(68:78,*),arra3(71:81,*),color=100,psym=3
oplot,arra2(79:89,*),arra3(82:92,*),color=120,psym=1
oplot,arra2(90:210,*),arra3(93:213,*),color=180,psym=1

fit=0
if fit eq 1 then begin
plot,[0,0],[0,0],color=0,psym=3,background=255,yrange=[1,14000],xrange=[50,1200],ystyle=1,xstyle=1,xtitle='D06',ytitle='A06',charsize=2
oplot,arra3(93:213,*),arra2(90:210,*),color=180,psym=1
oplot,arra3(16:815,*),arra2(13:812,*),color=0,psym=1
result=poly_fit(arra3(93:213,*),arra2(90:210,*),1)
print,result
b=-5404
a=105.15
x1=50
x2=90
y1=a*x1+b
y2=a*x2+b
oplot,[x1,x2],[y1,y2],color=0
b=-7097
a=133.779
x1=50
x2=90
y1=a*x1+b
y2=a*x2+b
oplot,[x1,x2],[y1,y2],color=20
endif
;x=fltarr(800*30)
;y=fltarr(800*30)
;ind=0
;for ip=0,799 do begin
;for jp=0,29 do begin
;ind=ind+1
;x(ind)=
;y=fltarr(800*30)
;endfor


plot,[0,0],[0,0],color=0,background=255,xrange=[10,10000],/xlog,xstyle=1,yrange=[1,100],ystyle=1
for it=0,29 do begin
;plot,arra1(13:112,it),z,color=0,background=255,xrange=[10,10000],/xlog,xstyle=1,yrange=[1,100],ystyle=1
oplot,arra1(13:2012,it),z,color=0
oplot,arra2(13:2012,it),z,color=20
oplot,10*(arra3(16:815,it)-50)-10,z(0:799),color=40
endfor
stop
arra1=fltarr(2013,200)
arra2=fltarr(2013,200)
arra3=fltarr(816,200)
plot,[0,0],[0,0],color=0,psym=3,background=255,xrange=[1,10000],yrange=[1,10000],/xlog,/ylog,xtitle='D04',ytitle='D06'

openr,1,dove+'10\100824_____D05.out',error=err1
openr,2,dove+'10\100824_____D04.out',error=err2
openr,3,dove+'10\100824_____A04.out',error=err2
;for i=0,30-1 do begin
readf,1,arra1,format='(i4,5i3,i4,f6.1,1x,f7.3,f4.1,3i5,2000i6)'
readf,2,arra2,format='(i4,5i3,i4,f6.1,1x,f7.3,f4.1,3i5,2000i6)'
readf,3,arra3 ;,format='(i4,5i3,i4,f6.1,1x,f7.3,f4.1,3i5,2000i6)'
oplot,arra1(13:31,*),arra2(13:31,*),color=0,psym=1
oplot,arra1(32:72,*),arra2(32:72,*),color=20,psym=1
oplot,arra1(72:112,*),arra2(72:112,*),color=40,psym=1
oplot,arra1(113:2012,*),arra2(113:2012,*),color=80,psym=1
oplot,arra1(20:27,*),arra2(20:27,*),color=120,psym=1
close,1
close,2
close,3

;plot,[0,0],[0,0],color=0,psym=3,background=255,xrange=[100,3000],yrange=[50,90],ystyle=1,xstyle=1,xtitle='D06',ytitle='A06',charsize=2
plot,[0,0],[0,0],color=0,psym=3,background=255,xrange=[100,9000],yrange=[50,600],ystyle=1,xstyle=1,xtitle='D06',ytitle='A06',charsize=2
oplot,arra2(13:23,*),arra3(16:26,*),color=0,psym=1
oplot,arra2(24:34,*),arra3(27:37,*),color=20,psym=1
oplot,arra2(35:45,*),arra3(38:48,*),color=40,psym=1
oplot,arra2(46:56,*),arra3(49:59,*),color=60,psym=3
oplot,arra2(57:67,*),arra3(60:70,*),color=80,psym=3
oplot,arra2(68:78,*),arra3(71:81,*),color=100,psym=3
oplot,arra2(79:89,*),arra3(82:92,*),color=120,psym=1
oplot,arra2(90:210,*),arra3(93:213,*),color=180,psym=1


plot,[0,0],[0,0],color=0,background=255,xrange=[10,10000],/xlog,xstyle=1,yrange=[1,100],ystyle=1
for it=0,199 do begin
oplot,arra1(13:2012,it),z,color=0
oplot,10*(arra3(16:815,it)-50)-10,z(0:799),color=40
oplot,arra2(13:2012,it),z,color=20
endfor

write_gif,'d:/COMPATTI/10/analogicon2low.gif',tvrd()
end