pro plt,dt121,dt122,dt123,dt124 ttl='crl1=150.550, crl2=190.505' ofl='test14_2075_150.550_190.505.ps' tek_color psout,011,ofl mx=MAX(dt121) & mn=MIN(dt123) plot,dt121,xrange=[0,14],yrange=[mn,mx],xsty=1,ysty=1,title=ttl,xtitle='CR Number from 2075 to 2089',ytitle='Time difference from crl1 to crl2 in days' oplot,dt122,line=2,color=2 oplot,dt123,line=3,color=3 ; dt124=fltarr(15) ; for i=0,14 do dt124(i)=-3.02718 oplot,dt124,line=4,color=4 oplot,[0.5,1.5],[-3.02,-3.02],color=0 xyouts,2.0,-3.02,'Xuepu' oplot,[0.5,1.5],[-3.022,-3.022],line=4,color=4 xyouts,2.0,-3.022,'Roger' oplot,[0.5,1.5],[-3.024,-3.024],line=3,color=3 xyouts,2.0,-3.024,'suninfo' oplot,[0.5,1.5],[-3.026,-3.026],line=2,color=2 xyouts,2.0,-3.026,'times' ddt=(dt121-dt124)*86400 mx=STRTRIM(MAX(ddt),2) & mn=STRTRIM(MIN(ddt),2) xyouts,6.5,-3.033,'Difference between black and blue:' xyouts,5.5,-3.034,'Maximum='+mx+'sec, Minimum='+mn+'sec' psdone print,ofl end