Page 353 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
332
APPENDIX C. 2D CARTESIAN CODE
C ******************************************* May 25, 2005 11:59
DO 1 J=1,JN
DO 1 I=1,IN
PP(I,J)=0.0
P(I,J)=0.0
U(I,J)=0.0
V(I,J)=0.0
W(I,J)=0.0
E(I,J)=0.0
D(I,J)=0.0
T(I,J)=0.0
QW(I,J)=0.0
VIS(I,J)=VISCOS
VIST(I,J)=0.0
RHO(I,J)=DENSIT
RHOO(I,J)=DENSIT
AW(I,J)=0.0
AE(I,J)=0.0
AS(I,J)=0.0
AN(I,J)=0.0
APU(I,J)=GREAT
APV(I,J)=GREAT
AP1(I,J)=GREAT
AP(I,J)=GREAT
NTAG(I,J)=0
NTAGW(I,J)=0
NTAGE(I,J)=0
NTAGS(I,J)=0
1 NTAGN(I,J)=0
RETURN
END
C *******************************************
SUBROUTINE TAG(CHAR1,CHAR2,IB,IL,JB,JL)
INCLUDE ’COM2D.FOR’
C *******************************************
CHARACTER*10 CHAR1,CHAR2
IF(CHAR2.EQ.’PERIOD’)NPERIOD=1
DO 1 J=JB,JL
DO 1 I=IB,IL
IF(CHAR1.EQ.’BLOCK’)THEN
NTAG(I,J)=1
GO TO 1