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 Internal flow with blockages
 MESG(
 MESG(
 MESG( This case illustrates the benefit of using 
 MESG( the GMRES feature of MIGAL when blockages
 MESG( obstruct the domain. This case is derived 
 MESG( from VR tutorial case V404. 
 MESG(
 MESG(
 MESG(
 MESG( Enter required solver ID :  1 - MIGAL-GMRES (default)
 MESG(                             2 - MIGAL 
 MESG(
 READVDU(SLV,INT,1)

 ************************************************************
 IRUNN   =       1 ;LIBREF =     302
 ************************************************************
  Group 1. Run Title
 IF(SLV.EQ.1) THEN       
 + TEXT(Internal flow with blockages MIGAL/GMRES)
 ELSE
 + TEXT(Internal flow with blockages MIGAL      )
 ENDIF
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,80,40,1)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(148) =PRPS ; NAME(150) =WDIS
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  )
    * Stored variables list
 STORE(WDIS,PRPS)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 TURMOD(LVEL)
 
 ************************************************************
  Group 8. Terms & Devices
 TERMS (P1  ,Y,Y,Y,Y,Y,Y)
 ************************************************************
  Group 9. Properties
 PRESS0  = 1.000000E+05 ;TEMP0  = 2.730000E+02
 SETPRPS(1,  0)
 DRH1DP  = 1.000000E-05
 DVO1DT  = 3.410000E-03
 EMISS   = 0.000000E+00 ;SCATT  = 0.000000E+00
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(PRPS) = -1.000000E+00 ;FIINIT(WDIS) =  1.000000E-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
   No PATCHes used for this Group
 
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =      100
 SELREF  =    F
 ************************************************************
  Group 16. Terminate Iterations
 LITER (LTLS) =     1000
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 RELAX(U1  ,FALSDT, 1.000000E+10)
 RELAX(V1  ,FALSDT, 1.000000E+10)
 ************************************************************
  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
 ASAP    =    T

 SPEDAT(SET,MIGAL,SOLVED1,C,HYDRO)
 SPEDAT(SET,MIGAL,LINRLX1,R,0.90000E+00)
 SPEDAT(SET,MIGAL,RELAX1,R,0.8)
 SPEDAT(SET,MIGAL,NBRELAX1,I,10)
 SPEDAT(SET,MIGAL,RESFAC1,R,1.00000E-02)

 IF(SLV.EQ.1) THEN       
 + LSWEEP  = 50
 + SPEDAT(SET,MIGAL,IGMRES1,I,5)
 ENDIF

 ************************************************************
  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(V1  ,Y,Y,Y,N,Y,Y)
 OUTPUT(WDIS,Y,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =      21 ;IYMON  =       4 ;IZMON  =       1
 NPRMON  =  100000
 NPRMNT  =       1
 TSTSWP  =      -1
 ITABL   = 3
 NPLT    = 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
 NOWIPE  =    T
 
 GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
 GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
 
> DOM,    SIZE,        2.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       1.008750E+00, 1.750000E-01, 5.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
 
> OBJ,    NAME,        N-WALL
> OBJ,    POSITION,    0.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        S-WALL
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        INLET
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube3
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    VELOCITY,      2.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    TEMPERATURE,   0.000000E+00
 
> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    2.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube12
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,  -1.026000E+04
> OBJ,    COEFFICIENT,   1.000000E+03
 
> OBJ,    NAME,        BLOCK-1
> OBJ,    POSITION,    3.600000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B6
> OBJ,    POSITION,    6.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B7
> OBJ,    POSITION,    8.400000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B8
> OBJ,    POSITION,    1.080000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B9
> OBJ,    POSITION,    1.320000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B10
> OBJ,    POSITION,    1.560000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B11
> OBJ,    POSITION,    3.600000E-01, 7.500000E-01, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B12
> OBJ,    POSITION,    6.000000E-01, 7.500000E-01, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B13
> OBJ,    POSITION,    8.400000E-01, 7.500000E-01, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B14
> OBJ,    POSITION,    1.080000E+00, 7.500000E-01, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B15
> OBJ,    POSITION,    1.320000E+00, 7.500000E-01, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        B16
> OBJ,    POSITION,    1.560000E+00, 7.500000E-01, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 2.500000E-01, 1.000000E+00
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198