Make sure not to apply the IsothermalDiskPotential amplitude twice, f…

…ixes #400
jobovy · Oct 28, 2019 · 050b6c0 · 050b6c0
1 parent bfa34da
commit 050b6c0
Showing 1 changed file with 1 addition and 0 deletions.
diff --git 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.):