def mod_from_vels(Vp, Vs, Rho):
    """
    Convenience function to compute bulk and shear moduli from
    Vp, Vp, and Rho values. Input units should be m/s and kg/m3.

    Output units will be in Pascals
    """

    K = Rho*(Vp**2.0 - 4.0/3.0*Vs**2)
    G = Rho*Vs**2.0

    return K, G