Gráfico da emissividade da radiação de corpo negro
Radiação de Corpo Negro, teoria disponível em: http://coral.ufsm.br/gef/Moderna/moderna02.pdf
Dados gerados: https://plot.ly/~zrhans/408/
Fortran
x
43
1
!---------------------------------------------------
2
! Calcula e gera arquivos com valores de Radiância
3
! em função do comprimento de onda para um
4
! Corpo Negro.
5
!---------------------
6
! Autor: Hans Zimermann
7
! Data: 09-10-2016 Alt:
8
!---------------------------------------------------
9
program curvas_de_radiancia
10
implicit none
11
12
real, parameter :: c1 = 3.74e-16 ! [W][m^2] ou [J/s][m^2]
13
real, parameter :: c2 = 1.44e-2 ! [m][K]
14
real :: lambda = 0. , T = 0.
15
real :: radiancia = 0.
16
integer :: i, N = 20 ! N numero de pontos da função
17
real :: dl
18
19
T =7000.
20
21
! Relação entre lambda e N contador incremento = (b - a) / n
22
dl = ( 2e-6 - 0.)/ N
23
24
open(unit=1, file="radiacia.txt", status="replace")
25
open(unit=2, file="radiacia-x.txt", status="replace")
26
do i = 1,N
27
lambda = lambda + dl
28
radiancia = (c1)*(lambda**(-5.0))*( 1 / (exp(c2/(lambda*T))-1.0) )
29
print 100,lambda, radiancia
30
write(1,100) radiancia
31
write(2,100) lambda
32
end do
33
100 format (es12.2)
34
35
close(1)
36
close(2)
37
stop('Programma finalizado!')
38
end program curvas_de_radiancia
39
40
! TO DO (para os alunos)
41
! otimizar a escrita em arquivos empregando arranjos
42
! usar laços para gerar radiância para faixas distintas de temperaturas
43
! Saída de dados deve ser semelhante a https://plot.ly/~zrhans/408/