TALK=T;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 3.6, Date 21/10/05
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
    DISPLAY
     This unsteady case considers the prediction of vortex
     shedding past a square cylinder at Re = 22000 by means
     of a 2D Large Eddy Simulation. The cylinder dimension
     is 0.4m and the free-stream water velocity is 0.535m/s.
 
     The computations are initiated from solution fields
     produced by a 2000-sweep, unconverged steady computation.
     The basic Smagorinsky model is used with a coefficient
     of 0.1 and no wall damping. The linear upwind scheme is
     used for the momentum equations, and the near-wall region
     is bridged by standard wall functions. The simulation
     is advanced in time for 7 seconds using a uniform time
     step of 1ms, corresponding to roughly 14 shedding cycles.
     The INFORM facility is used to reduce gradually with
     time the number of sweeps per time step.
 
     The main results are listed below:
                                             LES      Data
 
        Strouhal number                      0.15     0.13
        Mean Drag Coefficient                2.17     2.10
        RMS  Drag Coefficient                0.14     0.12
        RMS  Lift Coefficient                1.17     1.20
 
     This run takes 2.25 hours to complete the 7000 time steps
     on a 3GHZ PC with 2GB RAM.
    ENDDIS
 ************************************************************
 IRUNN   =         1 ;LIBREF =         0
 ************************************************************
  Group 1. Run Title
 TEXT(2D LES SQUARE CYLINDER- unsteady run    )
 ************************************************************
  Group 2. Transience
 STEADY=F
    * Set overall time and no. of steps
 RSET(U,0.000000E+00,7.000000E+00,7000)
    * Modify regions
 RSET(T,1,7000,1.000000E+00)
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,120,85,1)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(137) =FMU  ; NAME(138) =AREN
 NAME(139) =AREE ; NAME(140) =GEN1
 NAME(142) =SHRY ; NAME(143) =SHRX
 NAME(144) =STRS ; NAME(145) =SKIN
 NAME(146) =YPLS ; NAME(147) =WDIS
 NAME(148) =DEN1 ; NAME(149) =EL1
 NAME(150) =ENUT
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  )
    * Stored variables list
 STORE(ENUT,EL1 ,DEN1,WDIS,YPLS,SKIN,STRS,SHRX)
 STORE(SHRY,GEN1,AREE,AREN,FMU )
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 TURMOD(SGSMOD_NOWD)
 
 ************************************************************
  Group 8. Terms & Devices
 DIFCUT  = 0.000000E+00
 SCHEME(LUS   ,U1  ,V1  )
 ************************************************************
  Group 9. Properties
 PRESS0  = 1.000000E+05 ;TEMP0  = 2.730000E+02
    * Domain material index is  67 signifying:
    * WATER at 20. deg C
 SETPRPS(1, 67)
 EL1A    = 1.000000E-01 ;EL1B   = 0.000000E+00
 EL1C    = 0.000000E+00
 DVO1DT  = 1.180000E-04
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 RESTRT(ALL)
   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.000000E-03
 ************************************************************
  Echo InForm settings for Group 15
  inform15begin
(stored lsweep is 200 with IF(ISTEP.gt.3.and.ISTEP.le.6)!TSTSTR)
(stored lsweep is 100 with IF(ISTEP.GT.6.AND.ISTEP.LE.8)!TSTSTR)
(stored lsweep is  50 with IF(ISTEP.GT.8.AND.ISTEP.LE.10)!TSTSTR)
(stored lsweep is  40 with IF(ISTEP.GT.10.AND.ISTEP.LE.14)!TSTSTR)
(stored lsweep is 30 with IF(ISTEP.GT.14.AND.ISTEP.LE.20)!TSTSTR)
(stored lsweep is 20 with IF(ISTEP.GT.20)!TSTSTR)
  inform15end
 ************************************************************
  Group 16. Terminate Iterations
 LITER (P1  ) =      300
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 GENK    =    T
 ASAP    =    T
 PARSOL  =    F
 CONWIZ  =    T
 CALFOR  =    T
 IDISPB  =         1 ;IDISPC =      7000
 CG(  1) ='DRAG'
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 DISTIL  =    T  ;NULLPR =    F
 NDST    =         0
 DSTTOL  = 1.000000E-02
 EX    (P1  ) =  0.000000E+00 ;EX    (U1  ) =  0.000000E+00
 EX    (V1  ) =  0.000000E+00 ;EX    (FMU ) =  0.000000E+00
 EX    (AREN) =  0.000000E+00 ;EX    (AREE) =  0.000000E+00
 EX    (GEN1) =  0.000000E+00 ;EX    (LTLS) =  0.000000E+00
 EX    (SHRY) =  0.000000E+00 ;EX    (SHRX) =  0.000000E+00
 EX    (STRS) =  0.000000E+00 ;EX    (SKIN) =  0.000000E+00
 EX    (YPLS) =  0.000000E+00 ;EX    (WDIS) =  0.000000E+00
 EX    (DEN1) =  0.000000E+00 ;EX    (EL1 ) =  0.000000E+00
 EX    (ENUT) =  0.000000E+00
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(GEN1,Y,N,N,N,N,N)
 OUTPUT(WDIS,Y,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =        68 ;IYMON  =        43 ;IZMON  =         1
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 NTPRIN  =    100000 ;ISTPRF =         1 ;ISTPRL =     10000
 NXPRIN  =         1
 IXPRF   =        20 ;IXPRL  =        56
 NYPRIN  =         1
 IYPRF   =        25 ;IYPRL  =        60
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 IDISPA  =       500 ;IDISPB =         1 ;IDISPC =      7000
 CSG1    ='m'
 
 GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
 GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
 
> DOM,    SIZE,        8.000000E-01, 5.600000E-01, 4.000000E-01
> DOM,    MONIT,       2.607040E-01, 2.800000E-01, 2.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,     27, 9.000000E-01,G
> GRID,   RSET_X_2,     22, 1.000000E+00
> GRID,   RSET_X_3,     71, 1.040000E+00,G
> GRID,   RSET_Y_1,     30, 9.000000E-01,G
> GRID,   RSET_Y_2,     25, 1.000000E+00
> GRID,   RSET_Y_3,     30, 1.100000E+00,G
 
> OBJ,    NAME,        SQCL
> OBJ,    POSITION,    1.800000E-01, 2.600000E-01, 0.000000E+00
> OBJ,    SIZE,        4.000000E-02, 4.000000E-02, 4.000000E-01
> OBJ,    GEOMETRY,    cube
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
> OBJ,    TIME_LIMITS,   ALWAYS_ACTIVE
 
> OBJ,    NAME,        INLET
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 5.600000E-01, 4.000000E-01
> OBJ,    GEOMETRY,    cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    VELOCITY,     5.350000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    TEMPERATURE,  0.000000E+00
> OBJ,    TIME_LIMITS,   ALWAYS_ACTIVE
 
> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    8.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 5.600000E-01, 4.000000E-01
> OBJ,    GEOMETRY,    cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    TEMPERATURE,  SAME
> OBJ,    COEFFICIENT,  1.000000E+03
> OBJ,    TIME_LIMITS,   ALWAYS_ACTIVE
STOP
  ** These parameters are used in GROUND and define the inlet veloc$
     and cylinder dimension, respectively.
     RG(1)=0.535;RG(2)=0.04