Page 360 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
0 521 85326 5
CB908/Date
0521853265appc
339
APPENDIX C. 2D CARTESIAN CODE
IF(N.EQ.2)SU(I,J)=ASNOW*U(I,J-1)+SU(I,J) May 25, 2005 11:59
IF(N.EQ.3)SU(I,J)=ASNOW*V(I,J-1)+SU(I,J)
IF(N.EQ.4)SU(I,J)=ASNOW*W(I,J-1)+SU(I,J)
IF(N.EQ.5)SU(I,J)=ASNOW*E(I,J-1)+SU(I,J)
IF(N.EQ.6)SU(I,J)=ASNOW*D(I,J-1)+SU(I,J)
IF(N.EQ.7)SU(I,J)=ASNOW*T(I,J-1)+SU(I,J)
SP(I,J)=ASNOW+SP(I,J)
ENDIF
C SYMMETRY
IF(LS.EQ.32)THEN
IF(N.EQ.3)SP(I,J)=ASNOW+SP(I,J)
AS(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(LS.EQ.33.OR.LS.EQ.35) THEN
RATIO=(Y(J)-YC(J))/DYMI(J+1)
IF(LS.EQ.33)RATIO=0.0
AS(I,J)=0.0
IF(N.EQ.2)U(I,J-1)=U(I,J)-RATIO*(U(I,J+1)-U(I,J))
IF(N.EQ.3)V(I,J-1)=V(I,J)-RATIO*(V(I,J+1)-V(I,J))
IF(N.EQ.4)W(I,J-1)=W(I,J)-RATIO*(W(I,J+1)-W(I,J))
IF(N.EQ.5)E(I,J-1)=E(I,J)-RATIO*(E(I,J+1)-E(I,J))
IF(N.EQ.6)D(I,J-1)=D(I,J)-RATIO*(D(I,J+1)-D(I,J))
IF(N.EQ.7)T(I,J-1)=T(I,J)-RATIO*(T(I,J+1)-T(I,J))
ENDIF
C WALL
IF(LS.EQ.34.OR.LS.EQ.36) THEN
AS(I,J)=0.0
DELTA=Y(J)-YC(J)
AREA=RC(J)*DXP(I)
UWAL=U(I,J-1)
VWAL=V(I,J-1)
WWAL=W(I,J-1)
VISWAL=VIS(I,J-1)
ANG=ATAN(W(I,J)/(U(I,J)+SMALL))
VT=(U(I,J)-UWAL)*COS(ANG)+(W(I,J)-WWAL)*SIN(ANG)