TALK=T;RUN(1,1)
  DISPLAY
  The problem considered is plane turbulent couette flow in a
  channel with one moving wall, as simulated for library case T100
  using various high-Re turbulence models with wall functions. Here,
  the turbulence is simulated by use of various low-Re k-e and
  k-w models. The calculation integrates down to the wall and the
  solution is performed by use of the single-slab solver. A non-
  uniform grid is employed so as to concentrate cells very close to
  the walls. The Reynolds number is 1.E5 based on the channel height
  and the average velocity.

  For this flow, the shear stress is uniform across the flow and
  hence the axial pressure-gradient is zero. The velocity profile
  is S-shaped and symmetrical about the central plane, and so the
  average velocity is one half of the velocity of the moving wall.
  ENDDIS
   AUTOPLOT USE
   file
   phi 5

   da 1 w1;col9 1
   msg Velocity (W1) profile
   msg Press RETURN to continue
   pause
   clear
   da 1 ke;col9 1
   msg KE profile
   msg Press RETURN to continue
   pause
   clear
   da 1 enut;col9 1
   msg ENUT profile
   msg Press e to END
   ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(LAM-BRE KE_1D PLANE COUETTE FLOW  :T205 )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 205
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 1
 XULAST =1.
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 132
 YVLAST =0.1
 YFRAC(1)=2.55102E-04 ;YFRAC(2)=5.306122E-04
 YFRAC(3)=8.281632E-04 ;YFRAC(4)=1.149518E-03
 YFRAC(5)=1.496582E-03 ;YFRAC(6)=1.87141E-03
 YFRAC(7)=2.276225E-03 ;YFRAC(8)=2.713425E-03
 YFRAC(9)=3.185601E-03 ;YFRAC(10)=3.695552E-03
 YFRAC(11)=4.246298E-03 ;YFRAC(12)=4.841104E-03
 YFRAC(13)=5.483495E-03 ;YFRAC(14)=6.177276E-03
 YFRAC(15)=6.926561E-03 ;YFRAC(16)=7.735788E-03
 YFRAC(17)=8.609753E-03 ;YFRAC(18)=9.553635E-03
 YFRAC(19)=0.010573 ;YFRAC(20)=0.011674
 YFRAC(21)=0.012863 ;YFRAC(22)=0.014147
 YFRAC(23)=0.015534 ;YFRAC(24)=0.017032
 YFRAC(25)=0.018649 ;YFRAC(26)=0.020397
 YFRAC(27)=0.022283 ;YFRAC(28)=0.024321
 YFRAC(29)=0.026522 ;YFRAC(30)=0.028899
 YFRAC(31)=0.031466 ;YFRAC(32)=0.034238
 YFRAC(33)=0.037232 ;YFRAC(34)=0.040466
 YFRAC(35)=0.043958 ;YFRAC(36)=0.04773
 YFRAC(37)=0.051804 ;YFRAC(38)=0.056203
 YFRAC(39)=0.060954 ;YFRAC(40)=0.066086
 YFRAC(41)=0.071628 ;YFRAC(42)=0.077613
 YFRAC(43)=0.084077 ;YFRAC(44)=0.091059
 YFRAC(45)=0.098599 ;YFRAC(46)=0.106741
 YFRAC(47)=0.115536 ;YFRAC(48)=0.125034
 YFRAC(49)=0.135292 ;YFRAC(50)=0.14637
 YFRAC(51)=0.158335 ;YFRAC(52)=0.171257
 YFRAC(53)=0.185212 ;YFRAC(54)=0.200285
 YFRAC(55)=0.216562 ;YFRAC(56)=0.234143
 YFRAC(57)=0.253129 ;YFRAC(58)=0.273634
 YFRAC(59)=0.29578 ;YFRAC(60)=0.319698
 YFRAC(61)=0.345529 ;YFRAC(62)=0.373426
 YFRAC(63)=0.403556 ;YFRAC(64)=0.436095
 YFRAC(65)=0.471238 ;YFRAC(66)=0.5
 YFRAC(67)=0.528762 ;YFRAC(68)=0.563905
 YFRAC(69)=0.596444 ;YFRAC(70)=0.626574
 YFRAC(71)=0.654471 ;YFRAC(72)=0.680302
 YFRAC(73)=0.70422 ;YFRAC(74)=0.726366
 YFRAC(75)=0.746871 ;YFRAC(76)=0.765857
 YFRAC(77)=0.783438 ;YFRAC(78)=0.799715
 YFRAC(79)=0.814788 ;YFRAC(80)=0.828743
 YFRAC(81)=0.841665 ;YFRAC(82)=0.85363
 YFRAC(83)=0.864708 ;YFRAC(84)=0.874966
 YFRAC(85)=0.884464 ;YFRAC(86)=0.893259
 YFRAC(87)=0.901402 ;YFRAC(88)=0.908941
 YFRAC(89)=0.915923 ;YFRAC(90)=0.922387
 YFRAC(91)=0.928372 ;YFRAC(92)=0.933914
 YFRAC(93)=0.939046 ;YFRAC(94)=0.943797
 YFRAC(95)=0.948196 ;YFRAC(96)=0.95227
 YFRAC(97)=0.956042 ;YFRAC(98)=0.959534
 YFRAC(99)=0.962768 ;YFRAC(100)=0.965762
 YFRAC(101)=0.968534 ;YFRAC(102)=0.971101
 YFRAC(103)=0.973478 ;YFRAC(104)=0.975679
 YFRAC(105)=0.977717 ;YFRAC(106)=0.979603
 YFRAC(107)=0.981351 ;YFRAC(108)=0.982968
 YFRAC(109)=0.984466 ;YFRAC(110)=0.985853
 YFRAC(111)=0.987137 ;YFRAC(112)=0.988326
 YFRAC(113)=0.989427 ;YFRAC(114)=0.990446
 YFRAC(115)=0.99139 ;YFRAC(116)=0.992264
 YFRAC(117)=0.993073 ;YFRAC(118)=0.993823
 YFRAC(119)=0.994516 ;YFRAC(120)=0.995159
 YFRAC(121)=0.995754 ;YFRAC(122)=0.996304
 YFRAC(123)=0.996814 ;YFRAC(124)=0.997287
 YFRAC(125)=0.997724 ;YFRAC(126)=0.998129
 YFRAC(127)=0.998503 ;YFRAC(128)=0.99885
 YFRAC(129)=0.999172 ;YFRAC(130)=0.999469
 YFRAC(131)=0.999745 ;YFRAC(132)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =1.
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(7)=W1 ;NAME(12)=KE
 NAME(13)=EP ;NAME(140)=FTWO
 NAME(141)=FONE ;NAME(142)=REYN
 NAME(143)=REYT ;NAME(144)=FMU
 NAME(145)=LTLS ;NAME(146)=WDIS
 NAME(147)=CF ;NAME(148)=STRS
 NAME(149)=LEN1 ;NAME(150)=ENUT
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(W1,Y,Y,N,N,N,N)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(FTWO,Y,N,N,N,N,Y)
 SOLUTN(FONE,Y,N,N,N,N,Y)
 SOLUTN(REYN,Y,N,N,N,N,Y)
 SOLUTN(REYT,Y,N,N,N,N,Y)
 SOLUTN(FMU,Y,N,N,N,N,Y)
 SOLUTN(LTLS,Y,Y,Y,N,N,Y)
 SOLUTN(WDIS,Y,N,N,N,N,N)
 SOLUTN(CF,Y,N,N,N,N,N)
 SOLUTN(STRS,Y,N,N,N,N,Y)
 SOLUTN(LEN1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 VIST = 150
 LEN1 = 149
 ************************************************************
  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(KE,Y,N,Y,Y,Y,N)
 TERMS(EP,Y,N,Y,Y,Y,N)
 TERMS(LTLS,N,N,Y,N,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWENL = T
 NEWENT = 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 = GRND1 ;EL1 = GRND4
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 TMP1A =0. ;TMP1B =0. ;TMP1C =0.
 TMP2A =0. ;TMP2B =0. ;TMP2C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0. ;EL1C =0.
 ENUL = GRND1 ;ENUT = GRND3
 ENULA =5.0E-07 ;ENULB =0. ;ENULC =0. ;ENULD =0.
 ENULE =0. ;ENULF =0. ;ENULG =0.
 IENULA = 0 ;IENULB = 0
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 3
 PRNDTL(W1)=1. ;PRNDTL(KE)=1.
 PRNDTL(EP)=1. ;PRNDTL(LTLS)=1.
 PRT(W1)=1. ;PRT(KE)=1.
 PRT(EP)=1.314 ;PRT(LTLS)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(W1)=1.0E-10 ;FIINIT(KE)=1.280533E-03
 FIINIT(EP)=1.673059E-03 ;FIINIT(FTWO)=1.0E-10
 FIINIT(FONE)=1.0E-10 ;FIINIT(REYN)=1.0E-10
 FIINIT(REYT)=1.0E-10 ;FIINIT(FMU)=1.0E-10
 FIINIT(LTLS)=1.0E-10 ;FIINIT(WDIS)=1.0E-02
 FIINIT(CF)=1.0E-10 ;FIINIT(STRS)=1.0E-10
 FIINIT(LEN1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
 
 PATCH(ICOUF ,LINVLY, 1, 1, 1, 132, 1, 1, 1, 1)
 INIT(ICOUF ,W1 ,10. ,0. )
 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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(WALLN ,NWALL , 1, 1, 132, 132, 1, 1, 1, 1)
 COVAL(WALLN ,W1 , GRND2 ,1. )
 COVAL(WALLN ,KE ,1. ,0. )
 
 PATCH(WALLS ,SWALL , 1, 1, 1, 1, 1, 1, 1, 1)
 COVAL(WALLS ,W1 , GRND2 ,0. )
 COVAL(WALLS ,KE ,1. ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 80 ;ISWC1 = 1
 LITHYD = 6 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(W1)=10 ;LITER(KE)=20
 LITER(EP)=20 ;LITER(LTLS)=20
 ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
 ENDIT(EP)=1.0E-03 ;ENDIT(LTLS)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(W1,FALSDT,0.2)
 RELAX(KE,FALSDT,0.05)
 RELAX(EP,FALSDT,0.05)
 RELAX(FTWO,LINRLX,1.)
 RELAX(FONE,LINRLX,1.)
 RELAX(REYN,LINRLX,1.)
 RELAX(REYT,LINRLX,1.)
 RELAX(FMU,LINRLX,1.)
 RELAX(LTLS,FALSDT,1.0E+09)
 RELAX(WDIS,LINRLX,1.)
 RELAX(CF,LINRLX,1.)
 RELAX(STRS,LINRLX,1.)
 RELAX(LEN1,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 KELIN = 1
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=1.0E-10
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(FTWO)=1.0E+10 ;VARMIN(FTWO)=-1.0E+10
 VARMAX(FONE)=1.0E+10 ;VARMIN(FONE)=-1.0E+10
 VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10
 VARMAX(REYT)=1.0E+10 ;VARMIN(REYT)=-1.0E+10
 VARMAX(FMU)=1.0E+10 ;VARMIN(FMU)=-1.0E+10
 VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
 VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
 VARMAX(CF)=1.0E+10 ;VARMIN(CF)=-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
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 IENUTA = 3
 ISG62 = 1
 SPEDAT(SET,STORED,CF,C,=2.*STRS/(RHO1*0.5*0.5)!ZSLFIN)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(W1)=0.5 ;EX(KE)=1.589E-03
 EX(EP)=0.0169 ;EX(FTWO)=0.9653
 EX(FONE)=2.24 ;EX(REYN)=773.700012
 EX(REYT)=1833. ;EX(FMU)=0.7624
 EX(LTLS)=3.596E-04 ;EX(WDIS)=9.707E-03
 EX(CF)=5.737E-05 ;EX(STRS)=7.172E-06
 EX(LEN1)=3.769E-03 ;EX(ENUT)=8.077E-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(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(FTWO,Y,N,Y,N,N,N)
 OUTPUT(FONE,Y,N,Y,N,N,N)
 OUTPUT(REYN,Y,N,Y,N,N,N)
 OUTPUT(REYT,Y,N,Y,N,N,N)
 OUTPUT(FMU,Y,N,Y,N,N,N)
 OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
 OUTPUT(WDIS,Y,N,Y,N,N,N)
 OUTPUT(CF,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)
 WALPRN = T
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 2 ;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 = 1 ;IYPRF = 1 ;IYPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = 5
 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