#!/usr/bin/env python
# -----------------------------------------------------------------------------
# GENHERNQUIST.VCIRC
# Laura L Watkins [lauralwatkins@gmail.com]
# -----------------------------------------------------------------------------

from numpy import *
from astropy import units as u, constants as c
from scipy import special
from encmass import encmass


def vcirc(r, norm, rs, alpha, beta, gamma):
    
    """
    Enclosed mass profile of a generalised Hernquist model.
    
    INPUTS
      r     : radial variable (requires unit)
      norm  : normalisation (requires unit)
      rs    : scale radius of model (requires unit)
      alpha : sharpness of transition between inner and outer
      beta  : outer logarithmic slope
      gamma : inner logarithmic slope
    """
    
    mass = encmass(r, norm, rs, alpha, beta, gamma)
    
    vcirc = (sqrt(c.G*mass/r)).to("km/s")
    
    return vcirc

Add routine to calculate circular velocity.