c c CROSP.FOR (Cross-correlation of the Kundic c (main) event in 1995 and the main event in 1996) c QSO 0957+561A,B - April/May 2001 c implicit double precision (a-h),(o-z) dimension ta(13),a(13),ea(13),x(13),y(13) dimension tb(14),u(14),v(14),b(14),eb(14) dimension am(13),eam(13),bm(14),ebm(14) dimension siga(13),sigb(14),udcf(13,14),tlag(13,14) open(unit=1,file='crosp15.dat',status='new') open(unit=2,file='95p.dat',status='old') open(unit=3,file='96p.dat',status='old') read(2,*)(ta(i),a(i),ea(i),x(i),y(i),i=1,13) read(3,*)(tb(i),u(i),v(i),b(i),eb(i),i=1,14) do 10 i=1,13 if(i.eq.1)am(i)=a(i) if(i.gt.1)am(i)=a(i)+am(i-1) if(i.eq.1)eam(i)=ea(i) if(i.gt.1)eam(i)=ea(i)+eam(i-1) 10 continue amedio=am(13)/13.D0 errora2=(eam(13)/13.D0)**2.D0 do 20 i=1,14 if(i.eq.1)bm(i)=b(i) if(i.gt.1)bm(i)=b(i)+bm(i-1) if(i.eq.1)ebm(i)=eb(i) if(i.gt.1)ebm(i)=eb(i)+ebm(i-1) 20 continue bmedio=bm(14)/14.D0 errorb2=(ebm(14)/14.D0)**2.D0 do 30 i=1,13 if(i.eq.1)siga(i)=(a(i)-amedio)**2.D0 if(i.gt.1)siga(i)=(a(i)-amedio)**2.D0+siga(i-1) 30 continue sigmaa2=siga(13)/12.D0 do 40 i=1,14 if(i.eq.1)sigb(i)=(b(i)-bmedio)**2.D0 if(i.gt.1)sigb(i)=(b(i)-bmedio)**2.D0+sigb(i-1) 40 continue sigmab2=sigb(14)/13.D0 raiz=dsqrt((sigmaa2-errora2)*(sigmab2-errorb2)) do 50 i=1,13 do 60 j=1,14 udcf(i,j)=(a(i)-amedio)*(b(j)-bmedio)/raiz tlag(i,j)=tb(j)-ta(i) 60 continue 50 continue c *********************************************************** c BINS WITH SIZE OF 2*alfa c *********************************************************** alfa=7.5D0 c *********************************************************** do 70 k=1,51 tau=394.D0+k bdcf=0.D0 bdat=0.D0 do 80 i=1,13 do 90 j=1,14 if(tlag(i,j).ge.(tau-alfa).and.tlag(i,j).lt.(tau+alfa)) *bdcf=bdcf+udcf(i,j) if(tlag(i,j).ge.(tau-alfa).and.tlag(i,j).lt.(tau+alfa)) *bdat=bdat+1.D0 90 continue 80 continue if(bdat.ge.1.D0)write(1,200)tau,bdcf/bdat 70 continue stop 200 format(E10.3,4x,E10.3) end