From 050b6c0e62c76c91011068db4f23db6cc05547a6 Mon Sep 17 00:00:00 2001 From: Jo Bovy Date: Mon, 28 Oct 2019 12:52:30 -0400 Subject: [PATCH] Make sure not to apply the IsothermalDiskPotential amplitude twice, fixes #400 --- galpy/potential/IsothermalDiskPotential.py | 1 + 1 file changed, 1 insertion(+) diff --git a/galpy/potential/IsothermalDiskPotential.py b/galpy/potential/IsothermalDiskPotential.py index 795f4a1ce..48f2192d6 100644 --- a/galpy/potential/IsothermalDiskPotential.py +++ b/galpy/potential/IsothermalDiskPotential.py @@ -46,6 +46,7 @@ def __init__(self,amp=1.,sigma=0.1,ro=None,vo=None): sigma= sigma.to(units.km/units.s).value/self._vo self._sigma2= sigma**2. self._H= sigma/numpy.sqrt(8.*numpy.pi*self._amp) + self._amp= 1. # Need to manually set to 1, because amp is now contained in the combination of H and sigma^2 self.hasC= True def _evaluate(self,x,t=0.):