begin # Begin blocking entire file for automatic inclusion function blackbody(Tkelvin, descrip, LamLoHi, vec) begin global LamDomain LamLoHi = LamDomain # format double "f8.0" descrip = [ntoa(Tkelvin), " K blackbody"] NumData = LamLoHi(2) - LamLoHi(1) + 1 c1 = 3.74177107d-16 # watts/m^2 c2 = 1.4387752d7 # (m K) * (10^9 nm/m), SI value has 10^-2 vec = zeros(NumData) Index = 0 for lambda = LamLoHi(1) to LamLoHi(2) begin Index = Index+1 LamWork = lambda*1.0d-9 vec(Index) = c1*(LamWork^-5)/(exp(c2/(Tkelvin*lambda))-1) end # End for lambda = etc. # Normalization = max( vec ) Normalization = vec( 555 - LamLoHi(1) + 1 ) # format double "g15.4" vec = vec/(Normalization * 0.01) end # End function blackbody() end # End blocking entire file for automatic inclusion