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