TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2007, Date 16/01/07
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN   =         1 ;LIBREF =       128
 ************************************************************
  Group 1. Run Title
 TEXT(Test 23. Flow around small body.    )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
BOOLEAN(lAUTO)
mesg(Do you want to use automatic mode (y) or manual (n)?
readvdu(ans,char,n)
if(:ans:.eq.y)then
lAUTO=T
 RSET(M,10,1,10,1.000000E-04)
else
lAUTO=F
 RSET(M,30,1,28,1.000000E-04)
endif

 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 NAME(150) =PRPS
    * Solved variables list
 SOLVE(P1  ,U1  ,W1  )
    * Stored variables list
 STORE(PRPS)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 SOLUTN(U1  ,Y,Y,Y,N,N,Y)
 SOLUTN(W1  ,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
INTEGER(LEVEL)
LEVEL = 0
mesg(Level of refinement is :LEVEL:
mesg(Do you want to change it? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
 mesg(Enter level of refinement
 readvdu(LEVEL,int,LEVEL)
 mesg(New level of refinement is :LEVEL:
endif

USP    = T
UAUTO  = lAUTO
USPDBG = F
UTCPLT = T
USPVTK = T
USPIMB = F
MXLEV  = 0 
MYLEV  = 0
MZLEV  = 0
DOMAT  = -1
MINPRP = -1
MAXPRP = 100
CELLST = 10
FACEST = 10

UMXLEV  = LEVEL 
UNOLAY = 2
UIOLEV = 0
UWALEV = 0
UMXCRS = 16
UGRITR= 10

if(lAUTO)then
mesg(Do you want to use PARSOL=T (y) or PARSOL=F (n)?
readvdu(ans,char,n)
if(:ans:.eq.y)then
PARSOL=T
else
PARSOL=F
endif
else
PARSOL=F
endif

 ************************************************************
  Group 9. Properties
 RHO1    = 1.000000E+00
 ENUL    = 1.000000E-05
 CP1     = 1.000000E+00
 ENUT    = 0.000000E+00
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(PRPS) = -1.000000E+00
   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  =       800
 RESFAC  = 1.000000E-06
 ************************************************************
  Group 16. Terminate Iterations
 LITER (P1  ) =      200 ;LITER (U1  ) =       50
 LITER (W1  ) =       50
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 7.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, 0.3)
 RELAX(U1  ,FALSDT, 1.E-00)
 RELAX(W1  ,FALSDT, 1.E-00)
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
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-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
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =        26 ;IYMON  =         1 ;IZMON  =         9
 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.000000E+00,-9.999875E-01,4.999979E-03)
 GVIEW(UP,0.000000E+00,4.999979E-03,9.999875E-01)
 
> DOM,    SIZE,        1.000000E+00, 1.000000E-02, 5.000000E-01
> DOM,    MONIT,       5.701250E-01, 5.000000E-03, 2.263850E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
if(lAUTO)then
> GRID,   RSET_X_1,      10, 1.
> GRID,   RSET_Y_1,       1, 1.
> GRID,   RSET_Z_1,      10, 1.
else
> GRID,   RSET_X_1,      3, 5.000000E-01
> GRID,   RSET_X_2,      5, 7.000000E-01
> GRID,   RSET_X_3,     10, 1.000000E+00
> GRID,   RSET_X_4,      7, 1.300000E+00
> GRID,   RSET_X_5,      5, 1.500000E+00
> GRID,   RSET_Y_1,      1, 1.000000E+00,G
> GRID,   RSET_Z_1,      4, 5.000000E-01
> GRID,   RSET_Z_2,      5, 7.000000E-01
> GRID,   RSET_Z_3,     10, 1.000000E+00
> GRID,   RSET_Z_4,      5, 1.300000E+00
> GRID,   RSET_Z_5,      4, 1.500000E+00
endif

> OBJ,    NAME,        INLET
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E-02, 5.000000E-01
> OBJ,    GEOMETRY,    cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
if(lAUTO)then
> OBJ,    GRID,        NO
endif
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    VELOCITY,     1.000000E-02, 0.000000E+00, 0.000000E+00
> OBJ,    TEMPERATURE,  0.000000E+00
 
> OBJ,    NAME,        OUT
> OBJ,    POSITION,    1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E-02, 5.000000E-01
> OBJ,    GEOMETRY,    cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
if(lAUTO)then
> OBJ,    GRID,        NO
endif
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    TEMPERATURE,  SAME
> OBJ,    COEFFICIENT,  1.000000E+03
 
> OBJ,    NAME,        BODY
> OBJ,    POSITION,    4.000000E-01, 0.000000E+00, 2.300000E-01
> OBJ,    SIZE,        4.000000E-02, 1.000000E-02, 4.000000E-02
> OBJ,    GEOMETRY,    cylinder
> OBJ,    ROTATION24,        5
> OBJ,    TYPE,        BLOCKAGE
if(lAUTO)then
> OBJ,    GRID,        NO
endif
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
if(lAUTO)then
else 
> OBJ,    NAME,        NGRID
> OBJ,    POSITION,    3.500000E-01, 0.000000E+00, 1.800000E-01
> OBJ,    SIZE,        2.000000E-01, 1.000000E-02, 1.400000E-01
> OBJ,    GEOMETRY,    wirexyz
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        NULL
endif
STOP