Page 365 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
11:59
0 521 85326 5
May 25, 2005
0521853265appc
CB908/Date
344
X(IN)=XC(IN) APPENDIX C. 2D CARTESIAN CODE
ENDIF
C *** CALCULATE INTERPOLATION FACTORS
DXMI(1)=0.0
DO 4 I=2,IN
4 DXMI(I)=X(I)-X(I-1)
DYMI(1)=0.0
DO 5 J=2,JN
5 DYMI(J)=Y(J)-Y(J-1)
DO 6 J=1,JN
R(J)=1.0
RC(J)=1.0
IF(AXISYMM)THEN
R(J)=Y(J)
RC(J)=YC(J)
ENDIF
6 CONTINUE
C *** CALCULATE CELL VOLUME
SUMVOL=0.0
DO 7 J=2,JNM
DO 7 I=2,INM
VOL(I,J)=R(J)*(XC(I+1)-XC(I))*(YC(J+1)-YC(J))
SUMVOL=SUMVOL+VOL(I,J)
7 CONTINUE
WRITE(6,*)’ DOMAIN VOLUME = ’,SUMVOL
C *** CALCULATE AREAS
DO 9 I=1,INM
DXP(I)=(XC(I+1)-XC(I))
9 DXP(IN)=0.0
DO 10 J=1,JNM
DYP(J)=(YC(J+1)-YC(J))
10 DYP(JN)=0.0
RETURN
END
C *******************************************
SUBROUTINE COEF(NN,PRN,PRTN)
INCLUDE ’COM2D.FOR’
C *******************************************
N=NN
IF(N.EQ.1) GO TO 1000
C COEFFICIENTS OF TRANSPORT EQUATIONS
PRINV=1./PRN