parameter(n=28)
	real fra1(n),err1(n),fra2(n),err2(n),fra3(n),
     +	err3(n),fra4(n),err4(n),agn_zspec(n)
	real R200(n)

	open(1,file='../data/final_agn_gal.tbl',status='old')

	do 10 i=1,n 
	 read(1,101)agn_zspec(i),R200(i)
10	continue

101	format(92x,f5.3,52x,f6.3)

   
c************** begin to plot *******************

c	call pgbegin(0,'r_200_z.eps/vcps',1,1)
	call pgbegin(0,'?',1,1)
	call pgpaper(0.0,0.8)
	call pgscf(2)
	call pgsls(1)
	call pgslw(3)
	call pgsch(1.3) 
	call pgsci(1)

c************ plot the p > 0.1**********

	call pgenv(0.0,1.1,-0.1,1.0,0,0)
	call pglabel('\fiz\fr\dspec\u','\fiR/R\fr\d200\u','')
	call pgsch(1.8) 
	call pgpt(27,agn_zspec,R200,4)

	 call pgend
        
	close(1)
	end