PHOTON USE
  p;;;;
 
  GR OU Z 1
  msg   Y-X POLAR GRID  DOWNFLOW  V RAD. , U AZI.
  MSG          ---   FLOW DOMAIN   ---
  MSG Press RETURN to plot velocity vectors
  PAUSE;RED
  VEC Z 1 SH
  MSG Press E to end
  ENDUSE
 ************************************************************
  Group 1. Run Title
TEXT(POLAR GRID DOWNFLOW V RAD., U AZI.  :562
TITLE
 ************************************************************
  Groups 3, 4, 5  Grid Information
NX=6 ; NY=5
 ************************************************************
  Group 6. Body-Fitted coordinates
BFC=T; NONORT=T
REAL(PI,VIN); PI=3.14159; VIN=-1.0
INTEGER(IXFIN,IXLIN,IXFU,IXLU)
SETLIN(XC,LNJ*SIN(PI*LNI));SETLIN(YC,LNJ*COS(PI*LNI))
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
SOLVE(P1,U1,V1); STORE(UCRT,VCRT)
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=5.E-5; BFCA=RHO1
PATCH(IBFC,INIVAL,1,NX,1,NY,1,1,1,1)
INIT(IBFC,U1,0.0,GRND1); INIT(IBFC,V1,0.0,GRND1)
INIT(IBFC,UCRT,0.0,0.0); INIT(IBFC,VCRT,0.0,VIN)
 ************************************************************
  Group 13. Boundary & Special Sources
PATCH(BFCIN,NORTH,1,NX/2,NY,NY,1,NZ,1,1)
COVAL(BFCIN,P1,FIXFLU,GRND1); COVAL(BFCIN,V1,FIXVAL,GRND1)
COVAL(BFCIN,U1,FIXVAL,GRND1); COVAL(BFCIN,VCRT,0.0,VIN)
COVAL(BFCIN,UCRT,0.0,0.0)
PATCH(BFCOUT,NORTH,NX/2+1,NX,NY,NY,1,NZ,1,1)
COVAL(BFCOUT,P1,1.E8,0.0)
COVAL(BFCOUT,V1,0.0,0.0);COVAL(BFCOUT,U1,0.0,0.0)
 ************************************************************
  Group 15. Terminate Sweeps
LSWEEP=250;NPRINT=LSWEEP;NPLT=1;IPLTL=LSWEEP
TSTSWP=-1