;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;; PROCEDURE lUMDIST   ;;;
;;; By Darrin Casebeer  ;;;
;;;    Kevin Tubbs      ;;;
;;;    Summer 98        ;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;

PRO lumdist
COMMON COSMO,Omegamass,OmegaLambda,Ho,LIMITS
Ho=60.0
LIMITS=1000
deltaz=0.005
 
znow=DBLARR(LIMITS)
z=DBLARR(LIMITS)

holder=0.0

FOR I=0,LIMITS-1 DO BEGIN
  znow=0.00
  z(I)=holder
  holder=holder+deltaz
ENDFOR
z(0)=0.001

mu=DBLARR(LIMITS)
dm=DBLARR(LIMITS)
dl=DBLARR(LIMITS)

Omegamass=1.0
OmegaLambda=0.0
dm=QSIMP('dl',0.00,z)
Omegak=1-Omegamass-OmegaLambda
dm=1/Ho*dm
dl=dm*(1+z)*3e5
mu=5*ALOG10(dl)+25
PLOT,z,dl

Omegamass=0.5
OmegaLambda=0.5
dm=QSIMP('dl',0.00,z)
Omegak=1-Omegamass-OmegaLambda
dm=1/Ho*dm
dl=dm*(1+z)*3e5
mu=5*ALOG10(dl)+25
OPLOT,z,dl

Omegamass=0.1
OmegaLambda=0.9
dm=QSIMP('dl',0.00,z)
Omegak=1-Omegamass-OmegaLambda
dm=1/Ho*dm
dl=dm*(1+z)*3e5
mu=5*ALOG10(dl)+25
OPLOT,z,dl
END