Page 336 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
May 25, 2005
0 521 85326 5
CB908/Date
0521853265appc
315
APPENDIX C. 2D CARTESIAN CODE
WRITE(6,*)’*************************************************’ 11:59
WRITE(6,*)’ PROGRAM TO CALCULATE POROS MEDIUM FLOW ’
WRITE(6,*)’*************************************************’
RETURN
END
C **************************************
SUBROUTINE INIT
INCLUDE ’COM2D.FOR’
C **************************************
C INITIAL GUESS
PIN=RESIST
DO 1 I=1,IN
DO 1 J=1,JN
C P(I,J)=PIN*(1-X(I)/X(IN))
1 CONTINUE
P(1,2)=PIN
RETURN
END
C **************************************
SUBROUTINE BSPEC
INCLUDE ’COM2D.FOR’
C **************************************
C **** PROVIDE BOUNDARY & BLOCKED REGIONS
C
CHARACTER*10 BLOCK,WEST,EAST,SOUTH,NORTH
CHARACTER*10 INFLOW,EXIT1,SYMM,EXIT2,WALLT,WALLQ,PERIOD
DATA BLOCK,WEST,EAST,SOUTH,NORTH
1 /’BLOCK’,’WEST’,’EAST’,’SOUTH’,’NORTH’/
DATA INFLOW,EXIT1,SYMM,EXIT2,WALLT,WALLQ,PERIOD
1 /’INFLOW’,’EXIT1’,’SYMM’,’EXIT2’,’WALLT’,’WALLQ’,’PERIOD’/
C ***** BLOCKED REGIONS
C CALL TAG(BLOCK,BLOCK,2,IB1-1,JB2,JNM)
C CALL TAG(BLOCK,BLOCK,IB1,INM,2,JB1-1)
C ***** DEFINESW&E BOUNDARIES
CALL TAG(WEST,EXIT1,2,2, 2,JNM)
CALL TAG(EAST,EXIT1,INM,INM,2,JNM)
C ***** DEFINES N&S BOUNDARIES
CALL TAG(NORTH,SYMM,2,INM,JNM,JNM)
CALL TAG(SOUTH,SYMM,2,INM,2,2)
RETURN
END