TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2009, Date 17/02/09
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 1. Run Title
 TEXT(Flow in pipe. Effective viscosity.     )
 ************************************************************
  Group 2. Transience
 STEADY = T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
mesg(Do you want to use uniform grid (y) or AGG(n)?
readvdu(ans,char,n)
if(:ans:.eq.y)then
 RSET(M,24,48,18,2.0E-03)
 UAUTO = F
else
 RSET(M,8,16,6,2.0E-03)
 UAUTO = T
endif
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 SOLVE(P1,U1,V1,W1)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
 USP = T ; USPDBG= F
 UTCPLT = T ; USPIMB= F
 UGTGRD= F
 USPVTK = T ; USPTER= F
 USCANX = T ; USCANY= T
 USCANZ = T
 MXLEV = 0 ;MYLEV = 0 ;MZLEV = 0
 DOMAT = 67
 CELLST = 100 ; FACEST= 100
 MINPRP = -1 ; MAXPRP= 250
 UMXLEV = 3 ; UMXCLS= 200000
 UNOLAY = 2 ; UIOLEV= 1
 UWALEV = 1 ; UGRITR= 2
 UMXCRS = 16
 ************************************************************
  Group 9. Properties
 PRESS0 =1.0E+05 ;TEMP0 =273.
    * Domain material index is 198 signifying:
    * Solid with smooth-wall friction
 SETPRPS(1,198)
 ENUT =3.e-2
 RHO1=1000.0
 
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(V1)=-20.
   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 = 700
 RESFAC =1.0E-06
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,0.4 )
 RELAX(U1 ,FALSDT,1.E-02 )
 RELAX(V1 ,FALSDT,1.E-02 )
 RELAX(W1 ,FALSDT,1.E-02)
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, 0.2)
 RELAX(U1  ,FALSDT, 1.E-04)
 RELAX(V1  ,FALSDT, 1.E-04)
 RELAX(W1  ,FALSDT, 1.E-04)
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
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 5 ;IYMON = 8 ;IZMON = 6
 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
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,-0.939341,0.135782,0.314966)
 GVIEW(UP,0.337958,0.209726,0.917497)
 
> DOM,    SIZE,        2.000000E+00, 6.593945E+00, 1.889835E+00
> DOM,    MONIT,       1.162625E+00, 2.936366E+00, 1.667376E+00
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
if(UAUTO)then
  > GRID,   RSET_X_1,      4, 1.000000E+00
  > GRID,   RSET_X_2,      2, 1.000000E+00
  > GRID,   RSET_X_3,      2, 1.000000E+00
  > GRID,   RSET_Y_1,      3, 1.000000E+00
  > GRID,   RSET_Y_2,     13, 1.000000E+00
  > GRID,   RSET_Z_1,      3, 1.000000E+00
  > GRID,   RSET_Z_2,      1, 1.000000E+00
  > GRID,   RSET_Z_3,      2, 1.000000E+00
else
  > GRID,   RSET_X_1,      12, 1.000000E+00
  > GRID,   RSET_X_2,      6, 1.000000E+00
  > GRID,   RSET_X_3,      6, 1.000000E+00
  > GRID,   RSET_Y_1,      9, 1.000000E+00
  > GRID,   RSET_Y_2,     39, 1.000000E+00
  > GRID,   RSET_Z_1,      9, 1.000000E+00
  > GRID,   RSET_Z_2,      3, 1.000000E+00
  > GRID,   RSET_Z_3,      6, 1.000000E+00
endif
> OBJ,    NAME,        B2
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    cube12t
> OBJ,    TYPE,        OUTLET
> OBJ,    WIREFRAME,   YES
> OBJ,    PRESSURE,    0.
> OBJ,    TEMPERATURE, SAME
> OBJ,    COEFFICIENT, 1000.0
 
> OBJ,    NAME,        B3
> OBJ,    POSITION,    1.000000E+00, 6.593945E+00, 8.898350E-01
> OBJ,    SIZE,        1.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    cube3t
> OBJ,    TYPE,        INLET
> OBJ,    WIREFRAME,   YES
> OBJ,    PRESSURE,    0.
> OBJ,    VELOCITY,    0. ,-25. ,0.
> OBJ,    TEMPERATURE, 0.
 
> OBJ,    NAME,        B4
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.650500E+00, 6.593945E+00, 1.889835E+00
> OBJ,    GEOMETRY,    fluidreg
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    67, WATER at 20. deg C
STOP