Page 345 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
324
APPENDIX C. 2D CARTESIAN CODE
C ************************************** May 25, 2005 11:59
SUBROUTINE OMEGA
INCLUDE ’COM2D.FOR’
C **************************************
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,1.0,1.0,1.0,1.0,1.0,1.0,0.5,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,3.0 /
DATA PRT/1.0,1.0,1.0,1.0,1.0,1.3,0.9 /
DATA DENSIT,VISCOS,SPHEAT/1.0,0.5E-3,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 /.TRUE.,.TRUE., .FALSE., .TRUE., .TRUE., .TRUE./
DATA UPWIND,HYBRID,POWER/.TRUE.,.FALSE.,.FALSE./
DATA SLVE/3*.TRUE.,.FALSE.,2*.TRUE.,.FALSE./
DATA BSOR/9*.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./
DATA IN,JN/67,24/
DATA JB1/14/
DATA YC/0.0,0.0,0.02,0.05,0.08,0.1,0.15,0.2,0.25
1 ,0.30,0.35,0.40,0.45,0.5
1 ,0.55,0.6,0.65,0.70,0.75
1 ,0.8,0.85,0.9,0.95,1.0/
DATA XC/0.0,0.0,0.03,0.06,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45
1 ,0.50,0.55,0.6,0.65,0.7,0.75,0.8,0.85,0.9,0.95,1.0,1.05
1 ,1.12,1.2,1.3,1.4,1.5,1.6,1.75,1.9,2.05,2.2,2.4,2.6,2.8