Page 310 - Introduction to Computational Fluid Dynamics

P1: ICD
0 521 85326 5
May 11, 2005
CB908/Date
0521853265appb
289
APPENDIX B. 1D CONDUCTION CODE
COMMON/CTRAN/DELT,TIMEMX,MXSTEP,PSI,ITER,NTIME,TTIME 15:43
COMMON/CPROPS/CONDREF,RHOREF,SPHREF
COMMON/CDAT1/T1,TN,QB1,QBN,HB1,HBN,TINF1,TINFN,HPREF,TINF
COMMON/CDAT1/QB1O,QBNO,HB1O,HBNO,TINF1O,TINFNO,HPREFO,TINFO
COMMON/CDUM/DUM1(5000),DUM2(5000),DUM3(5000)
C ADDITIONAL PROBLEM-DEPENDENT VARIABLES
C VARIABLES FOR PROB2
COMMON/CP2/BREADTH,THICK
C VARIABLES FOR PROB3
COMMON/CRADS/R1,R2,R3

USER File for Problem 1 – Chapter 2
C *************************************************
PROGRAM ONED
INCLUDE ’COM1D.FOR’
C *************************************************
OPEN(6,FILE=’OO’)
WRITE(6,*)’ **********************************************’
WRITE(6,*)’ ADHESION OF PLASTIC SHEETS - PROB1-CHAPTER2’
WRITE(6,*)’ **********************************************’
DO 1 I=1,N
1 XCF(I)=XCF(I)*AL
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
TIN=30
DO 1 I=1,N
T(I)=30
IF(I.EQ.1.OR.I.EQ.N)T(I)=250
1 CONTINUE
RETURN
END
C *************************************************
SUBROUTINE NEWVAL
INCLUDE ’COM1D.FOR’
C *************************************************