def schomate(T):
	"""
	Shchomate equation to calculate specific heat of water vapor for a given temperature T.
	:param T: Temperature [K]
	:return: Heat Capacity [J/mol*K]
	"""
	t = T / 1000
	if 500 <= T < 1700:
		a, b, c, d, e = [30.092, 6.832514, 6.793425, -2.53448, 0.082139]
	elif T == 1700:
		return 2.7175
	elif 1700 < T <= 6000:
		a, b, c, d, e = [41.96426, 8.622053, -1.49978, 0.098119, -11.1576]
	else:
		raise ValueError('Input temperature outside valid domain')
	c_p = (a + (b * t) + (c * t**2) + (d * t**3) + (e / t**2)) / 18
	R = 0.4615
	c_v = c_p - R
	g = c_p / c_v
	return g