pro analogico_N2_scatter2 device,decomposed=1 loadct,38 dove='d:\compatti\' z=fltarr(2000) for iz=1,2000 do z(iz-1)=iz arra1=fltarr(2013,31) arra2=fltarr(2013,31) arra3=fltarr(816,31) !p.multi=[0,4,3] window,0,xsize=1100,ysize=750;,title=cchan(ich) nome='10\100914_satu' chan=['1','2','4','5','6','7'] for ich=0,5 do begin openr,2,dove+nome+'D0'+chan(ich)+'.out',error=err2 openr,3,dove+nome+'A0'+chan(ich)+'.out',error=err2 readf,2,arra2,format='(i4,5i3,i4,f6.1,1x,f7.3,f4.1,3i5,2000i6)' readf,3,arra3 close,2 close,3 amin=min(arra3(16:210,*)) amaxv=max(arra3(16:210,*)) symm=3 symz=0.3 plot,[0,0],[0,0],color=0,psym=3,background=255,xrange=[1,20000],yrange=[0.9*amin,1.1*amaxv],ystyle=1,xstyle=1,xtitle='D',ytitle='A',charsize=1.5,/xlog,/ylog,title=chan(ich) oplot,arra2(13:23,*),arra3(16:26,*),color=0,psym=symm,symsize=symz oplot,arra2(24:34,*),arra3(27:37,*),color=20,psym=symm,symsize=symz oplot,arra2(35:45,*),arra3(38:48,*),color=40,psym=symm,symsize=symz oplot,arra2(46:56,*),arra3(49:59,*),color=60,psym=symm,symsize=symz oplot,arra2(57:67,*),arra3(60:70,*),color=80,psym=symm,symsize=symz oplot,arra2(68:78,*),arra3(71:81,*),color=100,psym=symm,symsize=symz oplot,arra2(79:89,*),arra3(82:92,*),color=120,psym=symm,symsize=symz oplot,arra2(90:210,*),arra3(93:213,*),color=180,psym=symm,symsize=symz xmax=max(arra2(13:214,*)) amin=min(arra3(16:217,*)) amaxv=max(arra3(16:217,*)) amax=max(arra3(16:217,*)-amin) coeff=xmax/amax ;xmin=min(arra2(13:214,*)) plot,[0,0],[0,0],color=0,background=255,xrange=[1,1.1*xmax],/xlog,xstyle=1,yrange=[1,200],ystyle=1,charsize=1.5,/ylog,title=string(xmax)+string(amaxv) for it=0,30 do begin oplot,arra2(13:2012,it),z,color=20 ;oplot,10*(arra3(16:815,it)-50)-10,z(0:799),color=40 oplot,coeff*(arra3(16:815,it)-amin),z(0:799),color=40 endfor endfor write_gif,'d:/COMPATTI/10/100914_satu.gif',tvrd() end