TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2007, Date 16/01/07
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN   =         1 ;LIBREF =         0
 ************************************************************
  Group 1. Run Title
 TEXT(Test 20. Laminar flow around of cylinder.
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,20,1,30)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 NAME(148) =DEN1
    * Solved variables list
 SOLVE(P1  ,U1  ,W1  )
    * Stored variables list
 STORE(DEN1)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
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
else
lAUTO=F
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  =         0 ; UMXCLS=    100000
UNOLAY = 2
UIOLEV = 1
UWALEV = 3
UMXCRS = 16

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
 PRESS0  = 1.000000E+05 ;TEMP0  = 2.730000E+02
    * Domain material index is   0 signifying:
    * Air at 20 deg C, 1 atm, treated as incompressible
 SETPRPS(1,  0)
 ENUT    = 0.000000E+00
 DVO1DT  = 3.410000E-03
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
   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  =       500
 RESFAC  = 1.000000E-06
 ************************************************************
  Group 16. Terminate Iterations
            LITER (P1  ) =      200
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 7.000000E-01)
 RELAX(U1,FALSDT,1.E-00)
 RELAX(W1,FALSDT,1.E-00)
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, 3.E-01)
 RELAX(V1  ,FALSDT, 3.E-01)
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   =        10 ;IYMON  =         1 ;IZMON  =        24
 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.995500E-01,-2.999551E-02)
 GVIEW(UP,1.000000E+00,0.000000E+00,0.000000E+00)
 
> DOM,    SIZE,        2.000000E-01, 1.000000E-02, 3.000000E-01
> DOM,    MONIT,       9.500000E-02, 5.000000E-03, 2.650000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   AUTO,         F T F
 
> OBJ,    NAME,        CYL
> OBJ,    POSITION,    2.000000E-02, 0.000000E+00, 5.000000E-02
> OBJ,    SIZE,        1.600000E-01, 1.000000E-02, 1.600000E-01
> OBJ,    GEOMETRY,    cylinder
> OBJ,    ROTATION24,        7
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        INLET
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E-01, 1.000000E-02, 0.000000E+00
> OBJ,    GEOMETRY,    cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    VELOCITY,     0.000000E+00, 0.000000E+00, 1.000000E-02
> OBJ,    TEMPERATURE,  0.000000E+00
 
> OBJ,    NAME,        OUT
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 3.000000E-01
> OBJ,    SIZE,        2.000000E-01, 1.000000E-02, 0.000000E+00
> OBJ,    GEOMETRY,    cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    TEMPERATURE,  SAME
> OBJ,    COEFFICIENT,  1.000000E+03
STOP