c c ACCP15FIG.FOR c real x1(51),y1(51) real x2(51),y2(51) real x3(51) call pgbegin(0,'?',1,1) call pgslw(2) call pgscf(2) open(unit=1,file='crosp15.dat',status='old') open(unit=2,file='autop15.dat',status='old') read(1,*)(x1(i),y1(i),i=1,51) read(2,*)(x2(i),y2(i),i=1,51) do 10 i=1,51 x3(i)=x2(i)+418. 10 continue call pgenv(394.,446.,-1.5,1.5,0,1) call pglabel('Lag (days)', *'DCC vs. DAC(+ 418 days)','Main events (Bin = 15 days)') call pgpoint(51,x1,y1,17) call pgpoint(51,x3,y2,23) call pgend end