PROGRAM Compute_da USE cosmo USE angular_distance ! A sample program for computing the proper angular diameter distance ! in units of h^-1 Mpc ! August 23, 2008: E.Komatsu IMPLICIT none integer :: j double precision :: da,z external da ! specify four cosmological parameters ode0=0.723d0 om0=0.277d0 w=-1d0 ! tabulate da(z) by calling "setup_da" call setup_da ! now output da(z) as a function of z... open(1,file='redshift_da.txt') z=0d0 do j=1,10000 write(1,'(2F18.8)')z,da(z) z = z+0.01d0 enddo close(1) stop END PROGRAM Compute_da