Page 355 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
11:59
0 521 85326 5
May 25, 2005
CB908/Date
0521853265appc
334
DO 1 J=2,JNM APPENDIX C. 2D CARTESIAN CODE
DO 1 I=2,INM
VOLP=VOL(I,J)
RHOP=RHO(I,J)
C *** BLOCKED REGION
IF(NTAG(I,J).EQ.1)THEN
IF(N.EQ.2)SU(I,J)=GREAT*U(I,J)
IF(N.EQ.3)SU(I,J)=GREAT*V(I,J)
IF(N.EQ.4)SU(I,J)=GREAT*W(I,J)
IF(N.EQ.5)SU(I,J)=GREAT*E(I,J)
IF(N.EQ.6)SU(I,J)=GREAT*D(I,J)
IF(N.EQ.7)SU(I,J)=GREAT*T(I,J)
SP(I,J)=GREAT
GO TO 1
END IF
C *** WEST BOUNDARY
LW=NTAGW(I,J)
IF(LW.EQ.0)GO TO 100
AWNOW=AW(I,J)
C INLET
IF(LW.EQ.11.OR.LW.EQ.17)THEN
AW(I,J)=0.0
IF(N.EQ.2)SU(I,J)=AWNOW*U(I-1,J)+SU(I,J)
IF(N.EQ.3)SU(I,J)=AWNOW*V(I-1,J)+SU(I,J)
IF(N.EQ.4)SU(I,J)=AWNOW*W(I-1,J)+SU(I,J)
IF(N.EQ.5)SU(I,J)=AWNOW*E(I-1,J)+SU(I,J)
IF(N.EQ.6)SU(I,J)=AWNOW*D(I-1,J)+SU(I,J)
IF(N.EQ.7)SU(I,J)=AWNOW*T(I-1,J)+SU(I,J)
SP(I,J)=AWNOW+SP(I,J)
ENDIF
C SYMMETRY
IF(LW.EQ.12)THEN
IF(N.EQ.2)SP(I,J)=AWNOW+SP(I,J)
AW(I,J)=0.0
IF(N.EQ.2)U(I-1,J)=0.0
IF(N.EQ.3)V(I-1,J)=V(I,J)
IF(N.EQ.4)W(I-1,J)=W(I,J)
IF(N.EQ.5)E(I-1,J)=E(I,J)
IF(N.EQ.6)D(I-1,J)=D(I,J)
IF(N.EQ.7)T(I-1,J)=T(I,J)
ENDIF
C EXIT