TALK=T;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 3.5, Date 15/03/02
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************

INTEGER(SLV)
MESG(
MESG(
MESG(
MESG(    2D Flow-driven Cavity K-EPS Model
MESG(
MESG(
MESG( This case illustrates the speed-up that MIGAL 
MESG( also produces with turbulent flows. The k-esp
MESG( equations are solved with MIGAL. This case is
MESG( derived from library case 291. 
MESG(
MESG(
MESG(
MESG( Enter required solver ID :  1 - MIGAL (default)
MESG(                             2 - SIMPLEST 
MESG(
READVDU(SLV,INT,1)

 ************************************************************
 IRUNN   =       1 ;LIBREF =     291
 ************************************************************
  Group 1. Run Title
IF(SLV.EQ.1) THEN
+ TEXT(Flow-driven Cavity K-EPS MIGAL    )
ELSE
+ TEXT(Flow-driven Cavity K-EPS SIMPLEST )
ENDIF
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,1,120,120,1.000000E-04)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(149) =VPOR ; NAME(150) =ENUT
    * Solved variables list
 SOLVE(P1  ,V1  ,W1  )
    * Stored variables list
 STORE(ENUT,VPOR)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,N)
 TURMOD(KEMODL)
 
 ************************************************************
  Group 8. Terms & Devices
 DIFCUT  = 0.500000E+00
 ************************************************************
  Group 9. Properties
 RHO1    = 1.000000E+00
 ENUL    = 1.000000E-05
 CP1     = 1.005000E+03
 PRT   (EP  ) =  1.314000E+00
 EMISS   = 0.000000E+00 ;SCATT  = 0.000000E+00
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(W1  ) =  8.000000E+00 ;FIINIT(KE  ) =  2.080000E-01
 FIINIT(EP  ) =  6.920000E-01 ;FIINIT(VPOR) =  1.000000E+00
 FIINIT(ENUT) =  1.001000E-10
   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
 
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =      4300
 RESFAC  = 1.000000E-03
 SELREF=F
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 3.500000E-01)
 RELAX(V1  ,FALSDT, 3.000000E-02)
 RELAX(W1  ,FALSDT, 3.000000E-02)
 RELAX(KE  ,FALSDT, 3.000000E-02)
 RELAX(EP  ,FALSDT, 3.000000E-02)
 KELIN   =       0
 ************************************************************
  Group 18. Limits
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 GENK    =    T
 ASAP    =    T

  MIGAL parameters
  ----------------

 IF(SLV.EQ.1) THEN
 + LSWEEP  = 60
 + spedat(migal,solved1,c,hydro)
 + spedat(migal,linrlx1,R,.8)
 + spedat(migal,solved2,c,kemodl)
 + spedat(migal,linrlx2,R,1.)
 ENDIF

 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =       1 ;IYMON  =       54 ;IZMON  =       42
 NPRMON  =      10
 NPRMNT  =       1
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =  100000
 NYPRIN  =       4
 NZPRIN  =       4
 NPLT    =       1
 ITABL   = 3
 ISWPRF  =       1 ;ISWPRL =  100000
 
 PATCH (MAP     ,CONTUR,1,1,1,20,1,20,1,1)
 PLOT(MAP     ,P1  , 0.000000E+00, 1.000000E+01)
 PLOT(MAP     ,V1  , 0.000000E+00, 1.000000E+01)
 PLOT(MAP     ,W1  , 0.000000E+00, 1.000000E+01)
 PLOT(MAP     ,KE  , 0.000000E+00, 1.000000E+01)
 PLOT(MAP     ,EP  , 0.000000E+00, 1.000000E+01)
 PLOT(MAP     ,ENUT, 0.000000E+00, 1.000000E+01)
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,-1.000000E+00,0.000000E+00,0.000000E+00)
 GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
 
> DOM,    SIZE,        1.000000E+00, 2.500000E-01, 2.000000E-01
> DOM,    MONIT,       5.000000E-01, 7.083330E-02, 7.083330E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_Y_1,     72, 1.000000E+00
> GRID,   RSET_Y_2,     48, 1.000000E+00
> GRID,   RSET_Z_1,     24, 1.000000E+00
> GRID,   RSET_Z_2,     72, 1.000000E+00
> GRID,   RSET_Z_3,     24, 1.000000E+00
 
> OBJ,    NAME,        CMP1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 1.500000E-01
> OBJ,    SIZE,        1.000000E+00, 1.000000E-01, 4.999998E-02
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        CMP2
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 1.000000E-01, 5.000000E-02
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        WFUNSOUT
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 5.000000E-02
> OBJ,    SIZE,        1.000000E+00, 0.000000E+00, 1.000000E-01
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        INLET
> OBJ,    POSITION,    0.000000E+00, 1.000000E-01, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 1.500000E-01, 0.000000E+00
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,      1.000000E+00
> OBJ,    VELOCITY,      0.000000E+00, 0.000000E+00, 8.000000E+00
> OBJ,    TEMPERATURE,   0.000000E+00
> OBJ,    KE_IN,         2.080000E-01
> OBJ,    EP_IN,         6.920000E-01
 
> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    0.000000E+00, 1.000000E-01, 2.000000E-01
> OBJ,    SIZE,        1.000000E+00, 1.500000E-01, 0.000000E+00
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,  -1.026000E+04
> OBJ,    COEFFICIENT,   1.000000E+08
> OBJ,    TURBULENCE,    0.000000E+00, 0.000000E+00
 
> OBJ,    NAME,        STREAM
> OBJ,    POSITION,    0.000000E+00, 2.500000E-01, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 0.000000E+00, 2.000000E-01
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,  -1.026000E+04
> OBJ,    COEFFICIENT,   1.000000E+08
> OBJ,    VELOCITY,      0.000000E+00, 0.000000E+00, 8.000000E+00
> OBJ,    TURBULENCE,    1.000000E-10, 1.000000E-10