code: def jr(self,pot=None,**kwargs):
          
          if not pot is None: pot= flatten_potential(pot)
          _check_consistent_units(self,pot)
          self._orb._setupaA(pot=pot,**kwargs)
          if self._orb._aAType.lower() == 'isochroneapprox':
              return float(self._orb._aA(self(),use_physical=False)[0])
          else:
              return float(self._orb._aA(self,use_physical=False)[0])
docs: NAME:
  
             jr
  
          PURPOSE:
  
             calculate the radial action
  
          INPUT:
  
             pot - potential
  
             type= ('staeckel') type of actionAngle module to use
  
                1) 'adiabatic'
  
                2) 'staeckel'
  
                3) 'isochroneApprox'
  
                4) 'spherical'
                
             +actionAngle module setup kwargs
  
             ro= (Object-wide default) physical scale for distances to use to convert (can be Quantity)
  
             vo= (Object-wide default) physical scale for velocities to use to convert (can be Quantity)
  
             use_physical= use to override Object-wide default for using a physical scale for output
  
          OUTPUT:
  
             jr
  
          HISTORY:
  
             2010-11-30 - Written - Bovy (NYU)
  
             2013-11-27 - Re-written using new actionAngle modules - Bovy (IAS)