TALK=T;RUN(1,1) DISPLAY The case considered is fully-developed turbulent flow and heat transfer in a circular pipe at Re=1.E5 and Pr=3.0. The tube wall is held at a constant temperature, and the calculation integrates down to the wall through the viscous sublayer. A non-uniform grid is employed so as to concentrate cells very close to the wall. For this purpose a grid is generated which is a geometric progression with the property that the ratio of any two adjacent cell lengths is a constant. The turbulence is simulated using the low Reynolds version of the Spalart-Allmaras model. This model predicts the following results compared to experimental data: f vdp Nu SA 0.018 4.07 361 Data 0.018 3.75 392 where f =8.*(wstar/wbulk)**2 is the friction factor, vdp=(wcl-wbulk)/u* is the velocity defect parameter and Nu=h*d/k is the Nusselt number. In-Form is used to print the f, vdp and Nu to the PHIDA, RESULT and inforout files. ENDDIS ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(S-A Low Re_1D DEVELOPED PIPE FLOW ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = F NX = 1 XULAST =0.1 XFRAC(1)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 58 YVLAST =0.05 YFRAC(1)=0.038842 ;YFRAC(3)=0.206387 YFRAC(5)=0.344854 ;YFRAC(7)=0.459289 YFRAC(9)=0.553864 ;YFRAC(11)=0.632025 YFRAC(13)=0.696621 ;YFRAC(15)=0.750006 YFRAC(17)=0.794126 ;YFRAC(19)=0.830588 YFRAC(21)=0.860723 ;YFRAC(23)=0.885627 YFRAC(25)=0.906209 ;YFRAC(27)=0.92322 YFRAC(29)=0.937277 ;YFRAC(31)=0.948896 YFRAC(33)=0.958497 ;YFRAC(35)=0.966433 YFRAC(37)=0.972991 ;YFRAC(39)=0.978411 YFRAC(41)=0.98289 ;YFRAC(43)=0.986592 YFRAC(45)=0.989651 ;YFRAC(47)=0.99218 YFRAC(49)=0.994269 ;YFRAC(51)=0.995996 YFRAC(53)=0.997424 ;YFRAC(55)=0.998603 YFRAC(57)=0.999578 ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 1 ZWLAST =1.0E-02 ZFRAC(1)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(5)=V1 ;NAME(7)=W1 NAME(124)=NUSS ;NAME(125)=TB NAME(126)=ASUM ;NAME(127)=TSUM NAME(128)=YPL ;NAME(129)=UPL NAME(130)=VDP ;NAME(131)=USTR NAME(132)=UTAU ;NAME(133)=FRIC NAME(134)=PRPS ;NAME(135)=HTCO NAME(136)=STAN ;NAME(137)=AREH NAME(138)=KOND ;NAME(139)=EPKE NAME(140)=DWDY ;NAME(141)=LTLS NAME(142)=WDIS ;NAME(143)=VOR1 NAME(144)=SKIN ;NAME(145)=YPLS NAME(146)=STRS ;NAME(147)=LEN1 NAME(148)=ENUT ;NAME(149)=TEM1 NAME(150)=ENTI * Y in SOLUTN argument list denotes: * 1-stored 2-solved 3-whole-field * 4-point-by-point 5-explicit 6-harmonic averaging SOLUTN(V1,Y,N,N,N,N,Y) SOLUTN(W1,Y,Y,N,N,N,N) SOLUTN(NUSS,Y,N,N,N,N,N) SOLUTN(TB,Y,N,N,N,N,N) SOLUTN(ASUM,Y,N,N,N,N,N) SOLUTN(TSUM,Y,N,N,N,N,N) SOLUTN(YPL,Y,N,N,N,N,N) SOLUTN(UPL,Y,N,N,N,N,N) SOLUTN(VDP,Y,N,N,N,N,N) SOLUTN(USTR,Y,N,N,N,N,N) SOLUTN(UTAU,Y,N,N,N,N,N) SOLUTN(FRIC,Y,N,N,N,N,N) SOLUTN(PRPS,Y,N,N,N,N,Y) SOLUTN(HTCO,Y,N,N,N,N,Y) SOLUTN(STAN,Y,N,N,N,N,Y) SOLUTN(AREH,Y,N,N,N,N,Y) SOLUTN(KOND,Y,N,N,N,N,Y) SOLUTN(EPKE,Y,N,N,N,N,N) SOLUTN(DWDY,Y,N,N,N,N,N) SOLUTN(LTLS,Y,Y,Y,N,N,Y) SOLUTN(WDIS,Y,N,N,N,N,N) SOLUTN(VOR1,Y,N,N,N,N,N) SOLUTN(SKIN,Y,N,N,N,N,Y) SOLUTN(YPLS,Y,N,N,N,N,Y) SOLUTN(STRS,Y,N,N,N,N,Y) SOLUTN(LEN1,Y,N,N,N,N,Y) SOLUTN(ENUT,Y,N,N,N,N,Y) SOLUTN(TEM1,Y,Y,N,N,N,Y) SOLUTN(ENTI,Y,Y,Y,N,N,Y) VIST = 148 LEN1 = 147 PRPS = 134 ************************************************************ Group 8. Terms & Devices * Y in TERMS argument list denotes: * 1-built-in source 2-convection 3-diffusion 4-transient * 5-first phase variable 6-interphase transport TERMS(W1,N,N,Y,Y,Y,Y) TERMS(LTLS,N,N,Y,N,Y,Y) TERMS(TEM1,N,N,Y,Y,N,Y) TERMS(ENTI,N,Y,Y,Y,Y,Y) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWENT = F USOURC = T ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1 ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 =1. ;TMP1 =0. ;EL1 =0. TSURR =0. ;TEMP0 =0. ;PRESS0 =0. DVO1DT =0. ;DRH1DP =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. ENUL =1.0E-06 ;ENUT = GRND2 ENUTA =0. ;ENUTB =0. ;ENUTC =0. IENUTA = 22 PRNDTL(W1)=1. ;PRNDTL(LTLS)=1. PRNDTL(TEM1)=CONDFILE ;PRNDTL(ENTI)=0.666667 PRT(W1)=1. ;PRT(LTLS)=1. PRT(TEM1)=1. ;PRT(ENTI)=0.666667 CP1 =1. ;CP2 =1. * List of user-defined materials to be read by EARTH MATFLG=T;IMAT=1 * Name *Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr. *33 1. 1.E-6 1.0 3.333E-7 0 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10 FIINIT(NUSS)=1.0E-10 ;FIINIT(TB)=1.0E-10 FIINIT(ASUM)=1.0E-10 ;FIINIT(TSUM)=1.0E-10 FIINIT(YPL)=1.0E-10 ;FIINIT(UPL)=1.0E-10 FIINIT(VDP)=1.0E-10 ;FIINIT(USTR)=1.0E-10 FIINIT(UTAU)=1.0E-10 ;FIINIT(FRIC)=1.0E-10 FIINIT(PRPS)=33. ;FIINIT(HTCO)=1.0E-10 FIINIT(STAN)=1.0E-10 ;FIINIT(AREH)=1.0E-10 FIINIT(KOND)=1.0E-10 ;FIINIT(EPKE)=1.0E-10 FIINIT(DWDY)=1.0E-10 ;FIINIT(LTLS)=1.0E-10 FIINIT(WDIS)=1.0E-03 ;FIINIT(VOR1)=1.0E-10 FIINIT(SKIN)=1.0E-10 ;FIINIT(YPLS)=1.0E-10 FIINIT(STRS)=1.0E-10 ;FIINIT(LEN1)=1.0E-10 FIINIT(ENUT)=1.0E-10 ;FIINIT(TEM1)=5. FIINIT(ENTI)=1.4025E-04 PATCH(IN1 ,INIVAL, 1, 1, 1, 1, 1, 1, 1, 1) INIT(IN1 ,W1 ,0. ,1.144765 ) PATCH(IN2 ,INIVAL, 1, 1, 2, 2, 1, 1, 1, 1) INIT(IN2 ,W1 ,0. ,1.137051 ) PATCH(IN3 ,INIVAL, 1, 1, 3, 3, 1, 1, 1, 1) INIT(IN3 ,W1 ,0. ,1.12598 ) PATCH(IN4 ,INIVAL, 1, 1, 4, 4, 1, 1, 1, 1) INIT(IN4 ,W1 ,0. ,1.114904 ) PATCH(IN5 ,INIVAL, 1, 1, 5, 5, 1, 1, 1, 1) INIT(IN5 ,W1 ,0. ,1.103823 ) PATCH(IN6 ,INIVAL, 1, 1, 6, 6, 1, 1, 1, 1) INIT(IN6 ,W1 ,0. ,1.092735 ) PATCH(IN7 ,INIVAL, 1, 1, 7, 7, 1, 1, 1, 1) INIT(IN7 ,W1 ,0. ,1.081639 ) PATCH(IN8 ,INIVAL, 1, 1, 8, 8, 1, 1, 1, 1) INIT(IN8 ,W1 ,0. ,1.070536 ) PATCH(IN9 ,INIVAL, 1, 1, 9, 9, 1, 1, 1, 1) INIT(IN9 ,W1 ,0. ,1.059425 ) PATCH(IN10 ,INIVAL, 1, 1, 10, 10, 1, 1, 1, 1) INIT(IN10 ,W1 ,0. ,1.048303 ) PATCH(IN11 ,INIVAL, 1, 1, 11, 11, 1, 1, 1, 1) INIT(IN11 ,W1 ,0. ,1.037171 ) PATCH(IN12 ,INIVAL, 1, 1, 12, 12, 1, 1, 1, 1) INIT(IN12 ,W1 ,0. ,1.026028 ) PATCH(IN13 ,INIVAL, 1, 1, 13, 13, 1, 1, 1, 1) INIT(IN13 ,W1 ,0. ,1.014872 ) PATCH(IN14 ,INIVAL, 1, 1, 14, 14, 1, 1, 1, 1) INIT(IN14 ,W1 ,0. ,1.003701 ) PATCH(IN15 ,INIVAL, 1, 1, 15, 15, 1, 1, 1, 1) INIT(IN15 ,W1 ,0. ,0.992515 ) PATCH(IN16 ,INIVAL, 1, 1, 16, 16, 1, 1, 1, 1) INIT(IN16 ,W1 ,0. ,0.981312 ) PATCH(IN17 ,INIVAL, 1, 1, 17, 17, 1, 1, 1, 1) INIT(IN17 ,W1 ,0. ,0.97009 ) PATCH(IN18 ,INIVAL, 1, 1, 18, 18, 1, 1, 1, 1) INIT(IN18 ,W1 ,0. ,0.958847 ) PATCH(IN19 ,INIVAL, 1, 1, 19, 19, 1, 1, 1, 1) INIT(IN19 ,W1 ,0. ,0.947581 ) PATCH(IN20 ,INIVAL, 1, 1, 20, 20, 1, 1, 1, 1) INIT(IN20 ,W1 ,0. ,0.93629 ) PATCH(IN21 ,INIVAL, 1, 1, 21, 21, 1, 1, 1, 1) INIT(IN21 ,W1 ,0. ,0.92497 ) PATCH(IN22 ,INIVAL, 1, 1, 22, 22, 1, 1, 1, 1) INIT(IN22 ,W1 ,0. ,0.91362 ) PATCH(IN23 ,INIVAL, 1, 1, 23, 23, 1, 1, 1, 1) INIT(IN23 ,W1 ,0. ,0.902234 ) PATCH(IN24 ,INIVAL, 1, 1, 24, 24, 1, 1, 1, 1) INIT(IN24 ,W1 ,0. ,0.890811 ) PATCH(IN25 ,INIVAL, 1, 1, 25, 25, 1, 1, 1, 1) INIT(IN25 ,W1 ,0. ,0.879345 ) PATCH(IN26 ,INIVAL, 1, 1, 26, 26, 1, 1, 1, 1) INIT(IN26 ,W1 ,0. ,0.867833 ) PATCH(IN27 ,INIVAL, 1, 1, 27, 27, 1, 1, 1, 1) INIT(IN27 ,W1 ,0. ,0.856268 ) PATCH(IN28 ,INIVAL, 1, 1, 28, 28, 1, 1, 1, 1) INIT(IN28 ,W1 ,0. ,0.844646 ) PATCH(IN29 ,INIVAL, 1, 1, 29, 29, 1, 1, 1, 1) INIT(IN29 ,W1 ,0. ,0.832959 ) PATCH(IN30 ,INIVAL, 1, 1, 30, 30, 1, 1, 1, 1) INIT(IN30 ,W1 ,0. ,0.821202 ) PATCH(IN31 ,INIVAL, 1, 1, 31, 31, 1, 1, 1, 1) INIT(IN31 ,W1 ,0. ,0.809364 ) PATCH(IN32 ,INIVAL, 1, 1, 32, 32, 1, 1, 1, 1) INIT(IN32 ,W1 ,0. ,0.797438 ) PATCH(IN33 ,INIVAL, 1, 1, 33, 33, 1, 1, 1, 1) INIT(IN33 ,W1 ,0. ,0.785412 ) PATCH(IN34 ,INIVAL, 1, 1, 34, 34, 1, 1, 1, 1) INIT(IN34 ,W1 ,0. ,0.773276 ) PATCH(IN35 ,INIVAL, 1, 1, 35, 35, 1, 1, 1, 1) INIT(IN35 ,W1 ,0. ,0.761015 ) PATCH(IN36 ,INIVAL, 1, 1, 36, 36, 1, 1, 1, 1) INIT(IN36 ,W1 ,0. ,0.748613 ) PATCH(IN37 ,INIVAL, 1, 1, 37, 37, 1, 1, 1, 1) INIT(IN37 ,W1 ,0. ,0.736054 ) PATCH(IN38 ,INIVAL, 1, 1, 38, 38, 1, 1, 1, 1) INIT(IN38 ,W1 ,0. ,0.723317 ) PATCH(IN39 ,INIVAL, 1, 1, 39, 39, 1, 1, 1, 1) INIT(IN39 ,W1 ,0. ,0.710377 ) PATCH(IN40 ,INIVAL, 1, 1, 40, 40, 1, 1, 1, 1) INIT(IN40 ,W1 ,0. ,0.697208 ) PATCH(IN41 ,INIVAL, 1, 1, 41, 41, 1, 1, 1, 1) INIT(IN41 ,W1 ,0. ,0.683776 ) PATCH(IN42 ,INIVAL, 1, 1, 42, 42, 1, 1, 1, 1) INIT(IN42 ,W1 ,0. ,0.670044 ) PATCH(IN43 ,INIVAL, 1, 1, 43, 43, 1, 1, 1, 1) INIT(IN43 ,W1 ,0. ,0.655964 ) PATCH(IN44 ,INIVAL, 1, 1, 44, 44, 1, 1, 1, 1) INIT(IN44 ,W1 ,0. ,0.641482 ) PATCH(IN45 ,INIVAL, 1, 1, 45, 45, 1, 1, 1, 1) INIT(IN45 ,W1 ,0. ,0.626529 ) PATCH(IN46 ,INIVAL, 1, 1, 46, 46, 1, 1, 1, 1) INIT(IN46 ,W1 ,0. ,0.611022 ) PATCH(IN47 ,INIVAL, 1, 1, 47, 47, 1, 1, 1, 1) INIT(IN47 ,W1 ,0. ,0.594856 ) PATCH(IN48 ,INIVAL, 1, 1, 48, 48, 1, 1, 1, 1) INIT(IN48 ,W1 ,0. ,0.577897 ) PATCH(IN49 ,INIVAL, 1, 1, 49, 49, 1, 1, 1, 1) INIT(IN49 ,W1 ,0. ,0.55997 ) PATCH(IN50 ,INIVAL, 1, 1, 50, 50, 1, 1, 1, 1) INIT(IN50 ,W1 ,0. ,0.540842 ) PATCH(IN51 ,INIVAL, 1, 1, 51, 51, 1, 1, 1, 1) INIT(IN51 ,W1 ,0. ,0.495807 ) PATCH(IN52 ,INIVAL, 1, 1, 52, 52, 1, 1, 1, 1) INIT(IN52 ,W1 ,0. ,0.407648 ) PATCH(IN53 ,INIVAL, 1, 1, 53, 53, 1, 1, 1, 1) INIT(IN53 ,W1 ,0. ,0.327501 ) PATCH(IN54 ,INIVAL, 1, 1, 54, 54, 1, 1, 1, 1) INIT(IN54 ,W1 ,0. ,0.254644 ) PATCH(IN55 ,INIVAL, 1, 1, 55, 55, 1, 1, 1, 1) INIT(IN55 ,W1 ,0. ,0.188408 ) PATCH(IN56 ,INIVAL, 1, 1, 56, 56, 1, 1, 1, 1) INIT(IN56 ,W1 ,0. ,0.128203 ) PATCH(IN57 ,INIVAL, 1, 1, 57, 57, 1, 1, 1, 1) INIT(IN57 ,W1 ,0. ,0.073459 ) PATCH(IN58 ,INIVAL, 1, 1, 58, 58, 1, 1, 1, 1) INIT(IN58 ,W1 ,0. ,0.023697 ) INIADD = F FSWEEP = 1 NAMFI =CHAM ************************************************************ Group 12. Patchwise adjustment of terms Patches for this group are printed with those for Group 13. Their names begin either with GP12 or & ************************************************************ Group 13. Boundary & Special Sources PATCH(SASOPRO ,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(SASOPRO ,ENTI, FIXFLU , GRND4 ) PATCH(SASODIS ,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(SASODIS ,ENTI, GRND4 ,0. ) PATCH(SASODIF ,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(SASODIF ,ENTI, FIXFLU , GRND4 ) PATCH(WALLN ,NWALL , 1, 1, 58, 58, 1, 1, 1, 1) COVAL(WALLN ,W1 , GRND2 ,0. ) COVAL(WALLN ,TEM1, GRND2 ,10. ) PATCH(FDFW1DP ,VOLUME, 1, 1, 1, 58, 1, 1, 1, 1) COVAL(FDFW1DP ,W1 ,1.25E-04 , GRND1 ) PATCH(FDFCWT ,PHASEM, 1, 1, 1, 58, 1, 1, 1, 1) COVAL(FDFCWT ,TEM1,0.23795 ,10. ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 70 ;ISWC1 = 1 LITHYD = 6 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations LITER(W1)=10 ;LITER(LTLS)=20 LITER(TEM1)=20 ;LITER(ENTI)=20 ENDIT(W1)=1.0E-03 ;ENDIT(LTLS)=1.0E-03 ENDIT(TEM1)=1.0E-03 ;ENDIT(ENTI)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(V1,LINRLX,1.) RELAX(W1,FALSDT,0.05) RELAX(NUSS,LINRLX,1.) RELAX(TB,LINRLX,1.) RELAX(ASUM,LINRLX,1.) RELAX(TSUM,LINRLX,1.) RELAX(YPL,LINRLX,1.) RELAX(UPL,LINRLX,1.) RELAX(VDP,LINRLX,1.) RELAX(USTR,LINRLX,1.) RELAX(UTAU,LINRLX,1.) RELAX(FRIC,LINRLX,1.) RELAX(PRPS,LINRLX,1.) RELAX(HTCO,LINRLX,1.) RELAX(STAN,LINRLX,1.) RELAX(AREH,LINRLX,1.) RELAX(KOND,LINRLX,1.) RELAX(EPKE,LINRLX,1.) RELAX(DWDY,LINRLX,1.) RELAX(LTLS,FALSDT,1.0E+09) RELAX(WDIS,LINRLX,1.) RELAX(VOR1,LINRLX,1.) RELAX(SKIN,LINRLX,1.) RELAX(YPLS,LINRLX,1.) RELAX(STRS,LINRLX,1.) RELAX(LEN1,LINRLX,1.) RELAX(ENUT,LINRLX,1.) RELAX(TEM1,FALSDT,1.0E+09) RELAX(ENTI,LINRLX,0.3) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=1.0E-10 VARMAX(NUSS)=1.0E+10 ;VARMIN(NUSS)=-1.0E+10 VARMAX(TB)=1.0E+10 ;VARMIN(TB)=-1.0E+10 VARMAX(ASUM)=1.0E+10 ;VARMIN(ASUM)=-1.0E+10 VARMAX(TSUM)=1.0E+10 ;VARMIN(TSUM)=-1.0E+10 VARMAX(YPL)=1.0E+10 ;VARMIN(YPL)=-1.0E+10 VARMAX(UPL)=1.0E+10 ;VARMIN(UPL)=-1.0E+10 VARMAX(VDP)=1.0E+10 ;VARMIN(VDP)=-1.0E+10 VARMAX(USTR)=1.0E+10 ;VARMIN(USTR)=-1.0E+10 VARMAX(UTAU)=1.0E+10 ;VARMIN(UTAU)=-1.0E+10 VARMAX(FRIC)=1.0E+10 ;VARMIN(FRIC)=-1.0E+10 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(HTCO)=1.0E+10 ;VARMIN(HTCO)=-1.0E+10 VARMAX(STAN)=1.0E+10 ;VARMIN(STAN)=-1.0E+10 VARMAX(AREH)=1.0E+10 ;VARMIN(AREH)=-1.0E+10 VARMAX(KOND)=1.0E+10 ;VARMIN(KOND)=-1.0E+10 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10 VARMAX(DWDY)=1.0E+10 ;VARMIN(DWDY)=-1.0E+10 VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10 VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10 VARMAX(VOR1)=1.0E+10 ;VARMIN(VOR1)=-1.0E+10 VARMAX(SKIN)=1.0E+10 ;VARMIN(SKIN)=-1.0E+10 VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10 VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-1.0E+10 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10 VARMAX(TEM1)=1.0E+10 ;VARMIN(TEM1)=-1.0E+10 VARMAX(ENTI)=1.0E+10 ;VARMIN(ENTI)=1.0E-10 ************************************************************ Group 19. Data transmitted to GROUND GENK = T FDFSOL = T PARSOL = F IENUTA = 22 ISG62 = 1 CSG10 ='q1' SPEDAT(SET,STORED,FRIC,C,=8.*STRS/(1.*1.)!ZSLFIN) SPEDAT(SET,STORED,UTAU,C,=STRS^0.5!ZSLFIN) SPEDAT(SET,STORED,USTR,C,=UTAU[&NY&]!ZSLFIN) SPEDAT(SET,STORED,VDP,C,=(W1[&1&]-1.)/UTAU[&NY&]!ZSLFIN) SPEDAT(SET,STORED,UPL,C,=W1/USTR!ZSLFIN) SPEDAT(SET,STORED,YPL,C,=((YVLAST-YG)*USTR)/1.0E-06!ZSLFIN) SPEDAT(SET,STORED,AREH,C,=AHIGH!ZSLFIN) SPEDAT(SET,STORED,TSUM,C,=0.0!ZSLFIN) SPEDAT(SET,STORED,TSUM!FDFCWT,C,=SUM(W1*AREH*TEM1)!ZSLFIN) SPEDAT(SET,STORED,ASUM,C,=0.0!ZSLFIN) SPEDAT(SET,STORED,ASUM!FDFCWT,C,=SUM(W1*AREH)!ZSLFIN) SPEDAT(SET,STORED,TB!FDFCWT,C,=TSUM/ASUM!ZSLFIN) SPEDAT(SET,STORED,NUSS,C,=HTCO*0.1*(10.-TEM1[&NY&1])/(3.333333E-0$) SPEDAT(SET,STORED,NUSS,C,7*(10.-TB))!ZSLFIN) SPEDAT(SET,MAKE,FFAC,C,=0.0) SPEDAT(SET,MAKE,USTAR,C,=0.0) SPEDAT(SET,MAKE,VDUS,C,=0.0) SPEDAT(SET,MAKE,NUS,C,=0.0) SPEDAT(SET,STORED,FFAC!WALLN,C,=FRIC!ZSLFIN) SPEDAT(SET,STORED,USTAR!WALLN,C,=UTAU!ZSLFIN) SPEDAT(SET,STORED,VDUS!WALLN,C,=VDP!ZSLFIN) SPEDAT(SET,STORED,NUS!WALLN,C,=NUSS!ZSLFIN) SPEDAT(SET,PRINT,F,C,=FFAC) SPEDAT(SET,PRINT,U*,C,=USTAR) SPEDAT(SET,PRINT,(UCL-UB)/U*,C,=VDUS) SPEDAT(SET,PRINT,NU,C,=NUS) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,GXMONI,CLASSIC,L,T) SPEDAT(SET,MATERIAL,33,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(V1)=1.0E-10 ;EX(W1)=0.7849 EX(NUSS)=6.222 ;EX(TB)=5.055 EX(ASUM)=1.25E-04 ;EX(TSUM)=6.318E-04 EX(YPL)=439.5 ;EX(UPL)=16.68 EX(VDP)=4.074 ;EX(USTR)=0.04706 EX(UTAU)=8.212E-04 ;EX(FRIC)=3.055E-04 EX(PRPS)=33. ;EX(HTCO)=5.447E-04 EX(STAN)=0.02331 ;EX(AREH)=2.155E-06 EX(KOND)=3.333E-07 ;EX(EPKE)=1.0E-10 EX(DWDY)=354. ;EX(LTLS)=1.714E-04 EX(WDIS)=7.667E-03 ;EX(VOR1)=354. EX(SKIN)=0.06994 ;EX(YPLS)=8.56E-03 EX(STRS)=3.819E-05 ;EX(LEN1)=0. EX(ENUT)=7.414E-05 ;EX(TEM1)=5.893 EX(ENTI)=7.494E-05 ************************************************************ Group 21. Print-out of Variables INIFLD = F ;SUBWGR = F * Y in OUTPUT argument list denotes: * 1-field 2-correction-eq. monitor 3-selective dumping * 4-whole-field residual 5-spot-value table 6-residual table OUTPUT(V1,Y,N,Y,N,N,N) OUTPUT(W1,Y,N,Y,Y,Y,Y) OUTPUT(NUSS,Y,N,Y,N,N,N) OUTPUT(TB,Y,N,Y,N,N,N) OUTPUT(ASUM,Y,N,Y,N,N,N) OUTPUT(TSUM,Y,N,Y,N,N,N) OUTPUT(YPL,Y,N,Y,N,N,N) OUTPUT(UPL,Y,N,Y,N,N,N) OUTPUT(VDP,Y,N,Y,N,N,N) OUTPUT(USTR,Y,N,Y,N,N,N) OUTPUT(UTAU,Y,N,Y,N,N,N) OUTPUT(FRIC,Y,N,Y,N,N,N) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(HTCO,Y,N,Y,N,N,N) OUTPUT(STAN,Y,N,Y,N,N,N) OUTPUT(AREH,Y,N,Y,N,N,N) OUTPUT(KOND,Y,N,Y,N,N,N) OUTPUT(EPKE,Y,N,Y,N,N,N) OUTPUT(DWDY,Y,N,Y,N,N,N) OUTPUT(LTLS,Y,N,Y,Y,Y,Y) OUTPUT(WDIS,Y,N,Y,N,N,N) OUTPUT(VOR1,Y,N,Y,N,N,N) OUTPUT(SKIN,Y,N,Y,N,N,N) OUTPUT(YPLS,Y,N,Y,N,N,N) OUTPUT(STRS,Y,N,Y,N,N,N) OUTPUT(LEN1,Y,N,Y,N,N,N) OUTPUT(ENUT,Y,N,Y,N,N,N) OUTPUT(TEM1,Y,N,Y,Y,Y,Y) OUTPUT(ENTI,Y,N,Y,Y,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 56 ;IZMON = 1 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ;NUMCLS = 5 NYPRIN = 2 ;IYPRF = 1 ;IYPRL = 10000 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.4 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20 No PATCHes yet used for this Group ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM STOP