function Te = tempr(g,w)
% Local effective temperature for a rotating star.
%
% D. Holmgren holmgren@phobos.astro.uwo.ca

a = zeros(10,1);
a(1) = 1; a(4) = 0.7770854; a(5) = -0.5863339;
a(6) = -0.1162862; a(10) = 0.4589665;
Cn = 0;
for k = 1:10
 Cn = Cn + a(k)*w^(k-1);
end
% need stellar parameters here...
L = 7.0281e29; G = 6.67e-11; M = 1.5913e31; sigma = 5.67e-08;
Cw = Cn*L/(4*pi*G*M*sigma);
Te = (Cw*g)^0.25;