c c ACCS15FIG.FOR c real x1(71),y1(71) real x2(71),y2(71) real x3(71) call pgbegin(0,'?',1,1) call pgslw(2) call pgscf(2) open(unit=1,file='cross15.dat',status='old') open(unit=2,file='autos15.dat',status='old') read(1,*)(x1(i),y1(i),i=1,71) read(2,*)(x2(i),y2(i),i=1,71) do 10 i=1,71 x3(i)=x2(i)+432. 10 continue call pgenv(384.,456.,-1.5,1.5,0,1) call pglabel('Lag (days)', *'DCC vs. DAC(+ 432 days)','Secondary events (Bin = 15 days)') call pgpoint(71,x1,y1,17) call pgpoint(71,x3,y2,23) call pgend end