Page 317 - Introduction to Computational Fluid Dynamics

P1: ICD
May 11, 2005
15:43
CB908/Date
0 521 85326 5
0521853265appb
296
USER File for Problem 3 – Chapter 2 APPENDIX B. 1D CONDUCTION CODE
C *************************************************
C THIS IS USER FILE USER1D.FOR - A. W. DATE
C *************************************************
PROGRAM ONED
INCLUDE ’COM1D.FOR’
C *************************************************
OPEN(6,FILE=’OO’)
WRITE(6,*)’ **********************************************’
WRITE(6,*)’ ANNULAR COMPOSITE FIN - PROB3-CHAPTER2’
WRITE(6,*)’ SOLVE BY TDMA’
WRITE(6,*)’ **********************************************’
DX=(R3-R1)/FLOAT(N-2)
XCF(1)=0
XCF(2)=0.0
DO 1 I=3,N
1 XCF(I)=XCF(I-1) + DX
CALL MAIN
CALL RESULT
STOP
END
C *************************************************
SUBROUTINE INIT
INCLUDE ’COM1D.FOR’
C *************************************************
C GIVE INITIAL GUESS AT TIME=0.0 OR AT ITER=0 FOR STEADY STATE
T(1)=T1
T(N)=TN
RETURN
END
C *************************************************
SUBROUTINE NEWVAL
INCLUDE ’COM1D.FOR’
C *************************************************
C SET NEW VALUES OF HB1,HBN,QB1,QBN,TINF1,TINFN OR SOURCES
RETURN
END
C *************************************************
SUBROUTINE PROPS
INCLUDE ’COM1D.FOR’
C *************************************************