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