Page 326 - Introduction to Computational Fluid Dynamics

P1: ICD
CB908/Date
0521853265appb
305
APPENDIX B. 1D CONDUCTION CODE
QBNO=QBN 0 521 85326 5 May 11, 2005 15:43
HB1O=HB1
HBNO=HBN
TINF1O=TINF1
TINFNO=TINFN
HPREFO=HPREF
TINFO=TINF
RETURN
END
C *************************************************
SUBROUTINE PRINT(F,HEADER)
INCLUDE ’COM1D.FOR’
CHARACTER*20 HEADER
C *************************************************
DIMENSION F(IT)
WRITE(6,*)’*****************************************’
WRITE(6,*)’DISTRIBUTION OF ’,HEADER
IB=1
IE=IB+6
IF(IE.GT.N)IE=N
100 CONTINUE
WRITE(6,500)(F(I),I=IB,IE)
WRITE(6,600)(I,I=IB,IE)
IF(IE.LT.N) THEN
IB=IE+1
IE=IB+6
IF(IE.GT.N)IE=N
GO TO 100
ENDIF
WRITE(6,*)’*****************************************’
500 FORMAT(7E10.3)
600 FORMAT(4X,I3,6I10)
RETURN
END