parameter(n_agn=43)
        character Name_AGN*22
        integer XID,id_gal,id_group_zbest,S_Class
        real x_offset,y_offset,diff,ra_agn,dec_agn,
     +zspec_agn,zpho_agn,P_mem_zbest,sm, sm_err,
     +R_ratio_zbest,zpho_gal

        open(1,file='../data/cross_agn_gal.tbl',status='old')
        open(2,file='../data/sel_agn_gal.tbl')

       write(*,*)'XID,id_gal,id_group,zspec_agn-zpho_agn,'
       write(*,*)'zspec_agn-zpho_gal,(zspec_agn-zpho_gal)/(1.+zpho_gal)'
        ic1=0
        ic2=0
        do 100 i=1,n_agn
          read(1,101) 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

          diff1=abs(zspec_agn-zpho_gal)
          diff2=abs(zpho_agn-zpho_gal)
          diff3=abs(zspec_agn-zpho_gal)/(1.+zpho_gal)

          if(P_mem_zbest.gt.0.1.and.diff3.lt.0.02) then
            ic1=ic1+1
            write(2,101)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
          else
            ic2=ic2+1
            write(*,102)XID,id_gal,id_group_zbest,diff1,diff2,diff3
          endif

100     continue

        write(*,*) 'Number of AGN counterpart AGN is ',ic1
        write(*,*) 'Number of deviated is ',ic2

101     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)
102     format(i5,1x,i8,1x,i8,3(1x,f7.4))
103     format(i5,1x,i8,1x,i8)
        close(1);close(2)
        end