Page 311 - Introduction to Computational Fluid Dynamics

P1: ICD
0 521 85326 5
May 11, 2005
CB908/Date
0521853265appb
290
APPENDIX B. 1D CONDUCTION CODE
C SET NEW VALUES OF HB1,HBN,QB1,QBN,TINF1,TINFN OR SOURCES 15:43
RETURN
END
C *************************************************
SUBROUTINE PROPS
INCLUDE ’COM1D.FOR’
C *************************************************
C COND(I) AND SPH(I) ARE DEFINED AT NODE P
DO 1 I=1,N
RHO(I)=RHOREF
COND(I)=CONDREF
1 SPH(I)=SPHREF
RETURN
END
C *************************************************
SUBROUTINE SORCE
INCLUDE ’COM1D.FOR’
C *************************************************
C FORM PROBLEM DEPENDENT SOURCE TERM INCLUDING SU AND SP
DO 1 I=2,N-1
SU(I)=SU(I)+0.0
1 CONTINUE
RETURN
END
C *************************************************
SUBROUTINE INTPRI
INCLUDE ’COM1D.FOR’
CHARACTER*20 HEADER
C *************************************************
WRITE(6,*)’ TIMESTEP = ’,NTIME,’ TOTAL TIME = ’,TTIME
C PRINT TEMPERATURES AT THE CURRENT STEP
HEADER=’ TEMP ’
CALL PRINT(T,HEADER)
C STORE MID-POINT TEMPERATURE
DUM1(NTIME)=T(4)
C GIVE STOP CONDITION
IMID=4
IF(T(IMID).GT.140)ISTOP=1
RETURN
END