TALK=T;RUN(1,1)
  DISPLAY
  The case considered is 2D incompressible, turbulent flow
  past a thick flat rectangular plate with a sharp leading
  edge, as studied experimentally by Djilali,N. & Gartshore,I.S.,
  "Turbulent flow around a bluff rectangular plate, Part I:
  Experimental Investigation", Journal of Fluid Mechanics,
  Vol.113, pp51-59, (1991).

  The flow separates at the leading edge of the plate, and then
  reattaches further downstream to form a long separation zone on
  top of the plate. The thickness of the plate H is taken as 0.1m,
  and the working fluid is air at standard temperature and
  pressure.The flow Reynolds number based on the free-stream
  velocity and plate thickness H is 50,000.

  The inlet and outlet planes are located 10H upstream and
  downstream of the leading edge of the plate. The height of the
  solution domain is taken as 10H, and for simplicity a zero
  flux boundary is assumed here. Symmetry is exploited so that
  only one half of the flow is simulated. A fixed-pressure boundary
  condition is applied at the outlet, and uniform flow profiles are
  specified at the inlet. Scalable wall functions are used at
  the walls of the plate.

  The default simulation is made with the realisable k-e model, but
  the case can also be run with the standard k-e and k-w models, the
  revised Wilcox k-w, Menter k-w and k-w-SST model, and the followi$
