docs: 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