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