Page 352 - Introduction to Computational Fluid Dynamics

P2: IWV
P1: ICD/GKJ
0 521 85326 5
CB908/Date
0521853265appc
331
APPENDIX C. 2D CARTESIAN CODE
DATA IB1,IB2,JB1,IN,JN/16,21,24,37,37/ May 25, 2005 11:59
Ch=2ANDH=8
DATA YC/0.0,0.0,0.02,0.04,0.07,0.12,0.18,0.25,0.35,0.5
1 ,0.65,0.8,0.95,1.1,1.25,1.4,1.55,1.7,1.8,1.85,1.9
1 ,1.95,1.98,2.0
1 ,2.02,2.05,2.1,2.2,2.4,2.8,3.3,4.0,5.0,6.0,7.0,7.5,8.0/
C L/2 = 8, t = 0.1, l/2=0.5
DATA XC/0.0,0.0,0.02,0.04,0.07,0.1,0.15,0.2,0.25,0.3,0.35,0.4
1 ,0.43,0.46,0.48,0.5
1 ,0.52,0.54,0.56,0.58,0.6
1 ,0.62,0.65,0.7,0.8,0.9,1.0,1.2,1.5,2.0,3.0,4.0,5.0
1 ,6.0,7.0,7.5,8.0/
END

Library File LIB2D.FOR
C *******************************************
C THIS IS LIBRARY FILE LIB2D.FOR ---- A W DATE
C RESIDUALS ARE STORED FOR PLOTTING IN SUBROUTINE EQN
C *******************************************
SUBROUTINE MAINPR
INCLUDE ’COM2D.FOR’
C *******************************************
NPERIOD=0
CALL TITLE
CALL GRID
CALL INITIA
CALL BSPEC
CALL INIT
CALL PROPS
CALL INDATA
CALL INFLUX
IF(IREAD) CALL IPT
IF(UNSTDY)CALL UPDATE
CALL EQN
CALL BOUNDP
IF(IWRITE)CALL OPT
CALL RESULT
RETURN
END
C *******************************************
SUBROUTINE INITIA
INCLUDE ’COM2D.FOR’