Page 342 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
CB908/Date
0 521 85326 5
0521853265appc
321
APPENDIX C. 2D CARTESIAN CODE
DATA INFLOW,EXIT1,SYMM,EXIT2,WALLT,WALLQ,PERIOD May 25, 2005 11:59
1 /’INFLOW’,’EXIT1’,’SYMM’,’EXIT2’,’WALLT’,’WALLQ’,’PERIOD’/
C ***** BLOCKED REGIONS
C CALL TAG(BLOCK,BLOCK,IB1,IB2-1,JB1,JNM)
C ***** DEFINESW&E BOUNDARIES
CALL TAG(WEST,INFLOW,2,2, 2,JB1-1)
CALL TAG(WEST,WALLQ,2,2,JB1,JNM)
CALL TAG(EAST,EXIT2,INM,INM,2,JNM)
C ***** DEFINES N&S BOUNDARIES
CALL TAG(NORTH,WALLQ,2,INM, JNM,JNM)
CALL TAG(SOUTH,SYMM,2,INM,2,2)
RETURN
END
C **************************************
SUBROUTINE RESULT
INCLUDE ’COM2D.FOR’
C **************************************
CHARACTER*20 HEADER
DO 1 J=1,JN
DO 1 I=1,IN
1 VIST(I,J)=VIST(I,J)/VISCOS

JSTEP=-1
WRITE(6,*)’ NITER=’,NITER
HEADER=’ U-VEL’
CALL PRINTK(U,1,IN,1,JN,HEADER,JSTEP)
HEADER=’ V-VEL’
CALL PRINTK(V,1,IN,1,JN,HEADER,JSTEP)
HEADER=’ PRESS’
CALL PRINTK(P,1,IN,1,JN,HEADER,JSTEP)
IF(TURBUL)THEN
HEADER=’ K ENERGY’
CALL PRINTK(E,1,IN,1,JN,HEADER,JSTEP)
HEADER=’ DISS’
CALL PRINTK(D,1,IN,1,JN,HEADER,JSTEP)
HEADER=’ MUT’
CALL PRINTK(VIST,1,IN,1,JN,HEADER,JSTEP)
ENDIF
IF(SLVE(7))THEN
HEADER=’ TEMP’
CALL PRINTK(T,1,IN,1,JN,HEADER,JSTEP)
ENDIF