L(220)
  PHOTON USE
  p
 
 
 
 
 
  view z
  red
  msg Velocity vectors
  vec z 1 sh
  pause
  vec off;red
  msg circumferential-velocity contours
  con u1 z 1 fi;0.001
  pause
  con off;red
  msg radial-velocity contours
  con v1 z 1 fi;0.001
  pause
  msg streamlines
  stre 2d z 1;int 50
  g z 1
  ENDUSE
    GROUP 1. Run title and other preliminaries
TEXT(Rectilinear Flow In Polar Quadrant 
TITLE
  DISPLAY
 
  The flow is from west to east across the north-east quadrant.
  An initial field is inserted by way of an INIPOL patch which
  is identical with the exact solution. The numerical solution
  differs (of course) very little.
 
  This Q1 contains a PHOTON USE file.
 
  ENDDIS
 
    GROUP 3. X-direction grid specification
GRDPWR(X,12,1.5708,1.0)
    GROUP 4. Y-direction grid specification
RINNER=0.0;XCYCLE=F
PATCH(INIPOL,INIVAL,1,NX,1,NY,1,NZ,1,1)
INIT(INIPOL,U1,0.0,SETSPEED);INIT(INIPOL,V1,0.0,SETSPEED)
INIADD=F;FIINIT(U1)=-1.5708;FIINIT(V1)=1.0
  ** Inlet
INLET(IN,WEST,1,1,1,NY,1,1,1,1)
VALUE(IN,P1,1.0);VALUE(IN,U1,0.9914)
VALUE(IN,V1,SAME)
  ** Outlet
PATCH(OUTLET,NORTH,1,NX,NY,NY,1,1,1,1)
COVAL(OUTLET,P1,FIXVAL,0.0)
LSWEEP=100
 
NXPRIN=NX/5;NYPRIN=NY/5;NROWCO=40
PATCH(MAP,CONTUR,1,NX,1,NY,1,1,1,1)
PLOT(MAP,P1,1.0,10.0);PLOT(MAP,U1,1.0,10.0);PLOT(MAP,V1,1.0,10.0)
PATCH(OUTERY,PROFIL,1,NX,NY,NY,1,1,1,1)
PLOT(OUTERY,P1,0.0,0.0);PLOT(OUTERY,U1,0.0,0.0)
PATCH(INNERY,PROFIL,1,NX,1,1,1,1,1,1)
PLOT(INNERY,P1,0.0,0.0);PLOT(INNERY,U1,0.0,0.0)
INIFLD=T;TSTSWP=-1