TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2009, Date 17/02/09
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  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
 
  Uin=0.1*caseno
  ENDDIS
    
REAL(Uin)
INTEGER(caseno)
caseno=1
Uin=0.1*caseno
 ************************************************************
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 1. Run Title
 TEXT(USP labirinth with divided grid.        )
 ************************************************************
  Group 2. Transience
 STEADY = T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,10,10,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)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
 USP = T ; USPDBG= F
 UTCPLT = F ; USPIMB= F
 UAUTO = T ; UGTGRD= F
 USPVTK = T ; USPTER= F
 CELLST  =        10 ;FACEST =        10
 USCANX = T ; USCANY= T
 USCANZ = T
 MXLEV = 0 ;MYLEV = 0 ;MZLEV = 0
 DOMAT = -1
 MINPRP = -1 ; MAXPRP= 100
 UMXLEV = 5 ; UMXCLS= 10000
 UNOLAY = 1 ; UIOLEV= 0
 UWALEV = 0 ; UGRITR= 2
 UMXCRS = 16
 
 ************************************************************
  Group 9. Properties
 PRESS0 =1.0E+05 ;TEMP0 =273.
    * Domain material index is   0 signifying:
    * Air at 20 deg C, 1 atm, treated as incompressible
 SETPRPS(1, 0)
 ENUT =0.
 DVO1DT =3.41E-03
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(PRPS)=-1.
   No PATCHes used for this Group
 
 
 INIADD = F
 ************************************************************
  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.0E-06
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,0.7 )
 RELAX(U1  ,FALSDT, 1.000000E-01)
 RELAX(V1  ,FALSDT, 1.000000E-01)
mesg(Do you want to use collocated arrangement (y) or staggered one (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
 SPEDAT(SET,USP,METHOD,I,1)
 RELAX(P1  ,LINRLX, 3.000000E-01)
 RELAX(U1  ,FALSDT, 5.000000E-02)
 RELAX(V1  ,FALSDT, 5.000000E-02)
mesg(Do you want to use SIMPLEST (y) or SIMPLE (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
 SPEDAT(SET,USP,SIMPLEST,L,T)
endif

endif


 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 PARSOL = F
 ISG62 = 0
 ISG60 = 1
 SPEDAT(SET,USP,METHOD,I,1)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 5 ;IYMON = 4 ;IZMON = 1
 NPRMON = 100000
 NPRMNT = 1
 TSTSWP = -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000
 ISWPRF = 1 ;ISWPRL = 100000
   No PATCHes used for this Group

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


 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,-0.042817,0.054706,0.997584)
 GVIEW(UP,0.017703,0.998385,-0.05399)
 
> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       4.750000E-01, 3.750000E-01, 5.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
  > GRID,   RSET_X_1,     10, 1.000000E+00
  > GRID,   RSET_Y_1,     10, 1.000000E+00
  > GRID,   RSET_Z_1,      1, 1.000000E+00
 
> 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,    TYPE,        INLET
> OBJ,    GRID,        NO
> OBJ,    PRESSURE,    0.
> OBJ,    VELOCITY,    0.1 ,0. ,0.
 
> 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,    TYPE,        OUTLET
> OBJ,    GRID,        NO
> OBJ,    PRESSURE,    0.
> OBJ,    COEFFICIENT, 1000.
 
> OBJ,    NAME,        B3
> OBJ,    POSITION,    0.000000E+00, 8.500000E-01, 0.000000E+00
> OBJ,    SIZE,        TO_END,       TO_END,       TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        B4
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E-01, 8.000000E-01, 1.000000E+00
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        B5
> OBJ,    POSITION,    2.500000E-01, 6.000000E-01, 0.000000E+00
> OBJ,    SIZE,        TO_END,       2.500000E-01, 1.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        B6
> OBJ,    POSITION,    2.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.500000E-01, 5.500000E-01, TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        B7
> OBJ,    POSITION,    5.000000E-01, 4.000000E-01, 0.000000E+00
> OBJ,    SIZE,        TO_END,       2.000000E-01, 1.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        B8
> OBJ,    POSITION,    4.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E-01, 3.500000E-01, 1.000000E+00
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        B9
> OBJ,    POSITION,    7.000000E-01, 2.000000E-01, 0.000000E+00
> OBJ,    SIZE,        TO_END,       2.000000E-01, 1.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        B10
> OBJ,    POSITION,    6.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.500000E-01, 1.500000E-01, 1.000000E+00
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        B11
> OBJ,    POSITION,    8.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        TO_END,       2.000000E-01, 1.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
STOP