Page 362 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0521853265appc
341
APPENDIX C. 2D CARTESIAN CODE
ENDIF 0 521 85326 5 May 25, 2005 11:59
ENDIF
ENDIF
C *** NORTH BOUNDARY
300 LN=NTAGN(I,J)
IF(LN.EQ.0)GO TO 1
ANNOW=AN(I,J)
C INLET
IF(LN.EQ.41.OR.LN.EQ.47)THEN
AN(I,J)=0.0
IF(N.EQ.2)SU(I,J)=ANNOW*U(I,J+1)+SU(I,J)
IF(N.EQ.3)SU(I,J)=ANNOW*V(I,J+1)+SU(I,J)
IF(N.EQ.4)SU(I,J)=ANNOW*W(I,J+1)+SU(I,J)
IF(N.EQ.5)SU(I,J)=ANNOW*E(I,J+1)+SU(I,J)
IF(N.EQ.6)SU(I,J)=ANNOW*D(I,J+1)+SU(I,J)
IF(N.EQ.7)SU(I,J)=ANNOW*T(I,J+1)+SU(I,J)
SP(I,J)=ANNOW+SP(I,J)
ENDIF
C SYMMETRY
IF(LN.EQ.42)THEN
IF(N.EQ.3)SP(I,J)=ANNOW+SP(I,J)
AN(I,J)=0.0
IF(N.EQ.3)V(I,J+1)=0.0
IF(N.EQ.2)U(I,J+1)=U(I,J)
IF(N.EQ.4)W(I,J+1)=W(I,J)
IF(N.EQ.5)E(I,J+1)=E(I,J)
IF(N.EQ.6)D(I,J+1)=D(I,J)
IF(N.EQ.7)T(I,J+1)=T(I,J)
ENDIF
C EXIT
IF(LN.EQ.43.OR.LN.EQ.45) THEN
AN(I,J)=0.0
RATIO=(YC(J+1)-Y(J))/DYMI(J)
IF(LN.EQ.43)RATIO=0.0
IF(N.EQ.2)U(I,J+1)=U(I,J)+RATIO*(U(I,J)-U(I,J-1))
IF(N.EQ.3)V(I,J+1)=V(I,J)+RATIO*(V(I,J)-V(I,J-1))
IF(N.EQ.4)W(I,J+1)=W(I,J)+RATIO*(W(I,J)-W(I,J-1))
IF(N.EQ.5)E(I,J+1)=E(I,J)+RATIO*(E(I,J)-E(I,J-1))
IF(N.EQ.6)D(I,J+1)=D(I,J)+RATIO*(D(I,J)-D(I,J-1))
IF(N.EQ.7)T(I,J+1)=T(I,J)+RATIO*(T(I,J)-T(I,J-1))
ENDIF
C WALL