TALK=T;RUN( 1, 1) 
 ************************************************************
   Q1 created by VDI menu, Version 3.5, Date 22/09/03
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  DISPLAY
 
  
  SUPERSONIC FLOW OVER DIAMOND-SHAPE BODY                                 
 
  The supersonic flow over a 2D diamond-shaped body is solved
  here on Cartesian grid by "cut-off" technique of PARSOL.
 
  The case is aimed to demonstrate the PARSOL performance for
  compressible flow along with the relevant treatments and settings.
 
 
  The demonstration case considers the supersonic air flow over a
  diamond-shaped profile placed in a rectangular channel.
 
  In this case the pressure distribution and accompanying velocity
  field have to be calculated along with auxiliary field of density
  to get the Mach number contours.
 
  Symmetry of the flow allows to consider the half of the actual
  flow pattern.
 
  ENDDIS
 ************************************************************
 IRUNN   =         1 ;LIBREF =         0
 ************************************************************
  Group 1. Run Title
 TEXT(SUPERSONIC FLOW OVER DIAMOND-SHAPE BODY )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,60,30,1)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(146) =VABS ; NAME(147) =MACH
 NAME(148) =PTOT ; NAME(149) =SMOK
 NAME(150) =RHO1
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  ,SMOK)
    * Stored variables list
 STORE(RHO1,PTOT,MACH,VABS)
 
 ************************************************************
  Group 8. Terms & Devices
 TERMS (U1  ,Y,Y,N,Y,Y,Y)
 TERMS (V1  ,Y,Y,N,Y,Y,Y)
 NEWRH1  =    T
 ************************************************************
  Group 9. Properties
 PRESS0  = 1.000000E+00
 RHO1    = GRND3
 RHO1A   = 1.000000E+00 ;RHO1B  = 7.140000E-01
 RHO1C   = 0.000000E+00
 ENUL    = 1.000000E-05
 CP1     = 1.005000E+03
 ENUT    = 0.000000E+00
 DRH1DP  = GRND3
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1  ) =  0.000000E+00 ;FIINIT(U1  ) =  2.000000E+00
 FIINIT(SMOK) =  1.000000E+00 ;FIINIT(RHO1) =  1.000000E+00
   No PATCHes used for this Group
 INIADD  =    T
 ************************************************************
  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  =       500
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 0.3)
 RELAX(U1  ,FALSDT, 5.000000E-00)
 RELAX(V1  ,FALSDT, 5.000000E-00)
 RELAX(SMOK,LINRLX, 1.000000E+00)
 ************************************************************
  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
 PARSOL  =    T
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(SMOK,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =         4 ;IYMON  =        21 ;IZMON  =         1
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 NXPRIN  =         4
 NYPRIN  =         2
 IPLTF   =         1 ;IPLTL  =        40
 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,        5.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       2.916670E-01, 6.833330E-01, 5.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,     60, 1.000000E+00
> GRID,   RSET_Y_1,     30, 1.000000E+00
 
> OBJ,    NAME,        IN
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    VELOCITY,      4.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    TEMPERATURE,   0.000000E+00
 
> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    5.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,   SAME
> OBJ,    COEFFICIENT,   FIXVAL
 
> OBJ,    NAME,        B4
> OBJ,    POSITION,    0.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        B5
> OBJ,    POSITION,    1.333000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.700000E-01, 2.400000E-01, 1.000000E+00
> OBJ,    CLIPART,     wedge2
> OBJ,    ROTATION24,        9
> OBJ,    GRID,              2
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      199
 
> OBJ,    NAME,        B6
> OBJ,    POSITION,    7.600000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.740000E-01, 2.400000E-01, 1.000000E+00
> OBJ,    CLIPART,     wedge2
> OBJ,    ROTATION24,       18
> OBJ,    GRID,              2
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      199
STOP
distil=t
 EX(P1  )=   4.431E-01  
 EX(U1  )=   3.691E+00  
 EX(V1  )=   3.671E-01  
 EX(PRPS)=   5.268E+00  
 EX(VABS)=   3.733E+00  
 EX(MACH)=   3.173E+00  
 EX(PTOT)=   5.356E+01  
 EX(SMOK)=   5.637E-20  
 EX(RHO1)=   1.075E+00