/// Radii are in physical units here, so be careful void cosmology::modelNFWhalo(double m200,double z,double &Mvir, double &Rvir, double &cvir, double &R200, double &c200) { Mvir=getMvir(m200,z); Rvir=getRvirfromMvir(Mvir, z)/(1.+z); double Delta=200.; R200=getRDelfromMDel(m200, z, Delta)/(1.+z); cvir=conc(Mvir,z); c200=getc200(cvir,z); }
/// Radii are in physical units here, so be careful void cosmology::modelNFWhalo(double m200,double z,double &Mvir, double &Rvir, double &cvir, double &R200, double &c200) { Mvir=getMvir(m200,z); Rvir=0.169*pow(Mvir/1.e12,1./3.); Rvir*=pow(Delta_crit(z)/178.0,-1.0/3.0); Rvir*=pow(Eofz(z),-2./3.); double Delta=200.; R200=pow( 1.0e12/(4./3.*Delta*3e4*0.3/(8*gee)),1./3.); R200*=pow(m200/1.e12,1./3.); R200*=pow(Omega(z)/0.3,-1.0/3.0); R200*=pow(Eofz(z),-2./3.); cvir=conc(Mvir,z); c200=getc200(cvir,z); }
/// Radii are in comoving units here void cosmology::modelNFWhalo_com(double m200,double z,double &Mvir, double &Rvir, double &cvir, double &R200, double &c200) { Mvir=getMvir(m200,z); Rvir=getRvirfromMvir(Mvir, z); double Delta=200.; R200=getRDelfromMDel(m200, z, Delta); /* R200=pow( 1.0e12/(4./3.*Delta*3e4*0.3/(8*gee)),1./3.); R200*=pow(m200/1.e12,1./3.); R200*=pow(Omega(z)/0.3,-1.0/3.0); R200*=pow(Eofz(z),-2./3.); */ cvir=conc(Mvir,z); c200=getc200(cvir,z); }