~abenson/galacticus/v0.9.0

!% Contains a module which implements a isothermal halo spin distribution.

 

module Dark_Matter_Profiles_Concentrations_Gao20008

  !% Implements a isothermal halo spin distribution.

  use Tree_Nodes

  private

  public :: Dark_Matter_Concentrations_Gao20008_Initialize

 

contains

 

  !# <darkMatterConcentrationMethod>

  !#  <unitName>Dark_Matter_Concentrations_Gao20008_Initialize</unitName>

  !# </darkMatterConcentrationMethod>

  subroutine Dark_Matter_Concentrations_Gao20008_Initialize(darkMatterConcentrationMethod,Dark_Matter_Profile_Concentration_Get)

    !% Initializes the ``Gao 2008'' halo concentration module.

    use ISO_Varying_String

    implicit none

    type(varying_string),          intent(in)    :: darkMatterConcentrationMethod

    procedure(),          pointer, intent(inout) :: Dark_Matter_Profile_Concentration_Get

    

    if (darkMatterConcentrationMethod == 'Gao 2008') Dark_Matter_Profile_Concentration_Get => Dark_Matter_Profile_Concentration_Gao2008

  

    return

  end subroutine Dark_Matter_Concentrations_Gao20008_Initialize

 

  double precision function Dark_Matter_Profile_Concentration_Gao2008(thisNode)

    !% Returns the concentration of the dark matter profile of {\tt thisNode} using the method of \cite{gao_redshift_2008}. More

    !% specifically, we fit the redshift dependence of the parameters $A$ and $B$ in the fitting formula of

    !% \cite{gao_redshift_2008} and use these fits to find $A$ and $B$ at any given redshift, from which we then compute the

    !% concentration. Note that the fits of \cite{gao_redshift_2008} were computed using Einasto profile fits and utilizing

    !% $M_{200}$ and $c_{200}$.

    use Tree_Nodes

    use Tree_Node_Methods

    use Cosmology_Functions

    implicit none

    type(treeNode),   intent(inout), pointer :: thisNode

    double precision, parameter              :: littleHubbleConstantGao2008=0.73d0

    double precision                         :: logarithmHaloMass,logarithmExpansionFactor,parameterA,parameterB

  

    logarithmHaloMass                        =dlog10(littleHubbleConstantGao2008*Tree_Node_Mass(thisNode) )

    logarithmExpansionFactor                 =dlog10(Expansion_Factor(           Tree_Node_Time(thisNode)))

    parameterA                               =-0.140d0*dexp(-((logarithmExpansionFactor+0.05d0)/0.35d0)**2)

    parameterB                               = 2.646d0*dexp(-((logarithmExpansionFactor+0.00d0)/0.50d0)**2)

    Dark_Matter_Profile_Concentration_Gao2008=10.0d0**(parameterA*logarithmHaloMass+parameterB)

    return

  end function Dark_Matter_Profile_Concentration_Gao2008

  

end module Dark_Matter_Profiles_Concentrations_Gao20008