Page 314 - Introduction to Computational Fluid Dynamics

P1: ICD
0521853265appb
293
APPENDIX B. 1D CONDUCTION CODE
XCF(I)=XCF(I)*AL
1 CB908/Date 0 521 85326 5 May 11, 2005 15:43
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
BB=(TN-T1)/AL
DO 1 I=1,N
1 T(I)=T1+BB*X(I)
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 *************************************************
C COND(I) AND SPH(I) ARE DEFINED AT NODE P
DO 1 I=1,N
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
TERM=HPERI(I)*PERIM(I)*(XCF(I+1)-XCF(I))
SU(I)=SU(I)+TERM*TINF
SP(I)=SP(I)+TERM