TALK=T;RUN( 1, 1)
  DISPLAY
  
        
  Air flows through a porous plate with an area of 2m^2, a porosity of
  0.15, and a loss coefficient of 0.4. The expected pressure drop is:

        Dp = 0.5*K*Vel^2 = 0.5*0.4*(18/0.15)^2 = 2880

  The resistive force on the U1 should be dP*Area*porosity

        F = 2880*2*0.15 = 864

  These values are exactly those which are computed and printed in 
  RESULT.


  ENDDIS
 
 ************************************************************
   Q1 created by VDI menu, Version 3.3, Date 10/02/00
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN   =       1 ;LIBREF =       0
 ************************************************************
  Group 1. Run Title
 TEXT(Pressure Drop Across Porous Plate
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,10,3,3)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast
            name
 XSI= 1.000000E+01; YSI= 2.000000E+00; ZSI= 1.000000E+00
 RSET(D,CHAM    )
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  ,W1  )
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
 SETPRPS(1,  0)
 RHO1    = 1.189000E+00
 PRESS0  = 1.000000E+05
 TEMP0   = 2.730000E+02
 CP1     = 1.005000E+03
 ENUL    = 1.544000E-05 ;ENUT   = 0.000000E+00
 DVO1DT  = 3.410000E-03
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
   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
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 RELAX(U1  ,FALSDT, 1.000000E+00)
 RELAX(V1  ,FALSDT, 1.000000E+00)
 RELAX(W1  ,FALSDT, 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
 VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 ASAP    =    T
 SPEDAT(SET,SATLIT,TALK,L,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(V1  ,Y,Y,Y,N,Y,Y)
 OUTPUT(W1  ,Y,Y,Y,N,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =       1 ;IYMON  =       2 ;IZMON  =       2
 NPRMON  =  100000
 NPRMNT  =       1
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =  100000
 NXPRIN  =       1
 ISWPRF  =       1 ;ISWPRL =  100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 NOWIPE  =    T
 
 GVIEW(P,7.390102E-01,-6.736943E-01,0.000000E+00)
 GVIEW(UP,0.000000E+00,0.000000E+00,1.000000E+00)
 
> DOM,    SIZE,        1.000000E+01, 2.000000E+00, 1.000000E+00
> DOM,    MONIT,       5.000000E-01, 1.000000E+00, 5.000000E-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,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ1,   SIZE,        0.000000E+00, 2.000000E+00, 1.000000E+00
> OBJ1,   CLIPART,     cube3t
> OBJ1,   ROTATION24,        1
> OBJ1,   TYPE,        INLET
> OBJ1,   PRESSURE,      0.000000E+00
> OBJ1,   VELOCITY,      1.800000E+01, 0.000000E+00, 0.000000E+00
> OBJ1,   TEMPERATURE,   0.000000E+00
 
> OBJ2,   NAME,        B2
> OBJ2,   POSITION,    5.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ2,   SIZE,        0.000000E+00, 2.000000E+00, 1.000000E+00
> OBJ2,   CLIPART,     cube11
> OBJ2,   ROTATION24,        1
> OBJ2,   TYPE,        PLATE
> OBJ2,   POROSITY,      1.500000E-01
> OBJ2,   RESISTANCE,    4.000000E-01
> OBJ2,   SIDE,        BOTH
 
> OBJ3,   NAME,        B3
> OBJ3,   POSITION,    1.000000E+01, 0.000000E+00, 0.000000E+00
> OBJ3,   SIZE,        0.000000E+00, 2.000000E+00, 1.000000E+00
> OBJ3,   CLIPART,     cube12t
> OBJ3,   ROTATION24,        1
> OBJ3,   TYPE,        OUTLET
> OBJ3,   PRESSURE,      0.000000E+00
> OBJ3,   TEMPERATURE,  -1.026000E+04
> OBJ3,   COEFFICIENT,   1.000000E+03
 
M
STOP
 distil=t
 EX(P1  )=   2.277E+03 
 EX(U1  )=   2.933E+01 
 EX(V1  )=   2.097E-04 
 EX(W1  )=   1.461E-04 
 EX(EPOR)=   9.150E-01