```
TALK=F;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.6, Date 31/01/06
CPVNAM=VDI;SPPNAM=Core
************************************************************
title
DISPLAY

This case solves a two-dimensional steady hydrodynamics
problem about flow in labirinth channel.

The analytical solution is absent. The solution can be compared
with solution of Structured PHOENICS.

Cases changes inlet velocity as Uin=0.1*caseno

Domain material is air.
ENDDIS

************************************************************
IRUNN   =         1 ;LIBREF =         0
************************************************************
Group 1. Run Title
TEXT(USP labirinth.                          )
************************************************************
Group 2. Transience
************************************************************
Groups 3, 4, 5  Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
INTEGER(NCUBE,NCHAN)
NCUBE = 80
NCHAN = NCUBE/20
RSET(M,NCUBE,NCUBE,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
NAME(148) = PRPS
* Solved variables list
SOLVE(P1  ,U1  ,V1)
* Stored variables list
STORE(PRPS)
SOLUTN(P1  ,Y,Y,Y,N,N,Y)

************************************************************
Group 8. Terms & Devices
Usp related variables
USP    = T
UAUTO  = F
USPDBG = F
UTCPLT = F
USPVTK = T
USPIMB = F
MXLEV  = 0
MYLEV  = 0
MZLEV  = 0
DOMAT  = -1
MINPRP = -1
MAXPRP = 250
CELLST = 100
FACEST = 100
PARSOL = F
************************************************************
Group 9. Properties
PRESS0  = 1.000000E+05 ;TEMP0  = 2.730000E+02
* Domain material index is   0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
RHO1 = 1.
ENUL = 1.e-5
REAL(Uin)
Uin=0.1
REAL(REYNO)
REYNO = UIN*0.05/ENUL
mesg(REYNO = :REYNO:
DVO1DT  = 3.410000E-03
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
No PATCHes used for this Group
PATCH(B1,INIVAL, 4*NX/5+NCHAN+1,NX,1,NY,1,NZ,1,1)
COVAL(B1, PRPS, 0.0, 198)
PATCH(B2,INIVAL, 3*NX/5+NCHAN+1,4*NX/5+NCHAN,NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B2, PRPS, 0.0, 198)
PATCH(B3,INIVAL, 2*NX/5+NCHAN+1,3*NX/5+NCHAN,2*NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B3, PRPS, 0.0, 198)
PATCH(B4,INIVAL, NX/5+NCHAN+1,2*NX/5+NCHAN,3*NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B4, PRPS, 0.0, 198)
PATCH(B5,INIVAL, 1,NX/5+NCHAN,4*NY/5+NCHAN+1,NY,1,NZ,1,1)
COVAL(B5, PRPS, 0.0, 198)
PATCH(B6,INIVAL, 1,4*NX/5,1,NY/5,1,NZ,1,1)
COVAL(B6, PRPS, 0.0, 198)
PATCH(B7,INIVAL, 1,3*NX/5,NY/5+1,2*NY/5,1,NZ,1,1)
COVAL(B7, PRPS, 0.0, 198)
PATCH(B8,INIVAL, 1,2*NX/5,2*NY/5+1,3*NY/5,1,NZ,1,1)
COVAL(B8, PRPS, 0.0, 198)
PATCH(B9,INIVAL, 1,NX/5,3*NY/5+1,4*NY/5,1,NZ,1,1)
COVAL(B9, PRPS, 0.0, 198)
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
No PATCHes used for this Group

EGWF    =    T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP  =       1000
RESFAC  = 1.000000E-06
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(P1  ,LINRLX, 8.000000E-01)
RELAX(U1  ,FALSDT, 1.000000E+00)
RELAX(V1  ,FALSDT, 1.000000E+00)
mesg(Do you want to use collocated arrangement (y) or staggered one (n)? (y/n)
if(:ans:.eq.y)then
SPEDAT(SET,USP,METHOD,I,1)
RELAX(P1  ,LINRLX, 5.000000E-01)
RELAX(U1  ,FALSDT, 1.000000E-01)
RELAX(V1  ,FALSDT, 1.000000E-01)
mesg(Do you want to use SIMPLEST (y) or SIMPLE (n)? (y/n)
if(:ans:.eq.y)then
SPEDAT(SET,USP,SIMPLEST,L,T)
endif

endif

************************************************************
Group 18. Limits
VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
ASAP    =    T
PARSOL  =    F
************************************************************
Group 20. Preliminary Printout
ECHO    =    T
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON   =        50 ;IYMON  =        42 ;IZMON  =         1
NPRMON  =    100000
NPRMNT  =         1
TSTSWP  =        -1
mesg(Do you want to use max&min(y) on monitor instead of default(n)? (y/n)
if(:ans:.eq.y)then
#maxmin
endif

mesg(Do you want to view results in the centres of cells? (y/n)
if(:ans:.eq.y)then
SPEDAT(SET,USPIO,VERTCENT,L,F)
endif

************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT  =    100000
ISWPRF  =         1 ;ISWPRL =    100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts

GVIEW(P,0.000000E+00,0.000000E+00,-1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)

> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       4.950000E-01, 4.150000E-01, 5.000000E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02

> OBJ,    NAME,        IN
> OBJ,    POSITION,    0.000000E+00, 8.000000E-01, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 5.000000E-02, 1.000000E+00
> OBJ,    GEOMETRY,    cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    VELOCITY,     Uin, 0.000000E+00, 0.000000E+00
> OBJ,    TEMPERATURE,  0.000000E+00

> OBJ,    NAME,        OUT
> OBJ,    POSITION,    8.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 0.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    TEMPERATURE,  SAME
> OBJ,    COEFFICIENT,  1.000000E+07
STOP
```