Page 322 - Introduction to Computational Fluid Dynamics

P1: ICD
0 521 85326 5
CB908/Date
0521853265appb
301
APPENDIX B. 1D CONDUCTION CODE
IF(XCELL)THEN May 11, 2005 15:43
XCF(2)=XCF(1)
X(1)=XCF(1)
DO 1 I=2,N-1
1 X(I)=0.5*(XCF(I)+XCF(I+1))
X(N)=XCF(N)
ELSE
XCF(1)=X(1)
XCF(2)=X(1)
DO 2 I=3,N-1
2 XCF(I)=0.5*(X(I)+X(I-1))
XCF(N)=X(N)
ENDIF
C CALCULATE PERIMETER,CELL-FACE AREA AND CELL VOLUME
C AREA AND PERI ARE FUNCTION ROUTINES (USER FILE)
DO 3 I=1,N
ACF(I)=AREA(XCF(I))
PERIM(I)=PERI(X(I))
3 CONTINUE
DO 4 I=2,N-1
4 VOL(I)=AREA(X(I))*(XCF(I+1)-XCF(I))
HEADER=’ X(I) ’
CALL PRINT(X,HEADER)
HEADER=’ XCF(I) ’
CALL PRINT(XCF,HEADER)
HEADER=’ CELL FACE AREA ’
CALL PRINT(ACF,HEADER)
HEADER=’ PERIMETER ’
CALL PRINT(PERIM,HEADER)
HEADER=’ CELL-VOLUME ’
CALL PRINT(VOL,HEADER)
RETURN
END
C *************************************************
SUBROUTINE COEF
INCLUDE ’COM1D.FOR’
C *************************************************
DO 1 I=2,N-1
C INITIALISE SU ANS SP
STAB(I)=0.0
SU(I)=0.0
SP(I)=0.0