program Gibbs_1p3_main

! Cálculo de propriedades para reações químicas em equilíbrio e não-equilíbrio
! para o propelente H2/O2

! Implementação da versão atual (1.3):
!      Programador: Carlos Henrique Marchi
!      Local:       LENA/DEMEC/UFPR, Curitiba, PR
!      Período:     14 de janeiro a 24 de fevereiro de 2005
!      Última alteração: 3 de maio de 2005

! Implementação da versão atual (1.2):
!      Programador: Carlos Henrique Marchi
!      Local:       LENA/DEMEC/UFPR, Curitiba, PR
!      Período:     10 a 13 de janeiro de 2005
!      Última alteração: 13 de janeiro de 2005

! Implementação da versão atual (1.1):
!      Programador: Carlos Henrique Marchi
!      Local:       LENA/DEMEC/UFPR, Curitiba, PR
!      Período:     7 a 7 de janeiro de 2005
!      Última alteração: 7 de janeiro de 2005

! Implementação da versão inicial (1.0):
!      Programador: Carlos Henrique Marchi
!      Local:       LENA/DEMEC/UFPR, Curitiba, PR
!      Período:     4 a 6 de janeiro de 2005

!-------------------------------------------------------------------------------

! MODELO NUMÉRICO (resumo):
!   Até 8 espécies químicas disponíveis:
!     1 = H2O   2 = O2   3 = H2    4 = OH   
!     5 = O     6 = H    7 = HO2   8 = H2O2
!   Propriedades termodinâmicas e de trasnporte são calculadas conforme
!     NASA TM-4513 (1993)
!   O equilíbrio químico é resolvido com o método das constantes de	equilíbrio, 
!     que necessita definir as reações e as espécies envolvidas, e o método de
!     Newton simples para resolver as taxas de dissociação
!   A temperatura de combustão é resolvida com o método da bisseção
!   Linguagem Fortran
!   Aplicativo usado: Microsoft Developer Studio, Fortran PowerStation 4.0
!   Tipo de projeto: Console Application
!   O programa permite resolver 5 tipos de problemas:
!     1 = TERMO: calcular propriedades termodinâmicas e de transporte de cada
!                espécie química, na base mássica ou molar
!     2 = CONGELADO: calcular propriedades termodinâmicas e de transporte de
!                    uma mistura de espécies químicas
!     3 = EQUILIBRIO: calcular a composição de uma mistura de gases;
!                     há 9 modelos disponíveis, com 0 a 18 reações de
!                     dissociação, envolvendo 3 a 8 espécies químicas
!     4 = CÂMARA: calcular a temperatura de combustão e a composição de uma
!                 mistura gases; há 9 modelos disponíveis, com 0 a 18 reações de
!                 dissociação, envolvendo 3 a 8 espécies químicas
!     5 = TAXA FINITA: calcular a taxa de geração de massa por unidade de
!                      volume de cada espécie química;
!                      há 6 modelos disponíveis, com 4 a 18 reações de
!                      dissociação, envolvendo 6 ou 8 espécies químicas

! ARQUIVOS envolvidos no programa:
!   Gibbs_1p3_modulo.f90   = sub-rotinas funcionais
!   Gibbs_1p3_main.f90     = programa principal
!   Gibbs_coeficientes_cp_h_g.dat = coeficientes para propriedades termodinâmicas
!   Gibbs_coeficientes_ki_mi.dat  = coeficientes para propriedades de transporte
!   Gibbs_1p3.in              = arquivo de dados principal
!   Gibbs_1p3_termo.in        = arquivo de dados do bloco TERMO
!   Gibbs_1p3_congelado.in    = arquivo de dados do bloco CONGELADO
!   Gibbs_1p3_equilibrio.in   = arquivo de dados do bloco EQUILIBRIO
!   Gibbs_1p3_camara.in       = arquivo de dados do bloco CAMARA
!   Gibbs_1p3_taxa_finita.in  = arquivo de dados 1 do bloco TAXA FINITA
!   Gibbs_1p3_taxa_finita2.in = arquivo de dados 2 do bloco TAXA FINITA
!   *.*         = arquivo de saída do programa (nome definido pelo usuário)
!   Notepad.exe = editor de arquivos
!   Gibbs_1p3_LEIA_ME.txt = informações básicas sobre o programa
!   Gibbs_1p3.exe         = programa executável

!-------------------------------------------------------------------------------

Seqüência para compilação dos módulos pela primeira vez após criar o projeto:

1) Gibbs_1p3_modulo.f90
2) Gibbs_1p3_main.f90