Page 334 - Introduction to Computational Fluid Dynamics
P. 334

P2: IWV
            P1: ICD/GKJ
                            CB908/Date
                                         0 521 85326 5
            0521853265appc
                                                                                                    313
                        APPENDIX C. 2D CARTESIAN CODE
                        Subroutines PR1D This routine is used to print out 1D variables.  May 25, 2005  11:59
                        Function STAN This function routine is called from subroutine BOUND to im-
                        plement a wall-function boundary condition for temperature and mass-fraction
                        variables in a turbulent flow. The routine evaluates the Stanton number based on
                        specification of PF, which the USER may change if required. For mass-fraction
                        variables, the function is called from routine OMEGA.
                        Functions FINTW, FINTE, FINTS, and FINTN These function routines evalu-
                        ate variable   at cell faces, w, e, s, and n, respectively, using linear interpolation.


                        COMMON BLOCK COM2D.FOR
                        C*******************************
                        C THIS IS COM2D.FOR
                        C ******************************
                        C IT AND JT CHANGE WITH THE PROBLEM
                                 PARAMETER(IT=37,JT=37,MXGR=37)
                                 PARAMETER(GREAT=1.0E+20,SMALL=1.0E-20,PI=3.1415926)
                                 LOGICAL TURBUL,STEADY,UNSTDY,FTRAN,CONMAS,AXISYMM,BSOR
                                 LOGICAL UPWIND,HYBRID,POWER,SLVE,IREAD,IWRITE
                                 LOGICAL GRCELL,GRNODE
                        C
                                 COMMON/CFLOW/TURBUL,STEADY,UNSTDY,FTRAN,CONMAS,AXISYMM
                                 COMMON/SCHEME/UPWIND,HYBRID,POWER,SLVE(7),IREAD,IWRITE
                                 COMMON/CGRID/IN,JN,INM,JNM,IPREF,JPREF,CORP,NPERIOD,MFREQ
                              1        ,GRCELL,GRNODE
                                 COMMON/CONTR1/CC,MXIT,CCTM,MXSTEP,DELT,STIME,TTIME
                                 COMMON/CONTR2/RP(9),NSWEEP(7),NITER,RSDU(7),FDIF(7),RNORM(7)
                                 COMMON/CPROP/DENSIT,VISCOS,PR(7),PRT(7),RHO(IT,JT)
                              1         ,SPHEAT,SPH(IT,JT),VIS(IT,JT),VIST(IT,JT),GAMMA
                                 COMMON/CTURB/CD1,CD2,CMU,ELOG,CAPPA
                                 COMMON/CTAG/NTAG(IT,JT),NTAGW(IT,JT),NTAGE(IT,JT)
                              1             ,NTAGS(IT,JT),NTAGN(IT,JT)
                                 COMMON/COFV/AW(IT,JT),AE(IT,JT),AS(IT,JT),AN(IT,JT)
                                 COMMON/CSOR/SU(IT,JT),SP(IT,JT),BSOR(9)
                                 COMMON/CAP/AP1(IT,JT),AP(IT,JT),APU(IT,JT),APV(IT,JT)
                                 COMMON/CVAR/U(IT,JT),V(IT,JT),W(IT,JT),P(IT,JT)
                              1             ,E(IT,JT),D(IT,JT),T(IT,JT),PP(IT,JT),PSM(IT,JT)
                                COMMON/CVAO/UO(IT,JT),VO(IT,JT),WO(IT,JT),PO(IT,JT)
                              1            ,EO(IT,JT),DO(IT,JT),TO(IT,JT),RHOO(IT,JT)
                                 COMMON/CORD/X(IT),Y(JT),XC(IT),YC(JT),R(JT),RC(JT)
                              1           ,DXMI(IT),DYMI(JT),DXP(IT),DYP(JT),VOL(IT,JT)
                                 COMMON/CHEAT/QW(IT,JT),TAUW(IT,JT),AMW(IT,JT)
   329   330   331   332   333   334   335   336   337   338   339