Page 315 - Introduction to Computational Fluid Dynamics
P. 315

P1: ICD
                                                                                                 15:43
                                                                                   May 11, 2005
            0521853265appb
                     294
                                   CONTINUE
                            1 CB908/Date  0 521 85326 5           APPENDIX B. 1D CONDUCTION CODE
                                   RETURN
                                   END
                            C *************************************************
                                   SUBROUTINE INTPRI
                                   INCLUDE ’COM1D.FOR’
                            C *************************************************
                                   RETURN
                                   END
                            C *************************************************
                            C       FUNCTION ROUTINES
                            C *************************************************
                                   FUNCTION HPERI(II)
                                   INCLUDE ’COM1D.FOR’
                            C H AT PERIMETER
                                   I=II
                                   HPERI=HPREF+X(I)*0.0+T(I)*0.0
                                   RETURN
                                   END
                            C --------------------------------------------
                                   FUNCTION AREA(XX)
                                   INCLUDE ’COM1D.FOR’
                            C AREA OF CROSS-SECTION
                                   AREA=BREADTH*THICK+0.0*XX
                                   RETURN
                                   END
                            C --------------------------------------------
                                   FUNCTION PERI(XX)
                                   INCLUDE ’COM1D.FOR’
                            C PERIMETER
                                   PERI=2*BREADTH+0.0*XX
                                   RETURN
                                   END
                            C *************************************************
                                   SUBROUTINE RESULT
                                   INCLUDE ’COM1D.FOR’
                                   CHARACTER*20 HEADER
                            C *************************************************
                                   HEADER=’ FINAL-TEMP ’
                                   CALL PRINT(T,HEADER)
                            C EXTRACT PROBLEM DEPENDENT PARAMETERS IF ANY
   310   311   312   313   314   315   316   317   318   319   320