@@ -0,0 +1,918 @@
+#!/usr/bin/env python%0A# -----------------------------------------------------------------------------%0A# GENHERNQUIST.VCIRC%0A# Laura L Watkins %5Blauralwatkins@gmail.com%5D%0A# -----------------------------------------------------------------------------%0A%0Afrom numpy import *%0Afrom astropy import units as u, constants as c%0Afrom scipy import special%0Afrom encmass import encmass%0A%0A%0Adef vcirc(r, norm, rs, alpha, beta, gamma):%0A    %0A    %22%22%22%0A    Enclosed mass profile of a generalised Hernquist model.%0A    %0A    INPUTS%0A      r     : radial variable (requires unit)%0A      norm  : normalisation (requires unit)%0A      rs    : scale radius of model (requires unit)%0A      alpha : sharpness of transition between inner and outer%0A      beta  : outer logarithmic slope%0A      gamma : inner logarithmic slope%0A    %22%22%22%0A    %0A    mass = encmass(r, norm, rs, alpha, beta, gamma)%0A    %0A    vcirc = (sqrt(c.G*mass/r)).to(%22km/s%22)%0A    %0A    return vcirc%0A

Add routine to calculate circular velocity.