PHOTON USE
  P;;;;
 
  USE PATGEO
  VEC Z 1 SH
  MSG Velocity Vectors
  MSG Press RETURN to plot pressure contours
  PAUSE;VEC CL;RED
  CON P1 Z 1 FI;.001
  ENDUSE
 ************************************************************
  Group 1. Run Title
 TEXT(Flow through a circular domain      :561
 TITLE
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,30,10,1)
    * Set overall domain extent:
    *        xulast  yvlast  zwlast    name
 XSI= 6.283E+00;YSI= 1.000E+00;ZSI= 1.000E+00;RSET(D,CHAM    )
    * Set objects: x0     y0     z0
    *              dx     dy     dz    name
 XPO= 0.000E+00;YPO= 0.000E+00;ZPO= 0.000E+00
 XSI= 2.485E-01;YSI= 1.513E-01;ZSI= 0.000E+00;RSET(B,B1      )
 XPO= 5.402E+00;YPO= 1.000E+00;ZPO= 0.000E+00
 XSI= 8.813E-01;YSI= 0.000E+00;ZSI= 1.000E+00;RSET(B,B2      )
    * Cylindrical-polar grid
 CARTES=F
 ************************************************************
  Group 6. Body-Fitted coordinates
 BFC=T; NONORT=T
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME( 48) =VABS ; NAME( 49) =PRPS ; NAME( 50) =DEN1
    * Solved variables list
 SOLVE(P1  ,U1  ,V1  )
    * Stored variables list
 STORE(DEN1,PRPS,VABS)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,N)
 ************************************************************
  Group 8. Terms & Devices
 NEWRH1  =    T; NEWENL  =    T; NEWENT  =    T
 ************************************************************
  Group 9. Properties
 RHO1    = FILE; CP1     = 1.005E+03; ENUL    = FILE
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(PRPS) =  0.000E+00
 INIADD  =    F
 ************************************************************
  Group 13. Boundary & Special Sources
 INLET (B1,NORTH,1,3,#2,#2,#1,#1,1,1)
 VALUE (B1,P1  , 10.*1.189E+00); VALUE (B1,V1  ,-10.000E+00)
 
 PATCH (B2,NORTH ,15,17,#2,#2,#1,#1,1,1)
 COVAL (B2,P1  , 1.000E+00, 0.000E+00)
 
 PATCH (B3,NWALL  ,1,NX,#2,#2,#1,#1,1,1)
 COVAL (B3,U1  , 1.000E+01, 0.000E+00)
 
 PATCH (BUOYANCY,PHASEM,#1,#NREGX,#1,#NREGY,#1,#NREGZ,#1,#NREGT)
 COVAL (BUOYANCY,U1  , FIXFLU   , DENSITY)
 COVAL (BUOYANCY,V1  , FIXFLU   , DENSITY)
 
 BUOYA   = 0.000E+00 ; BUOYB =-9.810E+00 ; BUOYC = 0.000E+00
 EGWF    =    T
 XCYCLE=T
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =      2000; RESFAC  = 1.000E-02
 ************************************************************
  Group 16. Terminate Iterations
 LITER (U1  ) =  10 ;LITER (V1  ) =  10
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000E+00); RELAX(U1  ,FALSDT, 6.283E-03)
 RELAX(V1  ,FALSDT, 6.283E-03)
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000E+06 ;VARMIN(U1  ) =-1.000E+06
 VARMAX(V1  ) = 1.000E+06 ;VARMIN(V1  ) =-1.000E+06
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =      15 ;IYMON  =       5 ;IZMON  =       2
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NXPRIN  =       6; NPLT    =       1; ITABL   =       3