'''
using https://arxiv.org/pdf/1607.04275.pdf
convert stellar mass to velocity disperson
'''
import sys
import numpy as np
import ipdb as pdb
def main( stellarMass ):

    Mb = 10**(10.26)
    alpha1 = 0.403
    alpha2 = 0.293
    sigmaB = 10**2.073
    if stellarMass <= Mb:
        velocity = sigmaB*(np.float(stellarMass)/Mb)**alpha1
    else:
        velocity = sigmaB*(np.float(stellarMass)/Mb)**alpha2
    
    return velocity


if __name__ == '__main__':
    main(sys.argv[1])

d05ce99080229ca76e2d28d0f593504ac10b3e3a
<mask token>

<mask token>


def main(stellarMass):
    Mb = 10 ** 10.26
    alpha1 = 0.403
    alpha2 = 0.293
    sigmaB = 10 ** 2.073
    if stellarMass <= Mb:
        velocity = sigmaB * (np.float(stellarMass) / Mb) ** alpha1
    else:
        velocity = sigmaB * (np.float(stellarMass) / Mb) ** alpha2
    return velocity


<mask token>

<mask token>


def main(stellarMass):
    Mb = 10 ** 10.26
    alpha1 = 0.403
    alpha2 = 0.293
    sigmaB = 10 ** 2.073
    if stellarMass <= Mb:
        velocity = sigmaB * (np.float(stellarMass) / Mb) ** alpha1
    else:
        velocity = sigmaB * (np.float(stellarMass) / Mb) ** alpha2
    return velocity


if __name__ == '__main__':
    main(sys.argv[1])

<mask token>
import sys
import numpy as np
import ipdb as pdb


def main(stellarMass):
    Mb = 10 ** 10.26
    alpha1 = 0.403
    alpha2 = 0.293
    sigmaB = 10 ** 2.073
    if stellarMass <= Mb:
        velocity = sigmaB * (np.float(stellarMass) / Mb) ** alpha1
    else:
        velocity = sigmaB * (np.float(stellarMass) / Mb) ** alpha2
    return velocity


if __name__ == '__main__':
    main(sys.argv[1])

'''
using https://arxiv.org/pdf/1607.04275.pdf
convert stellar mass to velocity disperson
'''
import sys
import numpy as np
import ipdb as pdb
def main( stellarMass ):

    Mb = 10**(10.26)
    alpha1 = 0.403
    alpha2 = 0.293
    sigmaB = 10**2.073
    if stellarMass <= Mb:
        velocity = sigmaB*(np.float(stellarMass)/Mb)**alpha1
    else:
        velocity = sigmaB*(np.float(stellarMass)/Mb)**alpha2
    
    return velocity


if __name__ == '__main__':
    main(sys.argv[1])
