USE linearpk double precision :: k_ov_h,linear_pk character(len=128) :: filename integer :: n external linear_pk filename='wmap5baosn_max_likelihood_matterpower.dat' n=896 ! # of lines in the file CALL open_linearpk(filename,n) k_ov_h=0.01d0 ! h Mpc^-1 print*,'P(k) at k=',k_ov_h,' h Mpc^-1 is',linear_pk(k_ov_h),' h^-3 Mpc^3' CALL close_linearpk end