TALK=T;RUN(1,1)
  DISPLAY
  Numerical Schemes validation example:
  2-d x-y, Cartesian, steady, elliptic simulation
  The case considered is steady incompressible, laminar backward-
  facing-step flow, i.e. flow through a straight channel having a
  sudden asymmetric expansion. The flow is characterised by the
  step height H and the Reynolds number Re based on the bulk inlet
  inlet velocity and 2h, where h is the flow inlet height. The
  channel expansion ratio is 1.94, and the total length of the
  domain is 40H. A fully-developed parabolic velocity profile is
  prescribed at the inlet boundary located at the step. The
  calculation is made for Re=150 with the cubic upwind scheme, and
  the option exists to select the hybrid or van Albada scheme.
  ENDDIS
  PHOTON USE
   P


  0.20443E+04 0.15633E+04 CR
  gr ou z 1;mag gr 5
  0.28838E+03 0.17522E+04 CR
  msg Reynolds number = 150  Cubic upwind scheme
  vec z 1 x 1 40 y 1 m sh
  STREAM 2D Z 1 X 1 40 Y 1 M
  -.699 0. 5
  STREAM 2D Z 1 X 1 40 Y 1 M
  0. 12. 8
  msg press  to continue
  msg press  to end
  pause
  ENDUSE
  AUTOPLOT USE
  file
  phi 5

  d 1 u1 y 1;plot;div x .49 1
  scale x 0 6;level y 0
  msg Reynolds number = 150  Cubic upwind scheme
  msg horizontal (U1) velocity distribution along bottom wall
  msg reattachment point is where U1 crosses zero point
  msg press  to continue
  msg press  to end
  pause
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(YX LAMINAR BACKWARD-FACING-STEP FLOW    )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 200
 XULAST =19.6
 XFRAC(1)=4.166662E-03 ;XFRAC(2)=8.333324E-03
 XFRAC(3)=0.0125 ;XFRAC(4)=0.016667
 XFRAC(5)=0.020833 ;XFRAC(6)=0.025
 XFRAC(7)=0.029167 ;XFRAC(8)=0.033333
 XFRAC(9)=0.0375 ;XFRAC(10)=0.041667
 XFRAC(11)=0.045833 ;XFRAC(12)=0.05
 XFRAC(13)=0.054167 ;XFRAC(14)=0.058333
 XFRAC(15)=0.0625 ;XFRAC(16)=0.066667
 XFRAC(17)=0.070833 ;XFRAC(18)=0.075
 XFRAC(19)=0.079167 ;XFRAC(20)=0.083333
 XFRAC(21)=0.0875 ;XFRAC(22)=0.091667
 XFRAC(23)=0.095833 ;XFRAC(24)=0.1
 XFRAC(25)=0.104167 ;XFRAC(26)=0.108333
 XFRAC(27)=0.1125 ;XFRAC(28)=0.116667
 XFRAC(29)=0.120833 ;XFRAC(30)=0.125
 XFRAC(31)=0.129167 ;XFRAC(32)=0.133333
 XFRAC(33)=0.1375 ;XFRAC(34)=0.141667
 XFRAC(35)=0.145833 ;XFRAC(36)=0.15
 XFRAC(37)=0.154167 ;XFRAC(38)=0.158333
 XFRAC(39)=0.1625 ;XFRAC(40)=0.166666
 XFRAC(41)=0.170833 ;XFRAC(42)=0.175
 XFRAC(43)=0.179166 ;XFRAC(44)=0.183333
 XFRAC(45)=0.1875 ;XFRAC(46)=0.191666
 XFRAC(47)=0.195833 ;XFRAC(48)=0.2
 XFRAC(49)=0.204166 ;XFRAC(50)=0.208333
 XFRAC(51)=0.2125 ;XFRAC(52)=0.216666
 XFRAC(53)=0.220833 ;XFRAC(54)=0.225
 XFRAC(55)=0.229166 ;XFRAC(56)=0.233333
 XFRAC(57)=0.2375 ;XFRAC(58)=0.241666
 XFRAC(59)=0.245833 ;XFRAC(60)=0.25
 XFRAC(61)=0.254166 ;XFRAC(62)=0.258333
 XFRAC(63)=0.2625 ;XFRAC(64)=0.266666
 XFRAC(65)=0.270833 ;XFRAC(66)=0.275
 XFRAC(67)=0.279166 ;XFRAC(68)=0.283333
 XFRAC(69)=0.2875 ;XFRAC(70)=0.291666
 XFRAC(71)=0.295833 ;XFRAC(72)=0.3
 XFRAC(73)=0.304166 ;XFRAC(74)=0.308333
 XFRAC(75)=0.3125 ;XFRAC(76)=0.316666
 XFRAC(77)=0.320833 ;XFRAC(78)=0.325
 XFRAC(79)=0.329166 ;XFRAC(80)=0.333333
 XFRAC(81)=0.3375 ;XFRAC(82)=0.341666
 XFRAC(83)=0.345833 ;XFRAC(84)=0.35
 XFRAC(85)=0.354166 ;XFRAC(86)=0.358333
 XFRAC(87)=0.3625 ;XFRAC(88)=0.366666
 XFRAC(89)=0.370833 ;XFRAC(90)=0.375
 XFRAC(91)=0.379166 ;XFRAC(92)=0.383333
 XFRAC(93)=0.3875 ;XFRAC(94)=0.391666
 XFRAC(95)=0.395833 ;XFRAC(96)=0.4
 XFRAC(97)=0.404166 ;XFRAC(98)=0.408333
 XFRAC(99)=0.4125 ;XFRAC(100)=0.416666
 XFRAC(101)=0.420833 ;XFRAC(102)=0.425
 XFRAC(103)=0.429166 ;XFRAC(104)=0.433333
 XFRAC(105)=0.4375 ;XFRAC(106)=0.441666
 XFRAC(107)=0.445833 ;XFRAC(108)=0.45
 XFRAC(109)=0.454166 ;XFRAC(110)=0.458333
 XFRAC(111)=0.4625 ;XFRAC(112)=0.466666
 XFRAC(113)=0.470833 ;XFRAC(114)=0.474999
 XFRAC(115)=0.479166 ;XFRAC(116)=0.483333
 XFRAC(117)=0.487499 ;XFRAC(118)=0.491666
 XFRAC(119)=0.495833 ;XFRAC(120)=0.499999
 XFRAC(121)=0.504166 ;XFRAC(122)=0.508333
 XFRAC(123)=0.512499 ;XFRAC(124)=0.516666
 XFRAC(125)=0.520833 ;XFRAC(126)=0.524999
 XFRAC(127)=0.529166 ;XFRAC(128)=0.533333
 XFRAC(129)=0.537499 ;XFRAC(130)=0.541666
 XFRAC(131)=0.545833 ;XFRAC(132)=0.549999
 XFRAC(133)=0.554166 ;XFRAC(134)=0.558333
 XFRAC(135)=0.562499 ;XFRAC(136)=0.566666
 XFRAC(137)=0.570833 ;XFRAC(138)=0.574999
 XFRAC(139)=0.579166 ;XFRAC(140)=0.583333
 XFRAC(141)=0.587499 ;XFRAC(142)=0.591666
 XFRAC(143)=0.595833 ;XFRAC(144)=0.599999
 XFRAC(145)=0.604166 ;XFRAC(146)=0.608333
 XFRAC(147)=0.612499 ;XFRAC(148)=0.616666
 XFRAC(149)=0.620833 ;XFRAC(150)=0.624999
 XFRAC(151)=0.629166 ;XFRAC(152)=0.633333
 XFRAC(153)=0.637499 ;XFRAC(154)=0.641666
 XFRAC(155)=0.645833 ;XFRAC(156)=0.649999
 XFRAC(157)=0.654166 ;XFRAC(158)=0.658333
 XFRAC(159)=0.662499 ;XFRAC(160)=0.666666
 XFRAC(161)=0.670651 ;XFRAC(162)=0.675821
 XFRAC(163)=0.681558 ;XFRAC(164)=0.687698
 XFRAC(165)=0.694156 ;XFRAC(166)=0.700879
 XFRAC(167)=0.707831 ;XFRAC(168)=0.714985
 XFRAC(169)=0.72232 ;XFRAC(170)=0.729821
 XFRAC(171)=0.737473 ;XFRAC(172)=0.745266
 XFRAC(173)=0.753191 ;XFRAC(174)=0.761238
 XFRAC(175)=0.769401 ;XFRAC(176)=0.777673
 XFRAC(177)=0.78605 ;XFRAC(178)=0.794526
 XFRAC(179)=0.803097 ;XFRAC(180)=0.811758
 XFRAC(181)=0.820507 ;XFRAC(182)=0.829339
 XFRAC(183)=0.838252 ;XFRAC(184)=0.847242
 XFRAC(185)=0.856308 ;XFRAC(186)=0.865447
 XFRAC(187)=0.874657 ;XFRAC(188)=0.883935
 XFRAC(189)=0.893279 ;XFRAC(190)=0.902688
 XFRAC(191)=0.91216 ;XFRAC(192)=0.921694
 XFRAC(193)=0.931287 ;XFRAC(194)=0.940939
 XFRAC(195)=0.950647 ;XFRAC(196)=0.960411
 XFRAC(197)=0.97023 ;XFRAC(198)=0.980101
 XFRAC(199)=0.990025 ;XFRAC(200)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 32
 YVLAST =1.01
 YFRAC(1)=0.030322 ;YFRAC(2)=0.060644
 YFRAC(3)=0.090965 ;YFRAC(4)=0.121287
 YFRAC(5)=0.151609 ;YFRAC(6)=0.181931
 YFRAC(7)=0.212252 ;YFRAC(8)=0.242574
 YFRAC(9)=0.272896 ;YFRAC(10)=0.303218
 YFRAC(11)=0.33354 ;YFRAC(12)=0.363861
 YFRAC(13)=0.394183 ;YFRAC(14)=0.424505
 YFRAC(15)=0.454827 ;YFRAC(16)=0.485149
 YFRAC(17)=0.517327 ;YFRAC(18)=0.549505
 YFRAC(19)=0.581683 ;YFRAC(20)=0.613861
 YFRAC(21)=0.64604 ;YFRAC(22)=0.678218
 YFRAC(23)=0.710396 ;YFRAC(24)=0.742574
 YFRAC(25)=0.774752 ;YFRAC(26)=0.806931
 YFRAC(27)=0.839109 ;YFRAC(28)=0.871287
 YFRAC(29)=0.903465 ;YFRAC(30)=0.935644
 YFRAC(31)=0.967822 ;YFRAC(32)=1.
 ************************************************************
  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(1)=P1 ;NAME(3)=U1
 NAME(5)=V1
    * 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)
 ************************************************************
  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)
 DIFCUT =0. ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 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.2E-03 ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =0.16 ;ENUT =0.
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=11.538462
 FIINIT(V1)=1.0E-10
   No PATCHes yet used for this Group
 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(HOCS ,CELL , 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(HOCS ,U1 , FIXFLU , GRND1 )
 COVAL(HOCS ,V1 , FIXFLU , GRND1 )
 
 PATCH(IN1 ,WEST , 1, 1, 17, 17, 1, 1, 1, 1)
 COVAL(IN1 ,P1 , FIXFLU ,5.030048E-03 )
 COVAL(IN1 ,U1 ,0. ,4.191707 )
 COVAL(IN1 ,V1 ,0. ,0. )
 
 PATCH(IN2 ,WEST , 1, 1, 18, 18, 1, 1, 1, 1)
 COVAL(IN2 ,P1 , FIXFLU ,0.014117 )
 COVAL(IN2 ,U1 ,0. ,11.763823 )
 COVAL(IN2 ,V1 ,0. ,0. )
 
 PATCH(IN3 ,WEST , 1, 1, 19, 19, 1, 1, 1, 1)
 COVAL(IN3 ,P1 , FIXFLU ,0.021905 )
 COVAL(IN3 ,U1 ,0. ,18.254206 )
 COVAL(IN3 ,V1 ,0. ,0. )
 
 PATCH(IN4 ,WEST , 1, 1, 20, 20, 1, 1, 1, 1)
 COVAL(IN4 ,P1 , FIXFLU ,0.028395 )
 COVAL(IN4 ,U1 ,0. ,23.662863 )
 COVAL(IN4 ,V1 ,0. ,0. )
 
 PATCH(IN5 ,WEST , 1, 1, 21, 21, 1, 1, 1, 1)
 COVAL(IN5 ,P1 , FIXFLU ,0.033588 )
 COVAL(IN5 ,U1 ,0. ,27.989784 )
 COVAL(IN5 ,V1 ,0. ,0. )
 
 PATCH(IN6 ,WEST , 1, 1, 22, 22, 1, 1, 1, 1)
 COVAL(IN6 ,P1 , FIXFLU ,0.037482 )
 COVAL(IN6 ,U1 ,0. ,31.234976 )
 COVAL(IN6 ,V1 ,0. ,0. )
 
 PATCH(IN7 ,WEST , 1, 1, 23, 23, 1, 1, 1, 1)
 COVAL(IN7 ,P1 , FIXFLU ,0.040078 )
 COVAL(IN7 ,U1 ,0. ,33.398438 )
 COVAL(IN7 ,V1 ,0. ,0. )
 
 PATCH(IN8 ,WEST , 1, 1, 24, 24, 1, 1, 1, 1)
 COVAL(IN8 ,P1 , FIXFLU ,0.041376 )
 COVAL(IN8 ,U1 ,0. ,34.480171 )
 COVAL(IN8 ,V1 ,0. ,0. )
 
 PATCH(IN9 ,WEST , 1, 1, 25, 25, 1, 1, 1, 1)
 COVAL(IN9 ,P1 , FIXFLU ,0.041376 )
 COVAL(IN9 ,U1 ,0. ,34.480171 )
 COVAL(IN9 ,V1 ,0. ,0. )
 
 PATCH(IN10 ,WEST , 1, 1, 26, 26, 1, 1, 1, 1)
 COVAL(IN10 ,P1 , FIXFLU ,0.040078 )
 COVAL(IN10 ,U1 ,0. ,33.398441 )
 COVAL(IN10 ,V1 ,0. ,0. )
 
 PATCH(IN11 ,WEST , 1, 1, 27, 27, 1, 1, 1, 1)
 COVAL(IN11 ,P1 , FIXFLU ,0.037482 )
 COVAL(IN11 ,U1 ,0. ,31.234976 )
 COVAL(IN11 ,V1 ,0. ,0. )
 
 PATCH(IN12 ,WEST , 1, 1, 28, 28, 1, 1, 1, 1)
 COVAL(IN12 ,P1 , FIXFLU ,0.033588 )
 COVAL(IN12 ,U1 ,0. ,27.989786 )
 COVAL(IN12 ,V1 ,0. ,0. )
 
 PATCH(IN13 ,WEST , 1, 1, 29, 29, 1, 1, 1, 1)
 COVAL(IN13 ,P1 , FIXFLU ,0.028395 )
 COVAL(IN13 ,U1 ,0. ,23.662859 )
 COVAL(IN13 ,V1 ,0. ,0. )
 
 PATCH(IN14 ,WEST , 1, 1, 30, 30, 1, 1, 1, 1)
 COVAL(IN14 ,P1 , FIXFLU ,0.021905 )
 COVAL(IN14 ,U1 ,0. ,18.254211 )
 COVAL(IN14 ,V1 ,0. ,0. )
 
 PATCH(IN15 ,WEST , 1, 1, 31, 31, 1, 1, 1, 1)
 COVAL(IN15 ,P1 , FIXFLU ,0.014117 )
 COVAL(IN15 ,U1 ,0. ,11.76382 )
 COVAL(IN15 ,V1 ,0. ,0. )
 
 PATCH(IN16 ,WEST , 1, 1, 32, 32, 1, 1, 1, 1)
 COVAL(IN16 ,P1 , FIXFLU ,5.030053E-03 )
 COVAL(IN16 ,U1 ,0. ,4.19171 )
 COVAL(IN16 ,V1 ,0. ,0. )
 
 PATCH(TWALL ,NWALL , 1, 200, 32, 32, 1, 1, 1, 1)
 COVAL(TWALL ,U1 ,1. ,0. )
 
 PATCH(BWALL ,SWALL , 1, 200, 1, 1, 1, 1, 1, 1)
 COVAL(BWALL ,U1 ,1. ,0. )
 
 PATCH(STPWL ,WWALL , 1, 1, 1, 16, 1, 1, 1, 1)
 COVAL(STPWL ,V1 ,1. ,0. )
 
 PATCH(OUT ,EAST , 200, 200, 1, 32, 1, 1, 1, 1)
 COVAL(OUT ,P1 ,10. ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 100 ;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
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,4.246667E-03)
 RELAX(V1,FALSDT,4.246667E-03)
 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
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,SCHEME,INLCS003,I,3)
 SPEDAT(SET,SCHEME,INLCS005,I,3)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  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)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 130 ;IYMON = 29 ;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
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 IPLTF = 1 ;IPLTL = -1 ;NPLT = 20
 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