Page 339 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
11:59
0 521 85326 5
May 25, 2005
CB908/Date
0521853265appc
318
DO 91 I=1,IN APPENDIX C. 2D CARTESIAN CODE
UU(I,NITER)=U(I,2)
PPS(I,NITER)=PSM(I,2)/(P(I,2)+SMALL)
PRES(I,NITER)=P(I,2)
PPP(I,NITER)=PP(I,2)/(P(I,2)+SMALL)
91 CONTINUE
1000 CONTINUE
RETURN
END
C **************************************
SUBROUTINE OMEGA
INCLUDE ’COM2D.FOR’
C **************************************
C SPECIES EQUATION
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-06,3,2,55,0.5/
DATA CCTM,MXSTEP,DELT,STIME,MFREQ/1.0E-06,65,1.00,0.0,20/
C PP U V W E D T VIS P
DATA RP/1.0,0.95,0.5,0.5,0.5,1.0,1.0,1.0,0.95/
DATA NSWEEP/ 1,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/2*.TRUE.,5*.FALSE./
DATA BSOR/.TRUE.,.TRUE.,6*.FALSE.,.TRUE./
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/.FALSE.,.TRUE./