code: def plotEzJz(self,*args,**kwargs):
          
          if not kwargs.get('pot',None) is None: kwargs['pot']= flatten_potential(kwargs.get('pot'))
          return self._orb.plotEzJz(*args,**kwargs)
docs: NAME:
  
             plotEzJzt
  
          PURPOSE:
  
             plot E_z(t)/sqrt(dens(R)) / (E_z(0)/sqrt(dens(R(0)))) along the orbit (an approximation to the vertical action)
  
          INPUT:
  
             pot - Potential instance or list of instances in which the orbit was integrated
  
             d1= plot Ez vs d1: e.g., 't', 'z', 'R'
  
             +bovy_plot.bovy_plot inputs
  
          OUTPUT:
  
             figure to output device
  
          HISTORY:
  
             2010-08-08 - Written - Bovy (NYU)