ng
  k-e variants: Chen-Kim, RNG,  Kato-Launder, Murakami et al. The
  measured and predicted reattachment points are listed below:

           KE   RKE   CK   RNG  MMK   KL   KW   KWR  KWM SST EXPT

    Xr/H=  1.1  4.8  4.7   3.9  3.2  3.2  0.73  3.7 0.94 5.1  4.7

  These results are not grid independent. It can be seen that the
  standard k-e and k-w models seriously underestimate the separation
  length. This is because this model predicts excessive turbulence
  production at the stagnation zone on the front of the plate, and
  the high turbulence levels are then convected into the separation
  zone. The other models perform much better, with the revised
  Wilcox, Chen-Kim and Realisable models showing close agreement
  with the measurement.
  ENDDIS
   AUTOPLOT USE
   file
   phida 3

   da 1 u1 y 11
   divide x .1 1
   shift x -10 1
   colf 1
   level y 0;level x 0
   scale x 0 6
   redraw
   msg Velocity (U1) profile
   msg Press e to END
   ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Realisable_KE Blunt Flat Plate:T311     )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 311
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 115
 XULAST =2.
 XFRAC(1)=0.021251 ;XFRAC(24)=0.336972
 XFRAC(47)=0.465068 ;XFRAC(70)=0.531403
 XFRAC(93)=0.673179
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 75
 YVLAST =1.
 YFRAC(1)=4.164548E-03 ;YFRAC(16)=0.071363
 YFRAC(31)=0.152963 ;YFRAC(46)=0.29992
 YFRAC(61)=0.564581
 ************************************************************
  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(137)=PRPS
 NAME(138)=SKIN ;NAME(139)=STRS
 NAME(140)=YPLS ;NAME(141)=C1E
 NAME(142)=DVDY ;NAME(143)=DVDX
 NAME(144)=DUDY ;NAME(145)=DUDX
 NAME(146)=EPKE ;NAME(147)=CMU
 NAME(148)=DEN1 ;NAME(149)=EL1
 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,N,N,N,Y)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(PRPS,Y,N,N,N,N,Y)
 SOLUTN(SKIN,Y,N,N,N,N,Y)
 SOLUTN(STRS,Y,N,N,N,N,Y)
 SOLUTN(YPLS,Y,N,N,N,N,Y)
 SOLUTN(C1E,Y,N,N,N,N,Y)
 SOLUTN(DVDY,Y,N,N,N,N,N)
 SOLUTN(DVDX,Y,N,N,N,N,N)
 SOLUTN(DUDY,Y,N,N,N,N,N)
 SOLUTN(DUDX,Y,N,N,N,N,N)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(CMU,Y,N,N,N,N,N)
 SOLUTN(DEN1,Y,N,N,N,N,Y)
 SOLUTN(EL1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 DEN1 = 148
 VIST = 150
 LEN1 = 149
 PRPS = 137
 ************************************************************
  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)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWENT = T
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 DENPCO = T
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1.225 ;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.46E-05 ;ENUT = GRND5
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 14
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(KE)=1. ;PRT(EP)=1.2
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=7.3
 FIINIT(V1)=7.300001E-03 ;FIINIT(KE)=5.329E-03
 FIINIT(EP)=1.27831E-03 ;FIINIT(PRPS)=-1.
 FIINIT(SKIN)=1.0E-10 ;FIINIT(STRS)=1.0E-10
 FIINIT(YPLS)=1.0E-10 ;FIINIT(C1E)=1.0E-10
 FIINIT(DVDY)=1.0E-10 ;FIINIT(DVDX)=1.0E-10
 FIINIT(DUDY)=1.0E-10 ;FIINIT(DUDX)=1.0E-10
 FIINIT(EPKE)=1.0E-10 ;FIINIT(CMU)=0.09
 FIINIT(DEN1)=1.0E-10 ;FIINIT(EL1)=1.0E-10
 FIINIT(ENUT)=1.0E-10
 
 PATCH(PLATE ,INIVAL, 61, 115, 1, 10, 1, 1, 1, 1)
 INIT(PLATE ,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 )
 
 PATCH(REKESO ,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(REKESO ,EP , GRND4 , GRND4 )
 
 PATCH(INLET ,WEST , 1, 1, 1, 75, 1, 1, 1, 1)
 COVAL(INLET ,P1 , FIXFLU ,8.9425 )
 COVAL(INLET ,U1 ,0. ,7.3 )
 COVAL(INLET ,V1 ,0. ,0. )
 COVAL(INLET ,KE ,0. ,5.329E-03 )
 COVAL(INLET ,EP ,0. ,1.27831E-03 )
 
 PATCH(OUTLET ,EAST , 115, 115, 11, 75, 1, 1, 1, 1)
 COVAL(OUTLET ,P1 ,1000. ,0. )
 COVAL(OUTLET ,U1 ,0. ,0. )
 COVAL(OUTLET ,V1 ,0. ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 SCALWF = T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 1500 ;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)=10
 LITER(V1)=10 ;LITER(KE)=20
 LITER(EP)=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
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(KE,LINRLX,0.5)
 RELAX(EP,LINRLX,0.5)
 RELAX(PRPS,LINRLX,1.)
 RELAX(SKIN,LINRLX,1.)
 RELAX(STRS,LINRLX,1.)
 RELAX(YPLS,LINRLX,1.)
 RELAX(C1E,LINRLX,1.)
 RELAX(DVDY,LINRLX,1.)
 RELAX(DVDX,LINRLX,1.)
 RELAX(DUDY,LINRLX,1.)
 RELAX(DUDX,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(CMU,LINRLX,0.5)
 RELAX(DEN1,LINRLX,1.)
 RELAX(EL1,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(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(SKIN)=1.0E+10 ;VARMIN(SKIN)=-1.0E+10
 VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-1.0E+10
 VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10
 VARMAX(C1E)=1.0E+10 ;VARMIN(C1E)=-1.0E+10
 VARMAX(DVDY)=1.0E+10 ;VARMIN(DVDY)=-1.0E+10
 VARMAX(DVDX)=1.0E+10 ;VARMIN(DVDX)=-1.0E+10
 VARMAX(DUDY)=1.0E+10 ;VARMIN(DUDY)=-1.0E+10
 VARMAX(DUDX)=1.0E+10 ;VARMIN(DUDX)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(CMU)=1.0E+10 ;VARMIN(CMU)=-1.0E+10
 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
 VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10
 VARMAX(ENUT)=0.1 ;VARMIN(ENUT)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 CONWIZ = T
 IENUTA = 14
 ISG62 = 1
 SPEDAT(SET,KECONST,C2E,R,1.9)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,MATERIAL,198,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=5.421 ;EX(U1)=6.629
 EX(V1)=0.4262 ;EX(KE)=0.123
 EX(EP)=4.972 ;EX(PRPS)=0.9362
 EX(SKIN)=5.051E-05 ;EX(STRS)=2.497E-04
 EX(YPLS)=0.1502 ;EX(C1E)=0.5868
 EX(DVDY)=5.134 ;EX(DVDX)=3.958
 EX(DUDY)=18.360001 ;EX(DUDX)=5.133
 EX(EPKE)=5.302 ;EX(CMU)=0.07094
 EX(DEN1)=1.147 ;EX(EL1)=0.04079
 EX(ENUT)=1.602E-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(PRPS,Y,N,Y,N,N,N)
 OUTPUT(SKIN,Y,N,Y,N,N,N)
 OUTPUT(STRS,Y,N,Y,N,N,N)
 OUTPUT(YPLS,Y,N,Y,N,N,N)
 OUTPUT(C1E,Y,N,Y,N,Y,Y)
 OUTPUT(DVDY,Y,N,Y,N,N,N)
 OUTPUT(DVDX,Y,N,Y,N,N,N)
 OUTPUT(DUDY,Y,N,Y,N,N,N)
 OUTPUT(DUDX,Y,N,Y,N,N,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(CMU,Y,N,Y,N,Y,Y)
 OUTPUT(DEN1,Y,N,Y,N,N,N)
 OUTPUT(EL1,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 73 ;IYMON = 12 ;IZMON = 1
 NPRMON = 100 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NXPRIN = 23 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = 15 ;IYPRF = 1 ;IYPRL = 10000
 IPLTF = 1 ;IPLTL = 1500 ;NPLT = 10
 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
 IDISPA = 50 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='SW'
STOP