def equivalentConvection(h, eps, theta_inf, theta_s):
    """ Equivalent convection coefficient for radiation effect

	Parameters
	----------		
	h : float or (m,n) np.array
		convection coefficient [W/(m2K)]

	eps : float or (m,n) np.array
		emissivity/absorptivity [-]

	theta_inf : float or (m,n) np.array
		exposure temperature [°C]

	theta_s : float or (m,n) np.array
		surface temperature	

	Returns
	-------
	h_eq : float or (m,) np.array
		equivalent covection coefficient [W/(m2K)]

	Example use
	-----------
	>>> XXX	
	"""
    from scipy.constants import sigma as sig
    Tinf = theta_inf + 273.15
    Ts = theta_s + 273.15
    return h + eps * sig * (Tinf ** 2 + Ts ** 2) * (Tinf + Ts)
