12 .......... N:      número de volumes de controle
    5.0d+0 ..... L:      comprimento do domínio [m]
    1.0d-3 ..... mi:     viscosidade absoluta [Pa.s]
    1.0d+3 ..... ro:     massa específica [kg/m3]
    2.0d-2 ..... f:      fator de atrito de Darcy
    1.0d+1 ..... Uin:    velocidade de entrada [m/s]
    2.0d-2 ..... D0:     diâmetro de entrada [m]
    4.0d-3 ..... cd:     fator de expansão do diâmetro
    1.0d+0 ..... dt:     passo de tempo
    0.0d+0 ..... beta:   coeficiente do termo difusivo (funções de interpolação):  0.0d0 = UDS / 1.0d0 = CDS
    1  .......... itmax:  número máximo de iterações (geral)
    1 .......... itmaxm: número máximo de iterações (ciclo da massa)