TALK=T;RUN(1,1)
  DISPLAY
   The problem concerns 2d incompressible, turbulent flow and heat
   transfer over a backward-facing step. The case is similar to that
   described for library case T103, excepting that heat transfer is
   considered and the calculation is performed with any one of the
   following low-Re models: the Lam-Bremhorst k-e model, the Chen-
   Kim k-e model, the 2-layer k-e model, and the Wilcox k-omega
   model. The thermal field is calculated for a constant heat flux
   through the south wall downstream of the step and an adiabatic
   condition on all other walls. The laminar Prandtl number is taken
   as 1.0 and the turbulent Prandtl as 0.86.

  ENDDIS
   AUTOPLOT USE
   file
   phi 5

   da 1 u1 y 1
   divide x .0381 1
   shift x -4 1
   col9 1
   level y 0;level x 0
   msg Velocity (U1) profile
   msg Press e to END
   ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(LAM-BRE K-E BKWRD-FACING STEP Y-X  :T208)
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 208
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 80
 XULAST =0.762
 XFRAC(1)=0.027551 ;XFRAC(2)=0.053061
 XFRAC(3)=0.076531 ;XFRAC(4)=0.097959
 XFRAC(5)=0.117347 ;XFRAC(6)=0.134694
 XFRAC(7)=0.15 ;XFRAC(8)=0.163265
 XFRAC(9)=0.17449 ;XFRAC(10)=0.183673
 XFRAC(11)=0.190816 ;XFRAC(12)=0.195918
 XFRAC(13)=0.19898 ;XFRAC(14)=0.2
 XFRAC(15)=0.200124 ;XFRAC(16)=0.200498
 XFRAC(17)=0.20112 ;XFRAC(18)=0.201992
 XFRAC(19)=0.203112 ;XFRAC(20)=0.204482
 XFRAC(21)=0.2061 ;XFRAC(22)=0.207967
 XFRAC(23)=0.210084 ;XFRAC(24)=0.212449
 XFRAC(25)=0.215063 ;XFRAC(26)=0.217927
 XFRAC(27)=0.221039 ;XFRAC(28)=0.2244
 XFRAC(29)=0.22801 ;XFRAC(30)=0.23187
 XFRAC(31)=0.235978 ;XFRAC(32)=0.240335
 XFRAC(33)=0.244941 ;XFRAC(34)=0.249796
 XFRAC(35)=0.2549 ;XFRAC(36)=0.260254
 XFRAC(37)=0.265856 ;XFRAC(38)=0.271707
 XFRAC(39)=0.277807 ;XFRAC(40)=0.284156
 XFRAC(41)=0.290754 ;XFRAC(42)=0.297601
 XFRAC(43)=0.304697 ;XFRAC(44)=0.312042
 XFRAC(45)=0.319636 ;XFRAC(46)=0.327479
 XFRAC(47)=0.33557 ;XFRAC(48)=0.343911
 XFRAC(49)=0.352501 ;XFRAC(50)=0.36134
 XFRAC(51)=0.370428 ;XFRAC(52)=0.379765
 XFRAC(53)=0.38935 ;XFRAC(54)=0.399185
 XFRAC(55)=0.409269 ;XFRAC(56)=0.419602
 XFRAC(57)=0.430183 ;XFRAC(58)=0.441014
 XFRAC(59)=0.452094 ;XFRAC(60)=0.463422
 XFRAC(61)=0.475 ;XFRAC(62)=0.490334
 XFRAC(63)=0.510228 ;XFRAC(64)=0.532306
 XFRAC(65)=0.555933 ;XFRAC(66)=0.580784
 XFRAC(67)=0.606655 ;XFRAC(68)=0.633406
 XFRAC(69)=0.660936 ;XFRAC(70)=0.689164
 XFRAC(71)=0.718027 ;XFRAC(72)=0.747475
 XFRAC(73)=0.777463 ;XFRAC(74)=0.807956
 XFRAC(75)=0.838922 ;XFRAC(76)=0.870334
 XFRAC(77)=0.902168 ;XFRAC(78)=0.934403
 XFRAC(79)=0.967019 ;XFRAC(80)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 80
 YVLAST =0.1143
 YFRAC(1)=1.023523E-03 ;YFRAC(2)=3.325437E-03
 YFRAC(3)=6.625275E-03 ;YFRAC(4)=0.010804
 YFRAC(5)=0.015789 ;YFRAC(6)=0.021526
 YFRAC(7)=0.027975 ;YFRAC(8)=0.035104
 YFRAC(9)=0.042885 ;YFRAC(10)=0.051298
 YFRAC(11)=0.060321 ;YFRAC(12)=0.069937
 YFRAC(13)=0.080131 ;YFRAC(14)=0.09089
 YFRAC(15)=0.1022 ;YFRAC(16)=0.114052
 YFRAC(17)=0.126433 ;YFRAC(18)=0.139335
 YFRAC(19)=0.152749 ;YFRAC(20)=0.166667
 YFRAC(21)=0.180584 ;YFRAC(22)=0.193998
 YFRAC(23)=0.2069 ;YFRAC(24)=0.219282
 YFRAC(25)=0.231133 ;YFRAC(26)=0.242444
 YFRAC(27)=0.253202 ;YFRAC(28)=0.263397
 YFRAC(29)=0.273013 ;YFRAC(30)=0.282036
 YFRAC(31)=0.290448 ;YFRAC(32)=0.29823
 YFRAC(33)=0.305359 ;YFRAC(34)=0.311808
 YFRAC(35)=0.317545 ;YFRAC(36)=0.322529
 YFRAC(37)=0.326708 ;YFRAC(38)=0.330008
 YFRAC(39)=0.33231 ;YFRAC(40)=0.333333
 YFRAC(41)=0.334167 ;YFRAC(42)=0.336667
 YFRAC(43)=0.340833 ;YFRAC(44)=0.346667
 YFRAC(45)=0.354167 ;YFRAC(46)=0.363333
 YFRAC(47)=0.374167 ;YFRAC(48)=0.386667
 YFRAC(49)=0.400833 ;YFRAC(50)=0.416667
 YFRAC(51)=0.434167 ;YFRAC(52)=0.453333
 YFRAC(53)=0.474167 ;YFRAC(54)=0.496667
 YFRAC(55)=0.520833 ;YFRAC(56)=0.546667
 YFRAC(57)=0.574167 ;YFRAC(58)=0.603333
 YFRAC(59)=0.634167 ;YFRAC(60)=0.666667
 YFRAC(61)=0.699167 ;YFRAC(62)=0.73
 YFRAC(63)=0.759167 ;YFRAC(64)=0.786667
 YFRAC(65)=0.8125 ;YFRAC(66)=0.836667
 YFRAC(67)=0.859167 ;YFRAC(68)=0.88
 YFRAC(69)=0.899167 ;YFRAC(70)=0.916667
 YFRAC(71)=0.9325 ;YFRAC(72)=0.946667
 YFRAC(73)=0.959167 ;YFRAC(74)=0.97
 YFRAC(75)=0.979167 ;YFRAC(76)=0.986667
 YFRAC(77)=0.9925 ;YFRAC(78)=0.996667
 YFRAC(79)=0.999167 ;YFRAC(80)=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(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(12)=KE
 NAME(13)=EP ;NAME(14)=H1
 NAME(141)=EPKE ;NAME(143)=REYN
 NAME(144)=LTLS ;NAME(145)=WDIS
 NAME(146)=PRPS ;NAME(147)=FTWO
 NAME(148)=FMU ;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(P1,Y,Y,Y,N,N,N)
 SOLUTN(U1,Y,Y,N,N,N,N)
 SOLUTN(V1,Y,Y,N,N,N,N)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(H1,Y,Y,N,N,N,Y)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(REYN,Y,N,N,N,N,Y)
 SOLUTN(LTLS,Y,Y,Y,N,N,Y)
 SOLUTN(WDIS,Y,N,N,N,N,N)
 SOLUTN(PRPS,Y,N,N,N,N,Y)
 SOLUTN(FTWO,Y,N,N,N,N,Y)
 SOLUTN(FMU,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
 PRPS = 146
 ************************************************************
  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(P1,Y,Y,Y,N,Y,Y)
 TERMS(U1,Y,Y,Y,Y,Y,Y)
 TERMS(V1,Y,Y,Y,Y,Y,Y)
 TERMS(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 TERMS(H1,N,Y,Y,Y,Y,Y)
 TERMS(LTLS,N,N,Y,N,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 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 =0.
 EL1 = GRND4
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0.
 EL1C =0.
 ENUL =1.100667E-05 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0.
 ENUTC =0.
 IENUTA = 3
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRNDTL(H1)=1. ;PRNDTL(LTLS)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 PRT(H1)=0.86 ;PRT(LTLS)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.3E-04 ;FIINIT(U1)=1.3
 FIINIT(V1)=0. ;FIINIT(KE)=6.76
 FIINIT(EP)=31.777433 ;FIINIT(H1)=0.5
 FIINIT(EPKE)=1.0E-10 ;FIINIT(REYN)=1.0E-10
 FIINIT(LTLS)=1.0E-10 ;FIINIT(WDIS)=0.01143
 FIINIT(PRPS)=-1. ;FIINIT(FTWO)=1.0E-10
 FIINIT(FMU)=1.0E-10 ;FIINIT(LEN1)=1.0E-10
 FIINIT(ENUT)=1.0E-10
 
 PATCH(CMP1 ,INIVAL, 1, 14, 1, 40, 1, 1, 1, 1)
 INIT(CMP1 ,PRPS,0. ,198. )
 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(INLET ,WEST , 1, 1, 41, 80, 1, 1, 1, 1)
 COVAL(INLET ,P1 , FIXFLU ,13. )
 COVAL(INLET ,U1 ,0. ,13. )
 COVAL(INLET ,V1 ,0. ,0. )
 COVAL(INLET ,KE ,0. ,0.7605 )
 COVAL(INLET ,EP ,0. ,31.777433 )
 COVAL(INLET ,H1 ,0. ,0. )
 
 PATCH(OUTLET ,EAST , 80, 80, 1, 80, 1, 1, 1, 1)
 COVAL(OUTLET ,P1 ,1.0E+05 ,0. )
 COVAL(OUTLET ,H1 ,0. , SAME )
 
 PATCH(WALLN ,NWALL , 1, 80, 80, 80, 1, 1, 1, 1)
 COVAL(WALLN ,U1 , GRND2 ,0. )
 COVAL(WALLN ,KE ,1. ,0. )
 
 PATCH(WALLS ,SWALL , 15, 80, 1, 1, 1, 1, 1, 1)
 COVAL(WALLS ,U1 , GRND2 ,0. )
 COVAL(WALLS ,KE ,1. ,0. )
 
 PATCH(HEATIN ,SOUTH , 15, 80, 1, 1, 1, 1, 1, 1)
 COVAL(HEATIN ,H1 , FIXFLU ,1.625 )
 
 PATCH(CMP1-WW ,WWALL , 15, 15, 1, 40, 1, 1, 1, 1)
 COVAL(CMP1-WW ,V1 , GRND2 , SAME )
 COVAL(CMP1-WW ,KE ,1. ,0. )
 
 PATCH(CMP1-SW ,SWALL , 1, 14, 41, 41, 1, 1, 1, 1)
 COVAL(CMP1-SW ,U1 , GRND2 , SAME )
 COVAL(CMP1-SW ,KE ,1. ,0. )
 XCYCLE = F
 EGWF = F
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 1000 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=2
 LITER(V1)=2 ;LITER(KE)=5
 LITER(EP)=5 ;LITER(H1)=10
 LITER(LTLS)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(KE)=1.0E-03
 ENDIT(EP)=1.0E-03 ;ENDIT(H1)=1.0E-03
 ENDIT(LTLS)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.5)
 RELAX(U1,FALSDT,7.326923E-03)
 RELAX(V1,FALSDT,7.326923E-03)
 RELAX(KE,LINRLX,0.3)
 RELAX(EP,LINRLX,0.5)
 RELAX(H1,FALSDT,7.326923E-03)
 RELAX(EPKE,LINRLX,1.)
 RELAX(REYN,LINRLX,1.)
 RELAX(LTLS,FALSDT,1.0E+09)
 RELAX(WDIS,LINRLX,1.)
 RELAX(PRPS,LINRLX,1.)
 RELAX(FTWO,LINRLX,1.)
 RELAX(FMU,LINRLX,1.)
 RELAX(LEN1,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 KELIN = 3
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10
 VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
 VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(FTWO)=1.0E+10 ;VARMIN(FTWO)=-1.0E+10
 VARMAX(FMU)=1.0E+10 ;VARMIN(FMU)=-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,MATERIAL,198,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=19.76 ;EX(U1)=6.848
 EX(V1)=0.2304 ;EX(KE)=1.992
 EX(EP)=401.799988 ;EX(H1)=4.132
 EX(EPKE)=251. ;EX(REYN)=2129.
 EX(LTLS)=6.673E-04 ;EX(WDIS)=0.01796
 EX(PRPS)=0.9125 ;EX(FTWO)=0.9083
 EX(FMU)=0.825 ;EX(LEN1)=3.382E-03
 EX(ENUT)=2.477E-03
 ************************************************************
  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(P1,Y,N,Y,Y,Y,Y)
 OUTPUT(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(H1,Y,N,Y,Y,Y,Y)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(REYN,Y,N,Y,N,N,N)
 OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
 OUTPUT(WDIS,Y,N,N,N,N,N)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 OUTPUT(FTWO,Y,N,Y,N,N,N)
 OUTPUT(FMU,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 = 16 ;IYMON = 38 ;IZMON = 1
 NPRMON = 10000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 IPLTF = 1 ;IPLTL = 2000 ;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