parameter(n=28) parameter(m1=6181) parameter(m2=5285) parameter(m3=4348) parameter(m4=2923) real x(n),y1(m1),y2(m2),y3(m3),y4(m4) open(1,file='../data/agn_zspec.tbl',status='old') open(21,file='../data/gal_zpho_0.1.tbl',status='old') open(22,file='../data/gal_zpho_0.3.tbl',status='old') open(23,file='../data/gal_zpho_0.5.tbl',status='old') open(24,file='../data/gal_zpho_0.7.tbl',status='old') do 10 i=1,n read(1,'(f3.1)')x(i) 10 continue do 20 j=1,m1 read(21,'(e15.4)')y1(j) 20 continue do 30 j=1,m2 read(22,'(e15.4)')y2(j) 30 continue do 40 j=1,m3 read(23,'(e15.4)')y3(j) 40 continue do 50 j=1,m4 read(24,'(e15.4)')y4(j) 50 continue c call pgbegin(0,'N_z.eps/vcps',2,1) call pgbegin(0,'?',2,1) call pgpaper(0.0,0.4) call pgscf(2) call pgsls(1) call pgslw(3) call pgsch(2.0) ************ plot the N_AGN-zspec ********** call pgenv(0.0,1.2,0.,10.,0,0) call pglabel('\fiz\d\frspec\u','\fiN\d\frAGN\u','') call pgtext(1.08,1.0,'(a)') call pghist(n,x,0.,1.3,13,1) ************ plot the N_gal-zpho ********** call pgenv(-0.1,1.2,0.,1600.,0,0) call pglabel('\fiz\d\frpho\u','\fiN\d\frgal\u','') call pgtext(1.08,150.,'(b)') call pgsci(1) call pghist(m1,y1,0.,1.3,13,1) call pgsch(1.3) call pgmove(0.9,1450.) call pgdraw(1.1,1450.) call pgtext(0.6,1450.,'\fiP\d\frmem\u>0.1') call pgsci(2) call pghist(m2,y2,0.,1.3,13,1) call pgmove(0.9,1300.) call pgdraw(1.1,1300.) call pgtext(0.6,1300.,'\fiP\d\frmem\u>0.3') call pgsci(3) call pghist(m3,y3,0.,1.3,13,1) call pgmove(0.9,1150.) call pgdraw(1.1,1150.) call pgtext(0.6,1150.,'\fiP\d\frmem\u>0.5') call pgsci(4) call pghist(m4,y4,0.,1.3,13,1) call pgmove(0.9,1000.) call pgdraw(1.1,1000.) call pgtext(0.6,1000.,'\fiP\d\frmem\u>0.7') call pgend end