From 5fe91fa71e5b3214a1dd8128f10680ab586a4c3f Mon Sep 17 00:00:00 2001 From: Andrew Benson Date: Tue, 4 Feb 2025 11:16:49 -0800 Subject: [PATCH] =?UTF-8?q?fix:=20The=20central=20density=20of=20the=20?= =?UTF-8?q?=CE=B2-profile=20was=20evaluated=20at=20the=20wrong=20radius=20?= =?UTF-8?q?in=20cooling=20radius=20calculations?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit This was being evaluated at the outer radius instead of at zero radius. --- source/cooling.cooling_radius.beta_profile.F90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/source/cooling.cooling_radius.beta_profile.F90 b/source/cooling.cooling_radius.beta_profile.F90 index b6efc9c3c..bf19a8b87 100644 --- a/source/cooling.cooling_radius.beta_profile.F90 +++ b/source/cooling.cooling_radius.beta_profile.F90 @@ -18,8 +18,8 @@ !! along with Galacticus. If not, see . !!{ - Implementation of a cooling radius class for $\beta$-profile halos, assuming collisional ionization equilibrium such that cooling - time scales as inverse density. + Implementation of a cooling radius class for $\beta$-profile halos, assuming collisional ionization equilibrium such that + cooling time scales as inverse density. !!} use :: Cooling_Times , only : coolingTimeClass @@ -285,7 +285,7 @@ double precision function betaProfileRadiusGrowthRate(self,node) temperature = kinematicsDistribution_ %temperature(coordinates ) ! Compute density and cooling time at outer radius and zero radius. densityOuter = massDistribution_ %density (coordinates ) - coordinates = [outerRadius,0.0d0,0.0d0] + coordinates = [0.0d0 ,0.0d0,0.0d0] densityZero = massDistribution_ %density (coordinates ) coolingTimeZero = self %coolingTime_%time (node,temperature,densityZero ,hotAbundances,chemicalFractions*densityZero ,self%radiation) coolingTimeOuter = self %coolingTime_%time (node,temperature,densityOuter,hotAbundances,chemicalFractions*densityOuter,self%radiation) @@ -377,7 +377,7 @@ double precision function betaProfileRadius(self,node) temperature = kinematicsDistribution_ %temperature(coordinates ) ! Compute density and cooling time at outer radius and zero radius. densityOuter = massDistribution_ %density (coordinates ) - coordinates = [outerRadius,0.0d0,0.0d0] + coordinates = [0.0d0 ,0.0d0,0.0d0] densityZero = massDistribution_ %density (coordinates ) coolingTimeZero = self %coolingTime_%time (node,temperature,densityZero ,hotAbundances,chemicalFractions*densityZero ,self%radiation) coolingTimeOuter = self %coolingTime_%time (node,temperature,densityOuter,hotAbundances,chemicalFractions*densityOuter,self%radiation)