def p3sat_s(s):
    """function p3sat_s = p3sat_s(s)

    Section 4.2 Region 3. pSat_h  & pSat_s
    """
    Ii = [0, 1, 1, 4, 12, 12, 16, 24, 28, 32]
    Ji = [0, 1, 32, 7, 4, 14, 36, 10, 0, 18]
    ni = [0.639767553612785, -12.9727445396014, -2245951258484030.0, 
        1774667.41801846, 7170793495.71538, -3.78829107169011e+17, -
        9.55586736431328e+34, 1.87269814676188e+23, 119254746466.473, 
        1.10649277244882e+36]
    Sigma = s / 5.2
    Pi = 0
    for i in range(0, 10):
        Pi = Pi + ni[i] * (Sigma - 1.03) ** Ii[i] * (Sigma - 0.699) ** Ji[i]
    return Pi * 22
