pro pgal 

;Purpose: to determine the numerical distribution of galaxies by redshift.

;Histogram must be made in order to bin redshifts.  

;Calling Sequence: pgal

!p.charsize = 2

galaxy='~/Data/Galaxy25.csv'
e=read_csv(galaxy)

gra=e.field1
gdec=e.field2
gz=e.field7

window,0,xsize=800,ysize=600

bin=.01
min=0.0
max=1.0

h=histogram(gz,bin=bin,min=min,max=max)
n=n_elements(h)
x=findgen(n)*bin+min

plot,x,h,color=!red,xtit='Redshift Bins',ytit='# of Galaxies'

window,1,xsize=800,ysize=600

;a = area
;p = density
;h = # of galaxies

a = 25
t = total(h)
c = (t - total(h,/cumu))/a
p=h/a/bin

plot,x,c,color=!lightblue,psym=-7,xtit='Redshift Bins',ytit='# Galaxies per Deg^2'

get_lun,unit
file='~/data/pgal.txt'
openw,unit,file

for i=0L,n_elements(x)-1 do begin
   printf,unit,x[i],c[i],p[i]
endfor

free_lun,unit

return

end