Lignes permettant de calculer le Rayleigh et le Nusselt


C * ------------------- SECTION 7 ----

Finish of sweep.

c calcul du module de la vitesse

CONTINUE OPEN(50,FILE='nusselt')

L0H1=L0F(H1)

L0YG2D=L0F(YG2D)

L0DXGD=L0F(DXG2D)

FLUXDD=0.

IY=1

DO 1973 IX=1,NX

ICELL=IY+(IX-1)*NY

FLUXDD=FLUXDD+(F(L0H1+ICELL)-RG(5))*F(L0DXGD+ICELL)*RG(2)*

& RG(3)/(RG(4)*F(L0YG2D+ICELL))

1973 CONTINUE

ALAMBD=RG(2)*RG(3)*RG(6)/RG(4)

ANu=(FLUXDD/XULAST)/(ALAMBD*(RG(7)-RG(5))/RG(6))*YVLAST

WRITE(50,*)'Le nombre de Nussel est : ',ANu

WRITE(50,*)'Le nombre de Rayleigh est : ',RG(1)

CLOSE(50)

RETURN

198CONTINUE

C * ------------------- SECTION 8 ---- Finish of time step.