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 <http://www.gnu.org/licenses/>.
 
   !!{
-  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)