calculate static pressure at 300 K from Hugoniot data using the constq
    formulation

    :param v: unit-cell volume in A^3
    :param v0: unit-cell volume in A^3 at 1 bar
    :param c0: velocity at 1 bar in km/s
    :param s: slope of the velocity change
    :param gamma0: Gruneisen parameter at 1 bar
    :param q: logarithmic derivative of Gruneisen parameter
    :param theta0: Debye temperature in K
    :param n: number of elements in a chemical formula
    :param z: number of formula unit in a unit cell
    :param mass: molar mass in gram
    :param c_v: heat capacity
    :param three_r: 3 times gas constant.
        Jamieson modified this value to compensate for mismatches
    :param t_ref: reference temperature, 300 K
    :return: static pressure in GPa
    :note: 2017/05/18 I am unsure if this is actually being used in pytheos

def jamieson_pst(v, v0, c0, s, gamma0, q, theta0, n, z, mass, c_v,
                 three_r=3. * constants.R, t_ref=300.):
    """
    calculate static pressure at 300 K from Hugoniot data using the constq
    formulation

    :param v: unit-cell volume in A^3
    :param v0: unit-cell volume in A^3 at 1 bar
    :param c0: velocity at 1 bar in km/s
    :param s: slope of the velocity change
    :param gamma0: Gruneisen parameter at 1 bar
    :param q: logarithmic derivative of Gruneisen parameter
    :param theta0: Debye temperature in K
    :param n: number of elements in a chemical formula
    :param z: number of formula unit in a unit cell
    :param mass: molar mass in gram
    :param c_v: heat capacity
    :param three_r: 3 times gas constant.
        Jamieson modified this value to compensate for mismatches
    :param t_ref: reference temperature, 300 K
    :return: static pressure in GPa
    :note: 2017/05/18 I am unsure if this is actually being used in pytheos
    """
    rho = mass / vol_uc2mol(v, z) * 1.e-6
    rho0 = mass / vol_uc2mol(v0, z) * 1.e-6
    p_h = hugoniot_p(rho, rho0, c0, s)
    p_th_h = jamieson_pth(v, v0, c0, s, gamma0, q, theta0, n, z, mass, c_v,
                          three_r=three_r, t_ref=t_ref)
    p_st = p_h - p_th_h
    return p_st