code: def estimateTdisrupt(self,deltaAngle):
          
          return deltaAngle\
              /numpy.sqrt(numpy.sum(self._dsigomeanProg**2.))
docs: NAME:
  
             estimateTdisrupt
  
          PURPOSE:
  
             estimate the time of disruption
  
          INPUT:
  
             deltaAngle- spread in angle since disruption
  
          OUTPUT:
  
             time in natural units
  
          HISTORY:
  
             2013-11-27 - Written - Bovy (IAS)