	function gr_goes_tem, f4, f8, flux17 = flux17

; Temperature in 10^6 K, Emission Measure in 10^45 cm-3

R = double(f4/f8)

T = 3.15 + 77.2*R + 164*R^2 + 205*R^3

b8 = -3.86 + 1.17*T - 1.31d-2*T^2 + 1.78d-4*T^3

EM = F8/(b8)*1d55

flux17 = 3.1e7*em/1d55/sqrt(t)

return, [[float(T)], [float(EM/1d45)]]

	end