Page scientifique personnelle de Patrick Valageas

Voici le programme Fortran 90 qui est utilisé dans l'article astro-ph/0403593 pour calculer la distribution de probabilité de la densité. Il se compose de deux modules:

• Module "phizetaP.f90"
  Ce module calcule les fonctions génératrices des cumulants ainsi que la distribution de probabilité P(\rho) où \rho=1+\delta est la sur-densité dans une sphère de rayon physique R.
  
• Module "dens.f90"
  Ce module calcule la courbe P(\rho) pour un certain redshift z et une certaine échelle physique R en utilisant le module "phizetaP.f90".
  

Le programme principal, contenu dans le fichier "dens.f90", appelle la subroutine "calPrho(z,R)" qui calcule P(\rho) à partir du modèle décris dans astro-ph/0403593 ainsi que la prédiction Pln(\rho) du modèle lognormal. Le résultat est écris dans le fichier "PrhoR1z0.dat".

Avant de calculer P(\rho) pour diverses valeurs de \rho, la subroutine "calPrho(z,R)" appelle "nspecRz(R,z)" qui initialise divers paramètres pour le redshift z et l'échelle physique R. La subroutine "nspecRz(R,z)", située au début du module "phizetaP.f90", nécessite trois fonctions qui doivent etre fournies par l'utilisateur:

• "nsig(R,z)": donne la pente locale "n" du spectre de puissance linéaire à l'échelle physique R. Elle est définie à partir de la dérivée logarithmique de la variance linéaire \sigma^2.
  
• "Xir(R,z)": donne la variance non-linéaire <\delta^2> à l'échelle physique R et au redshift z.
  
• "XirL(R,z)": donne la variance linéaire <\delta_L^2> à l'échelle physique R et au redshift z (i.e. \sigma^2).
  

Ceci fait, la subroutine "calPrho(z,R)" calcule P(\rho) pour diverses valeurs de \rho en appelant la fonction "Prho(rho)" contenue dans le module "phizetaP.f90". Dans les fichiers à télécharger les appels aux trois fonctions "nsig(R,z),Xir(R,z),XirL(R,z)" sont mis en commentaire et le programme calcule P(\rho) pour le cas l=1,z=0 étudié dans l'article astro-ph/0403593.

En cas de problème ou pour des suggestions, commentaires,..., contactez-moi par e-mail.