Page 313 - Introduction to Computational Fluid Dynamics

P1: ICD
May 11, 2005
15:43
CB908/Date
0 521 85326 5
0521853265appb
292
TT=FLOAT(NTIME-1)*DELT APPENDIX B. 1D CONDUCTION CODE
TIME=(140-TOLD)/(TNOW-TOLD)*DELT+TT
WRITE(6,*)’ TIME FOR ADHESION = ’,TIME
RETURN
END
C *************************************************
BLOCK DATA
INCLUDE ’COM1D.FOR’
C *************************************************
C LOGICAL DECLARATIONS
DATA STEADY,UNSTEADY,GAUSS,THOMAS/.FALSE.,.TRUE.,.TRUE.,.FALSE./
C --------------------------------------------
C CONTROL PARAMETERS
C FULLY IMPLICIT (PSI=1),FULLY EXPLICIT (PSI=0),SEMI IMPLICIT (0<PSI<1)
DATA PSI,DELT,MXSTEP,ITERMX,RP,CC/0.0,10,10000,500,1.0,1E-5/
C --------------------------------------------
C BOUNDARY SPECIFICATION
DATA T1SPEC,Q1SPEC,H1SPEC/.TRUE.,2*.FALSE./
DATA TNSPEC,QNSPEC,HNSPEC/.TRUE.,2*.FALSE./
DATA T1,TN,QB1,QBN,HB1,HBN/250.0,250.0,0.0,0.0,0.0,0.0/
C DATA TINF,TINF1,TINFN,HPREF/25,150,250,12.0/
DATA CONDREF,RHOREF,SPHREF/0.25,1300,2000.0/
C --------------------------------------------
C GRID SPECIFICATION
DATA XCELL,XNODE/.TRUE.,.FALSE./
DATA N,AL/7,0.01/
DATA XCF/0.0,0.0,0.2,0.4,0.6,0.8,1.0,43*1.0/
END

USER File for Problem 2 – Chapter 2
C *************************************************
PROGRAM ONED
INCLUDE ’COM1D.FOR’
C *************************************************
OPEN(6,FILE=’OO’)
WRITE(6,*)’ **********************************************’
WRITE(6,*)’ RECTANGULAR FIN - PROB2-CHAPTER2’
WRITE(6,*)’ SOLVE BY GS AND TDMA’
WRITE(6,*)’ **********************************************’
DO 1 I=1,N