Page 335 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
314
APPENDIX C. 2D CARTESIAN CODE
COMMON/CPERIOD/DP1,DP2,IPERIOD,JPERIOD May 25, 2005 11:59
COMMON/CDUMT/DUM1(IT,JT),DUM2(IT,JT),DUM3(IT,JT)
COMMON/CRES/RESIU(5000),RESIV(5000),RESIM(5000)
C ADDITIONAL PROBLEM-DEPENDENT COMMON STATEMENTS
COMMON/CEVAP/IB1,IB2,JB1,GRM,SC,OBR,OWT,O(IT,JT),HH(IT,JT)
COMMON/CKRAL/D2,D1,U2,U1
COMMON/CPOROS/UU(IT,50),PRES(IT,50),PPP(IT,50),PPS(IT,50)
1 ,EPSI,RESIST,PIN

1D Porous Body Problem – Chapter 5
The USER file that follows shows how a fixed-pressure boundary condition (i.e.,

p = 0) is implemented in the ADSORB subroutine.
C **************************************
C THIS IS USER FILE POROS.FOR - PFIX BOUNDARY CONDITION
C **************************************
PROGRAM MAIN
INCLUDE ’COM2D.FOR’
OPEN(UNIT=6,FILE=’OO’)
WRITE(*,*)’-------- output is in OO file --------------’
C **** INITIAL DATA
WRITE(*,*)’GIVE ----- MXIT,IREAD,GAMMA’
READ(*,*)MXIT,IREAD,GAMMA
DX=1/FLOAT(IN-1)
X(1)=0.0
DO 1 I=2,IN
1 X(I)=X(I-1)+DX
Y(1)=0
Y(2)=0.5
Y(3)=1.0
INM=IN-1
JNM=JN-1
VISCOS=VISCOS/EPSI
DENSIT=DENSIT/EPSI**2
C
CALL MAINPR
STOP
END
C **************************************
SUBROUTINE TITLE
INCLUDE ’COM2D.FOR’
C **************************************