photon use
  p;;;;;
 
  use patgeo; msg Velocity vectors. Press RETURN for pressure
  vec z 1 sh;pause
  msg Pressure contours. Press RETURN for streamlines
  con p1 z 1 fi;0.001
  pause;  msg streamlines
  stre 2d z 1 ;int 30
  enduse
 ************************************************************
  Group 1. Run Title
 TEXT(Simple Labyrinth; Q1 Made By Menu  
   originated by pks/dbs 11.01.94
 TITLE
  DISPLAY
 
  This Q1, which sets up a steady 2D flow through a passage created
  by blockages, has been set up by means of the general input menu.
 
  Enter the menu (Type M, return, R) to see the geometry and
  settings, and to change them if desired.
 
  On re-entering command mode, another means of interpreting the
  settings can be activated by typing: L(20). This is always
  available.
 
  PHOTON USE commands have been included in the Q1.
 
  ENDDIS
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,30,25,1)
    * Overall domain extent, RSET(D,name,XULAST,YVLAST,ZWLAST)
 RSET(D,CHAM,2.500E+00,1.500E+00,1.000E+00)
    * Set objects: name  x0       y0        z0
    *                    dx       dy        dz
 RSET(B,B1      , 0.000E+00, 0.000E+00, 0.000E+00                , $
7.500E-01, 1.000E+00, 1.000E+00)
 RSET(B,B2      , 1.150E+00, 7.500E-01, 0.000E+00                , $
2.000E-01, 7.500E-01, 1.000E+00)
 RSET(B,B3      , 1.750E+00, 0.000E+00, 0.000E+00                , $
7.500E-01, 1.000E+00, 1.000E+00)
 RSET(B,IN1     , 0.000E+00, 1.000E+00, 0.000E+00                , $
0.000E+00, 5.000E-01, 1.000E+00)
 RSET(B,OUT1    , 2.500E+00, 1.000E+00, 0.000E+00                , $
0.000E+00, 5.000E-01, 1.000E+00)
 RSET(B,BWALL   , 0.000E+00, 0.000E+00, 0.000E+00                , $
2.500E+00, 0.000E+00, 1.000E+00)
 RSET(B,TWALL   , 0.000E+00, 1.500E+00, 0.000E+00                , $
2.500E+00, 0.000E+00, 1.000E+00)
    * Modify default grid
 RSET(X,2,5,1.000E+00)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(49) =DEN1 ; NAME(50) =VPOR
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  )
    * Stored variables list
 STORE(DEN1)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,N)
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
 RHO1    = 1.161E+00
 ENUL    = 1.589E-05 ;ENUT   = 3.178E-03
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(VPOR) =  1.000E+00
 
 CONPOR(B1      ,  0.00,VOLUME,#1,#1,#1,#2,#1,#1)
 
 CONPOR(B2      ,  0.00,VOLUME,#3,#3,#2,#3,#1,#1)
 
 CONPOR(B3      ,  0.00,VOLUME,#5,#5,#1,#2,#1,#1)
 
 INIADD  =    F
 ************************************************************
  Group 12. Convection and diffusion adjustments
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH (KESOURCE,PHASEM,1,30,1,25,1,1,1,1)
 
 INLET (IN1     ,WEST  ,#1,#1,#3,#3,#1,#1,#1,#1)
 VALUE (IN1     ,P1  , 5.805E+00)
 VALUE (IN1     ,U1  , 5.000E+00)
 
 PATCH (OUT1    ,EAST  ,#5,#5,#3,#3,#1,#1,#1,#1)
 COVAL (OUT1    ,P1  , 1.000E+00, 0.000E+00)
 
 PATCH (TWALL   ,NWALL ,#1,#5,#3,#3,#1,#1,#1,#1)
 COVAL (TWALL   ,U1  ,LOGLAW , 0.000E+00)
 
 PATCH (BWALL   ,SWALL ,#1,#5,#1,#1,#1,#1,#1,#1)
 COVAL (BWALL   ,U1  ,LOGLAW , 0.000E+00)
 
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =     100
 SELREF  =    T
 RESFAC  = 1.000E-02
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000E+00)
 RELAX(U1  ,FALSDT, 6.000E+00)
 RELAX(V1  ,FALSDT, 6.000E+00)
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 GENK    =    T
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =      19 ;IYMON  =      23 ;IZMON  =       1
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 ************************************************************
 MENSAV(S,RELX,DEF,6.0000E-02,1,100)
 MENSAV(S,PHSPROP,DEF,200,0,1.1610E+00,1.5890E-05)
 MENSAV(S,FLPRP,DEF,CONSTANT_EFFECTIVE,CONSTANT,AIR-CONSTANT)