integer ic_gal(4),ic_agn(4),ic(4) real ratio(4),d_ratio(4),complete(4),purity(4) data ic/4*0/ data complete/0.97,0.94,0.92,0.80/ data purity/0.63,0.67,0.69,0.78/ open(1,file='../data/pmem_agn_gal.tbl',status='old') ic1=0 ic2=0 ic3=0 ic4=0 c read(1,*) do 100 i=1,26 read(1,101) id_agn,id_gal,id_group, p_mem, ic_gal,ic_agn if(id_group.eq.124) then ic_agn(1)=2*ic_agn(1) ic_agn(2)=2*ic_agn(2) ic_agn(3)=2*ic_agn(3) ic_agn(4)=2*ic_agn(4) end if if(p_mem.gt.0.1) ic1=ic_gal(1)+ic1 if(p_mem.gt.0.3) ic2=ic_gal(2)+ic2 if(p_mem.gt.0.5) ic3=ic_gal(3)+ic3 if(p_mem.gt.0.7) ic4=ic_gal(4)+ic4 do 200 j=1,4 ratio(j)=1.0*ic_agn(j)/ic_gal(j) if(ic_gal(j).eq.0) print*,i if(ic_agn(j).eq.0) then d_ratio(j)=0.0 else d_ratio(j)=ratio(j)*sqrt(1.0/ic_agn(j)+1.0/ic_gal(j))*sqrt(complete(j)/purity(j)) end if 200 continue write(*,'(8f5.2)') ratio,d_ratio 100 continue 101 format(i5,i9,i9,86x,f7.5,50x,4i4,1x,4i2) write(*,*) '-------------------------------' write(*,*) 'P_mem > 0.1',ic1 write(*,*) 'P_mem > 0.3',ic2 write(*,*) 'P_mem > 0.5',ic3 write(*,*) 'P_mem > 0.7',ic4 close(1) end