************************************************************
   Q1 created by VDI menu, Version 3.2, Date 08/01/99
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN   =       1 ;LIBREF =     302
 ************************************************************
  Group 1. Run Title
 TEXT(Flow In A Closed Cavity          
 ************************************************************
  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,20)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
 XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00
 RSET(D,CHAM    )
    * Set objects: x0     y0     z0
    *              dx     dy     dz
                   name
 XPO= 1.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
 XSI= 0.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00
 RSET(B,B1      )
 XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
 XSI= 0.000000E+00; YSI= 1.000000E+00; ZSI= 1.000000E+00
 RSET(B,B2      )
 XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 1.000000E+00
 XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 0.000000E+00
 RSET(B,B3      )
 XPO= 0.000000E+00; YPO= 0.000000E+00; ZPO= 0.000000E+00
 XSI= 1.000000E+00; YSI= 1.000000E+00; ZSI= 0.000000E+00
 RSET(B,B4      )
 ************************************************************
  Group 6. Body-Fitted coordinates
    * X-cyclic boundaries switched
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(149) =EPKE ; NAME(150) =TEM1
    * Solved variables list
 SOLVE(P1  ,U1  ,W1  ,TEM1)
    * Stored variables list
 STORE(EPKE)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)
 TURMOD(KEMODL)
 
 ************************************************************
  Group 8. Terms & Devices
 TERMS (TEM1,Y,Y,Y,Y,Y,Y)
 NEWRH1  =    T
 NEWENL  =    T
 ************************************************************
  Group 9. Properties
 SETPRPS(1,  2)
 RHO1    = GRND5
 RHO1A   = 0.000E+00 ;RHO1B  = 3.484E-03 ;RHO1C  = 7.143E-01
 PRESS0  = 1.000E+05
 TEMP0   = 2.730E+02
 CP1     = 1.004E+03
 ENUL    = GRND2
 ENULA   =-4.947E-06 ;ENULB  = 4.584E-08 ;ENULC  = 8.097E-11
 DRH1DP  = GRND5
 DVO1DT  = 3.330E-03
 PRNDTL(TEM1) = -2.630E-02
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(EPKE) =  1.000E+00 ;FIINIT(TEM1) =  2.000E+01
   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
 
 PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1)
 
 BUOYA   = 0.000E+00 ; BUOYB = 0.000E+00 ; BUOYC =-9.810E+00
 BUOYD   = 1.276E+00
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =     200
 RESFAC  = 1.000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 RELAX(U1  ,FALSDT, 1.000000E+00)
 RELAX(W1  ,FALSDT, 1.000000E+00)
 RELAX(KE  ,LINRLX, 5.000000E-01)
 RELAX(EP  ,LINRLX, 5.000000E-01)
 RELAX(TEM1,FALSDT, 1.000000E+02)
 KELIN   =       3
 ************************************************************
  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
 GENK    =    T
 ASAP    =    T
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(P1  ,Y,Y,Y,N,Y,Y)
 OUTPUT(U1  ,Y,Y,Y,N,Y,Y)
 OUTPUT(W1  ,Y,Y,Y,N,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =      10 ;IYMON  =       1 ;IZMON  =      10
 NPRMNT  =       1
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 NOWIPE  =    T
 
 GVIEW(P,0.000000E+00,-1.000000E+00,0.000000E+00)
 GVIEW(UP,0.000000E+00,0.000000E+00,1.000000E+00)
 
> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       4.750000E-01, 5.000000E-01, 4.750000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> DOM,    RELAX,       5.000000E-01
 
> OBJ1,   NAME,        B1
> OBJ1,   POSITION,    1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ1,   SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ1,   CLIPART,     cube13
> OBJ1,   ROTATION,          1
> OBJ1,   TYPE,        PLATE
> OBJ1,   SURF_TEMP,     0.000000E+00, 2.000000E+01
 
> OBJ2,   NAME,        B2
> OBJ2,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2,   SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ2,   CLIPART,     cube13
> OBJ2,   ROTATION,          1
> OBJ2,   TYPE,        PLATE
> OBJ2,   SURF_TEMP,     0.000000E+00, 2.000000E+01
 
> OBJ3,   NAME,        B3
> OBJ3,   POSITION,    0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ3,   SIZE,        1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ3,   CLIPART,     cube13
> OBJ3,   ROTATION,          1
> OBJ3,   TYPE,        PLATE
> OBJ3,   SURF_TEMP,     0.000000E+00, 3.000000E+01
> OBJ3,   VELOCITY,      2.000000E+00, 0.000000E+00, 0.000000E+00
 
> OBJ4,   NAME,        B4
> OBJ4,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ4,   SIZE,        1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ4,   CLIPART,     cube13
> OBJ4,   ROTATION,          1
> OBJ4,   TYPE,        PLATE
> OBJ4,   SURF_TEMP,     0.000000E+00, 2.000000E+01
 
> OBJ5,   NAME,        B5
> OBJ5,   POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ5,   SIZE,        5.000000E-02, 1.000000E+00, 5.000000E-02
> OBJ5,   CLIPART,     cubet
> OBJ5,   ROTATION,          1
> OBJ5,   TYPE,        PRESSURE_RELIEF
> OBJ5,   PRES_RELIEF,   0.000000E+00, 0.000000E+00
STOP
DISTIL=T
 EX(P1  )=   4.938E-02 
 EX(U1  )=   9.619E-02 
 EX(W1  )=   9.743E-02 
 EX(KE  )=   2.825E-03 
 EX(EP  )=   5.043E-03 
 EX(EPKE)=   4.549E-01 
 EX(TEM1)=   2.544E+01