TALK=T;RUN(1,1)
  DISPLAY
  Numerical Schemes validation example:
  2-d x-y, Cartesian, steady, elliptic simulation
  The case considered is that of steady, incompressible, laminar
  recirculating flow inside a lid-driven cavity. The geometry is a
  1m square cavity with no inflow or outflow. The flow is driven
  by as moving wall of -1m/s at the top of cavity. The Reynolds
  number based on cavity height is 1000. The calculation may be
  performed with any of the 16 convection schemes in PHOENICS,
  although here the case is set up to be run with one of 3 linear
  schemes or 2 non-linear schemes.
  ENDDIS
  PHOTON USE
   p


  0.20443E+04 0.15633E+04 CR
  gr ou z 1;vec z 1 sh
  stream 2d z 1
  -.113 0 10
  stream 2d z 1
  0. .235e-2 10
  msg streamlines and vectors
  msg press  to end
  pause
  ENDUSE
  AUTOPLOT USE
  file
  phi 5

  d 1 u1 x 11;plot;swap;scale;level x 0
  msg horizontal velocity profile at x/h=0.5
  msg press  to end
  pause
  ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(YX LAMINAR WALL-DRIVEN CAVITY           )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 21
 XULAST =1.
 XFRAC(1)=0.018591 ;XFRAC(2)=0.049062
 XFRAC(3)=0.086552 ;XFRAC(4)=0.129476
 XFRAC(5)=0.176955 ;XFRAC(6)=0.228411
 XFRAC(7)=0.283428 ;XFRAC(8)=0.341689
 XFRAC(9)=0.402944 ;XFRAC(10)=0.466987
 XFRAC(11)=0.533013 ;XFRAC(12)=0.597056
 XFRAC(13)=0.658311 ;XFRAC(14)=0.716572
 XFRAC(15)=0.771589 ;XFRAC(16)=0.823045
 XFRAC(17)=0.870524 ;XFRAC(18)=0.913448
 XFRAC(19)=0.950938 ;XFRAC(20)=0.981409
 XFRAC(21)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 21
 YVLAST =1.
 YFRAC(1)=0.018591 ;YFRAC(2)=0.049062
 YFRAC(3)=0.086552 ;YFRAC(4)=0.129476
 YFRAC(5)=0.176955 ;YFRAC(6)=0.228411
 YFRAC(7)=0.283428 ;YFRAC(8)=0.341689
 YFRAC(9)=0.402944 ;YFRAC(10)=0.466987
 YFRAC(11)=0.533013 ;YFRAC(12)=0.597056
 YFRAC(13)=0.658311 ;YFRAC(14)=0.716572
 YFRAC(15)=0.771589 ;YFRAC(16)=0.823045
 YFRAC(17)=0.870524 ;YFRAC(18)=0.913448
 YFRAC(19)=0.950938 ;YFRAC(20)=0.981409
 YFRAC(21)=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
    * 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,N)
 SOLUTN(V1,Y,Y,N,N,N,N)
 ************************************************************
  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. ;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-03 ;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)=1.0E-10
 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(SOUTH ,SWALL , 1, 21, 1, 1, 1, 1, 1, 1)
 COVAL(SOUTH ,U1 ,1. ,0. )
 
 PATCH(MOVE ,NWALL , 1, 21, 21, 21, 1, 1, 1, 1)
 COVAL(MOVE ,U1 ,1. ,-1. )
 
 PATCH(WEST ,WWALL , 1, 1, 1, 21, 1, 1, 1, 1)
 COVAL(WEST ,V1 ,1. ,0. )
 
 PATCH(EAST ,EWALL , 21, 21, 1, 21, 1, 1, 1, 1)
 COVAL(EAST ,V1 ,1. ,0. )
 
 PATCH(FIXPRESS,CELL , 10, 10, 10, 10, 1, 1, 1, 1)
 COVAL(FIXPRESS,P1 ,1. ,0. )
 COVAL(FIXPRESS,U1 ,0. ,0. )
 COVAL(FIXPRESS,V1 ,0. ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 3000 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=5
 LITER(V1)=5
 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,0.238095)
 RELAX(V1,FALSDT,0.238095)
 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)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  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 = 10 ;IYMON = 20 ;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 = 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
STOP