PHOTON USE
  p
 
 
 
 
  use patgeo
  MSG Velocity vectors
  vec x 1 sh
  msg
  msg Press return to plot pressure contours
  pause
  gr x 1
  cont p1 x 1 fil;.001
  msg
  msg Type e to End
  ENDUSE
 
    GROUP 1. Run title
TEXT(PLANE DUCT + ANGLED SPLITTER PLATE: B513
TITLE
  DISPLAY
  As in the preceeding cases, the flow is in the y-z plane. The grid
  is distorted to follow the angle of the splitter plate. The flow
  blockage is represented by porosities in Group 11, and wall
  friction on both sides of the plate is activated in Group 13.
  ENDDIS
 
    GROUP 6. Body-fitted coordinates or grid distortion
BFC=T;NONORT=T
GSET(D,1,5,5,1.0,.05,.2)
GSET(C,J3,F,J3,1,NX,NZ,NZ,+,0.0,0.0075,0.0)
GSET(C,J2,F,J2,1,NX,NZ,NZ,+,0.0,0.00375,0.0)
GSET(T,J6,F,J3,1,NX,NZ,NZ,1.0)
GSET(T,K5,F,K2,1,NX,1,NY,1.0)
    GROUP 7. Variables stored, solved & named
   ** Solve for pressure (whole-field) and velocity.
SOLVE(P1,V1,W1);SOLUTN(P1,Y,Y,Y,N,N,N)
    GROUP 9. Properties of the medium (or media)
RHO1=1.0;ENUL=1.0E-04
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=1.0
   ** Block the areas covered by the plate
CONPOR(SPLT,0.0,NORTH,-#1,-#NREGX,-#2,-#2,-#2,-#2)
    GROUP 13. Boundary conditions and special sources
   ** Inlet
INLET(BFCIN,LOW,#1,#NREGX,#1,#NREGY,#1,#1,1,1)
VALUE(BFCIN,P1,GRND1);VALUE(BFCIN,W1,GRND1)
VALUE(BFCIN,WCRT,1.0)
  *  Transfer density for GXBFC subroutine
BFCA=RHO1
   ** Outlet
PATCH(OUTLET,HIGH,#1,#NREGX,#1,#NREGY,#NREGZ,#NREGZ,1,1)
COVAL(OUTLET,P1,1.0E5,0.0)
COVAL(OUTLET,V1,ONLYMS,0.0);COVAL(OUTLET,W1,ONLYMS,0.0)
   ** Wall function
WALL (WALL,NORTH,#1,#NREGX,#NREGY,#NREGY,#1,#NREGZ,1,1)
COVAL(WALL,W1,1.0,0.0)
   ** Wall function - South wall
WALL (WALL3,SOUTH,#1,#NREGX,#1,#1,#1,#NREGZ,1,1)
COVAL(WALL3,W1,1.0,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=20
    GROUP 22. Spot-value print-out
IYMON=3;IZMON=3;TSTSWP=-1
    GROUP 23. Field print-out and plot control
PATCH(YZ,CONTUR,1,1,1,NY,1,NZ,1,1)
PLOT(YZ,P1,0.0,20.0)