Page 351 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
11:59
May 25, 2005
CB908/Date
0521853265appc
330
ENDIF 0 521 85326 5 APPENDIX C. 2D CARTESIAN CODE
3 CONTINUE
RPO=1.0
DO 4 J=2,JNM
DO 4 I=2,INM
SUM=AW(I,J)+AE(I,J)+AS(I,J)+AN(I,J)
4 AP(I,J)=(SUM+SP(I,J))/RPO
C SOLVE THE O EQN
CALL SOLVE(O,RPO,RSU)
DO 5 J=1,JN
DO 5 I=1,IN
IF(O(I,J).GT.1.0)O(I,J)=1.0
5 IF(O(I,J).LT.0.0)O(I,J)=0.0
RETURN
END
C **************************************
BLOCK DATA
INCLUDE ’COM2D.FOR’
C **************************************
C **** INITIAL DATA(make sure that IN,JN equal IT,JT)

DATA CC,IPREF,JPREF,MXIT,GAMMA/1.0E-05,10,10,55,0.5/
DATA CCTM,MXSTEP,DELT,STIME,MFREQ/1.0E-05,65,1.00,0.0,20/
C PP U V W E D T VIS P
DATA RP/1.0,0.5,0.5,1.0,1.0,1.0,1.0,1.0,0.1/
DATA NSWEEP/ 10,1,1,1,1,1,1 /
DATA PR/1.0,1.0,1.0,1.0,1.0,1.0,0.7 /
DATA PRT/1.0,1.0,1.0,1.0,1.0,1.3,0.9 /
DATA DENSIT,VISCOS,SPHEAT/1.0,1.0,1.0/
DATA CD1,CD2,CMU,ELOG,CAPPA/1.44,1.92,0.09,9.793,0.4187/
C **** LOGICAL DATA
DATA TURBUL, STEADY, UNSTDY, FTRAN , CONMAS, AXISYMM
1 /.FALSE.,.TRUE., .FALSE., .FALSE., .FALSE., .FALSE./
DATA UPWIND,HYBRID,POWER/.TRUE.,.FALSE.,.FALSE./
DATA SLVE/3*.TRUE.,4*.FALSE./
DATA BSOR/2*.FALSE.,.TRUE.,6*.FALSE./
DATA IREAD,IWRITE/.FALSE. ,.TRUE./
C PERIODIC BC
DATA IPERIOD,JPERIOD/0,0/
DATA DP1,DP2/0.0,0.0/
C **** READ GRID DATA
DATA GRCELL,GRNODE/.TRUE.,.FALSE./