Page 347 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
May 25, 2005
0 521 85326 5
CB908/Date
0521853265appc
326
APPENDIX C. 2D CARTESIAN CODE
WRITE(6,*)’***************************************’ 11:59
WRITE(6,*)’ WATER EVAPORATION BY NATURAL CONVECTION ’
WRITE(6,*)’***************************************’
RETURN
END
C **************************************
SUBROUTINE INIT
INCLUDE ’COM2D.FOR’
C **************************************
C INITIAL GUESS (ONLY HEAT FLUX NEEDS TO BE SPECIFIED)
DO 10 I=1,IN
DO 10 J=1,JN
O(I,J)=0.0
IF(I.LT.IB1)O(I,1)=1.0
10 CONTINUE
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
CALL TAG(BLOCK,BLOCK,IB1,IB2-1,2,JB1-1)
C ***** DEFINESW&E BOUNDARIES
CALL TAG(WEST,SYMM,2,2, 2,JNM)
CALL TAG(WEST,WALLQ,IB2,IB2, 2,JB1-1)
CALL TAG(EAST,WALLQ,INM,INM,2,JNM)
CALL TAG(EAST,WALLQ,IB1-1,IB1-1,2,JB1-1)
C ***** DEFINES N&S BOUNDARIES
CALL TAG(NORTH,WALLQ,2,INM, JNM,JNM)
CALL TAG(SOUTH,WALLQ,2,IB1-1,2,2)
CALL TAG(SOUTH,WALLT,IB2,INM,2,2)
CALL TAG(SOUTH,WALLT,IB1,IB2-1,JB1,JB1)