O programa computacional PROG6_CFD1 resolve
 * convecção de calor bidimensional em regime permanente
   com o campo de velocidades conhecido (prescrito)
 * condições de contorno de Dirichlet
 * propriedade constante: difusividade térmica
 * malha uniforme em cada direção (x,y)
 * a variável T média não deve ser considerada
 * os perfis de T são corretos apenas quando
   se entra com número de volumes ímpares em cada
   direção
 * método numérico: volumes finitos
 * funções de interpolação: 
      CDS para os termos difusivos
      UDS ou CDS (com correção adiada) para os termos advectivos
 * solvers: Jacobi, Gauss-Seidel, SOR, ADI e MSI-5
 * precisão: dupla
 * com volumes fictícios
 * campo de velocidades é prescrito (fixo)
 * domínio de cálculo: Lx = Ly = 1
 * a solução analítica refere-se apenas à condução de calor
 * nos dados, usar sempre "número de malhas" = 1

Autor: Carlos Henrique Marchi
       DEMEC/UFPR, Curitiba, PR
       Setembro de 2003

Seqüência para compilação dos módulos
pela primeira vez após criar o projeto:
1) VARIA.f90
2) SOLVERS.f90
3) DADOS.f90
4) COEF.f90
5) RESULT.f90
6) PROG6_CFD1.f90