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(Cubical Building, K-E model, uniform Flow          )
 ************************************************************
  Group 2. Transience
 STEADY = T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,20,20,10)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
    * Non-default variable names
 NAME(145)=GEN1 ;NAME(146)=VPOR
 NAME(147)=EPKE ;NAME(148)=DEN1
 NAME(149)=EL1 ;NAME(150)=ENUT
    * Solved variables list
 SOLVE(P1,U1,V1,W1)
    * Stored variables list
 STORE(ENUT,EL1,DEN1,EPKE,VPOR,GEN1)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 TURMOD(KEMODL)
 
 ************************************************************
  Group 8. Terms & Devices
 USP = T ; USPDBG= F
 UTCPLT = T ; USPIMB= F
 UAUTO = T ; UGTGRD= F
 USPVTK = T ; USPTER= F
 USCANX = T ; USCANY= T
 USCANZ = T
 MXLEV = 0 ;MYLEV = 0 ;MZLEV = 0
 DOMAT = -1
 CELLST = 10 ; FACEST= 10
 MINPRP = -1 ; MAXPRP= 250
 UMXLEV = 2 ; UMXCLS= 10000
 UNOLAY = 3 ; UIOLEV= 1
 UWALEV = 1 ; 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)
 RHO1 =1.25
 ENUL =1.41E-05
 TMP1 = GRND1
 TMP1A =10. ;TMP1B =0.
 TMP1C =0.
 DVO1DT =3.41E-03
 PRT(EP)=1.314
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(VPOR)=1. ;FIINIT(EPKE)=9.
   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.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=500
 ************************************************************
  Group 17. Relaxation
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.4)
 RELAX(U1  ,FALSDT, 1.E-01)
 RELAX(V1  ,FALSDT, 1.E-01)
 RELAX(W1  ,FALSDT, 1.E-01)
 RELAX(KE ,LINRLX,0.5 )
 RELAX(EP ,LINRLX,0.5 )
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

else
 RELAX(P1 ,LINRLX,1. )
 RELAX(U1 ,FALSDT,10. )
 RELAX(V1 ,FALSDT,10. )
 RELAX(V1 ,FALSDT,10. )
 RELAX(KE ,LINRLX,0.5 )
 RELAX(EP ,LINRLX,0.5 )
endif

 KELIN = 3
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 GENK = T
 PARSOL = F
 CONWIZ = T
      RSG11 =0.2
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(GEN1,Y,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 10 ;IYMON = 12 ;IZMON = 3
 NPRMON = 100000
 NPRMNT = 1
 TSTSWP = -1
 
      #maxmin
 
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000
 ISWPRF = 1 ;ISWPRL = 100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 ************************************************************
  Group 25. DEBUG Group
 GVIEW(P,-0.987856,-0.09368,0.123949)
 GVIEW(UP,0.130697,-0.069691,0.98897)
 
> DOM,    SIZE,        3.500000E+02, 2.700000E+02, 9.000000E+01
> DOM,    MONIT,       1.735000E+02, 1.750000E+02, 1.950000E+01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID,   RSET_X_1,     20, 1.000000E+00
> GRID,   RSET_Y_1,     20,-1.300000E+00
> GRID,   RSET_Z_1,     10, 1.000000E+00
 
> OBJ,    NAME,        WIND
> OBJ,    POSITION,    0.000000E+00, 2.700000E+02, 0.000000E+00
> OBJ,    SIZE,        3.500000E+02, 0.000000E+00, 9.000000E+01
> OBJ,    GEOMETRY,    cube3t
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,    0.
> OBJ,    VELOCITY,    0. ,-5. ,0.
> OBJ,    TURB-INTENS, 5.
 
> OBJ,    NAME,        CUBUS
> OBJ,    POSITION,    1.600000E+02, 9.000000E+01, 0.000000E+00
> OBJ,    SIZE,        3.000000E+01, 3.000000E+01, 3.000000E+01
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    COLOR-MODE,  USER
> OBJ,    COLOR-VAL,   32
> OBJ,    MATERIAL,    198,Solid with smooth-wall friction
 
> OBJ,    NAME,        LEE
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        3.500000E+02, 0.000000E+00, 9.000000E+01
> OBJ,    GEOMETRY,    cube12t
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,    0.
> OBJ,    COEFFICIENT, 1000.
> OBJ,    TURBULENCE,  SAME , SAME
STOP