% Surface temperature from blackbody emissive flux
function T = InvEblk(Eb) % Eb=W/m^2 and T=K
  sig=5.6704e-8;
  T = (Eb/sig).^0.25;
end
