real ra_gal_sec,dec_gal_sec,ra_agn,dec_agn,ra_1,dec_1
	real ra_agn_sec,dec_agn_sec,ra_2,dec_2
	real zspec_agn,zpho_agn
	
	parameter(n_agn=730,n_xg=7984)
	parameter(pi=3.1415926, box=1.0)

	integer i,j
	integer id_gal,mmggs,id_group, Group_flag,mmggs_zbest,id_group_zbest, Group_flag_zbest
	real auto_mag,ra_gal,dec_gal,sm,sm_err,zpho_gal,zp_1,zp_2, P_mem, R_ratio,P_mem_zbest,R_ratio_zbest,color
	
	character Name_AGN*22
	integer XID,CFlag,chid,IFlag,CID,IID,S_Class,r_zspec
	real rmag,Imag,Kmag,mag36,mag45,mag58,mag80,mag24
	real SFlux,HFlux,UHFlux,HR

 	open(21,file='../COSMOS/table/xAGN_BNL_lim.tbl',status='old')
 	open(22,file='../data/xg_gal_p0.tbl',status='old')
 	open(1,file='../data/cross_agn_gal.tbl')

	itotal=0
 	do 100 i=1,n_agn
c	do 100 i=1,1
	read(21,101) Name_AGN,XID,ra_agn,dec_agn,SFlux,HFlux,UHFlux,CFlag,HR,chid,
     +       IFlag,CID,IID,rmag,Imag,Kmag,mag36,mag45,mag58,mag80,mag24, 
     +       zspec_agn,S_Class,r_zspec, zpho_agn

c	write(*,101)Name_AGN,XID,ra_agn,dec_agn,SFlux,HFlux,UHFlux,CFlag,HR,chid,
c     +       IFlag,CID,IID,rmag,Imag,Kmag,mag36,mag45,mag58,mag80,mag24, 
c     +       zspec_agn,S_Class,r_zspec, zpho_agn

	ra_agn_sec=ra_agn*3600.0   !degree --> seconds
	dec_agn_sec=dec_agn*3600.0
c	print*, ra_agn_sec,dec_agn_sec
c********searching box defined
	dec0=dec_agn*pi/180.0
	ra_1=ra_agn_sec-box/cos(dec0)
	ra_2=ra_agn_sec+box/cos(dec0)
	dec_1=dec_agn_sec-box
	dec_2=dec_agn_sec+box
c	print*, 'box=',ra_1,ra_2,dec_1,dec_2

	icount=0
 	

 	do 200 j=1,n_xg

c  	do 200 j=1,1

	read(22,102) id_gal,auto_mag,ra_gal,dec_gal,sm,sm_err,zpho_gal,zp_1,zp_2,
     +     mmggs, P_mem, id_group, Group_flag, R_ratio,
     +     mmggs_zbest,P_mem_zbest, id_group_zbest, Group_flag_zbest,
     +     R_ratio_zbest,color

c	write(*,103) XID,id_gal,id_group_zbest,x_offset,y_offset,diff,
c     +        Name_AGN,ra_agn,dec_agn,
c     +        zspec_agn,zpho_agn,P_mem_zbest,sm, sm_err,
c     +        R_ratio_zbest,zpho_gal,S_Class

	ra_gal_sec=ra_gal*3600.
	dec_gal_sec=dec_gal*3600.
	
	if(ra_gal_sec.lt.ra_1) goto 200
	if(ra_gal_sec.gt.ra_2) goto 200
	if(dec_gal_sec.lt.dec_1) goto 200
	if(dec_gal_sec.gt.dec_2) goto 200
	icount=icount+1
	x_offset=(ra_gal_sec-ra_agn_sec)*cos(dec0)
	y_offset=dec_gal_sec-dec_agn_sec

	diff=sqrt(x_offset**2.0+ y_offset**2.0)
	

	write(1,103) XID,id_gal,id_group_zbest,x_offset,y_offset,diff,
     +        Name_AGN,ra_agn,dec_agn,
     +        zspec_agn,zpho_agn,P_mem_zbest,sm, sm_err,
     +        R_ratio_zbest,zpho_gal,S_Class

	itotal=itotal+1
	if(icount.gt.1) write(*,*) 'got it again!','AGN_ID= ',XID

200     continue

c	if(icount.eq.0)write(*,*) '  no found!','COSMOS number= ',num

	rewind(22)
100	continue
	
	write(*,*) 'TOTAL NUMBER OF ASSOCIATED GALAXIES IS ',itotal
101	format(a22,x,i5,x,f11.7,x,f9.7,3(x,e9.3),x,i1,x,f5.2,x,i5,x,i2,x,i7,x,i7,
     +        x,3(f6.2,x),4(f7.2,x),f6.2,x,f7.3,x,i3,x,i3,x,f6.2)
102	format(9x,i6,4x,e15.9,4x,e23.17,4x,e23.17,4x,5(e15.9,4x),5x,
     +          i1,4x,e15.9,7x,i4,12x,i1,4x,e15.9,13x,
     +         i1,4x,e15.9,14x,i3,18x,i1,4x,e15.9,4x,e15.9)

1022	format(8x,i7,23x,e23.17,4x,e23.17,4x,e15.9,23x,e15.9,52x,e15.9, 
     +        8x,i3,16x,e16.9)

103	format(i5,1x,i8,1x,i8,3f7.2,1x,a22,1x,f11.7,1x,f9.7,1x,f7.3,1x,f6.2,1x,f11.5,1x,
     +          2e15.4,1x,f7.3,1x,f6.2,1x,i3)


	close(21)
	close(22)
	close(1)
	goto 998
999	write(*,*) 'readerr at i= ',i,' and j= ',j
998	stop
	end