Page 312 - Introduction to Computational Fluid Dynamics

P1: ICD
0 521 85326 5
CB908/Date
0521853265appb
291
APPENDIX B. 1D CONDUCTION CODE
C ************************************************* May 11, 2005 15:43
C FUNCTION ROUTINES
C *************************************************
FUNCTION HPERI(II)
INCLUDE ’COM1D.FOR’
C H AT PERIMETER
I=II
HPERI=HPREF*0.0+X(I)*0.0+T(I)*0.0
RETURN
END
C --------------------------------------------
FUNCTION AREA(XX)
INCLUDE ’COM1D.FOR’
C AREA OF CROSS-SECTION
AREA=1.0+0.0*XX
RETURN
END
C --------------------------------------------
FUNCTION PERI(XX)
INCLUDE ’COM1D.FOR’
C PERIMETER
PERI=0*XX
RETURN
END
C *************************************************
SUBROUTINE RESULT
INCLUDE ’COM1D.FOR’
CHARACTER*20 HEADER
C *************************************************
HEADER=’ FINAL-TEMP ’
CALL PRINT(T,HEADER)
HEADER=’ X(I) ’
CALL PRINT(X,HEADER)
HEADER=’ XCF(I) ’
CALL PRINT(XCF,HEADER)
C EXTRACT PROBLEM DEPENDENT PARAMETERS IF ANY
WRITE(6,*)’ PRINT MID-POINT TEMPERATURE’
DO 1 I=1,NTIME
TT=FLOAT(I)*DELT
1 WRITE(6,*)TT,DUM1(I)
TNOW=DUM1(NTIME)
TOLD=DUM1(NTIME-1)