program summary

C  Reads pulsar catalog data from bincat.dat (on lu 23) and provides
C  summary information on various pulsar catagories.

	include 'rdascat.h'
	parameter (NMAX=40)
	character*16 label(NMAX)
	integer ns(NMAX)

C	call rdascat(23,npsr)		!Read the ASCII catalog
	call rdbincat(23,npsr)		!Read the binary catalog (faster)

	do 10 j=1,NMAX
10	ns(j)=0

	ns(1)=npsr
	label(1)='Pulsars'

	do 20 i=1,npsr
	nt=ntype(i)
	j=2
	label(j)='Pdot > 0        '
	if(pdot(i).gt.0.d0) ns(j)=ns(j)+1
	j=j+1
	label(j)='Pdot > +Pdote   '
	if(pdot(i).gt.pdote(i)) ns(j)=ns(j)+1
	j=j+1
	label(j)='Pdot < -Pdote   '
	if(pdot(i).lt.-pdote(i)) ns(j)=ns(j)+1
	j=j+1
	label(j)='f2              '
	if(f2e(i).ne.0.d0) ns(j)=ns(j)+1
	j=j+1
	label(j)='w50             '
	if(w50(i).gt.0.d0) ns(j)=ns(j)+1
	j=j+1
	label(j)='S400            '
	if(s400(i).gt.0.d0) ns(j)=ns(j)+1
	j=j+1
	label(j)='S1400           '
	if(s1400(i).gt.0.d0) ns(j)=ns(j)+1
	j=j+1
	label(j)='Spectral Index  '
	if(s1400(i).gt.0.d0 .and. 
     +    (s400(i).gt.0.d0.or.s600(i).gt.0.d0)) ns(j)=ns(j)+1
	j=j+1
	label(j)='DM              '
	if(dm(i).gt.0.d0) ns(j)=ns(j)+1
	j=j+1
	label(j)='RM              '
	if(rme(i).ne.0.d0) ns(j)=ns(j)+1
	j=j+1
	label(j)='tau_s           '
	if(tau(i).ne.0.d0) ns(j)=ns(j)+1
	j=j+1
	label(j)='Binary          '
	if(iand(nt,8).ne.0) ns(j)=ns(j)+1		!Binary
	j=j+1
	label(j)='Cluster         '
	if(iand(nt,1).ne.0) ns(j)=ns(j)+1		!Cluster
	j=j+1
	label(j)='Glitches        '
	if(iand(nt,4).ne.0) ns(j)=ns(j)+1		!Glitches
	j=j+1
	label(j)='X or gamma-ray  '
	if(iand(nt,128).ne.0) ns(j)=ns(j)+1		!X or gamma ray
	j=j+1
	label(j)='Interpulse      '
	if(iand(nt,64).ne.0) ns(j)=ns(j)+1		!Interpulse
	j=j+1
	label(j)='P < 25 ms       '
	if(iand(nt,16).ne.0) ns(j)=ns(j)+1		!P < 25 ms
	j=j+1
	label(j)='Recycled        '
	if(iand(nt,32).ne.0) ns(j)=ns(j)+1		!Recycled
	j=j+1
	label(j)='SNR             '
	if(iand(nt,2).ne.0) ns(j)=ns(j)+1		!SNR
	j=j+1
	label(j)='Extragalactic   '
	if(iand(nt,256).ne.0) ns(j)=ns(j)+1		!Extragalactic
	j=j+1
	label(j)='Molonglo 1'
	if(iand(nscode(i),1).ne.0) ns(j)=ns(j)+1	!Molonglo 1 (lv71, etc)
	j=j+1
	label(j)='Jodrell Bank 1'
	if(iand(nscode(i),2).ne.0) ns(j)=ns(j)+1	!Jodrell Bank 1 (dls77)
	j=j+1
	label(j)='Arecibo 1'
	if(iand(nscode(i),4).ne.0) ns(j)=ns(j)+1	!Arecibo 1 (ht75a)
	j=j+1
	label(j)='Molonglo 2'
	if(iand(nscode(i),8).ne.0) ns(j)=ns(j)+1	!Molonglo 2 (mlt+78)
	j=j+1
	label(j)='Green Bank 1'
	if(iand(nscode(i),16).ne.0) ns(j)=ns(j)+1	!Green Bank 1 (dth78)
	j=j+1
	label(j)='Green Bank 2'
	if(iand(nscode(i),32).ne.0) ns(j)=ns(j)+1	!Green Bank 2 (dtws85)
	j=j+1
	label(j)='Jodrell Bank 2'
	if(iand(nscode(i),64).ne.0) ns(j)=ns(j)+1	!Jodrell Bank 2 (cl86)
	j=j+1
	label(j)='Green Bank 3'
	if(iand(nscode(i),128).ne.0) ns(j)=ns(j)+1	!Green Bank 3 (stwd85)
	j=j+1
	label(j)='Arecibo 2'
	if(iand(nscode(i),256).ne.0) ns(j)=ns(j)+1	!Arecibo 2 (fru89,nic92)
	j=j+1
	label(j)='Parkes 1'
	if(iand(nscode(i),512).ne.0) ns(j)=ns(j)+1	!Parkes 1 (jlm+92)
	j=j+1
	label(j)='Arecibo 3'
	if(iand(nscode(i),1024).ne.0) ns(j)=ns(j)+1	!Arecibo 3 (nic92)
20	continue 
	jz=j
	write(*,1020) (label(j),ns(j),j=1,jz)
1020	format(1x,a16,i4)
 
	end