L(224)
    GROUP 1. Run title
TEXT(2D Flow Over A Cylinder Polar Grid
TITLE
    Viscous flow over a circular cylinder in a semi-infinite domain
  is considered. The first example is laminar with Re=20.
REAL(REL,RET,EPSIN,TKEIN,AA,DY);INTEGER(NXF);VIN=2.0
INTEGER(JJ1)
 
    GROUP 3. X-direction grid specification
XCYCLE=F;NX=21;NXF=9
XFRAC(1)=-NXF;XFRAC(2)=0.5/NXF;XFRAC(3)=NXF;XFRAC(4)=0.4/NXF
XFRAC(5)=3.0;XFRAC(6)=0.1/3.0;XULAST=PI
 
    GROUP 4. Y-direction grid specification
NY=21;AA=1.32;RINNER=0.24;YVLAST=2.4-RINNER
   **First grid line interval is 1.25% of cylinder radius
     intervals increase in geometric progression
DY=RINNER*.0125/YVLAST
 
YFRAC(1)=DY
DO JJ=2,19
+ JJ1=JJ-1;YFRAC(JJ)=YFRAC(JJ1)+DY;DY=DY*AA
ENDDO
YFRAC(20)=0.92;YFRAC(21)=1.0
 
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1);SOLUTN( P1,Y,Y,Y,N,N,N )
 
    GROUP 9. Properties of the medium (or media)
REL=20.0;RHO1=1000;ENUL=2.0*RINNER/REL
 
    GROUP 13. Boundary conditions and special sources
PATCH(PINF,CELL,1,NX,NY,NY,1,NZ,1,1);COVAL(PINF,P1,1.0E6,0.0)
COVAL(PINF,U1,ONLYMS,0.0);COVAL(PINF,V1,ONLYMS,0.0)
PATCH(UPOL,CELL,1,NXF,NY,NY,1,NZ,1,1)
COVAL(UPOL,U1,FIXVAL,SETSPEED)
PATCH(VPOL,CELL,1,NXF,NY-1,NY-1,1,NZ,1,1)
COVAL(VPOL,V1,FIXVAL,SETSPEED)
WALL (WALL,SOUTH,1,NX,1,1,1,NZ,1,1)
COVAL(WALL,U1,1.0,0.0)
 
    GROUP 16. Termination of iterations
LITER(P1) = 12;LITER(U1)=1;LITER(V1)=1
 
    GROUP 17. Under-relaxation devices
RELAX(U1,FALSDT,0.02);RELAX(V1,FALSDT,0.02)
 
    GROUP 19. Special calls from EARTH to GROUND
   ** V1=-VIN at x=0
POLRA=-VIN
 
    GROUP 22. Spot-value print-out
IXMON = 2;IYMON=12;IZMON=1
 
    GROUP 23. Field print-out and plot control
NPLT=2;NYPRIN=NY/5
PATCH(MAP,CONTUR,1,NX,1,NY,1,1,1,1)